Changeset 4501


Ignore:
Timestamp:
2014-02-20T14:21:03+01:00 (7 years ago)
Author:
flavoni
Message:

add snow thickness plots, see ticket #724

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  
    2222    domdef, 20, 380, 50, 90 
    2323    vmin = 0. 
    24     vmax = 2. 
     24    vmax = 1. 
    2525    vint = 0.1 
    2626    ;fmt = '(i2)' 
     
    3636    domdef, 20, 380, -90, -50 
    3737    vmin = 0. 
    38     vmax = 2. 
     38    vmax = 1. 
    3939    vint = 0.1 
    4040    fmt = '(f4.1)' 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro

    r4395 r4501  
    209209  ENDELSE 
    210210; 
    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 
    233232; 
    234233  jpt = 1 
     
    365364  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: SEPT' 
    366365  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 
    380378;  
    381379; 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  
    7272  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_AMOC, masknp, POSTSCRIPT = postscript, _extra = ex 
    7373; 
    74   cnt = cnt+1   &   blabla = 'Sea-Ice cover' 
     74  cnt = cnt+1   &   blabla = 'Sea-Ice Cover' 
    7575  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE, masknp, POSTSCRIPT = postscript, _extra = ex 
    7676; 
     
    7878  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vol, masknp, POSTSCRIPT = postscript, _extra = ex 
    7979 
    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 = ex 
     80;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 
    8282; 
    8383  cnt = cnt+1   &   blabla = 'Sea-Ice Volume Export at Fram Strait' 
Note: See TracChangeset for help on using the changeset viewer.