Changeset 4501
- Timestamp:
- 2014-02-20T14:21:03+01:00 (10 years ago)
- Location:
- branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
- Files:
-
- 1 deleted
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_SnowThick.pro
r4320 r4501 22 22 domdef, 20, 380, 50, 90 23 23 vmin = 0. 24 vmax = 2.24 vmax = 1. 25 25 vint = 0.1 26 26 ;fmt = '(i2)' … … 36 36 domdef, 20, 380, -90, -50 37 37 vmin = 0. 38 vmax = 2.38 vmax = 1. 39 39 vint = 0.1 40 40 fmt = '(f4.1)' -
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 -
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro
r4394 r4501 72 72 IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_AMOC, masknp, POSTSCRIPT = postscript, _extra = ex 73 73 ; 74 cnt = cnt+1 & blabla = 'Sea-Ice cover'74 cnt = cnt+1 & blabla = 'Sea-Ice Cover' 75 75 IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE, masknp, POSTSCRIPT = postscript, _extra = ex 76 76 ; … … 78 78 IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vol, masknp, POSTSCRIPT = postscript, _extra = ex 79 79 80 ;SF cnt = cnt+1 & blabla = ' SNOW Volume'81 ;SF IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ SNOW_Vol, masknp, POSTSCRIPT = postscript, _extra = ex80 ;SF cnt = cnt+1 & blabla = 'ICE Velocity ' 81 ;SF IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex 82 82 ; 83 83 cnt = cnt+1 & blabla = 'Sea-Ice Volume Export at Fram Strait'
Note: See TracChangeset
for help on using the changeset viewer.