Changeset 4320 for branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro
- Timestamp:
- 2013-11-28T12:41:27+01:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro
r2755 r4320 15 15 PRINT, '' 16 16 ; 17 std_iodir_data = isadirectory(getenv('DIR_DATA'), title = 'path of data in NetCdf format')18 std_iodir_climato = isadirectory(getenv('DIR_CLIMATO'), title = 'path of climatological data')19 std_iodir_mask = isadirectory(getenv('DIR_MASK'), title = 'path of mask files (ex: subbasins)')17 std_iodir_data = isadirectory(getenv('DIR_DATA'), title = 'path of data in NetCdf format') 18 std_iodir_climato = isadirectory(getenv('DIR_CLIMATO'), title = 'path of climatological data') 19 std_iodir_mask = isadirectory(getenv('DIR_MASK'), title = 'path of mask files (ex: subbasins)') 20 20 ; meshmask 21 std_file_mesh = isafile(getenv('FILE_MESH_MASK'), title = 'mesh_mask', iodir = std_iodir_mask)21 std_file_mesh = isafile(getenv('FILE_MESH_MASK'), title = 'mesh_mask', iodir = std_iodir_mask) 22 22 std_file_msksub = isafile(getenv('FILE_MASK_SUBDOMAIN'), title = 'sub-bassin masks', iodir = std_iodir_mask) 23 23 24 ; Levitus 98 25 std_file_Levitus_T = isafile(getenv('FILE_TEMP_3D'), title = 'Levitus_T', iodir = std_iodir_climato) 26 std_file_Levitus_S = isafile(getenv('FILE_SAL_3D'), title = 'Levitus_S', iodir = std_iodir_climato) 27 std_file_reynolds = isafile(getenv('FILE_SST'), title = 'Reynolds', iodir = std_iodir_climato) 28 std_file_oaflux = isafile(getenv('FILE_FLUX'), title = 'oaflux', iodir = std_iodir_climato) 29 std_file_mld = isafile(getenv('FILE_MLD'), title = 'Mixed layer depth', iodir = std_iodir_climato) 24 ; climatologies 25 std_file_Levitus_T = isafile(getenv('FILE_TEMP_3D'), title = 'Levitus_T', iodir = std_iodir_climato) 26 std_file_Levitus_S = isafile(getenv('FILE_SAL_3D'), title = 'Levitus_S', iodir = std_iodir_climato) 27 std_file_reynolds = isafile(getenv('FILE_SST'), title = 'Reynolds', iodir = std_iodir_climato) 28 std_file_oaflux = isafile(getenv('FILE_FLUX'), title = 'oaflux', iodir = std_iodir_climato) 29 std_file_mld = isafile(getenv('FILE_MLD'), title = 'Mixed layer depth', iodir = std_iodir_climato) 30 std_file_ice = isafile(getenv('FILE_ICE'), title = 'ICE', iodir = std_iodir_climato) 31 std_file_snow_arc = isafile(getenv('FILE_SNOW_ARC'), title = 'SNOW_ARC', iodir = std_iodir_climato) 32 std_file_snow_ant = isafile(getenv('FILE_SNOW_ANT'), title = 'SNOW_ANT', iodir = std_iodir_climato) 30 33 31 34 IF strlowcase(getenv('FILE_GEOHEAT')) EQ 'no' THEN std_file_geoheat = 'no' $ … … 155 158 Ithi_1 = read_ncdf(getenv('VAR1_Ithick'), allrecords = allrec, filename = std_file1_I ) 156 159 caldat, time, mm 157 march = where(mm EQ 3, cnt)158 Ithi_ march_1 = {arr:1./float(cnt) * total(Ithi_1.arr[*, *, temporary(march)], 3), unit:Ithi_1.unit}160 april = where(mm EQ 4, cnt) 161 Ithi_april_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 159 162 sept = where(mm EQ 9, cnt) 160 Ithi_sept_1 = {arr:1./float(cnt) * total( Ithi_1.arr[*, *, temporary(sept)], 3), unit:Ithi_1.unit}163 Ithi_sept_1 = {arr:1./float(cnt) * total(reform(Ithi_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_1.unit} 161 164 undefine, Ithi_1 162 165 ; 163 166 IF std_file2_I NE std_file1_I THEN BEGIN 164 Ithi_2 = read_ncdf(getenv('VAR2_Ithick'), allrecords = allrec, filename = std_file2_I )165 caldat, time, mm166 march = where(mm EQ 3, cnt)167 Ithi_march_2 = {arr:1./float(cnt) * total(Ithi_2.arr[*, *, temporary(march)], 3), unit:Ithi_2.unit}168 sept = where(mm EQ 9, cnt)169 Ithi_sept_2 = {arr:1./float(cnt) * total(Ithi_2.arr[*, *, temporary(sept)], 3), unit:Ithi_2.unit}170 undefine, Ithi_2167 Ithi_2 = read_ncdf(getenv('VAR2_Ithick'), allrecords = allrec, filename = std_file2_I ) 168 caldat, time, mm 169 april = where(mm EQ 4, cnt) 170 Ithi_april_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Ithi_2.unit} 171 sept = where(mm EQ 9, cnt) 172 Ithi_sept_2 = {arr:1./float(cnt) * total(reform(Ithi_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ithi_2.unit} 173 undefine, Ithi_2 171 174 ENDIF ELSE BEGIN 172 Ithi_march_2 = {arr:-1}173 Ithi_sept_2 = {arr:-1}175 Ithi_april_2 = {arr:-1} 176 Ithi_sept_2 = {arr:-1} 174 177 ENDELSE 175 178 ; 176 179 Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 180 help, Ifra_1.arr, jpi,jpj 177 181 caldat, time, mm 178 182 march = where(mm EQ 3, cnt) 179 Ifra_march_1 = {arr:1./float(cnt) * total( Ifra_1.arr[*, *, temporary(march)], 3), unit:Ifra_1.unit}183 Ifra_march_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} 180 184 sept = where(mm EQ 9, cnt) 181 Ifra_sept_1 = {arr:1./float(cnt) * total( Ifra_1.arr[*, *, temporary(sept)], 3), unit:Ifra_1.unit}185 Ifra_sept_1 = {arr:1./float(cnt) * total(reform(Ifra_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ifra_1.unit} 182 186 undefine, Ifra_1 183 187 ; … … 186 190 caldat, time, mm 187 191 march = where(mm EQ 3, cnt) 188 Ifra_march_2 = {arr:1./float(cnt) * total( Ifra_2.arr[*, *, temporary(march)], 3), unit:Ifra_2.unit}192 Ifra_march_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 189 193 sept = where(mm EQ 9, cnt) 190 Ifra_sept_2 = {arr:1./float(cnt) * total( Ifra_2.arr[*, *, temporary(sept)], 3), unit:Ifra_2.unit}194 Ifra_sept_2 = {arr:1./float(cnt) * total(reform(Ifra_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Ifra_2.unit} 191 195 undefine, Ifra_2 192 196 ENDIF ELSE BEGIN … … 194 198 Ifra_sept_2 = {arr:-1} 195 199 ENDELSE 200 ; 201 ;; Sf commentato perché non funziona snow 202 ;; Sf Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 203 ;; Sf help, Isnow_1.arr, jpi,jpj 204 ;; Sf caldat, time, mm 205 ;; Sf april = where(mm EQ 4, cnt) 206 ;; Sf Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 207 ;; Sf sept = where(mm EQ 9, cnt) 208 ;; Sf Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 209 ;; Sf ;undefine, Isnow_1 210 ;; Sf ; 211 ;; Sf IF std_file2_I NE std_file1_I THEN BEGIN 212 ;; Sf Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 213 ;; Sf caldat, time, mm 214 ;; Sf april = where(mm EQ 4, cnt) 215 ;; Sf Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 216 ;; Sf sept = where(mm EQ 9, cnt) 217 ;; Sf Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 218 ;; Sf ; undefine, Isnow_2 219 ;; Sf ENDIF ELSE BEGIN 220 ;; Sf Isnow_april_2 = {arr:-1} 221 ;; Sf Isnow_sept_2 = {arr:-1} 222 ;; Sf ENDELSE 196 223 ; 197 224 jpt = 1 … … 299 326 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Pac', _extra = ex 300 327 301 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: MARCH' 302 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 303 328 ; cnt = cnt+1 & blabla = 'Arctic Ice Thickness: MARCH' 329 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 330 331 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: APRIL' 332 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 333 304 334 cnt = cnt+1 & blabla = 'Arctic Ice Thickness: SEPT' 305 335 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 306 336 307 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: MARCH' 308 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 337 ; cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: MARCH' 338 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 339 340 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: APRIL' 341 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex 309 342 310 343 cnt = cnt+1 & blabla = 'Antarctic Ice Thickness: SEPT' … … 322 355 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: SEPT' 323 356 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 357 ;; SF non funziona, dice che deve essere una struttura e non e!!!! 358 ;;SF cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: APRIL' 359 ;;SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex 360 ;; SF 361 ;; SF cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: SEPT' 362 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 363 ;; SF 364 ;; SF cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: APRIL' 365 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex 366 ;; SF 367 ;; SF cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT' 368 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 369 ; 370 ; cnt = cnt+1 & blabla = 'Arctic Ice Extend: ' 371 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC, _extra = ex 372 373 ; cnt = cnt+1 & blabla = 'Antarctic Ice Extend: ' 374 ; IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT, _extra = ex 324 375 325 376 cnt = cnt+1 & blabla = 'Meridional stream Function: Global (no Med)'
Note: See TracChangeset
for help on using the changeset viewer.