New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 4514 for branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM – NEMO

Ignore:
Timestamp:
2014-02-25T12:41:51+01:00 (10 years ago)
Author:
flavoni
Message:

add ice pots & time series, see ticket #724

Location:
branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts
Files:
3 added
5 deleted
12 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceSal.pro

    r4507 r4514  
    1 pro std_plot_IceSal, Isal1, Isal2, ARC = arc, ANT = ant, APRIL = april, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
     1pro std_plot_IceSal, Isal1, Isal2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
    22 
    33  compile_opt idl2, strictarrsubs 
     
    1212  IF keyword_set(arc) THEN var = var+'_Arc_' 
    1313  IF keyword_set(ant) THEN var = var+'_Ant_' 
    14   IF keyword_set(april) THEN var = var+'April' 
     14  IF keyword_set(febr) THEN var = var+'Febr' 
     15  IF keyword_set(march) THEN var = var+'March' 
    1516  IF keyword_set(sept) THEN var = var+'Sept' 
    1617 
     
    3132  ; 
    3233   title = var+'!C'+std_file1_I 
    33    plt, Isal1.arr, MIN = .1, MAX = 20.1, INTER = 1., /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $ 
     34   plt, Isal1.arr, MIN = 0., MAX = 12., INTER = 1., /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $ 
    3435        , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $ 
    3536        , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex 
     
    3738   if std_file1_I NE std_file2_I then begin 
    3839     title = title + std_file2_I 
    39      plt, Isal1.arr - Isal2.arr, MIN = -4, MAX = 4, INTER = .2, STYLE = 'so0so', format = '(f4.1)' $ 
     40     plt, Isal1.arr - Isal2.arr, MIN = -4, MAX = 4, INTER = 1., STYLE = 'so0so', format = '(f4.1)' $ 
    4041          , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $ 
    4142          , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_SnowThick.pro

    r4501 r4514  
    1 pro std_plot_SnowThick, Isnow1, Isnow2, ARC = arc, ANT = ant, APRIL = april, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
     1pro std_plot_SnowThick, Isnow1, Isnow2, ARC = arc, ANT = ant, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex 
    22 
    33  compile_opt idl2, strictarrsubs 
     
    1212  IF keyword_set(arc) THEN var = var+'_Arc_' 
    1313  IF keyword_set(ant) THEN var = var+'_Ant_' 
    14   IF keyword_set(april) THEN var = var+'April' 
     14  IF keyword_set(march) THEN var = var+'March' 
    1515  IF keyword_set(sept) THEN var = var+'Sept' 
    1616 
     
    2222    domdef, 20, 380, 50, 90 
    2323    vmin = 0. 
    24     vmax = 1. 
     24    vmax = .6 
    2525    vint = 0.1 
    2626    ;fmt = '(i2)' 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro

    r4507 r4514  
    187187  ENDELSE 
    188188; 
     189  Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I )  
     190  caldat, time, mm 
     191  march = where(mm EQ 3, cnt) 
     192  Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
     193  sept = where(mm EQ 9, cnt) 
     194  Iage_sept_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_1.unit}  
     195  undefine, Iage_1 
     196; 
     197  IF std_file2_I NE std_file1_I THEN BEGIN 
     198     Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I )  
     199     caldat, time, mm 
     200     march = where(mm EQ 3, cnt) 
     201     Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
     202     sept = where(mm EQ 9, cnt) 
     203     Iage_sept_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
     204     undefine, Iage_2 
     205  ENDIF ELSE BEGIN  
     206     Iage_march_2 = {arr:-1} 
     207     Iage_sept_2 = {arr:-1} 
     208  ENDELSE 
     209; 
    189210  Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )  
    190   ;help, Ifra_1.arr, jpi,jpj 
    191211  caldat, time, mm 
    192212  march = where(mm EQ 3, cnt) 
     
    210230; 
    211231   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I )  
    212    ;help, Isnow_1.arr, jpi,jpj 
    213232   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}  
     233   march = where(mm EQ 3, cnt) 
     234   Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}  
    216235   sept = where(mm EQ 9, cnt) 
    217236   Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 
     
    221240     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I )  
    222241     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} 
     242     march = where(mm EQ 3, cnt) 
     243     Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
    225244     sept = where(mm EQ 9, cnt) 
    226245     Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
    227246   ;  undefine, Isnow_2 
    228247   ENDIF ELSE BEGIN  
    229      Isnow_april_2 = {arr:-1} 
     248     Isnow_march_2 = {arr:-1} 
    230249     Isnow_sept_2 = {arr:-1} 
    231250   ENDELSE 
    232251; 
    233252   Isal_1 = read_ncdf(getenv('VAR1_Isal'), allrecords = allrec, filename = std_file1_I )  
    234    ;help, Isal_1.arr, jpi,jpj 
     253   Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 
     254   msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations  
     255   help, Ifra_1 
     256   help, mask 
    235257   caldat, time, mm 
    236    april = where(mm EQ 4, cnt) 
    237    Isal_april_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
     258   march = where(mm EQ 3, cnt) 
     259   febr = where(mm EQ 2, cnt) 
     260   ;SF  Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
     261   Isal_1.arr = Isal_1.arr * msk 
     262   help, Isal_1 
     263   Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
     264   help, Isal_march_1 
     265   Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
    238266   sept = where(mm EQ 9, cnt) 
    239267   Isal_sept_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 
    240268   ;undefine, Isal_1 
    241269 ; 
     270   Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I )  
     271   Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I ) 
     272   msk = Ifra_2.arr gt 0.15 ; remove 0.15% for observations  
    242273   IF std_file2_I NE std_file1_I THEN BEGIN 
    243274     Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I )  
     275     Isal_2.arr = Isal_2.arr * msk 
    244276     caldat, time, mm 
    245      april = where(mm EQ 4, cnt) 
    246      Isal_april_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
     277     march = where(mm EQ 3, cnt) 
     278     febr = where(mm EQ 2, cnt) 
     279     Isal_march_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
     280     Isal_febr_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
    247281     sept = where(mm EQ 9, cnt) 
    248282     Isal_sept_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
    249283   ;  undefine, Isal_2 
    250284   ENDIF ELSE BEGIN  
    251      Isal_april_2 = {arr:-1} 
     285     Isal_febr_2 = {arr:-1} 
     286     Isal_march_2 = {arr:-1} 
    252287     Isal_sept_2 = {arr:-1} 
    253288   ENDELSE 
     
    375410  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ANT, /SEPT, _extra = ex 
    376411; 
     412  cnt = cnt+1   &   blabla = 'Arctic Ice Age: MARCH' 
     413  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ARC, /MARCH, _extra = ex 
     414 
     415  cnt = cnt+1   &   blabla = 'Arctic Ice Age: SEPT' 
     416  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ARC, /SEPT, _extra = ex 
     417 
     418  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: MARCH' 
     419  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ANT, /MARCH, _extra = ex 
     420; 
     421  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: SEPT' 
     422  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ANT, /SEPT, _extra = ex 
     423; 
    377424  cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: MARCH' 
    378425  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 
     
    387434  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 
    388435; 
    389   cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: APRIL' 
    390   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex 
     436  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: MARCH' 
     437  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ARC, /MARCH, _extra = ex 
    391438;  
    392439  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT' 
    393440  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 
    394441 
    395   cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: APRIL' 
    396   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex 
     442  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: MARCH' 
     443  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex 
    397444;  
    398445  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT' 
    399446  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 
    400447;  
    401   cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: APRIL' 
    402   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_april_1, Isal_april_2, /ARC, /APRIL, _extra = ex 
     448  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: MARCH' 
     449  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_march_1, Isal_march_2, /ARC, /MARCH, _extra = ex 
    403450;  
    404451  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: SEPT' 
    405452  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ARC, /SEPT, _extra = ex 
    406453 
    407   cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: APRIL' 
    408   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_april_1, Isal_april_2, /ANT, /APRIL, _extra = ex 
     454  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: FEBRUARY' 
     455  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_febr_1, Isal_febr_2, /ANT, /FEBR, _extra = ex 
    409456;  
    410457  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: SEPT' 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_1Y.pro

    r4320 r4514  
    115115  EMP1 = read_ncdf(getenv('VAR1_EMP'), allrecords = allrec, direc = 't', filename = std_file1_T ) 
    116116  EMP1 = {arr:EMP1.arr * 86400., unit:'mm/day', grid:'T'} 
    117   IF std_file2_T NE std_file1_T THEN BEGIN 
    118     EMP2 = read_ncdf(getenv('VAR2_EMP'), allrecords = allrec, direc = 't', filename = std_file2_T ) 
    119     EMP2 = {arr:EMP2.arr * 86400., unit:'mm/day', grid:'T'} 
    120   ENDIF ELSE EMP2 = {arr:-1} 
    121    
    122 ;mixed layer depth 
    123   MLD1 = read_ncdf(getenv('VAR1_MLD'), allrecords = allrec, direc = 't', filename = std_file1_T ) ; 10 m 
    124   IF std_file2_T NE std_file1_T THEN BEGIN 
    125     MLD2 = read_ncdf(getenv('VAR2_MLD'), allrecords = allrec, direc = 't', filename = std_file2_T ) ; 10 m 
    126   ENDIF ELSE MLD2 = {arr:-1} 
    127 ;climatology  
    128   MLD = read_ncdf(getenv('VAR_MLD'), filename = std_file_mld ) 
    129  
    130 ; velocities 
    131   U1 = read_ncdf(getenv('VAR1_U'), allrecords = allrec, direc = 't', filename = std_file1_U ) 
    132   IF strlowcase(getenv('VAR1_U')) EQ 'uocetr_eff' THEN BEGIN 
    133     U1.arr = U1.arr / e3u_3d(/e2) * umask() 
    134     U1.unit = 'm/s' 
    135   ENDIF  
    136   IF std_file2_U NE std_file1_U THEN BEGIN 
    137     U2 = read_ncdf(getenv('VAR2_U'), allrecords = allrec, direc = 't', filename = std_file2_U ) 
    138     IF strlowcase(getenv('VAR2_U')) EQ 'uocetr_eff' THEN BEGIN 
    139       U2.arr = U2.arr / e3u_3d(/e2) * umask() 
    140       U2.unit = 'm/s' 
    141     ENDIF  
    142   ENDIF ELSE U2 = {arr:-1} 
    143 ; 
    144   V1 = read_ncdf(getenv('VAR1_V'), allrecords = allrec, direc = 't', filename = std_file1_V ) 
    145   IF strlowcase(getenv('VAR1_V')) EQ 'vocetr_eff' THEN BEGIN 
    146     V1.arr = V1.arr / e3v_3d(/e1) * vmask() 
    147     V1.unit = 'm/s' 
    148   ENDIF  
    149   IF std_file2_V NE std_file1_V THEN BEGIN 
    150     V2 = read_ncdf(getenv('VAR2_V'), allrecords = allrec, direc = 't', filename = std_file2_V ) 
    151     IF strlowcase(getenv('VAR2_V')) EQ 'vocetr_eff' THEN BEGIN 
    152       V2.arr = V2.arr / e3v_3d(/e1) * vmask() 
    153       V2.unit = 'm/s' 
    154     ENDIF  
    155   ENDIF ELSE V2 = {arr:-1} 
     117   IF std_file2_T NE std_file1_T THEN BEGIN 
     118     EMP2 = read_ncdf(getenv('VAR2_EMP'), allrecords = allrec, direc = 't', filename = std_file2_T ) 
     119     EMP2 = {arr:EMP2.arr * 86400., unit:'mm/day', grid:'T'} 
     120   ENDIF ELSE EMP2 = {arr:-1} 
     121    
     122 ;mixed layer depth 
     123   MLD1 = read_ncdf(getenv('VAR1_MLD'), allrecords = allrec, direc = 't', filename = std_file1_T ) ; 10 m 
     124   IF std_file2_T NE std_file1_T THEN BEGIN 
     125     MLD2 = read_ncdf(getenv('VAR2_MLD'), allrecords = allrec, direc = 't', filename = std_file2_T ) ; 10 m 
     126   ENDIF ELSE MLD2 = {arr:-1} 
     127 ;climatology  
     128   MLD = read_ncdf(getenv('VAR_MLD'), filename = std_file_mld ) 
     129 
     130 ; velocities 
     131   U1 = read_ncdf(getenv('VAR1_U'), allrecords = allrec, direc = 't', filename = std_file1_U ) 
     132   ; old formulation: we tested variable name 
     133   ; IF strlowcase(getenv('VAR1_U')) EQ 'uocetr_eff' OR strlowcase(getenv('VAR1_U')) EQ 'vozoeftr' THEN BEGIN 
     134   IF strlowcase(U1.unit) EQ 'm3/s' THEN BEGIN 
     135     ;IF it is a transport it is transofrmed in velocity 
     136     U1.arr = U1.arr / e3u_3d(/e2) * umask() 
     137     U1.unit = 'm/s' 
     138   ENDIF  
     139   IF std_file2_U NE std_file1_U THEN BEGIN 
     140     U2 = read_ncdf(getenv('VAR2_U'), allrecords = allrec, direc = 't', filename = std_file2_U ) 
     141     ; old formulation: we tested variable name 
     142     ; IF strlowcase(getenv('VAR2_U')) EQ 'uocetr_eff' OR strlowcase(getenv('VAR2_U')) EQ 'vozoeftr' THEN BEGIN 
     143     IF strlowcase(U2.unit) EQ 'm3/s' THEN BEGIN 
     144       U2.arr = U2.arr / e3u_3d(/e2) * umask() 
     145       U2.unit = 'm/s' 
     146     ENDIF  
     147   ENDIF ELSE U2 = {arr:-1} 
     148; 
     149   V1 = read_ncdf(getenv('VAR1_V'), allrecords = allrec, direc = 't', filename = std_file1_V ) 
     150   ; old formulation: we tested variable name 
     151   ; IF strlowcase(getenv('VAR1_V')) EQ 'vocetr_eff' OR strlowcase(getenv('VAR1_V')) EQ 'vomeeftr' THEN BEGIN 
     152   IF strlowcase(V1.unit) EQ 'm3/s' THEN BEGIN ; test on unit to understand if it is a transport or velocity 
     153     ;IF it is a transport it is transofrmed in velocity 
     154     V1.arr = V1.arr / e3v_3d(/e1) * vmask() 
     155     V1.unit = 'm/s' 
     156   ENDIF  
     157   IF std_file2_V NE std_file1_V THEN BEGIN 
     158     V2 = read_ncdf(getenv('VAR2_V'), allrecords = allrec, direc = 't', filename = std_file2_V ) 
     159     ; old formulation 
     160     ; IF strlowcase(getenv('VAR2_V')) EQ 'vocetr_eff' OR strlowcase(getenv('VAR2_V')) EQ 'vozoeftr' THEN BEGIN 
     161     IF strlowcase(V2.unit) EQ 'm3/s' THEN BEGIN  
     162       V2.arr = V2.arr / e3v_3d(/e1) * vmask() 
     163       V2.unit = 'm/s' 
     164     ENDIF  
     165   ENDIF ELSE V2 = {arr:-1} 
    156166 
    157167; ice 
     
    177187  ENDELSE 
    178188; 
     189  Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I )  
     190  caldat, time, mm 
     191  march = where(mm EQ 3, cnt) 
     192  Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
     193  sept = where(mm EQ 9, cnt) 
     194  Iage_sept_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_1.unit}  
     195  undefine, Iage_1 
     196; 
     197  IF std_file2_I NE std_file1_I THEN BEGIN 
     198     Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I )  
     199     caldat, time, mm 
     200     march = where(mm EQ 3, cnt) 
     201     Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
     202     sept = where(mm EQ 9, cnt) 
     203     Iage_sept_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
     204     undefine, Iage_2 
     205  ENDIF ELSE BEGIN  
     206     Iage_march_2 = {arr:-1} 
     207     Iage_sept_2 = {arr:-1} 
     208  ENDELSE 
     209; 
    179210  Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )  
    180   help, Ifra_1.arr, jpi,jpj 
    181211  caldat, time, mm 
    182212  march = where(mm EQ 3, cnt) 
     
    199229  ENDELSE 
    200230; 
    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 
     231   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I )  
     232   caldat, time, mm 
     233   march = where(mm EQ 3, cnt) 
     234   Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}  
     235   sept = where(mm EQ 9, cnt) 
     236   Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 
     237   ;undefine, Isnow_1 
     238 ; 
     239   IF std_file2_I NE std_file1_I THEN BEGIN 
     240     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I )  
     241     caldat, time, mm 
     242     march = where(mm EQ 3, cnt) 
     243     Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
     244     sept = where(mm EQ 9, cnt) 
     245     Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
     246   ;  undefine, Isnow_2 
     247   ENDIF ELSE BEGIN  
     248     Isnow_march_2 = {arr:-1} 
     249     Isnow_sept_2 = {arr:-1} 
     250   ENDELSE 
     251; 
     252   Isal_1 = read_ncdf(getenv('VAR1_Isal'), allrecords = allrec, filename = std_file1_I )  
     253   Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 
     254   msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations  
     255   help, Ifra_1 
     256   help, mask 
     257   caldat, time, mm 
     258   march = where(mm EQ 3, cnt) 
     259   febr = where(mm EQ 2, cnt) 
     260   ;SF  Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
     261   Isal_1.arr = Isal_1.arr * msk 
     262   help, Isal_1 
     263   Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
     264   help, Isal_march_1 
     265   Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
     266   sept = where(mm EQ 9, cnt) 
     267   Isal_sept_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 
     268   ;undefine, Isal_1 
     269 ; 
     270   Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I )  
     271   Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I ) 
     272   msk = Ifra_2.arr gt 0.15 ; remove 0.15% for observations  
     273   IF std_file2_I NE std_file1_I THEN BEGIN 
     274     Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I )  
     275     Isal_2.arr = Isal_2.arr * msk 
     276     caldat, time, mm 
     277     march = where(mm EQ 3, cnt) 
     278     febr = where(mm EQ 2, cnt) 
     279     Isal_march_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
     280     Isal_febr_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
     281     sept = where(mm EQ 9, cnt) 
     282     Isal_sept_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
     283   ;  undefine, Isal_2 
     284   ENDIF ELSE BEGIN  
     285     Isal_febr_2 = {arr:-1} 
     286     Isal_march_2 = {arr:-1} 
     287     Isal_sept_2 = {arr:-1} 
     288   ENDELSE 
    223289; 
    224290  jpt = 1 
     
    256322  cnt = cnt+1   &   blabla = 'Erp salinity damping term' 
    257323  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_erp, ERP1, ERP2, _extra = ex 
    258  
     324; 
    259325  cnt = cnt+1   &   blabla = 'Evaporation - Precipitation - Runoff term' 
    260326  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_emp, EMP1, EMP2, _extra = ex 
    261  
     327; 
    262328  cnt = cnt+1   &   blabla = 'Net heat flux' 
    263329  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_qnet, Q1, Q2, QNET, _extra = ex 
    264  
     330; 
    265331  cnt = cnt+1   &   blabla = 'Meridionnal Heat Transport' 
    266332  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mht, Q1.arr+geo, Q2.arr+geo, masknp, std_file_msksub, _extra = ex 
    267  
     333; 
    268334  cnt = cnt+1   &   blabla = 'Global Barotropic stream Function' 
    269335  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_bsf, U1, U2, _extra = ex 
    270  
     336; 
    271337  cnt = cnt+1   &   blabla = 'mean Temperature diff with New Reynolds' 
    272338  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sst, T1, T2, TRey, _extra = ex 
    273  
     339; 
    274340  cnt = cnt+1   &   blabla = 'mean Salinity diff with Levitus' 
    275341  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sss, S1, S2, SLev, _extra = ex 
    276  
     342; 
    277343  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus' 
    278344  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, _extra = ex        
    279  
     345; 
    280346  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus and exp2' 
    281347  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, _extra = ex 
    282  
     348; 
    283349  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus at z=100 meters' 
    284350  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, /z100, _extra = ex        
    285   
     351;  
    286352  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus and exp2 at z=100 meters' 
    287353  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, /z100, _extra = ex 
    288  
     354; 
    289355  cnt = cnt+1   &   blabla = 'mean Temperature diff with Levitus at z=100 meters' 
    290356  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_T100m, T1, T2, Tlev, _extra = ex 
    291    
     357;   
    292358  cnt = cnt+1   &   blabla = 'mean Salinity diff with Levitus at z=100 meters' 
    293359  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_S100m, S1, S2, SLev, _extra = ex 
    294    
     360;   
    295361  cnt = cnt+1   &   blabla = 'Mixed layer depth' 
    296362  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD, _extra = ex 
    297  
     363; 
    298364  cnt = cnt+1   &   blabla = 'Mixed layer depth differences' 
    299365  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD2, MLD, _extra = ex 
    300  
     366; 
    301367  cnt = cnt+1   &   blabla = 'Zonal mean Mixed layer depth' 
    302368  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ZonMld, MLD1, MLD2, MLD, _extra = ex 
    303    
     369;   
    304370  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Global' 
    305371  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, _extra = ex 
    306  
     372; 
    307373  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Atlantic' 
    308374  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Atl', _extra = ex 
    309    
     375;   
    310376  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Indian' 
    311377  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Ind', _extra = ex 
    312    
     378;   
    313379  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Pacific' 
    314380  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Pac', _extra = ex 
    315  
     381; 
    316382  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Global' 
    317383  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, _extra = ex 
    318  
     384; 
    319385  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Atlantic' 
    320386  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Atl', _extra = ex 
    321  
     387; 
    322388  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Indian' 
    323389  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Ind', _extra = ex  
    324  
     390; 
    325391  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Pacific' 
    326392  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Pac', _extra = ex  
    327    
     393;   
    328394;  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: MARCH' 
    329395;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 
    330  
     396; 
    331397  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: APRIL' 
    332398  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 
    333  
     399; 
    334400  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: SEPT' 
    335401  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 
    336    
     402;   
    337403;  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: MARCH' 
    338404;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 
    339  
     405;  
    340406  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: APRIL' 
    341407  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex 
    342  
     408; 
    343409  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: SEPT' 
    344410  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ANT, /SEPT, _extra = ex 
    345  
     411; 
     412  cnt = cnt+1   &   blabla = 'Arctic Ice Age: MARCH' 
     413  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ARC, /MARCH, _extra = ex 
     414 
     415  cnt = cnt+1   &   blabla = 'Arctic Ice Age: SEPT' 
     416  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ARC, /SEPT, _extra = ex 
     417 
     418  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: MARCH' 
     419  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ANT, /MARCH, _extra = ex 
     420; 
     421  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: SEPT' 
     422  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ANT, /SEPT, _extra = ex 
     423; 
    346424  cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: MARCH' 
    347425  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 
    348    
     426;   
    349427  cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: SEPT' 
    350428  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ARC, /SEPT, _extra = ex 
    351    
     429;   
    352430  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: MARCH' 
    353431  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex 
    354  
     432; 
    355433  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: SEPT' 
    356434  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 
    375  
     435; 
     436  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: MARCH' 
     437  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ARC, /MARCH, _extra = ex 
     438;  
     439  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT' 
     440  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 
     441 
     442  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: MARCH' 
     443  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex 
     444;  
     445  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT' 
     446  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 
     447;  
     448  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: MARCH' 
     449  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_march_1, Isal_march_2, /ARC, /MARCH, _extra = ex 
     450;  
     451  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: SEPT' 
     452  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ARC, /SEPT, _extra = ex 
     453 
     454  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: FEBRUARY' 
     455  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_febr_1, Isal_febr_2, /ANT, /FEBR, _extra = ex 
     456;  
     457  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: SEPT' 
     458  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ANT, /SEPT, _extra = ex 
     459;  
    376460  cnt = cnt+1   &   blabla = 'Meridional stream Function: Global (no Med)' 
    377461  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'GloNoMed', _extra = ex 
    378  
     462 
    379463  cnt = cnt+1   &   blabla = 'Meridional stream Function: Atlantic' 
    380464  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Atl', _extra = ex 
    381  
     465 ; 
    382466  cnt = cnt+1   &   blabla = 'Meridional stream Function: Indian' 
    383467  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Ind', _extra = ex 
    384  
     468 ; 
    385469  cnt = cnt+1   &   blabla = 'Meridional stream Function: Indo-Pacific' 
    386470  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'IndoPac', _extra = ex 
    387  
     471; 
    388472  cnt = cnt+1   &   blabla = 'Equatorial Temperature' 
    389473  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqT, T1, T2, Tlev, _extra = ex 
    390  
     474; 
    391475  cnt = cnt+1   &   blabla = 'Equatorial Salinity' 
    392476  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqS, S1, S2, SLev, _extra = ex 
    393  
     477; 
    394478  cnt = cnt+1   &   blabla = 'Equatorial zonal velocity' 
    395479  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqU, U1, U2, _extra = ex 
    396  
     480; 
    397481  cnt = cnt+1   &   blabla = 'Mediterranean salt tongue at depth=700' 
    398482  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 700, _extra = ex 
    399  
     483; 
    400484  cnt = cnt+1   &   blabla = 'Mediterranean salt tongue at depth=1000' 
    401485  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 1000, _extra = ex 
    402  
     486; 
    403487  cnt = cnt+1   &   blabla = 'Mediterranean water at lat=40°N' 
    404488  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 40, _extra = ex 
    405  
     489; 
    406490  cnt = cnt+1   &   blabla = 'Mediterranean water at lat=38°N' 
    407491  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex 
    408  
     492; 
    409493  cnt = cnt+1   &   blabla = 'Vertical Global mean T & S' 
    410494  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_GlobMeanTS, T1, T2, TLev, S1, S2, SLev, _extra = ex 
    411  
     495; 
    412496  IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_plot_html_body.txt', htmltxt[1:*] 
    413497   
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_last_year_5D.pro

    r4320 r4514  
    203203  ENDELSE 
    204204; 
     205  Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I )  
     206  caldat, time, mm 
     207  march = where(mm EQ 3, cnt) 
     208  Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
     209  sept = where(mm EQ 9, cnt) 
     210  Iage_sept_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_1.unit}  
     211  undefine, Iage_1 
     212; 
     213  IF std_file2_I NE std_file1_I THEN BEGIN 
     214     Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I )  
     215     caldat, time, mm 
     216     march = where(mm EQ 3, cnt) 
     217     Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
     218     sept = where(mm EQ 9, cnt) 
     219     Iage_sept_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_2.unit}  
     220     undefine, Iage_2 
     221  ENDIF ELSE BEGIN  
     222     Iage_march_2 = {arr:-1} 
     223     Iage_sept_2 = {arr:-1} 
     224  ENDELSE 
     225; 
    205226  Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )  
    206227  ;; help, Ifra_1.arr, jpi,jpj 
     
    225246  ENDELSE 
    226247; 
    227 ;; SF commentato perché non funziona snow 
    228 ;; SF   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I )  
    229 ;; SF   help, Isnow_1.arr, jpi,jpj 
    230 ;; SF   caldat, time, mm 
    231 ;; SF   april = where(mm EQ 4, cnt) 
    232 ;; SF   Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}  
    233 ;; SF   sept = where(mm EQ 9, cnt) 
    234 ;; SF   Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 
    235 ;; SF   ;undefine, Isnow_1 
    236 ;; SF ; 
    237 ;; SF   IF std_file2_I NE std_file1_I THEN BEGIN 
    238 ;; SF     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I )  
    239 ;; SF     caldat, time, mm 
    240 ;; SF     april = where(mm EQ 4, cnt) 
    241 ;; SF     Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
    242 ;; SF     sept = where(mm EQ 9, cnt) 
    243 ;; SF     Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
    244 ;; SF   ;  undefine, Isnow_2 
    245 ;; SF   ENDIF ELSE BEGIN  
    246 ;; SF     Isnow_april_2 = {arr:-1} 
    247 ;; SF     Isnow_sept_2 = {arr:-1} 
    248 ;; SF   ENDELSE 
     248   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I )  
     249   caldat, time, mm 
     250   march = where(mm EQ 3, cnt) 
     251   Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}  
     252   sept = where(mm EQ 9, cnt) 
     253   Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 
     254   ;undefine, Isnow_1 
     255 ; 
     256   IF std_file2_I NE std_file1_I THEN BEGIN 
     257     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I )  
     258     caldat, time, mm 
     259     march = where(mm EQ 3, cnt) 
     260     Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
     261     sept = where(mm EQ 9, cnt) 
     262     Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit} 
     263   ;  undefine, Isnow_2 
     264   ENDIF ELSE BEGIN  
     265     Isnow_march_2 = {arr:-1} 
     266     Isnow_sept_2 = {arr:-1} 
     267   ENDELSE 
     268; 
     269   Isal_1 = read_ncdf(getenv('VAR1_Isal'), allrecords = allrec, filename = std_file1_I )  
     270   Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 
     271   msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations  
     272   help, Ifra_1 
     273   help, mask 
     274   caldat, time, mm 
     275   march = where(mm EQ 3, cnt) 
     276   febr = where(mm EQ 2, cnt) 
     277   ;SF  Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
     278   Isal_1.arr = Isal_1.arr * msk 
     279   help, Isal_1 
     280   Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
     281   help, Isal_march_1 
     282   Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit}  
     283   sept = where(mm EQ 9, cnt) 
     284   Isal_sept_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 
     285   ;undefine, Isal_1 
     286 ; 
     287   Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I )  
     288   Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I ) 
     289   msk = Ifra_2.arr gt 0.15 ; remove 0.15% for observations  
     290   IF std_file2_I NE std_file1_I THEN BEGIN 
     291     Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I )  
     292     Isal_2.arr = Isal_2.arr * msk 
     293     caldat, time, mm 
     294     march = where(mm EQ 3, cnt) 
     295     febr = where(mm EQ 2, cnt) 
     296     Isal_march_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
     297     Isal_febr_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
     298     sept = where(mm EQ 9, cnt) 
     299     Isal_sept_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_2.unit} 
     300   ;  undefine, Isal_2 
     301   ENDIF ELSE BEGIN  
     302     Isal_febr_2 = {arr:-1} 
     303     Isal_march_2 = {arr:-1} 
     304     Isal_sept_2 = {arr:-1} 
     305   ENDELSE 
    249306; 
    250307  jpt = 1 
     
    282339  cnt = cnt+1   &   blabla = 'Erp salinity damping term' 
    283340  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_erp, ERP1, ERP2, _extra = ex 
    284  
     341; 
    285342  cnt = cnt+1   &   blabla = 'Evaporation - Precipitation - Runoff term' 
    286343  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_emp, EMP1, EMP2, _extra = ex 
    287  
     344; 
    288345  cnt = cnt+1   &   blabla = 'Net heat flux' 
    289346  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_qnet, Q1, Q2, QNET, _extra = ex 
    290  
     347; 
    291348  cnt = cnt+1   &   blabla = 'Meridionnal Heat Transport' 
    292349  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mht, Q1.arr+geo, Q2.arr+geo, masknp, std_file_msksub, _extra = ex 
    293  
     350; 
    294351  cnt = cnt+1   &   blabla = 'Global Barotropic stream Function' 
    295352  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_bsf, U1, U2, _extra = ex 
    296  
     353; 
    297354  cnt = cnt+1   &   blabla = 'mean Temperature diff with New Reynolds' 
    298355  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sst, T1, T2, TRey, _extra = ex 
    299  
     356; 
    300357  cnt = cnt+1   &   blabla = 'mean Salinity diff with Levitus' 
    301358  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sss, S1, S2, SLev, _extra = ex 
    302  
     359; 
    303360  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus' 
    304361  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, _extra = ex        
    305  
     362; 
    306363  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus and exp2' 
    307364  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, _extra = ex 
    308  
     365; 
    309366  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus at z=100 meters' 
    310367  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, /z100, _extra = ex        
    311   
     368;  
    312369  cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus and exp2 at z=100 meters' 
    313370  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, /z100, _extra = ex 
    314  
     371; 
    315372  cnt = cnt+1   &   blabla = 'mean Temperature diff with Levitus at z=100 meters' 
    316373  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_T100m, T1, T2, Tlev, _extra = ex 
    317    
     374;   
    318375  cnt = cnt+1   &   blabla = 'mean Salinity diff with Levitus at z=100 meters' 
    319376  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_S100m, S1, S2, SLev, _extra = ex 
    320   
    321 ;;SF need monthly climatology to plot difference of MAX of MLD   
    322   cnt = cnt+1   &   blabla = 'Mixed layer depth & difference with DeBoyer' 
    323   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD_1, MLD, _extra = ex 
    324     
     377 
     378  cnt = cnt+1   &   blabla = 'Mixed layer depth' 
     379  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD, _extra = ex 
     380; 
    325381  cnt = cnt+1   &   blabla = 'Mixed layer depth differences' 
    326382  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD_1, MLD_2, MLD, _extra = ex 
     
    339395  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Global' 
    340396  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, _extra = ex 
    341  
     397; 
    342398  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Atlantic' 
    343399  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Atl', _extra = ex 
    344    
     400;   
    345401  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Indian' 
    346402  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Ind', _extra = ex 
    347    
     403;   
    348404  cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Pacific' 
    349405  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Pac', _extra = ex 
    350  
     406; 
    351407  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Global' 
    352408  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, _extra = ex 
    353  
     409; 
    354410  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Atlantic' 
    355411  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Atl', _extra = ex 
    356  
     412; 
    357413  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Indian' 
    358414  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Ind', _extra = ex  
    359  
     415; 
    360416  cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Pacific' 
    361417  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Pac', _extra = ex  
    362    
     418;   
    363419;  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: MARCH' 
    364420;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex 
    365  
     421; 
    366422  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: APRIL' 
    367423  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex 
    368  
     424; 
    369425  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: SEPT' 
    370426  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex 
    371    
     427;   
    372428;  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: MARCH' 
    373429;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex 
    374  
     430;  
    375431  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: APRIL' 
    376432  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex 
    377  
     433; 
    378434  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: SEPT' 
    379435  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ANT, /SEPT, _extra = ex 
    380  
     436; 
     437  cnt = cnt+1   &   blabla = 'Arctic Ice Age: MARCH' 
     438  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ARC, /MARCH, _extra = ex 
     439 
     440  cnt = cnt+1   &   blabla = 'Arctic Ice Age: SEPT' 
     441  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ARC, /SEPT, _extra = ex 
     442 
     443  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: MARCH' 
     444  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ANT, /MARCH, _extra = ex 
     445; 
     446  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: SEPT' 
     447  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ANT, /SEPT, _extra = ex 
     448; 
    381449  cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: MARCH' 
    382450  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex 
    383    
     451;   
    384452  cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: SEPT' 
    385453  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ARC, /SEPT, _extra = ex 
    386    
     454;   
    387455  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: MARCH' 
    388456  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex 
    389  
     457; 
    390458  cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: SEPT' 
    391459  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex 
    392  
    393 ;; SF non funziona, dice che deve essere una struttura e non e!!!! 
    394 ;;SF  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: APRIL' 
    395 ;;SF  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex 
    396 ;; SF  
    397 ;; SF   cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT' 
    398 ;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 
    399 ;; SF   
    400 ;; SF   cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: APRIL' 
    401 ;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex 
    402 ;; SF   
    403 ;; SF   cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT' 
    404 ;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 
    405 ;  
    406 ;  cnt = cnt+1   &   blabla = 'Arctic Ice Extend: ' 
    407 ;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC,  _extra = ex 
    408    
    409 ;  cnt = cnt+1   &   blabla = 'Antarctic Ice Extend: ' 
    410 ;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT,  _extra = ex 
    411  
     460; 
     461  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: MARCH' 
     462  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ARC, /MARCH, _extra = ex 
     463;  
     464  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT' 
     465  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex 
     466 
     467  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: MARCH' 
     468  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex 
     469;  
     470  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT' 
     471  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex 
     472;  
     473  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: MARCH' 
     474  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_march_1, Isal_march_2, /ARC, /MARCH, _extra = ex 
     475;  
     476  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: SEPT' 
     477  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ARC, /SEPT, _extra = ex 
     478 
     479  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: FEBRUARY' 
     480  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_febr_1, Isal_febr_2, /ANT, /FEBR, _extra = ex 
     481;  
     482  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: SEPT' 
     483  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ANT, /SEPT, _extra = ex 
     484;  
    412485  cnt = cnt+1   &   blabla = 'Meridional stream Function: Global (no Med)' 
    413486  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'GloNoMed', _extra = ex 
    414  
     487 
    415488  cnt = cnt+1   &   blabla = 'Meridional stream Function: Atlantic' 
    416489  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Atl', _extra = ex 
    417  
     490 ; 
    418491  cnt = cnt+1   &   blabla = 'Meridional stream Function: Indian' 
    419492  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Ind', _extra = ex 
    420  
     493 ; 
    421494  cnt = cnt+1   &   blabla = 'Meridional stream Function: Indo-Pacific' 
    422495  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'IndoPac', _extra = ex 
    423  
     496; 
    424497  cnt = cnt+1   &   blabla = 'Equatorial Temperature' 
    425498  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqT, T1, T2, Tlev, _extra = ex 
    426  
     499; 
    427500  cnt = cnt+1   &   blabla = 'Equatorial Salinity' 
    428501  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqS, S1, S2, SLev, _extra = ex 
    429  
     502; 
    430503  cnt = cnt+1   &   blabla = 'Equatorial zonal velocity' 
    431504  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqU, U1, U2, _extra = ex 
    432  
     505; 
    433506  cnt = cnt+1   &   blabla = 'Mediterranean salt tongue at depth=700' 
    434507  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 700, _extra = ex 
    435  
     508; 
    436509  cnt = cnt+1   &   blabla = 'Mediterranean salt tongue at depth=1000' 
    437510  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 1000, _extra = ex 
    438  
     511; 
    439512  cnt = cnt+1   &   blabla = 'Mediterranean water at lat=40°N' 
    440513  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 40, _extra = ex 
    441  
     514; 
    442515  cnt = cnt+1   &   blabla = 'Mediterranean water at lat=38°N' 
    443516  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex 
    444  
     517; 
    445518  cnt = cnt+1   &   blabla = 'Vertical Global mean T & S' 
    446519  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_GlobMeanTS, T1, T2, TLev, S1, S2, SLev, _extra = ex 
    447  
     520; 
    448521  IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_plot_html_body.txt', htmltxt[1:*] 
    449522   
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh

    r4507 r4514  
    6969VAR1_Ithick=sithic 
    7070VAR1_Ifrac=siconc 
    71 VAR1_Isnow=sndept 
     71VAR1_Isnow=snvolu 
    7272VAR1_Isal=sisali 
     73VAR1_Iage=siages 
    7374# 
    7475#===================== EXP2 ===================== 
     
    8889VAR2_Ithick=sithic 
    8990VAR2_Ifrac=siconc 
    90 VAR2_Isnow=sndept 
     91VAR2_Isnow=snvolu 
    9192VAR2_Isal=sisali 
     93VAR2_Iage=siages 
    9294# 
    9395######################### Export Variables ############################### 
     
    113115export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    114116export VAR1_U      VAR1_V 
    115 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 
     117export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 
    116118#===================== EXP2 ===================== 
    117119export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    118120export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    119121export VAR2_U      VAR2_V 
    120 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 
     122export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 
    121123# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.last_year_5D

    r4507 r4514  
    2626SAXO_DIR=/Users/sflod/SAXO_DIR 
    2727# 
    28 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format 
    29 DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
    30 DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
     28DIR_DATA=~/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format 
     29DIR_CLIMATO=~/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
     30DIR_MASK=~//idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
    3131# 
    3232#===================== Model GRID ===================== 
     
    5555#===================== EXP1 ===================== 
    5656# 
    57 #FILE1_T=l3fwb178_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file 
    58 #FILE1_U=l3fwb178_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file 
    59 #FILE1_V=l3fwb178_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file 
    60 #FILE1_I=l3fwb178_24910101_25001231_1M_icemod.nc   # exp1 ice    input file 
    6157FILE1_T=l3fwb178_25010101_25011231_5D_grid_T.nc   # exp1 grid T input file 
    6258FILE1_U=l3fwb178_25010101_25011231_5D_grid_U.nc   # exp1 grid U input file 
    6359FILE1_V=l3fwb178_25010101_25011231_5D_grid_V.nc   # exp1 grid V input file 
    6460FILE1_I=l3fwb178_25010101_25011231_5D_icemod.nc   # exp1 ice    input file 
    65 VAR1_T=votemper          
    66 VAR1_S=vosaline          
     61VAR1_T=votemper 
     62VAR1_S=vosaline 
    6763VAR1_QNET=sohflxdo 
    68 VAR1_ERP=sowafldp       
    69 VAR1_EMP=sowaflup        
    70 #VAR1_MLD=somxl010 
     64VAR1_ERP=sowafldp 
     65VAR1_EMP=sowaflup 
    7166VAR1_MLD=max5dmld 
    72 #VAR1_U=uocetr_eff 
    73 #VAR1_U=uoce 
    74 #VAR1_U=vozocrtx 
    7567VAR1_U=vozoeftr 
    76 #VAR1_V=vocetr_eff 
    77 #VAR1_V=voce 
    7868VAR1_V=vomeeftr 
    79 #VAR1_V=vomecrty 
    80 VAR1_Ithick=iicethic     
     69VAR1_Ithick=iicethic 
    8170VAR1_Ifrac=iiceconc 
    82 VAR1_Isnow=isnowthi 
     71VAR1_Isnow=isnowvol 
    8372VAR1_Isal=iicesali 
     73VAR1_Iage=iiceages 
    8474# 
    8575#===================== EXP2 ===================== 
    8676# 
    87 #FILE2_T=l3fwb178_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file 
    88 #FILE2_U=l3fwb178_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file 
    89 #FILE2_V=l3fwb178_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file 
    90 #FILE2_I=l3fwb178_24910101_25001231_1M_icemod.nc   # exp1 ice    input file 
    9177FILE2_T=l3fwb178_25010101_25011231_5D_grid_T.nc   # exp1 grid T input file 
    9278FILE2_U=l3fwb178_25010101_25011231_5D_grid_U.nc   # exp1 grid U input file 
    9379FILE2_V=l3fwb178_25010101_25011231_5D_grid_V.nc   # exp1 grid V input file 
    9480FILE2_I=l3fwb178_25010101_25011231_5D_icemod.nc   # exp1 ice    input file 
    95 #FILE2_T=l378sal6_25010101_25011231_5D_grid_T.nc   # exp1 grid T input file 
    96 #FILE2_U=l378sal6_25010101_25011231_5D_grid_U.nc   # exp1 grid U input file 
    97 #FILE2_V=l378sal6_25010101_25011231_5D_grid_V.nc   # exp1 grid V input file 
    98 #FILE2_I=l378sal6_25010101_25011231_5D_icemod.nc   # exp1 ice    input file 
    9981VAR2_T=votemper 
    10082VAR2_S=vosaline 
     
    10284VAR2_ERP=sowafldp 
    10385VAR2_EMP=sowaflup 
    104 #VAR2_MLD=somxl010 
    10586VAR2_MLD=max5dmld 
    106 #VAR2_U=uoce_eff 
    107 #VAR2_U=uoce 
    10887VAR2_U=vozoeftr 
    109 #VAR2_V=vocetr_eff 
    110 #VAR2_V=voce 
    11188VAR2_V=vomeeftr 
    11289VAR2_Ithick=iicethic 
     
    11491VAR2_Isnow=isnowthi 
    11592VAR2_Isal=iicesali 
     93VAR2_Iage=iiceages 
    11694# 
    11795######################### Export Variables ############################### 
     
    137115export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    138116export VAR1_U      VAR1_V 
    139 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 
     117export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 
    140118#===================== EXP2 ===================== 
    141119export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    142120export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    143121export VAR2_U      VAR2_V 
    144 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 
     122export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 
    145123# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.new_names

    r4507 r4514  
    6969VAR1_Ithick=sithic 
    7070VAR1_Ifrac=siconc 
    71 VAR1_Isnow=sndept 
     71VAR1_Isnow=snvolu 
    7272VAR1_Isal=sisali 
     73VAR1_Iage=siages 
    7374# 
    7475#===================== EXP2 ===================== 
     
    8889VAR2_Ithick=sithic 
    8990VAR2_Ifrac=siconc 
    90 VAR2_Isnow=sndept 
     91VAR2_Isnow=snvolu 
    9192VAR2_Isal=sisali 
     93VAR2_Iage=siages 
    9294# 
    9395######################### Export Variables ############################### 
     
    113115export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    114116export VAR1_U      VAR1_V 
    115 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 
     117export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 
    116118#===================== EXP2 ===================== 
    117119export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    118120export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    119121export VAR2_U      VAR2_V 
    120 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 
     122export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 
    121123# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example1

    r4507 r4514  
    2626SAXO_DIR=/Users/sflod/SAXO_DIR 
    2727# 
    28 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada      # path of data in NetCDF format 
    29 DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
    30 DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
     28DIR_DATA=~/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format 
     29DIR_CLIMATO=~/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
     30DIR_MASK=~//idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
    3131# 
    3232#===================== Model GRID ===================== 
     
    5555#===================== EXP1 ===================== 
    5656# 
    57 FILE1_T=tr33beta_00910101_01001231_1Y_grid_T.nc   # exp1 grid T input file 
    58 FILE1_U=tr33beta_00910101_01001231_1Y_grid_U.nc   # exp1 grid U input file 
    59 FILE1_V=tr33beta_00910101_01001231_1Y_grid_V.nc   # exp1 grid V input file 
    60 FILE1_I=tr33beta_00910101_01001231_1M_icemod.nc   # exp1 ice    input file 
    61 VAR1_T=thetao 
    62 VAR1_S=so 
    63 VAR1_QNET=qt 
    64 VAR1_ERP=wfcorr 
    65 VAR1_EMP=wfo 
    66 VAR1_MLD=mldr10_1 
    67 VAR1_U=uocetr_eff 
    68 VAR1_V=vocetr_eff 
    69 VAR1_Ithick=sit 
    70 VAR1_Ifrac=sic 
    71 VAR1_Isnow=isnowthi 
     57FILE1_T=500yfwb0_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file 
     58FILE1_U=500yfwb0_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file 
     59FILE1_V=500yfwb0_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file 
     60FILE1_I=500yfwb0_24910101_25001231_1M_icemod.nc   # exp1 ice    input file 
     61VAR1_T=votemper 
     62VAR1_S=vosaline 
     63VAR1_QNET=sohflxdo 
     64VAR1_ERP=sowafldp 
     65VAR1_EMP=sowaflup 
     66VAR1_MLD=somxl010 
     67VAR1_U=vozoeftr 
     68VAR1_V=vomeeftr 
     69VAR1_Ithick=iicethic 
     70VAR1_Ifrac=iiceconc 
     71VAR1_Isnow=isnowvol 
    7272VAR1_Isal=iicesali 
     73VAR1_Iage=iiceages 
    7374# 
    7475#===================== EXP2 ===================== 
     
    9091VAR2_Isnow=isnowthi 
    9192VAR2_Isal=iicesali 
     93VAR2_Iage=iiceages 
    9294# 
    9395######################### Export Variables ############################### 
     
    113115export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    114116export VAR1_U      VAR1_V 
    115 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 
     117export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 
    116118#===================== EXP2 ===================== 
    117119export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    118120export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    119121export VAR2_U      VAR2_V 
    120 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 
     122export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 
    121123# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example2

    r4507 r4514  
    2121#idl_command=idl71 
    2222# 
    23 PS_DIR=$( pwd )/500yfwb0_ps_plot_500y 
    24 PDF_DIR=$( pwd )/500yfwb0_pdf_plot_500y 
     23PS_DIR=$( pwd )/ps_plot 
     24PDF_DIR=$( pwd )/pdf_plot 
    2525HTML_DIR=$( pwd )/html_plot 
    2626SAXO_DIR=/Users/sflod/SAXO_DIR 
    2727# 
    28 DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format 
    29 DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
    30 DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
     28DIR_DATA=~/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format 
     29DIR_CLIMATO=~/idl_PLOTS/CLIMATOLOGIES    # path of climatological data 
     30DIR_MASK=~//idl_PLOTS/MASK  # path of mask files (ex: subbasins) 
    3131# 
    3232#===================== Model GRID ===================== 
     
    5555#===================== EXP1 ===================== 
    5656# 
    57 FILE1_T=tag33_20910101_21001231_1Y_grid_T.nc   # exp1 grid T input file 
    58 FILE1_U=tag33_20910101_21001231_1Y_grid_U.nc   # exp1 grid U input file 
    59 FILE1_V=tag33_20910101_21001231_1Y_grid_V.nc   # exp1 grid V input file 
    60 FILE1_I=tag33_20910101_21001231_1M_icemod.nc   # exp1 ice    input file 
    61 VAR1_T=thetao 
    62 VAR1_S=so 
    63 VAR1_QNET=qt 
    64 VAR1_ERP=wfcorr 
    65 VAR1_EMP=wfo 
    66 VAR1_MLD=mldr10_1 
    67 VAR1_U=uocetr_eff 
    68 VAR1_V=vocetr_eff 
    69 VAR1_Ithick=sit 
    70 VAR1_Ifrac=sic 
    71 VAR1_Isnow=isnowthi 
     57FILE1_T=500yfwb0_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file 
     58FILE1_U=500yfwb0_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file 
     59FILE1_V=500yfwb0_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file 
     60FILE1_I=500yfwb0_24910101_25001231_1M_icemod.nc   # exp1 ice    input file 
     61VAR1_T=votemper 
     62VAR1_S=vosaline 
     63VAR1_QNET=sohflxdo 
     64VAR1_ERP=sowafldp 
     65VAR1_EMP=sowaflup 
     66VAR1_MLD=somxl010 
     67VAR1_U=vozoeftr 
     68VAR1_V=vomeeftr 
     69VAR1_Ithick=iicethic 
     70VAR1_Ifrac=iiceconc 
     71VAR1_Isnow=isnowvol 
    7272VAR1_Isal=iicesali 
     73VAR1_Iage=iiceages 
    7374# 
    7475#===================== EXP2 ===================== 
    7576# 
    76 FILE2_T=tmx_2000_00910101_01001231_1Y_grid_T.nc   # exp1 grid T input file 
    77 FILE2_U=tmx_2000_00910101_01001231_1Y_grid_U.nc   # exp1 grid U input file 
    78 FILE2_V=tmx_2000_00910101_01001231_1Y_grid_V.nc   # exp1 grid V input file 
    79 FILE2_I=tmx_2000_00910101_01001231_1M_icemod.nc   # exp1 ice    input file 
    80 VAR2_T=thetao 
    81 VAR2_S=so 
    82 VAR2_QNET=qns+qsr 
    83 VAR2_ERP=wfcorr 
    84 VAR2_EMP=wfo 
    85 VAR2_MLD=mldr10_1 
    86 VAR2_U=uoce_eff 
    87 VAR2_V=voce_eff 
    88 VAR2_Ithick=sit 
    89 VAR2_Ifrac=sic 
    90 VAR2_Isnow=isnowthi 
    91 VAR2_Isal=iicesali 
     77FILE1_T=500yfwb0_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file 
     78FILE1_U=500yfwb0_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file 
     79FILE1_V=500yfwb0_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file 
     80FILE1_I=500yfwb0_24910101_25001231_1M_icemod.nc   # exp1 ice    input file 
     81VAR1_T=votemper 
     82VAR1_S=vosaline 
     83VAR1_QNET=sohflxdo 
     84VAR1_ERP=sowafldp 
     85VAR1_EMP=sowaflup 
     86VAR1_MLD=somxl010 
     87VAR1_U=vozoeftr 
     88VAR1_V=vomeeftr 
     89VAR1_Ithick=iicethic 
     90VAR1_Ifrac=iiceconc 
     91VAR1_Isnow=isnowvol 
     92VAR1_Isal=iicesali 
     93VAR1_Iage=iiceages 
    9294# 
    9395######################### Export Variables ############################### 
     
    113115export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD  
    114116export VAR1_U      VAR1_V 
    115 export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal 
     117export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage 
    116118#===================== EXP2 ===================== 
    117119export FILE2_T     FILE2_U    FILE2_V   FILE2_I 
    118120export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD  
    119121export VAR2_U      VAR2_V 
    120 export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal 
     122export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage 
    121123# 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_FRAM.pro

    r4504 r4514  
    4949 ; 
    5050 ;OBSERVATIONS : mean seasonal cycle/month 
    51   vol_obs = [0.261625, 0.230750, 0.325375, 0.252000, 0.172500, 0.0805000, 0.0805000, 0.0805000, 0.0805000, 0.176500, 0.148500, 0.235000] 
     51 ; vol_obs = [0.261625, 0.230750, 0.325375, 0.252000, 0.172500, 0.0805000, 0.0805000, 0.0805000, 0.0805000, 0.176500, 0.148500, 0.235000] 
     52  vol_obs = [261.625, 230.750, 325.375, 252.000, 172.500, 80.5000, 80.5000, 80.5000, 80.5000, 176.500, 148.500, 235.000] 
    5253  area_obs = [0.103292, 0.0997500, 0.107625, 0.0944167, 0.0612083, 0.0262500, 0.0262500, 0.0262500, 0.0262500, 0.0843750, 0.0914583, 0.104083] 
    5354 ; 
     
    8081  ICE_vol_export = reform(ICE_vol_export, 12, nyr) 
    8182  ICE_vol_export = total(ICE_vol_export,2)/nyr 
    82   ; ICE_vol_export = {arr:ICE_vol_export * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'}      ; annual mean 
    83   ICE_vol_export = {arr:ICE_vol_export * 1.e-12 * 86400 * 30 , unit : '10^3 Km^3/month'}      ; monthly mean 
    84  
    85   ; 
    86   ;title = 'Northern Hemisphere'+'!C'+prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)' 
     83  ; ICE_vol_export = {arr:ICE_vol_export * 1.e-9 * 86400 * 365 , unit : '10^3 Km^3/year'}      ; annual mean 
     84  ICE_vol_export = {arr:ICE_vol_export * 1.e-9 * 86400 * 30 , unit : '10^3 Km^3/month'}      ; monthly mean 
     85 
     86  ; 
    8787  title = 'Fram Strait Areal Export'+'!C'+prefix+' '+d1_d2 
    8888  jpt=12 
     
    9595  tot_area_expo = total(ICE_area_export.arr) 
    9696  ; 
    97   xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 0.851 million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
    98   xyouts, julday(6,15,1900), 0.46, 'Tot. Annual Export Model = '+strtrim(tot_area_expo, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
    99   xyouts, julday(7,15,1900), 0.42, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
    100   ;title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)' 
     97  xyouts, julday(5,15,1900), 0.49, 'Tot. Annual Export OBS = 0.851 million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     98  xyouts, julday(5,15,1900), 0.46, 'Tot. Annual Export Model = '+strtrim(tot_area_expo, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     99  xyouts, julday(7,15,1900), 0.41, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
     100  ; 
    101101  title = 'Fram Strait Volume Export'+'!C'+prefix+' '+d1_d2 
    102   pltt, ICE_vol_export, 't',  /REMPLI, MIN = 0., MAX = .6 , /NOERASE, XGRIDSTYLE = 1 $ 
     102  pltt, ICE_vol_export, 't',  /REMPLI, MIN = 60., MAX = 500. , /NOERASE, XGRIDSTYLE = 1 $ 
    103103       , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 
    104104  pltt, vol_obs, 't', /REMPLI, /NOERASE, psym = 4, THICK = 4  $    ; light blue 
     
    107107  tot_vol_expo = total(ICE_vol_export.arr) 
    108108  tot_vol_expo_Sv = tot_vol_expo * 1.e06 * 1/86400 * 1/365 ; annual mean in Sverdrup 
    109   xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 2124 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
    110   xyouts, julday(6,15,1900), 0.46, 'Tot. Annual Export Model = '+strtrim(tot_vol_expo, 1)+' Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
    111   xyouts, julday(9,15,1900), 0.44, 'in Sv = '+strtrim(tot_vol_expo_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
    112   xyouts, julday(7,15,1900), 0.38, 'Data from Kwok et al.(2004), 1992-1998 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
     109  xyouts, julday(5,15,1900), 490, 'Tot. Annual Export OBS = 2124 10^3 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     110  xyouts, julday(5,15,1900), 460, 'Tot. Annual Export Model = '+strtrim(tot_vol_expo, 1)+' 10^3 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     111  xyouts, julday(9,15,1900), 440, 'in Sv = '+strtrim(tot_vol_expo_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     112  xyouts, julday(7,15,1900), 360, 'Data from Kwok et al.(2004), 1992-1998 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
    113113 
    114114  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ] 
     
    147147    ICE_vol_export_2 = total(ICE_vol_export_2,2)/nyr 
    148148    ; ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'}     ; annual mean 
    149     ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-12 * 86400 * 30 , unit : '10^3 Km^3/month'}       ; monthly mean 
     149    ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-9 * 86400 * 30 , unit : '10^3 Km^3/month'}       ; monthly mean 
    150150    ; 
    151151    ; 
     
    165165    tot_area_expo_2 = total(ICE_area_export_2.arr) 
    166166    ; 
    167     xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 0.851 million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
    168     xyouts, julday(6,15,1900), 0.46, 'Tot. Annual Export Model 1= '+strtrim(tot_area_expo, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
    169     xyouts, julday(6,15,1900), 0.44, 'Tot. Annual Export Model 2= '+strtrim(tot_area_expo_2, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
    170     xyouts, julday(7,15,1900), 0.42, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
     167    xyouts, julday(5,15,1900), 0.49, 'Tot. Annual Export OBS = 0.851 million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     168    xyouts, julday(5,15,1900), 0.46, 'Tot. Annual Export Model 1= '+strtrim(tot_area_expo, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     169    xyouts, julday(5,15,1900), 0.43, 'Tot. Annual Export Model 2= '+strtrim(tot_area_expo_2, 1)+' million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     170    xyouts, julday(7,15,1900), 0.41, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
     171    ; 
    171172    title = 'Fram Strait Volume Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 
    172     pltt, ICE_vol_export, 't', /REMPLI,  MIN = 0., MAX = .6, /NOERASE, XGRIDSTYLE = 1 $ 
     173    pltt, ICE_vol_export, 't', /REMPLI,  MIN = 60., MAX = 500., /NOERASE, XGRIDSTYLE = 1 $ 
    173174        , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex 
    174175    pltt, ICE_vol_export_2, 't', /REMPLI, /NOERASE  $ 
     
    181182    tot_vol_expo_Sv = tot_vol_expo * 1.e06 * 1/86400 * 1/365 ; annual mean in Sverdrup 
    182183    tot_vol_expo_2_Sv = tot_vol_expo_2 * 1.e06 * 1/86400 * 1/365 ; annual mean in Sverdrup 
    183     xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 2124 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
    184     xyouts, julday(6,15,1900), 0.46, 'Tot. Annual Export Model 1 = '+strtrim(tot_vol_expo, 1)+' Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
    185     xyouts, julday(9,15,1900), 0.44, 'in Sv = '+strtrim(tot_vol_expo_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
    186     xyouts, julday(6,15,1900), 0.42, 'Tot. Annual Export Model 2 = '+strtrim(tot_vol_expo_2, 1)+' Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
    187     xyouts, julday(9,15,1900), 0.40, 'in Sv = '+strtrim(tot_vol_expo_2_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
    188     xyouts, julday(7,15,1900), 0.38, 'Data from Kwok et al.(2004), 1992-1998 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
     184    xyouts, julday(5,15,1900), 490, 'Tot. Annual Export OBS = 2124 10^3 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     185    xyouts, julday(5,15,1900), 460, 'Tot. Annual Export Model 1 = '+strtrim(tot_vol_expo, 1)+' 10^3 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     186    xyouts, julday(9,15,1900), 440, 'in Sv = '+strtrim(tot_vol_expo_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     187    xyouts, julday(5,15,1900), 410, 'Tot. Annual Export Model 2 = '+strtrim(tot_vol_expo_2, 1)+' 10^3 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     188    xyouts, julday(9,15,1900), 390, 'in Sv = '+strtrim(tot_vol_expo_2_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2 
     189    xyouts, julday(7,15,1900), 360, 'Data from Kwok et al.(2004), 1992-1998 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2 
    189190 
    190191    htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png  />  ' ] 
  • branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro

    r4504 r4514  
    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 = 'ICE Velocity ' 
    81 ;SF  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, 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.