Changeset 4501 for branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro
- Timestamp:
- 2014-02-20T14:21:03+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro
r4395 r4501 209 209 ENDELSE 210 210 ; 211 ;; Sf commentato perché non funziona snow 212 ;; Sf Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 213 ;; Sf help, Isnow_1.arr, jpi,jpj 214 ;; Sf caldat, time, mm 215 ;; Sf april = where(mm EQ 4, cnt) 216 ;; Sf Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 217 ;; Sf sept = where(mm EQ 9, cnt) 218 ;; Sf Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 219 ;; Sf ;undefine, Isnow_1 220 ;; Sf ; 221 ;; Sf IF std_file2_I NE std_file1_I THEN BEGIN 222 ;; Sf Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 223 ;; Sf caldat, time, mm 224 ;; Sf april = where(mm EQ 4, cnt) 225 ;; Sf Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 226 ;; Sf sept = where(mm EQ 9, cnt) 227 ;; Sf Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 228 ;; Sf ; undefine, Isnow_2 229 ;; Sf ENDIF ELSE BEGIN 230 ;; Sf Isnow_april_2 = {arr:-1} 231 ;; Sf Isnow_sept_2 = {arr:-1} 232 ;; Sf ENDELSE 211 Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 212 help, Isnow_1.arr, jpi,jpj 213 caldat, time, mm 214 april = where(mm EQ 4, cnt) 215 Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 216 sept = where(mm EQ 9, cnt) 217 Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 218 ;undefine, Isnow_1 219 ; 220 IF std_file2_I NE std_file1_I THEN BEGIN 221 Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 222 caldat, time, mm 223 april = where(mm EQ 4, cnt) 224 Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 225 sept = where(mm EQ 9, cnt) 226 Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 227 ; undefine, Isnow_2 228 ENDIF ELSE BEGIN 229 Isnow_april_2 = {arr:-1} 230 Isnow_sept_2 = {arr:-1} 231 ENDELSE 233 232 ; 234 233 jpt = 1 … … 365 364 cnt = cnt+1 & blabla = 'Antarctic Ice Fraction: SEPT' 366 365 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 367 ;; SF need to be reviewed, it does not work, it need a structure!!!! 368 ;; SF 369 ;; SF cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: APRIL' 370 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex 371 ;; SF 372 ;; SF cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: SEPT' 373 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 374 ;; SF 375 ;; SF cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: APRIL' 376 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex 377 ;; SF 378 ;; SF cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT' 379 ;; SF IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 366 ; 367 cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: APRIL' 368 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex 369 ; 370 cnt = cnt+1 & blabla = 'Arctic SNOW Thickness: SEPT' 371 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 372 ; 373 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: APRIL' 374 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex 375 ; 376 cnt = cnt+1 & blabla = 'Antarctic SNOW Thickness: SEPT' 377 IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 380 378 ; 381 379 ; NOT NEEDED, only for TIME SERIES NOW ICE EXTEND
Note: See TracChangeset
for help on using the changeset viewer.