Index: /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceAge.pro
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceAge.pro	(revision 4514)
+++ /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceAge.pro	(revision 4514)
@@ -0,0 +1,58 @@
+pro std_plot_IceAge, Iage1, Iage2, ARC = arc, ANT = ant, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex
+
+  compile_opt idl2, strictarrsubs
+
+@common  
+@std_common
+
+  cdti3 = string(cnt, format = '(i3.3)')
+  print, cdti3 + ') ' + blabla
+
+  var = 'IceAge'
+  IF keyword_set(arc) THEN var = var+'_Arc_'
+  IF keyword_set(ant) THEN var = var+'_Ant_'
+  IF keyword_set(march) THEN var = var+'March'
+  IF keyword_set(sept) THEN var = var+'Sept'
+
+  filename = cdti3 + '_'+var+'_'+std_file1_I
+  IF std_file1_I NE std_file2_I then filename = filename + '_'+std_file2_I
+  IF KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1
+  
+  IF keyword_set(arc) THEN BEGIN 
+    domdef, 20, 380, 50, 90
+    vmin = 0.
+    vmax = 7.
+    vint = 1. 
+    fmtd = '(f4.1)'
+    map = [90, 0, 0]
+  ENDIF
+  IF keyword_set(ant) THEN BEGIN 
+    domdef, 20, 380, -90, -50
+    vmin = 0.
+    vmax = 2.
+    vint = 0.5
+    fmt = '(f4.1)'
+    map = [-90, 0, 0]
+  ENDIF
+;
+  varunit = Iage1.unit
+;
+  title = var+'!C'+std_file1_I
+  plt, Iage1.arr, MIN = vmin, MAX = vmax, INTER = vint, /STRICTFILL, CELL_FILL = 2, format = fmt $
+       , small = [1, 2, 1], COAST_THICK = 2, TITLE = title, divisions = div  $
+       , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex
+;                                ;
+  if std_file1_I NE std_file2_I then begin            
+    title = title + std_file2_I
+    plt, Iage1.arr - Iage2.arr, MIN = vmind, MAX = vmaxd, INTER = vintd, STYLE = 'so0so', format = fmtd $
+         , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title  $
+         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex
+  endif                        
+  
+  domdef
+
+  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
+  if KEYWORD_SET(postscript) then closeps
+
+  return
+end
Index: /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceSal.pro
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceSal.pro	(revision 4513)
+++ /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_IceSal.pro	(revision 4514)
@@ -1,3 +1,3 @@
-pro std_plot_IceSal, Isal1, Isal2, ARC = arc, ANT = ant, APRIL = april, SEPT = sept, POSTSCRIPT = postscript, _extra = ex
+pro std_plot_IceSal, Isal1, Isal2, ARC = arc, ANT = ant, FEBR = febr, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex
 
   compile_opt idl2, strictarrsubs
@@ -12,5 +12,6 @@
   IF keyword_set(arc) THEN var = var+'_Arc_'
   IF keyword_set(ant) THEN var = var+'_Ant_'
-  IF keyword_set(april) THEN var = var+'April'
+  IF keyword_set(febr) THEN var = var+'Febr'
+  IF keyword_set(march) THEN var = var+'March'
   IF keyword_set(sept) THEN var = var+'Sept'
 
@@ -31,5 +32,5 @@
   ;
    title = var+'!C'+std_file1_I
-   plt, Isal1.arr, MIN = .1, MAX = 20.1, INTER = 1., /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $
+   plt, Isal1.arr, MIN = 0., MAX = 12., INTER = 1., /STRICTFILL, CELL_FILL = 2, format = '(f4.1)' $
         , small = [1, 2, 1], COAST_THICK = 2, TITLE = title $
         , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /PORTRAIT, _extra = ex
@@ -37,5 +38,5 @@
    if std_file1_I NE std_file2_I then begin
      title = title + std_file2_I
-     plt, Isal1.arr - Isal2.arr, MIN = -4, MAX = 4, INTER = .2, STYLE = 'so0so', format = '(f4.1)' $
+     plt, Isal1.arr - Isal2.arr, MIN = -4, MAX = 4, INTER = 1., STYLE = 'so0so', format = '(f4.1)' $
           , small = [1, 2, 2], COAST_THICK = 2, CELL_FILL = 2, TITLE = title $
           , CHARSIZE = 1.05, GLINETHICK = 2., /ORTHO, MAP = map, /NOERASE, _extra = ex
Index: /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_SnowThick.pro
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_SnowThick.pro	(revision 4513)
+++ /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_SnowThick.pro	(revision 4514)
@@ -1,3 +1,3 @@
-pro std_plot_SnowThick, Isnow1, Isnow2, ARC = arc, ANT = ant, APRIL = april, SEPT = sept, POSTSCRIPT = postscript, _extra = ex
+pro std_plot_SnowThick, Isnow1, Isnow2, ARC = arc, ANT = ant, MARCH = march, SEPT = sept, POSTSCRIPT = postscript, _extra = ex
 
   compile_opt idl2, strictarrsubs
@@ -12,5 +12,5 @@
   IF keyword_set(arc) THEN var = var+'_Arc_'
   IF keyword_set(ant) THEN var = var+'_Ant_'
-  IF keyword_set(april) THEN var = var+'April'
+  IF keyword_set(march) THEN var = var+'March'
   IF keyword_set(sept) THEN var = var+'Sept'
 
@@ -22,5 +22,5 @@
     domdef, 20, 380, 50, 90
     vmin = 0.
-    vmax = 1.
+    vmax = .6
     vint = 0.1
     ;fmt = '(i2)'
Index: /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro	(revision 4513)
+++ /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all.pro	(revision 4514)
@@ -187,6 +187,26 @@
   ENDELSE
 ;
+  Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I ) 
+  caldat, time, mm
+  march = where(mm EQ 3, cnt)
+  Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit}
+  sept = where(mm EQ 9, cnt)
+  Iage_sept_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
+  undefine, Iage_1
+;
+  IF std_file2_I NE std_file1_I THEN BEGIN
+     Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I ) 
+     caldat, time, mm
+     march = where(mm EQ 3, cnt)
+     Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 
+     sept = where(mm EQ 9, cnt)
+     Iage_sept_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 
+     undefine, Iage_2
+  ENDIF ELSE BEGIN 
+     Iage_march_2 = {arr:-1}
+     Iage_sept_2 = {arr:-1}
+  ENDELSE
+;
   Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 
-  ;help, Ifra_1.arr, jpi,jpj
   caldat, time, mm
   march = where(mm EQ 3, cnt)
@@ -210,8 +230,7 @@
 ;
    Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 
-   ;help, Isnow_1.arr, jpi,jpj
    caldat, time, mm
-   april = where(mm EQ 4, cnt)
-   Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 
+   march = where(mm EQ 3, cnt)
+   Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 
    sept = where(mm EQ 9, cnt)
    Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}
@@ -221,33 +240,49 @@
      Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 
      caldat, time, mm
-     april = where(mm EQ 4, cnt)
-     Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit}
+     march = where(mm EQ 3, cnt)
+     Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit}
      sept = where(mm EQ 9, cnt)
      Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit}
    ;  undefine, Isnow_2
    ENDIF ELSE BEGIN 
-     Isnow_april_2 = {arr:-1}
+     Isnow_march_2 = {arr:-1}
      Isnow_sept_2 = {arr:-1}
    ENDELSE
 ;
    Isal_1 = read_ncdf(getenv('VAR1_Isal'), allrecords = allrec, filename = std_file1_I ) 
-   ;help, Isal_1.arr, jpi,jpj
+   Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )
+   msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations 
+   help, Ifra_1
+   help, mask
    caldat, time, mm
-   april = where(mm EQ 4, cnt)
-   Isal_april_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 
+   march = where(mm EQ 3, cnt)
+   febr = where(mm EQ 2, cnt)
+   ;SF  Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 
+   Isal_1.arr = Isal_1.arr * msk
+   help, Isal_1
+   Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 
+   help, Isal_march_1
+   Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 
    sept = where(mm EQ 9, cnt)
    Isal_sept_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_1.unit}
    ;undefine, Isal_1
  ;
+   Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I ) 
+   Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I )
+   msk = Ifra_2.arr gt 0.15 ; remove 0.15% for observations 
    IF std_file2_I NE std_file1_I THEN BEGIN
      Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I ) 
+     Isal_2.arr = Isal_2.arr * msk
      caldat, time, mm
-     april = where(mm EQ 4, cnt)
-     Isal_april_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isal_2.unit}
+     march = where(mm EQ 3, cnt)
+     febr = where(mm EQ 2, cnt)
+     Isal_march_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_2.unit}
+     Isal_febr_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_2.unit}
      sept = where(mm EQ 9, cnt)
      Isal_sept_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_2.unit}
    ;  undefine, Isal_2
    ENDIF ELSE BEGIN 
-     Isal_april_2 = {arr:-1}
+     Isal_febr_2 = {arr:-1}
+     Isal_march_2 = {arr:-1}
      Isal_sept_2 = {arr:-1}
    ENDELSE
@@ -375,4 +410,16 @@
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ANT, /SEPT, _extra = ex
 ;
+  cnt = cnt+1   &   blabla = 'Arctic Ice Age: MARCH'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ARC, /MARCH, _extra = ex
+;  
+  cnt = cnt+1   &   blabla = 'Arctic Ice Age: SEPT'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ARC, /SEPT, _extra = ex
+;  
+  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: MARCH'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ANT, /MARCH, _extra = ex
+;
+  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: SEPT'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ANT, /SEPT, _extra = ex
+;
   cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: MARCH'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex
@@ -387,24 +434,24 @@
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex
 ;
-  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: APRIL'
-  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex
+  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: MARCH'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ARC, /MARCH, _extra = ex
 ; 
   cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex
 ;  
-  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: APRIL'
-  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex
+  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: MARCH'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex
 ; 
   cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex
 ; 
-  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: APRIL'
-  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_april_1, Isal_april_2, /ARC, /APRIL, _extra = ex
+  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: MARCH'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_march_1, Isal_march_2, /ARC, /MARCH, _extra = ex
 ; 
   cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: SEPT'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ARC, /SEPT, _extra = ex
 ;  
-  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: APRIL'
-  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_april_1, Isal_april_2, /ANT, /APRIL, _extra = ex
+  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: FEBRUARY'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_febr_1, Isal_febr_2, /ANT, /FEBR, _extra = ex
 ; 
   cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: SEPT'
Index: /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_1Y.pro
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_1Y.pro	(revision 4513)
+++ /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_1Y.pro	(revision 4514)
@@ -115,43 +115,53 @@
   EMP1 = read_ncdf(getenv('VAR1_EMP'), allrecords = allrec, direc = 't', filename = std_file1_T )
   EMP1 = {arr:EMP1.arr * 86400., unit:'mm/day', grid:'T'}
-  IF std_file2_T NE std_file1_T THEN BEGIN
-    EMP2 = read_ncdf(getenv('VAR2_EMP'), allrecords = allrec, direc = 't', filename = std_file2_T )
-    EMP2 = {arr:EMP2.arr * 86400., unit:'mm/day', grid:'T'}
-  ENDIF ELSE EMP2 = {arr:-1}
-  
-;mixed layer depth
-  MLD1 = read_ncdf(getenv('VAR1_MLD'), allrecords = allrec, direc = 't', filename = std_file1_T ) ; 10 m
-  IF std_file2_T NE std_file1_T THEN BEGIN
-    MLD2 = read_ncdf(getenv('VAR2_MLD'), allrecords = allrec, direc = 't', filename = std_file2_T ) ; 10 m
-  ENDIF ELSE MLD2 = {arr:-1}
-;climatology 
-  MLD = read_ncdf(getenv('VAR_MLD'), filename = std_file_mld )
-
-; velocities
-  U1 = read_ncdf(getenv('VAR1_U'), allrecords = allrec, direc = 't', filename = std_file1_U )
-  IF strlowcase(getenv('VAR1_U')) EQ 'uocetr_eff' THEN BEGIN
-    U1.arr = U1.arr / e3u_3d(/e2) * umask()
-    U1.unit = 'm/s'
-  ENDIF 
-  IF std_file2_U NE std_file1_U THEN BEGIN
-    U2 = read_ncdf(getenv('VAR2_U'), allrecords = allrec, direc = 't', filename = std_file2_U )
-    IF strlowcase(getenv('VAR2_U')) EQ 'uocetr_eff' THEN BEGIN
-      U2.arr = U2.arr / e3u_3d(/e2) * umask()
-      U2.unit = 'm/s'
-    ENDIF 
-  ENDIF ELSE U2 = {arr:-1}
-;
-  V1 = read_ncdf(getenv('VAR1_V'), allrecords = allrec, direc = 't', filename = std_file1_V )
-  IF strlowcase(getenv('VAR1_V')) EQ 'vocetr_eff' THEN BEGIN
-    V1.arr = V1.arr / e3v_3d(/e1) * vmask()
-    V1.unit = 'm/s'
-  ENDIF 
-  IF std_file2_V NE std_file1_V THEN BEGIN
-    V2 = read_ncdf(getenv('VAR2_V'), allrecords = allrec, direc = 't', filename = std_file2_V )
-    IF strlowcase(getenv('VAR2_V')) EQ 'vocetr_eff' THEN BEGIN
-      V2.arr = V2.arr / e3v_3d(/e1) * vmask()
-      V2.unit = 'm/s'
-    ENDIF 
-  ENDIF ELSE V2 = {arr:-1}
+   IF std_file2_T NE std_file1_T THEN BEGIN
+     EMP2 = read_ncdf(getenv('VAR2_EMP'), allrecords = allrec, direc = 't', filename = std_file2_T )
+     EMP2 = {arr:EMP2.arr * 86400., unit:'mm/day', grid:'T'}
+   ENDIF ELSE EMP2 = {arr:-1}
+   
+ ;mixed layer depth
+   MLD1 = read_ncdf(getenv('VAR1_MLD'), allrecords = allrec, direc = 't', filename = std_file1_T ) ; 10 m
+   IF std_file2_T NE std_file1_T THEN BEGIN
+     MLD2 = read_ncdf(getenv('VAR2_MLD'), allrecords = allrec, direc = 't', filename = std_file2_T ) ; 10 m
+   ENDIF ELSE MLD2 = {arr:-1}
+ ;climatology 
+   MLD = read_ncdf(getenv('VAR_MLD'), filename = std_file_mld )
+
+ ; velocities
+   U1 = read_ncdf(getenv('VAR1_U'), allrecords = allrec, direc = 't', filename = std_file1_U )
+   ; old formulation: we tested variable name
+   ; IF strlowcase(getenv('VAR1_U')) EQ 'uocetr_eff' OR strlowcase(getenv('VAR1_U')) EQ 'vozoeftr' THEN BEGIN
+   IF strlowcase(U1.unit) EQ 'm3/s' THEN BEGIN
+     ;IF it is a transport it is transofrmed in velocity
+     U1.arr = U1.arr / e3u_3d(/e2) * umask()
+     U1.unit = 'm/s'
+   ENDIF 
+   IF std_file2_U NE std_file1_U THEN BEGIN
+     U2 = read_ncdf(getenv('VAR2_U'), allrecords = allrec, direc = 't', filename = std_file2_U )
+     ; old formulation: we tested variable name
+     ; IF strlowcase(getenv('VAR2_U')) EQ 'uocetr_eff' OR strlowcase(getenv('VAR2_U')) EQ 'vozoeftr' THEN BEGIN
+     IF strlowcase(U2.unit) EQ 'm3/s' THEN BEGIN
+       U2.arr = U2.arr / e3u_3d(/e2) * umask()
+       U2.unit = 'm/s'
+     ENDIF 
+   ENDIF ELSE U2 = {arr:-1}
+;
+   V1 = read_ncdf(getenv('VAR1_V'), allrecords = allrec, direc = 't', filename = std_file1_V )
+   ; old formulation: we tested variable name
+   ; IF strlowcase(getenv('VAR1_V')) EQ 'vocetr_eff' OR strlowcase(getenv('VAR1_V')) EQ 'vomeeftr' THEN BEGIN
+   IF strlowcase(V1.unit) EQ 'm3/s' THEN BEGIN ; test on unit to understand if it is a transport or velocity
+     ;IF it is a transport it is transofrmed in velocity
+     V1.arr = V1.arr / e3v_3d(/e1) * vmask()
+     V1.unit = 'm/s'
+   ENDIF 
+   IF std_file2_V NE std_file1_V THEN BEGIN
+     V2 = read_ncdf(getenv('VAR2_V'), allrecords = allrec, direc = 't', filename = std_file2_V )
+     ; old formulation
+     ; IF strlowcase(getenv('VAR2_V')) EQ 'vocetr_eff' OR strlowcase(getenv('VAR2_V')) EQ 'vozoeftr' THEN BEGIN
+     IF strlowcase(V2.unit) EQ 'm3/s' THEN BEGIN 
+       V2.arr = V2.arr / e3v_3d(/e1) * vmask()
+       V2.unit = 'm/s'
+     ENDIF 
+   ENDIF ELSE V2 = {arr:-1}
 
 ; ice
@@ -177,6 +187,26 @@
   ENDELSE
 ;
+  Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I ) 
+  caldat, time, mm
+  march = where(mm EQ 3, cnt)
+  Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit}
+  sept = where(mm EQ 9, cnt)
+  Iage_sept_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
+  undefine, Iage_1
+;
+  IF std_file2_I NE std_file1_I THEN BEGIN
+     Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I ) 
+     caldat, time, mm
+     march = where(mm EQ 3, cnt)
+     Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 
+     sept = where(mm EQ 9, cnt)
+     Iage_sept_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 
+     undefine, Iage_2
+  ENDIF ELSE BEGIN 
+     Iage_march_2 = {arr:-1}
+     Iage_sept_2 = {arr:-1}
+  ENDELSE
+;
   Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 
-  help, Ifra_1.arr, jpi,jpj
   caldat, time, mm
   march = where(mm EQ 3, cnt)
@@ -199,26 +229,62 @@
   ENDELSE
 ;
-;; Sf commentato perché non funziona snow
-;; Sf   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 
-;; Sf   help, Isnow_1.arr, jpi,jpj
-;; Sf   caldat, time, mm
-;; Sf   april = where(mm EQ 4, cnt)
-;; Sf   Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 
-;; Sf   sept = where(mm EQ 9, cnt)
-;; Sf   Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}
-;; Sf   ;undefine, Isnow_1
-;; Sf ;
-;; Sf   IF std_file2_I NE std_file1_I THEN BEGIN
-;; Sf     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 
-;; Sf     caldat, time, mm
-;; Sf     april = where(mm EQ 4, cnt)
-;; Sf     Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit}
-;; Sf     sept = where(mm EQ 9, cnt)
-;; Sf     Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit}
-;; Sf   ;  undefine, Isnow_2
-;; Sf   ENDIF ELSE BEGIN 
-;; Sf     Isnow_april_2 = {arr:-1}
-;; Sf     Isnow_sept_2 = {arr:-1}
-;; Sf   ENDELSE
+   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 
+   caldat, time, mm
+   march = where(mm EQ 3, cnt)
+   Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 
+   sept = where(mm EQ 9, cnt)
+   Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}
+   ;undefine, Isnow_1
+ ;
+   IF std_file2_I NE std_file1_I THEN BEGIN
+     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 
+     caldat, time, mm
+     march = where(mm EQ 3, cnt)
+     Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit}
+     sept = where(mm EQ 9, cnt)
+     Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit}
+   ;  undefine, Isnow_2
+   ENDIF ELSE BEGIN 
+     Isnow_march_2 = {arr:-1}
+     Isnow_sept_2 = {arr:-1}
+   ENDELSE
+;
+   Isal_1 = read_ncdf(getenv('VAR1_Isal'), allrecords = allrec, filename = std_file1_I ) 
+   Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )
+   msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations 
+   help, Ifra_1
+   help, mask
+   caldat, time, mm
+   march = where(mm EQ 3, cnt)
+   febr = where(mm EQ 2, cnt)
+   ;SF  Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 
+   Isal_1.arr = Isal_1.arr * msk
+   help, Isal_1
+   Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 
+   help, Isal_march_1
+   Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 
+   sept = where(mm EQ 9, cnt)
+   Isal_sept_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_1.unit}
+   ;undefine, Isal_1
+ ;
+   Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I ) 
+   Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I )
+   msk = Ifra_2.arr gt 0.15 ; remove 0.15% for observations 
+   IF std_file2_I NE std_file1_I THEN BEGIN
+     Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I ) 
+     Isal_2.arr = Isal_2.arr * msk
+     caldat, time, mm
+     march = where(mm EQ 3, cnt)
+     febr = where(mm EQ 2, cnt)
+     Isal_march_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_2.unit}
+     Isal_febr_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_2.unit}
+     sept = where(mm EQ 9, cnt)
+     Isal_sept_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_2.unit}
+   ;  undefine, Isal_2
+   ENDIF ELSE BEGIN 
+     Isal_febr_2 = {arr:-1}
+     Isal_march_2 = {arr:-1}
+     Isal_sept_2 = {arr:-1}
+   ENDELSE
 ;
   jpt = 1
@@ -256,158 +322,176 @@
   cnt = cnt+1   &   blabla = 'Erp salinity damping term'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_erp, ERP1, ERP2, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Evaporation - Precipitation - Runoff term'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_emp, EMP1, EMP2, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Net heat flux'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_qnet, Q1, Q2, QNET, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Meridionnal Heat Transport'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mht, Q1.arr+geo, Q2.arr+geo, masknp, std_file_msksub, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Global Barotropic stream Function'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_bsf, U1, U2, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'mean Temperature diff with New Reynolds'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sst, T1, T2, TRey, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'mean Salinity diff with Levitus'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sss, S1, S2, SLev, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, _extra = ex       
-
+;
   cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus and exp2'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus at z=100 meters'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, /z100, _extra = ex       
- 
+; 
   cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus and exp2 at z=100 meters'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, /z100, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'mean Temperature diff with Levitus at z=100 meters'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_T100m, T1, T2, Tlev, _extra = ex
-  
+;  
   cnt = cnt+1   &   blabla = 'mean Salinity diff with Levitus at z=100 meters'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_S100m, S1, S2, SLev, _extra = ex
-  
+;  
   cnt = cnt+1   &   blabla = 'Mixed layer depth'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Mixed layer depth differences'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD2, MLD, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Zonal mean Mixed layer depth'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ZonMld, MLD1, MLD2, MLD, _extra = ex
-  
+;  
   cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Global'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Atlantic'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Atl', _extra = ex
-  
+;  
   cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Indian'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Ind', _extra = ex
-  
+;  
   cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Pacific'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Pac', _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Global'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Atlantic'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Atl', _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Indian'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Ind', _extra = ex 
-
+;
   cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Pacific'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Pac', _extra = ex 
-  
+;  
 ;  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: MARCH'
 ;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: APRIL'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: SEPT'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex
-  
+;  
 ;  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: MARCH'
 ;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex
-
+; 
   cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: APRIL'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: SEPT'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ANT, /SEPT, _extra = ex
-
+;
+  cnt = cnt+1   &   blabla = 'Arctic Ice Age: MARCH'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ARC, /MARCH, _extra = ex
+;  
+  cnt = cnt+1   &   blabla = 'Arctic Ice Age: SEPT'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ARC, /SEPT, _extra = ex
+;  
+  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: MARCH'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ANT, /MARCH, _extra = ex
+;
+  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: SEPT'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ANT, /SEPT, _extra = ex
+;
   cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: MARCH'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex
-  
+;  
   cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: SEPT'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ARC, /SEPT, _extra = ex
-  
+;  
   cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: MARCH'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: SEPT'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex
-;; SF non funziona, dice che deve essere una struttura e non e!!!!
-;;SF  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: APRIL'
-;;SF  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex
-;; SF 
-;; SF   cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT'
-;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex
-;; SF  
-;; SF   cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: APRIL'
-;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex
-;; SF  
-;; SF   cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT'
-;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex
-; 
-;  cnt = cnt+1   &   blabla = 'Arctic Ice Extend: '
-;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC,  _extra = ex
-  
-;  cnt = cnt+1   &   blabla = 'Antarctic Ice Extend: '
-;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT,  _extra = ex
-
+;
+  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: MARCH'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ARC, /MARCH, _extra = ex
+; 
+  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex
+;  
+  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: MARCH'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex
+; 
+  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex
+; 
+  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: MARCH'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_march_1, Isal_march_2, /ARC, /MARCH, _extra = ex
+; 
+  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: SEPT'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ARC, /SEPT, _extra = ex
+;  
+  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: FEBRUARY'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_febr_1, Isal_febr_2, /ANT, /FEBR, _extra = ex
+; 
+  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: SEPT'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ANT, /SEPT, _extra = ex
+; 
   cnt = cnt+1   &   blabla = 'Meridional stream Function: Global (no Med)'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'GloNoMed', _extra = ex
-
+;  
   cnt = cnt+1   &   blabla = 'Meridional stream Function: Atlantic'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Atl', _extra = ex
-
+ ;
   cnt = cnt+1   &   blabla = 'Meridional stream Function: Indian'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Ind', _extra = ex
-
+ ;
   cnt = cnt+1   &   blabla = 'Meridional stream Function: Indo-Pacific'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'IndoPac', _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Equatorial Temperature'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqT, T1, T2, Tlev, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Equatorial Salinity'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqS, S1, S2, SLev, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Equatorial zonal velocity'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqU, U1, U2, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Mediterranean salt tongue at depth=700'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 700, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Mediterranean salt tongue at depth=1000'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 1000, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Mediterranean water at lat=40°N'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 40, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Mediterranean water at lat=38°N'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Vertical Global mean T & S'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_GlobMeanTS, T1, T2, TLev, S1, S2, SLev, _extra = ex
-
+;
   IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_plot_html_body.txt', htmltxt[1:*]
   
Index: /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_last_year_5D.pro
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_last_year_5D.pro	(revision 4513)
+++ /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_all_last_year_5D.pro	(revision 4514)
@@ -203,4 +203,25 @@
   ENDELSE
 ;
+  Iage_1 = read_ncdf(getenv('VAR1_Iage'), allrecords = allrec, filename = std_file1_I ) 
+  caldat, time, mm
+  march = where(mm EQ 3, cnt)
+  Iage_march_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Iage_1.unit}
+  sept = where(mm EQ 9, cnt)
+  Iage_sept_1 = {arr:1./float(cnt) * total(reform(Iage_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_1.unit} 
+  undefine, Iage_1
+;
+  IF std_file2_I NE std_file1_I THEN BEGIN
+     Iage_2 = read_ncdf(getenv('VAR2_Iage'), allrecords = allrec, filename = std_file2_I ) 
+     caldat, time, mm
+     march = where(mm EQ 3, cnt)
+     Iage_march_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 
+     sept = where(mm EQ 9, cnt)
+     Iage_sept_2 = {arr:1./float(cnt) * total(reform(Iage_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Iage_2.unit} 
+     undefine, Iage_2
+  ENDIF ELSE BEGIN 
+     Iage_march_2 = {arr:-1}
+     Iage_sept_2 = {arr:-1}
+  ENDELSE
+;
   Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I ) 
   ;; help, Ifra_1.arr, jpi,jpj
@@ -225,26 +246,62 @@
   ENDELSE
 ;
-;; SF commentato perché non funziona snow
-;; SF   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 
-;; SF   help, Isnow_1.arr, jpi,jpj
-;; SF   caldat, time, mm
-;; SF   april = where(mm EQ 4, cnt)
-;; SF   Isnow_april_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 
-;; SF   sept = where(mm EQ 9, cnt)
-;; SF   Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}
-;; SF   ;undefine, Isnow_1
-;; SF ;
-;; SF   IF std_file2_I NE std_file1_I THEN BEGIN
-;; SF     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 
-;; SF     caldat, time, mm
-;; SF     april = where(mm EQ 4, cnt)
-;; SF     Isnow_april_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(april)],nxt,nyt,cnt), 3), unit:Isnow_2.unit}
-;; SF     sept = where(mm EQ 9, cnt)
-;; SF     Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit}
-;; SF   ;  undefine, Isnow_2
-;; SF   ENDIF ELSE BEGIN 
-;; SF     Isnow_april_2 = {arr:-1}
-;; SF     Isnow_sept_2 = {arr:-1}
-;; SF   ENDELSE
+   Isnow_1 = read_ncdf(getenv('VAR1_Isnow'), allrecords = allrec, filename = std_file1_I ) 
+   caldat, time, mm
+   march = where(mm EQ 3, cnt)
+   Isnow_march_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_1.unit} 
+   sept = where(mm EQ 9, cnt)
+   Isnow_sept_1 = {arr:1./float(cnt) * total(reform(Isnow_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_1.unit}
+   ;undefine, Isnow_1
+ ;
+   IF std_file2_I NE std_file1_I THEN BEGIN
+     Isnow_2 = read_ncdf(getenv('VAR2_Isnow'), allrecords = allrec, filename = std_file2_I ) 
+     caldat, time, mm
+     march = where(mm EQ 3, cnt)
+     Isnow_march_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isnow_2.unit}
+     sept = where(mm EQ 9, cnt)
+     Isnow_sept_2 = {arr:1./float(cnt) * total(reform(Isnow_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isnow_2.unit}
+   ;  undefine, Isnow_2
+   ENDIF ELSE BEGIN 
+     Isnow_march_2 = {arr:-1}
+     Isnow_sept_2 = {arr:-1}
+   ENDELSE
+;
+   Isal_1 = read_ncdf(getenv('VAR1_Isal'), allrecords = allrec, filename = std_file1_I ) 
+   Ifra_1 = read_ncdf(getenv('VAR1_Ifrac'), allrecords = allrec, filename = std_file1_I )
+   msk = Ifra_1.arr gt 0.15 ; remove 0.15% for observations 
+   help, Ifra_1
+   help, mask
+   caldat, time, mm
+   march = where(mm EQ 3, cnt)
+   febr = where(mm EQ 2, cnt)
+   ;SF  Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 
+   Isal_1.arr = Isal_1.arr * msk
+   help, Isal_1
+   Isal_march_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 
+   help, Isal_march_1
+   Isal_febr_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_1.unit} 
+   sept = where(mm EQ 9, cnt)
+   Isal_sept_1 = {arr:1./float(cnt) * total(reform(Isal_1.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_1.unit}
+   ;undefine, Isal_1
+ ;
+   Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I ) 
+   Ifra_2 = read_ncdf(getenv('VAR2_Ifrac'), allrecords = allrec, filename = std_file2_I )
+   msk = Ifra_2.arr gt 0.15 ; remove 0.15% for observations 
+   IF std_file2_I NE std_file1_I THEN BEGIN
+     Isal_2 = read_ncdf(getenv('VAR2_Isal'), allrecords = allrec, filename = std_file2_I ) 
+     Isal_2.arr = Isal_2.arr * msk
+     caldat, time, mm
+     march = where(mm EQ 3, cnt)
+     febr = where(mm EQ 2, cnt)
+     Isal_march_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(march)],nxt,nyt,cnt), 3), unit:Isal_2.unit}
+     Isal_febr_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(febr)],nxt,nyt,cnt), 3), unit:Isal_2.unit}
+     sept = where(mm EQ 9, cnt)
+     Isal_sept_2 = {arr:1./float(cnt) * total(reform(Isal_2.arr[*, *, temporary(sept)],nxt,nyt,cnt), 3), unit:Isal_2.unit}
+   ;  undefine, Isal_2
+   ENDIF ELSE BEGIN 
+     Isal_febr_2 = {arr:-1}
+     Isal_march_2 = {arr:-1}
+     Isal_sept_2 = {arr:-1}
+   ENDELSE
 ;
   jpt = 1
@@ -282,45 +339,44 @@
   cnt = cnt+1   &   blabla = 'Erp salinity damping term'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_erp, ERP1, ERP2, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Evaporation - Precipitation - Runoff term'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_emp, EMP1, EMP2, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Net heat flux'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_qnet, Q1, Q2, QNET, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Meridionnal Heat Transport'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mht, Q1.arr+geo, Q2.arr+geo, masknp, std_file_msksub, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Global Barotropic stream Function'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_bsf, U1, U2, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'mean Temperature diff with New Reynolds'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sst, T1, T2, TRey, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'mean Salinity diff with Levitus'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_sss, S1, S2, SLev, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, _extra = ex       
-
+;
   cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus and exp2'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus at z=100 meters'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, SLev, /z100, _extra = ex       
- 
+; 
   cnt = cnt+1   &   blabla = 'Arctic mean Salinity diff with Levitus and exp2 at z=100 meters'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_ArcSal, S1, S2, SLev, /z100, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'mean Temperature diff with Levitus at z=100 meters'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_T100m, T1, T2, Tlev, _extra = ex
-  
+;  
   cnt = cnt+1   &   blabla = 'mean Salinity diff with Levitus at z=100 meters'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_S100m, S1, S2, SLev, _extra = ex
- 
-;;SF need monthly climatology to plot difference of MAX of MLD  
-  cnt = cnt+1   &   blabla = 'Mixed layer depth & difference with DeBoyer'
-  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD_1, MLD, _extra = ex
-   
+;  
+  cnt = cnt+1   &   blabla = 'Mixed layer depth'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD1, MLD, _extra = ex
+;
   cnt = cnt+1   &   blabla = 'Mixed layer depth differences'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_mld, MLD_1, MLD_2, MLD, _extra = ex
@@ -339,111 +395,128 @@
   cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Global'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Atlantic'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Atl', _extra = ex
-  
+;  
   cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Indian'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Ind', _extra = ex
-  
+;  
   cnt = cnt+1   &   blabla = 'Zonal mean Temperature diff with Levitus: Pacific'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_T, T1, T2, TLev, SUBBASIN = 'Pac', _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Global'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Atlantic'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Atl', _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Indian'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Ind', _extra = ex 
-
+;
   cnt = cnt+1   &   blabla = 'Zonal mean Salinity diff with Levitus: Pacific'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_zonal_S, S1, S2, SLev, SUBBASIN = 'Pac', _extra = ex 
-  
+;  
 ;  cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: MARCH'
 ;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ARC, /MARCH, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: APRIL'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ARC, /APRIL, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Arctic Ice Thickness: SEPT'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ARC, /SEPT, _extra = ex
-  
+;  
 ;  cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: MARCH'
 ;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_march_1, Ithi_march_2, /ANT, /MARCH, _extra = ex
-
+; 
   cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: APRIL'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_april_1, Ithi_april_2, /ANT, /APRIL, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Antarctic Ice Thickness: SEPT'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceThick, Ithi_sept_1, Ithi_sept_2, /ANT, /SEPT, _extra = ex
-
+;
+  cnt = cnt+1   &   blabla = 'Arctic Ice Age: MARCH'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ARC, /MARCH, _extra = ex
+;  
+  cnt = cnt+1   &   blabla = 'Arctic Ice Age: SEPT'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ARC, /SEPT, _extra = ex
+;  
+  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: MARCH'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_march_1, Iage_march_2, /ANT, /MARCH, _extra = ex
+;
+  cnt = cnt+1   &   blabla = 'Antarctic Ice Age: SEPT'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceAge, Iage_sept_1, Iage_sept_2, /ANT, /SEPT, _extra = ex
+;
   cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: MARCH'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ARC, /MARCH, _extra = ex
-  
+;  
   cnt = cnt+1   &   blabla = 'Arctic Ice Fraction: SEPT'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ARC, /SEPT, _extra = ex
-  
+;  
   cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: MARCH'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_march_1, Ifra_march_2, /ANT, /MARCH, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Antarctic Ice Fraction: SEPT'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceFrac, Ifra_sept_1, Ifra_sept_2, /ANT, /SEPT, _extra = ex
-
-;; SF non funziona, dice che deve essere una struttura e non e!!!!
-;;SF  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: APRIL'
-;;SF  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ARC, /APRIL, _extra = ex
-;; SF 
-;; SF   cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT'
-;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex
-;; SF  
-;; SF   cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: APRIL'
-;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_april_1, Isnow_april_2, /ANT, /APRIL, _extra = ex
-;; SF  
-;; SF   cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT'
-;; SF   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex
-; 
-;  cnt = cnt+1   &   blabla = 'Arctic Ice Extend: '
-;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ARC,  _extra = ex
-  
-;  cnt = cnt+1   &   blabla = 'Antarctic Ice Extend: '
-;  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceExt, Ifra_1, Ifra_2, masknp, /ANT,  _extra = ex
-
+;
+  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: MARCH'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ARC, /MARCH, _extra = ex
+; 
+  cnt = cnt+1   &   blabla = 'Arctic SNOW Thickness: SEPT'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ARC, /SEPT, _extra = ex
+;  
+  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: MARCH'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_march_1, Isnow_march_2, /ANT, /MARCH, _extra = ex
+; 
+  cnt = cnt+1   &   blabla = 'Antarctic SNOW Thickness: SEPT'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_SnowThick, Isnow_sept_1, Isnow_sept_2, /ANT, /SEPT, _extra = ex
+; 
+  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: MARCH'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_march_1, Isal_march_2, /ARC, /MARCH, _extra = ex
+; 
+  cnt = cnt+1   &   blabla = 'Arctic Ice Salinity: SEPT'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ARC, /SEPT, _extra = ex
+;  
+  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: FEBRUARY'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_febr_1, Isal_febr_2, /ANT, /FEBR, _extra = ex
+; 
+  cnt = cnt+1   &   blabla = 'Antarctic Ice Salinity: SEPT'
+  IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_IceSal, Isal_sept_1, Isal_sept_2, /ANT, /SEPT, _extra = ex
+; 
   cnt = cnt+1   &   blabla = 'Meridional stream Function: Global (no Med)'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'GloNoMed', _extra = ex
-
+;  
   cnt = cnt+1   &   blabla = 'Meridional stream Function: Atlantic'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Atl', _extra = ex
-
+ ;
   cnt = cnt+1   &   blabla = 'Meridional stream Function: Indian'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'Ind', _extra = ex
-
+ ;
   cnt = cnt+1   &   blabla = 'Meridional stream Function: Indo-Pacific'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_msf, V1, V2, SUBBASIN = 'IndoPac', _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Equatorial Temperature'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqT, T1, T2, Tlev, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Equatorial Salinity'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqS, S1, S2, SLev, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Equatorial zonal velocity'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_EqU, U1, U2, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Mediterranean salt tongue at depth=700'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 700, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Mediterranean salt tongue at depth=1000'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sspread, S1, S2, SLev, 1000, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Mediterranean water at lat=40°N'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 40, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Mediterranean water at lat=38°N'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_Med_Sdepth, S1, S2, SLev, 38, _extra = ex
-
+;
   cnt = cnt+1   &   blabla = 'Vertical Global mean T & S'
   IF doplot EQ cnt OR doplot EQ 0 THEN std_plot_GlobMeanTS, T1, T2, TLev, S1, S2, SLev, _extra = ex
-
+;
   IF n_elements(htmltxt) GT 1 THEN putfile, psdir+'std_plot_html_body.txt', htmltxt[1:*]
   
Index: /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh	(revision 4513)
+++ /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh	(revision 4514)
@@ -69,6 +69,7 @@
 VAR1_Ithick=sithic
 VAR1_Ifrac=siconc
-VAR1_Isnow=sndept
+VAR1_Isnow=snvolu
 VAR1_Isal=sisali
+VAR1_Iage=siages
 #
 #===================== EXP2 =====================
@@ -88,6 +89,7 @@
 VAR2_Ithick=sithic
 VAR2_Ifrac=siconc
-VAR2_Isnow=sndept
+VAR2_Isnow=snvolu
 VAR2_Isal=sisali
+VAR2_Iage=siages
 #
 ######################### Export Variables ###############################
@@ -113,9 +115,9 @@
 export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 
 export VAR1_U      VAR1_V
-export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal
+export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage
 #===================== EXP2 =====================
 export FILE2_T     FILE2_U    FILE2_V   FILE2_I
 export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 
 export VAR2_U      VAR2_V
-export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal
+export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage
 #
Index: anches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.500yfwb0_100y
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.500yfwb0_100y	(revision 4513)
+++ 	(revision )
@@ -1,135 +1,0 @@
-#!/bin/sh 
-#
-# AUTHOR - date
-# ===========
-# Sebastien Masson - 04/2011 - LOCEAN
-#
-# DESCRIPTION
-# ===========
-# define all varibles needed by std_plot.sh and std_plot.pro
-#
-# EXAMPLES
-# ========
-# $ . ./std_plot_vardef.sh
-#
-#
-#===================== User PATHS =====================
-#
-#idl_command=/Applications/itt/idl64/bin/idl
-idl_command=/Applications/itt/idl71/bin/idl
-#idl_command=/usr/local_linux/idl/idl_6.4/idl64/bin/idl
-#idl_command=idl71
-#
-PS_DIR=$( pwd )/500yfwb0_ps_plot
-PDF_DIR=$( pwd )/500yfwb0_pdf_plot
-HTML_DIR=$( pwd )/html_plot
-SAXO_DIR=/Users/sflod/SAXO_DIR
-#
-DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format
-DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data
-DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins)
-#
-#===================== Model GRID =====================
-#
-FILE_MESH_MASK=/Users/sflod/idl_PLOTS/MASK/ORL2PISV35_mesh_mask.nc # meshmask
-FILE_MASK_SUBDOMAIN=subbasins_orca21_nored.nc              # sub-bassin masks
-#
-#===================== DATA =====================
-#
-VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=potT_annual_mean.nc                         # PHC3
-VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc                    # PHC3
-VAR_SST=sst              ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds
-VAR_FLUX=qnet            ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux
-VAR_MLD=mld              ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth
-VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent North Emisphere
-VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent South Emisphere
-VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area North Emisphere
-VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area South Emisphere
-VAR_SNOW_NH=SnowDepth    ;   FILE_SNOW_ARC=Warren99_1954_1991_arctic_snowdepth.nc    # Snow arctic 
-VAR_SNOW_SH=sno_dept_mm  ;   FILE_SNOW_ANT=SMMR_SSMI_antarc_sno_dep_mm1979_2004.nc   # Snow antartic
-#
-# Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing
-#                    -> define VAR_GEOHEAT to a constant if geothermal heating is constant over the domain
-VAR_GEOHEAT=heatflow   ;   FILE_GEOHEAT=geothermal_heating.nc 
-#
-#===================== EXP1 =====================
-#
-#FILE1_T=l3fwb178_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file
-#FILE1_U=l3fwb178_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file
-#FILE1_V=l3fwb178_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file
-#FILE1_I=l3fwb178_24910101_25001231_1M_icemod.nc   # exp1 ice    input file
-FILE1_T=500yfwb0_20010101_21001231_1Y_grid_T.nc   # exp1 grid T input file
-FILE1_U=500yfwb0_20010101_21001231_1Y_grid_U.nc   # exp1 grid U input file
-FILE1_V=500yfwb0_20010101_21001231_1Y_grid_V.nc   # exp1 grid V input file
-FILE1_I=500yfwb0_20010101_21001231_1M_icemod.nc   # exp1 ice    input file
-VAR1_T=votemper         
-VAR1_S=vosaline         
-VAR1_QNET=sohflxdo
-VAR1_ERP=sowafldp      
-VAR1_EMP=sowaflup       
-VAR1_MLD=somxl010
-#VAR1_U=uoce
-#VAR1_U=vozocrtx
-VAR1_U=vozoeftr
-#VAR1_V=voce
-VAR1_V=vomeeftr
-#VAR1_V=vomecrty
-VAR1_Ithick=iicethic    
-VAR1_Ifrac=iiceconc
-VAR1_Isnow=isnowthi
-VAR1_Isal=iicesali
-#
-#===================== EXP2 =====================
-#
-#FILE2_T=l3fwb178_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file
-#FILE2_U=l3fwb178_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file
-#FILE2_V=l3fwb178_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file
-#FILE2_I=l3fwb178_24910101_25001231_1M_icemod.nc   # exp1 ice    input file
-FILE2_T=500yfwb0_20010101_21001231_1Y_grid_T.nc   # exp1 grid T input file
-FILE2_U=500yfwb0_20010101_21001231_1Y_grid_U.nc   # exp1 grid U input file
-FILE2_V=500yfwb0_20010101_21001231_1Y_grid_V.nc   # exp1 grid V input file
-FILE2_I=500yfwb0_20010101_21001231_1M_icemod.nc   # exp1 ice    input file
-VAR2_T=votemper
-VAR2_S=vosaline
-VAR2_QNET=sohflxdo
-VAR2_ERP=sowafldp
-VAR2_EMP=sowaflup
-VAR2_MLD=somxl010
-#VAR2_U=uoce
-VAR2_U=vozoeftr
-#VAR2_V=voce
-VAR2_V=vomeeftr
-VAR2_Ithick=iicethic
-VAR2_Ifrac=iiceconc
-VAR2_Isnow=isnowthi
-VAR2_Isal=iicesali
-#
-######################### Export Variables ###############################
-#
-#===================== User PATHS =====================
-export PS_DIR PDF_DIR HTML_DIR SAXO_DIR
-export DIR_DATA DIR_CLIMATO DIR_MASK
-#===================== Model GRID =====================
-export FILE_MESH_MASK FILE_MASK_SUBDOMAIN
-#===================== DATA =====================
-export FILE_TEMP_3D   VAR_TEMP_3D
-export FILE_SAL_3D    VAR_SAL_3D
-export FILE_SST       VAR_SST
-export FILE_GEOHEAT   VAR_GEOHEAT
-export FILE_FLUX      VAR_FLUX
-export FILE_MLD       VAR_MLD
-export FILE_ICE       VAR_ICE_EXT_NH  VAR_ICE_EXT_SH VAR_ICE_area_NH VAR_ICE_area_SH
-export FILE_SNOW_ARC  VAR_SNOW_NH
-export FILE_SNOW_ANT  VAR_SNOW_SH
-#===================== EXP1 =====================
-export READ_ONLY_FIRST_RECORD
-export FILE1_T     FILE1_U  FILE1_V   FILE1_I
-export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 
-export VAR1_U      VAR1_V
-export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal
-#===================== EXP2 =====================
-export FILE2_T     FILE2_U    FILE2_V   FILE2_I
-export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 
-export VAR2_U      VAR2_V
-export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal
-#
Index: anches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.l378sal6
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.l378sal6	(revision 4513)
+++ 	(revision )
@@ -1,135 +1,0 @@
-#!/bin/sh 
-#
-# AUTHOR - date
-# ===========
-# Sebastien Masson - 04/2011 - LOCEAN
-#
-# DESCRIPTION
-# ===========
-# define all varibles needed by std_plot.sh and std_plot.pro
-#
-# EXAMPLES
-# ========
-# $ . ./std_plot_vardef.sh
-#
-#
-#===================== User PATHS =====================
-#
-#idl_command=/Applications/itt/idl64/bin/idl
-idl_command=/Applications/itt/idl71/bin/idl
-#idl_command=/usr/local_linux/idl/idl_6.4/idl64/bin/idl
-#idl_command=idl71
-#
-PS_DIR=$( pwd )/l378sal6_ps_plot
-PDF_DIR=$( pwd )/l378sal6_pdf_plot
-HTML_DIR=$( pwd )/html_plot
-SAXO_DIR=/Users/sflod/SAXO_DIR
-#
-DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format
-DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data
-DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins)
-#
-#===================== Model GRID =====================
-#
-FILE_MESH_MASK=/Users/sflod/idl_PLOTS/MASK/ORL2PISV35_mesh_mask.nc # meshmask
-FILE_MASK_SUBDOMAIN=subbasins_orca21_nored.nc              # sub-bassin masks
-#
-#===================== DATA =====================
-#
-VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=potT_annual_mean.nc                         # PHC3
-VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc                    # PHC3
-VAR_SST=sst              ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds
-VAR_FLUX=qnet            ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux
-VAR_MLD=mld              ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth
-VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent North Emisphere
-VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent South Emisphere
-VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area North Emisphere
-VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area South Emisphere
-VAR_SNOW_NH=SnowDepth    ;   FILE_SNOW_ARC=Warren99_1954_1991_arctic_snowdepth.nc    # Snow arctic 
-VAR_SNOW_SH=sno_dept_mm  ;   FILE_SNOW_ANT=SMMR_SSMI_antarc_sno_dep_mm1979_2004.nc   # Snow antartic
-#
-# Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing
-#                    -> define VAR_GEOHEAT to a constant if geothermal heating is constant over the domain
-VAR_GEOHEAT=heatflow   ;   FILE_GEOHEAT=geothermal_heating.nc 
-#
-#===================== EXP1 =====================
-#
-#FILE1_T=l378sal6_20010101_21001231_1Y_grid_T.nc   # exp1 grid T input file
-#FILE1_U=l378sal6_20010101_21001231_1Y_grid_U.nc   # exp1 grid U input file
-#FILE1_V=l378sal6_20010101_21001231_1Y_grid_V.nc   # exp1 grid V input file
-#FILE1_I=l378sal6_20010101_21001231_1M_icemod.nc   # exp1 ice    input file
-FILE1_T=l378sal6_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file
-FILE1_U=l378sal6_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file
-FILE1_V=l378sal6_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file
-FILE1_I=l378sal6_24910101_25001231_1M_icemod.nc   # exp1 ice    input file
-VAR1_T=votemper         
-VAR1_S=vosaline         
-VAR1_QNET=sohflxdo
-VAR1_ERP=sowafldp      
-VAR1_EMP=sowaflup       
-VAR1_MLD=somxl010
-#VAR1_U=uoce
-#VAR1_U=vozocrtx
-VAR1_U=vozoeftr
-#VAR1_V=voce
-VAR1_V=vomeeftr
-#VAR1_V=vomecrty
-VAR1_Ithick=iicethic    
-VAR1_Ifrac=iiceconc
-VAR1_Isnow=isnowthi
-VAR1_Isal=iicesali
-#
-#===================== EXP2 =====================
-#
-#FILE2_T=l378sal6_20010101_21001231_1Y_grid_T.nc   # exp1 grid T input file
-#FILE2_U=l378sal6_20010101_21001231_1Y_grid_U.nc   # exp1 grid U input file
-#FILE2_V=l378sal6_20010101_21001231_1Y_grid_V.nc   # exp1 grid V input file
-#FILE2_I=l378sal6_20010101_21001231_1M_icemod.nc   # exp1 ice    input file
-FILE2_T=l378sal6_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file
-FILE2_U=l378sal6_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file
-FILE2_V=l378sal6_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file
-FILE2_I=l378sal6_24910101_25001231_1M_icemod.nc   # exp1 ice    input file
-VAR2_T=votemper
-VAR2_S=vosaline
-VAR2_QNET=sohflxdo
-VAR2_ERP=sowafldp
-VAR2_EMP=sowaflup
-VAR2_MLD=somxl010
-#VAR2_U=uoce
-VAR2_U=vozoeftr
-#VAR2_V=voce
-VAR2_V=vomeeftr
-VAR2_Ithick=iicethic
-VAR2_Ifrac=iiceconc
-VAR2_Isnow=isnowthi
-VAR2_Isal=iicesali
-#
-######################### Export Variables ###############################
-#
-#===================== User PATHS =====================
-export PS_DIR PDF_DIR HTML_DIR SAXO_DIR
-export DIR_DATA DIR_CLIMATO DIR_MASK
-#===================== Model GRID =====================
-export FILE_MESH_MASK FILE_MASK_SUBDOMAIN
-#===================== DATA =====================
-export FILE_TEMP_3D   VAR_TEMP_3D
-export FILE_SAL_3D    VAR_SAL_3D
-export FILE_SST       VAR_SST
-export FILE_GEOHEAT   VAR_GEOHEAT
-export FILE_FLUX      VAR_FLUX
-export FILE_MLD       VAR_MLD
-export FILE_ICE       VAR_ICE_EXT_NH  VAR_ICE_EXT_SH VAR_ICE_area_NH VAR_ICE_area_SH
-export FILE_SNOW_ARC  VAR_SNOW_NH
-export FILE_SNOW_ANT  VAR_SNOW_SH
-#===================== EXP1 =====================
-export READ_ONLY_FIRST_RECORD
-export FILE1_T     FILE1_U  FILE1_V   FILE1_I
-export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 
-export VAR1_U      VAR1_V
-export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal
-#===================== EXP2 =====================
-export FILE2_T     FILE2_U    FILE2_V   FILE2_I
-export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 
-export VAR2_U      VAR2_V
-export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal
-#
Index: anches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.l3fwb178-500yfwb0
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.l3fwb178-500yfwb0	(revision 4513)
+++ 	(revision )
@@ -1,131 +1,0 @@
-#!/bin/sh 
-#
-# AUTHOR - date
-# ===========
-# Sebastien Masson - 04/2011 - LOCEAN
-#
-# DESCRIPTION
-# ===========
-# define all varibles needed by std_plot.sh and std_plot.pro
-#
-# EXAMPLES
-# ========
-# $ . ./std_plot_vardef.sh
-#
-#
-#===================== User PATHS =====================
-#
-#idl_command=/Applications/itt/idl64/bin/idl
-idl_command=/Applications/itt/idl71/bin/idl
-#idl_command=/usr/local_linux/idl/idl_6.4/idl64/bin/idl
-#idl_command=idl71
-#
-PS_DIR=$( pwd )/l3fwb178-500yfwb0_ps_plot
-PDF_DIR=$( pwd )/l3fwb178-500yfwb0_pdf_plot
-HTML_DIR=$( pwd )/html_plot
-SAXO_DIR=/Users/sflod/SAXO_DIR
-#
-DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format
-DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data
-DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins)
-#
-#===================== Model GRID =====================
-#
-FILE_MESH_MASK=/Users/sflod/idl_PLOTS/MASK/ORL2PISV35_mesh_mask.nc # meshmask
-FILE_MASK_SUBDOMAIN=subbasins_orca21_nored.nc              # sub-bassin masks
-#
-#===================== DATA =====================
-#
-VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=potT_annual_mean.nc                         # PHC3
-VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc                    # PHC3
-VAR_SST=sst              ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds
-VAR_FLUX=qnet            ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux
-VAR_MLD=mld              ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth
-VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent North Emisphere
-VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent South Emisphere
-VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area North Emisphere
-VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area South Emisphere
-VAR_SNOW_NH=SnowDepth    ;   FILE_SNOW_ARC=Warren99_1954_1991_arctic_snowdepth.nc    # Snow arctic 
-VAR_SNOW_SH=sno_dept_mm  ;   FILE_SNOW_ANT=SMMR_SSMI_antarc_sno_dep_mm1979_2004.nc   # Snow antartic
-#
-# Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing
-#                    -> define VAR_GEOHEAT to a constant if geothermal heating is constant over the domain
-VAR_GEOHEAT=heatflow   ;   FILE_GEOHEAT=geothermal_heating.nc 
-#
-#===================== EXP1 =====================
-#
-FILE1_T=l3fwb178_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file
-FILE1_U=l3fwb178_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file
-FILE1_V=l3fwb178_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file
-FILE1_I=l3fwb178_24910101_25001231_1M_icemod.nc   # exp1 ice    input file
-#FILE1_T=l3fwb178_20010101_21001231_1Y_grid_T.nc   # exp1 grid T input file
-#FILE1_U=l3fwb178_20010101_21001231_1Y_grid_U.nc   # exp1 grid U input file
-#FILE1_V=l3fwb178_20010101_21001231_1Y_grid_V.nc   # exp1 grid V input file
-#FILE1_I=l3fwb178_20010101_21001231_1M_icemod.nc   # exp1 ice    input file
-VAR1_T=votemper         
-VAR1_S=vosaline         
-VAR1_QNET=sohflxdo
-VAR1_ERP=sowafldp      
-VAR1_EMP=sowaflup       
-VAR1_MLD=somxl010
-#VAR1_U=uoce
-#VAR1_U=vozocrtx
-VAR1_U=vozoeftr
-#VAR1_V=voce
-VAR1_V=vomeeftr
-#VAR1_V=vomecrty
-VAR1_Ithick=iicethic    
-VAR1_Ifrac=iiceconc
-VAR1_Isnow=isnowthi
-VAR1_Isal=iicesali
-#
-#===================== EXP2 =====================
-#
-FILE2_T=500yfwb0_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file
-FILE2_U=500yfwb0_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file
-FILE2_V=500yfwb0_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file
-FILE2_I=500yfwb0_24910101_25001231_1M_icemod.nc   # exp1 ice    input file
-VAR2_T=votemper
-VAR2_S=vosaline
-VAR2_QNET=sohflxdo
-VAR2_ERP=sowafldp
-VAR2_EMP=sowaflup
-VAR2_MLD=somxl010
-#VAR2_U=uoce
-VAR2_U=vozoeftr
-#VAR2_V=voce
-VAR2_V=vomeeftr
-VAR2_Ithick=iicethic
-VAR2_Ifrac=iiceconc
-VAR2_Isnow=isnowthi
-VAR2_Isal=iicesali
-#
-######################### Export Variables ###############################
-#
-#===================== User PATHS =====================
-export PS_DIR PDF_DIR HTML_DIR SAXO_DIR
-export DIR_DATA DIR_CLIMATO DIR_MASK
-#===================== Model GRID =====================
-export FILE_MESH_MASK FILE_MASK_SUBDOMAIN
-#===================== DATA =====================
-export FILE_TEMP_3D   VAR_TEMP_3D
-export FILE_SAL_3D    VAR_SAL_3D
-export FILE_SST       VAR_SST
-export FILE_GEOHEAT   VAR_GEOHEAT
-export FILE_FLUX      VAR_FLUX
-export FILE_MLD       VAR_MLD
-export FILE_ICE       VAR_ICE_EXT_NH  VAR_ICE_EXT_SH VAR_ICE_area_NH VAR_ICE_area_SH
-export FILE_SNOW_ARC  VAR_SNOW_NH
-export FILE_SNOW_ANT  VAR_SNOW_SH
-#===================== EXP1 =====================
-export READ_ONLY_FIRST_RECORD
-export FILE1_T     FILE1_U  FILE1_V   FILE1_I
-export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 
-export VAR1_U      VAR1_V
-export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal
-#===================== EXP2 =====================
-export FILE2_T     FILE2_U    FILE2_V   FILE2_I
-export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 
-export VAR2_U      VAR2_V
-export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal
-#
Index: /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.last_year_5D
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.last_year_5D	(revision 4513)
+++ /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.last_year_5D	(revision 4514)
@@ -26,7 +26,7 @@
 SAXO_DIR=/Users/sflod/SAXO_DIR
 #
-DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format
-DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data
-DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins)
+DIR_DATA=~/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format
+DIR_CLIMATO=~/idl_PLOTS/CLIMATOLOGIES    # path of climatological data
+DIR_MASK=~//idl_PLOTS/MASK  # path of mask files (ex: subbasins)
 #
 #===================== Model GRID =====================
@@ -55,46 +55,28 @@
 #===================== EXP1 =====================
 #
-#FILE1_T=l3fwb178_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file
-#FILE1_U=l3fwb178_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file
-#FILE1_V=l3fwb178_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file
-#FILE1_I=l3fwb178_24910101_25001231_1M_icemod.nc   # exp1 ice    input file
 FILE1_T=l3fwb178_25010101_25011231_5D_grid_T.nc   # exp1 grid T input file
 FILE1_U=l3fwb178_25010101_25011231_5D_grid_U.nc   # exp1 grid U input file
 FILE1_V=l3fwb178_25010101_25011231_5D_grid_V.nc   # exp1 grid V input file
 FILE1_I=l3fwb178_25010101_25011231_5D_icemod.nc   # exp1 ice    input file
-VAR1_T=votemper         
-VAR1_S=vosaline         
+VAR1_T=votemper
+VAR1_S=vosaline
 VAR1_QNET=sohflxdo
-VAR1_ERP=sowafldp      
-VAR1_EMP=sowaflup       
-#VAR1_MLD=somxl010
+VAR1_ERP=sowafldp
+VAR1_EMP=sowaflup
 VAR1_MLD=max5dmld
-#VAR1_U=uocetr_eff
-#VAR1_U=uoce
-#VAR1_U=vozocrtx
 VAR1_U=vozoeftr
-#VAR1_V=vocetr_eff
-#VAR1_V=voce
 VAR1_V=vomeeftr
-#VAR1_V=vomecrty
-VAR1_Ithick=iicethic    
+VAR1_Ithick=iicethic
 VAR1_Ifrac=iiceconc
-VAR1_Isnow=isnowthi
+VAR1_Isnow=isnowvol
 VAR1_Isal=iicesali
+VAR1_Iage=iiceages
 #
 #===================== EXP2 =====================
 #
-#FILE2_T=l3fwb178_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file
-#FILE2_U=l3fwb178_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file
-#FILE2_V=l3fwb178_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file
-#FILE2_I=l3fwb178_24910101_25001231_1M_icemod.nc   # exp1 ice    input file
 FILE2_T=l3fwb178_25010101_25011231_5D_grid_T.nc   # exp1 grid T input file
 FILE2_U=l3fwb178_25010101_25011231_5D_grid_U.nc   # exp1 grid U input file
 FILE2_V=l3fwb178_25010101_25011231_5D_grid_V.nc   # exp1 grid V input file
 FILE2_I=l3fwb178_25010101_25011231_5D_icemod.nc   # exp1 ice    input file
-#FILE2_T=l378sal6_25010101_25011231_5D_grid_T.nc   # exp1 grid T input file
-#FILE2_U=l378sal6_25010101_25011231_5D_grid_U.nc   # exp1 grid U input file
-#FILE2_V=l378sal6_25010101_25011231_5D_grid_V.nc   # exp1 grid V input file
-#FILE2_I=l378sal6_25010101_25011231_5D_icemod.nc   # exp1 ice    input file
 VAR2_T=votemper
 VAR2_S=vosaline
@@ -102,11 +84,6 @@
 VAR2_ERP=sowafldp
 VAR2_EMP=sowaflup
-#VAR2_MLD=somxl010
 VAR2_MLD=max5dmld
-#VAR2_U=uoce_eff
-#VAR2_U=uoce
 VAR2_U=vozoeftr
-#VAR2_V=vocetr_eff
-#VAR2_V=voce
 VAR2_V=vomeeftr
 VAR2_Ithick=iicethic
@@ -114,4 +91,5 @@
 VAR2_Isnow=isnowthi
 VAR2_Isal=iicesali
+VAR2_Iage=iiceages
 #
 ######################### Export Variables ###############################
@@ -137,9 +115,9 @@
 export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 
 export VAR1_U      VAR1_V
-export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal
+export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage
 #===================== EXP2 =====================
 export FILE2_T     FILE2_U    FILE2_V   FILE2_I
 export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 
 export VAR2_U      VAR2_V
-export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal
+export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage
 #
Index: /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.new_names
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.new_names	(revision 4513)
+++ /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.new_names	(revision 4514)
@@ -69,6 +69,7 @@
 VAR1_Ithick=sithic
 VAR1_Ifrac=siconc
-VAR1_Isnow=sndept
+VAR1_Isnow=snvolu
 VAR1_Isal=sisali
+VAR1_Iage=siages
 #
 #===================== EXP2 =====================
@@ -88,6 +89,7 @@
 VAR2_Ithick=sithic
 VAR2_Ifrac=siconc
-VAR2_Isnow=sndept
+VAR2_Isnow=snvolu
 VAR2_Isal=sisali
+VAR2_Iage=siages
 #
 ######################### Export Variables ###############################
@@ -113,9 +115,9 @@
 export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 
 export VAR1_U      VAR1_V
-export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal
+export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage
 #===================== EXP2 =====================
 export FILE2_T     FILE2_U    FILE2_V   FILE2_I
 export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 
 export VAR2_U      VAR2_V
-export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal
+export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage
 #
Index: /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.old_names
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.old_names	(revision 4514)
+++ /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.old_names	(revision 4514)
@@ -0,0 +1,125 @@
+#!/bin/sh 
+#
+# AUTHOR - date
+# ===========
+# Sebastien Masson - 04/2011 - LOCEAN
+#
+# DESCRIPTION
+# ===========
+# define all varibles needed by std_plot.sh and std_plot.pro
+#
+# EXAMPLES
+# ========
+# $ . ./std_plot_vardef.sh
+#
+#
+#===================== User PATHS =====================
+#
+#idl_command=/Applications/itt/idl64/bin/idl
+idl_command=/Applications/itt/idl71/bin/idl
+#idl_command=/usr/local_linux/idl/idl_6.4/idl64/bin/idl
+#idl_command=idl71
+#
+PS_DIR=$( pwd )/500yfwb0_ps_plot_500y
+PDF_DIR=$( pwd )/500yfwb0_pdf_plot_500y
+HTML_DIR=$( pwd )/html_plot
+SAXO_DIR=/Users/sflod/SAXO_DIR
+#
+DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format
+DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data
+DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins)
+#
+#===================== Model GRID =====================
+#
+FILE_MESH_MASK=/Users/sflod/idl_PLOTS/MASK/ORL2PISV35_mesh_mask.nc # meshmask
+FILE_MASK_SUBDOMAIN=subbasins_orca21_nored.nc              # sub-bassin masks
+#
+#===================== DATA =====================
+#
+VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=potT_annual_mean.nc                         # PHC3
+VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc                    # PHC3
+VAR_SST=sst              ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds
+VAR_FLUX=qnet            ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux
+VAR_MLD=mld              ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth
+VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent North Emisphere
+VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent South Emisphere
+VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area North Emisphere
+VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area South Emisphere
+VAR_SNOW_NH=SnowDepth    ;   FILE_SNOW_ARC=Warren99_1954_1991_arctic_snowdepth.nc    # Snow arctic 
+VAR_SNOW_SH=sno_dept_mm  ;   FILE_SNOW_ANT=SMMR_SSMI_antarc_sno_dep_mm1979_2004.nc   # Snow antartic
+#
+# Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing
+#                    -> define VAR_GEOHEAT to a constant if geothermal heating is constant over the domain
+VAR_GEOHEAT=heatflow   ;   FILE_GEOHEAT=geothermal_heating.nc 
+#
+#===================== EXP1 =====================
+#
+FILE1_T=500yfwb0_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file
+FILE1_U=500yfwb0_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file
+FILE1_V=500yfwb0_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file
+FILE1_I=500yfwb0_24910101_25001231_1M_icemod.nc   # exp1 ice    input file
+VAR1_T=votemper
+VAR1_S=vosaline
+VAR1_QNET=sohflxdo
+VAR1_ERP=sowafldp
+VAR1_EMP=sowaflup
+VAR1_MLD=somxl010
+VAR1_U=vozoeftr
+VAR1_V=vomeeftr
+VAR1_Ithick=iicethic
+VAR1_Ifrac=iiceconc
+VAR1_Isnow=isnowvol
+VAR1_Isal=iicesali
+VAR1_Iage=iiceages
+#
+#===================== EXP2 =====================
+#
+FILE2_T=500yfwb0_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file
+FILE2_U=500yfwb0_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file
+FILE2_V=500yfwb0_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file
+FILE2_I=500yfwb0_24910101_25001231_1M_icemod.nc   # exp1 ice    input file
+VAR2_T=votemper
+VAR2_S=vosaline
+VAR2_QNET=sohflxdo
+VAR2_ERP=sowafldp
+VAR2_EMP=sowaflup
+VAR2_MLD=somxl010
+#VAR2_U=uoce
+VAR2_U=vozoeftr
+#VAR2_V=voce
+VAR2_V=vomeeftr
+VAR2_Ithick=iicethic
+VAR2_Ifrac=iiceconc
+VAR2_Isnow=isnowthi
+VAR2_Isal=iicesali
+VAR2_Iage=iiceages
+#
+######################### Export Variables ###############################
+#
+#===================== User PATHS =====================
+export PS_DIR PDF_DIR HTML_DIR SAXO_DIR
+export DIR_DATA DIR_CLIMATO DIR_MASK
+#===================== Model GRID =====================
+export FILE_MESH_MASK FILE_MASK_SUBDOMAIN
+#===================== DATA =====================
+export FILE_TEMP_3D   VAR_TEMP_3D
+export FILE_SAL_3D    VAR_SAL_3D
+export FILE_SST       VAR_SST
+export FILE_GEOHEAT   VAR_GEOHEAT
+export FILE_FLUX      VAR_FLUX
+export FILE_MLD       VAR_MLD
+export FILE_ICE       VAR_ICE_EXT_NH  VAR_ICE_EXT_SH VAR_ICE_area_NH VAR_ICE_area_SH
+export FILE_SNOW_ARC  VAR_SNOW_NH
+export FILE_SNOW_ANT  VAR_SNOW_SH
+#===================== EXP1 =====================
+export READ_ONLY_FIRST_RECORD
+export FILE1_T     FILE1_U  FILE1_V   FILE1_I
+export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 
+export VAR1_U      VAR1_V
+export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage
+#===================== EXP2 =====================
+export FILE2_T     FILE2_U    FILE2_V   FILE2_I
+export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 
+export VAR2_U      VAR2_V
+export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage
+#
Index: anches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.qsriceok
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.qsriceok	(revision 4513)
+++ 	(revision )
@@ -1,135 +1,0 @@
-#!/bin/sh 
-#
-# AUTHOR - date
-# ===========
-# Sebastien Masson - 04/2011 - LOCEAN
-#
-# DESCRIPTION
-# ===========
-# define all varibles needed by std_plot.sh and std_plot.pro
-#
-# EXAMPLES
-# ========
-# $ . ./std_plot_vardef.sh
-#
-#
-#===================== User PATHS =====================
-#
-#idl_command=/Applications/itt/idl64/bin/idl
-idl_command=/Applications/itt/idl71/bin/idl
-#idl_command=/usr/local_linux/idl/idl_6.4/idl64/bin/idl
-#idl_command=idl71
-#
-PS_DIR=$( pwd )/qsriceok_ps_plot
-PDF_DIR=$( pwd )/qsriceok_pdf_plot
-HTML_DIR=$( pwd )/html_plot
-SAXO_DIR=/Users/sflod/SAXO_DIR
-#
-DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format
-DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data
-DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins)
-#
-#===================== Model GRID =====================
-#
-FILE_MESH_MASK=/Users/sflod/idl_PLOTS/MASK/ORL2PISV35_mesh_mask.nc # meshmask
-FILE_MASK_SUBDOMAIN=subbasins_orca21_nored.nc              # sub-bassin masks
-#
-#===================== DATA =====================
-#
-VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=potT_annual_mean.nc                         # PHC3
-VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc                    # PHC3
-VAR_SST=sst              ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds
-VAR_FLUX=qnet            ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux
-VAR_MLD=mld              ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth
-VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent North Emisphere
-VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent South Emisphere
-VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area North Emisphere
-VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area South Emisphere
-VAR_SNOW_NH=SnowDepth    ;   FILE_SNOW_ARC=Warren99_1954_1991_arctic_snowdepth.nc    # Snow arctic 
-VAR_SNOW_SH=sno_dept_mm  ;   FILE_SNOW_ANT=SMMR_SSMI_antarc_sno_dep_mm1979_2004.nc   # Snow antartic
-#
-# Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing
-#                    -> define VAR_GEOHEAT to a constant if geothermal heating is constant over the domain
-VAR_GEOHEAT=heatflow   ;   FILE_GEOHEAT=geothermal_heating.nc 
-#
-#===================== EXP1 =====================
-#
-#FILE1_T=qsriceok_20010101_21001231_1Y_grid_T.nc   # exp1 grid T input file
-#FILE1_U=qsriceok_20010101_21001231_1Y_grid_U.nc   # exp1 grid U input file
-#FILE1_V=qsriceok_20010101_21001231_1Y_grid_V.nc   # exp1 grid V input file
-#FILE1_I=qsriceok_20010101_21001231_1M_icemod.nc   # exp1 ice    input file
-FILE1_T=qsriceok_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file
-FILE1_U=qsriceok_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file
-FILE1_V=qsriceok_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file
-FILE1_I=qsriceok_24910101_25001231_1M_icemod.nc   # exp1 ice    input file
-VAR1_T=votemper         
-VAR1_S=vosaline         
-VAR1_QNET=sohflxdo
-VAR1_ERP=sowafldp      
-VAR1_EMP=sowaflup       
-VAR1_MLD=somxl010
-#VAR1_U=uoce
-#VAR1_U=vozocrtx
-VAR1_U=vozoeftr
-#VAR1_V=voce
-VAR1_V=vomeeftr
-#VAR1_V=vomecrty
-VAR1_Ithick=iicethic    
-VAR1_Ifrac=iiceconc
-VAR1_Isnow=isnowthi
-VAR1_Isal=iicesali
-#
-#===================== EXP2 =====================
-#
-#FILE2_T=qsriceok_20010101_21001231_1Y_grid_T.nc   # exp1 grid T input file
-#FILE2_U=qsriceok_20010101_21001231_1Y_grid_U.nc   # exp1 grid U input file
-#FILE2_V=qsriceok_20010101_21001231_1Y_grid_V.nc   # exp1 grid V input file
-#FILE2_I=qsriceok_20010101_21001231_1M_icemod.nc   # exp1 ice    input file
-FILE2_T=qsriceok_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file
-FILE2_U=qsriceok_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file
-FILE2_V=qsriceok_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file
-FILE2_I=qsriceok_24910101_25001231_1M_icemod.nc   # exp1 ice    input file
-VAR2_T=votemper
-VAR2_S=vosaline
-VAR2_QNET=sohflxdo
-VAR2_ERP=sowafldp
-VAR2_EMP=sowaflup
-VAR2_MLD=somxl010
-#VAR2_U=uoce
-VAR2_U=vozoeftr
-#VAR2_V=voce
-VAR2_V=vomeeftr
-VAR2_Ithick=iicethic
-VAR2_Ifrac=iiceconc
-VAR2_Isnow=isnowthi
-VAR2_Isal=iicesali
-#
-######################### Export Variables ###############################
-#
-#===================== User PATHS =====================
-export PS_DIR PDF_DIR HTML_DIR SAXO_DIR
-export DIR_DATA DIR_CLIMATO DIR_MASK
-#===================== Model GRID =====================
-export FILE_MESH_MASK FILE_MASK_SUBDOMAIN
-#===================== DATA =====================
-export FILE_TEMP_3D   VAR_TEMP_3D
-export FILE_SAL_3D    VAR_SAL_3D
-export FILE_SST       VAR_SST
-export FILE_GEOHEAT   VAR_GEOHEAT
-export FILE_FLUX      VAR_FLUX
-export FILE_MLD       VAR_MLD
-export FILE_ICE       VAR_ICE_EXT_NH  VAR_ICE_EXT_SH VAR_ICE_area_NH VAR_ICE_area_SH
-export FILE_SNOW_ARC  VAR_SNOW_NH
-export FILE_SNOW_ANT  VAR_SNOW_SH
-#===================== EXP1 =====================
-export READ_ONLY_FIRST_RECORD
-export FILE1_T     FILE1_U  FILE1_V   FILE1_I
-export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 
-export VAR1_U      VAR1_V
-export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal
-#===================== EXP2 =====================
-export FILE2_T     FILE2_U    FILE2_V   FILE2_I
-export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 
-export VAR2_U      VAR2_V
-export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal
-#
Index: anches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.sbcmodMV-l334plus
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh.sbcmodMV-l334plus	(revision 4513)
+++ 	(revision )
@@ -1,121 +1,0 @@
-#!/bin/sh 
-#
-# AUTHOR - date
-# ===========
-# Sebastien Masson - 04/2011 - LOCEAN
-#
-# DESCRIPTION
-# ===========
-# define all varibles needed by std_plot.sh and std_plot.pro
-#
-# EXAMPLES
-# ========
-# $ . ./std_plot_vardef.sh
-#
-#
-#===================== User PATHS =====================
-#
-#idl_command=/Applications/itt/idl64/bin/idl
-idl_command=/Applications/itt/idl71/bin/idl
-#idl_command=/usr/local_linux/idl/idl_6.4/idl64/bin/idl
-#idl_command=idl71
-#
-PS_DIR=$( pwd )/sbcmodMV-l334plus_ps_plot
-PDF_DIR=$( pwd )/sbcmodMV-l334plus_pdf_plot
-HTML_DIR=$( pwd )/html_plot
-SAXO_DIR=/Users/sflod/SAXO_DIR
-#
-DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format
-DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data
-DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins)
-#
-#===================== Model GRID =====================
-#
-FILE_MESH_MASK=/Users/sflod/idl_PLOTS/MASK/ORL2PISV35_mesh_mask.nc # meshmask
-FILE_MASK_SUBDOMAIN=subbasins_orca21_nored.nc              # sub-bassin masks
-#
-#===================== DATA =====================
-#
-VAR_TEMP_3D=votemper     ;   FILE_TEMP_3D=potT_annual_mean.nc                         # PHC3
-VAR_SAL_3D=vosaline      ;   FILE_SAL_3D=Salt_1y_corr_PHC3WOA09.nc                    # PHC3
-VAR_SST=sst              ;   FILE_SST=NewREY_ORCA2_1991_2000_1y.nc                   # Reynolds
-VAR_FLUX=qnet            ;   FILE_FLUX=OAFlux_1my_01_12_1984_2004_orca2_qnet.nc      # flux
-VAR_MLD=mld              ;   FILE_MLD=mld_DR003_c1m_ORCA2_1y.nc                      # Mixed layer depth
-VAR_ICE_EXT_NH=extt_NH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent North Emisphere
-VAR_ICE_EXT_SH=extt_SH   ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Extent South Emisphere
-VAR_ICE_area_NH=area_NH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area North Emisphere
-VAR_ICE_area_SH=area_SH  ;   FILE_ICE=sea_ice_index_2000.nc                               # Ice Area South Emisphere
-VAR_SNOW_NH=SnowDepth    ;   FILE_SNOW_ARC=Warren99_1954_1991_arctic_snowdepth.nc    # Snow arctic 
-VAR_SNOW_SH=sno_dept_mm  ;   FILE_SNOW_ANT=SMMR_SSMI_antarc_sno_dep_mm1979_2004.nc   # Snow antartic
-#
-# Geothermal heating -> define FILE_GEOHEAT to 'NO' if there is not such forcing
-#                    -> define VAR_GEOHEAT to a constant if geothermal heating is constant over the domain
-VAR_GEOHEAT=heatflow   ;   FILE_GEOHEAT=geothermal_heating.nc 
-#
-#===================== EXP1 =====================
-#
-FILE1_T=sbcmodMV_20410101_20501231_1Y_grid_T.nc   # exp1 grid T input file
-FILE1_U=sbcmodMV_20410101_20501231_1Y_grid_U.nc   # exp1 grid U input file
-FILE1_V=sbcmodMV_20410101_20501231_1Y_grid_V.nc   # exp1 grid V input file
-FILE1_I=sbcmodMV_20410101_20501231_1M_icemod.nc   # exp1 ice    input file
-VAR1_T=thetao
-VAR1_S=so        
-VAR1_QNET=tohfls
-VAR1_ERP=sowafldp      
-VAR1_EMP=wfo
-VAR1_MLD=mldr10_1
-VAR1_U=uocetr_eff
-VAR1_V=vocetr_eff
-VAR1_Ithick=sithic
-VAR1_Ifrac=siconc
-VAR1_Isnow=sndept
-VAR1_Isal=sisali
-#
-#===================== EXP2 =====================
-#
-FILE2_T=l334plus_20410101_20501231_1Y_grid_T.nc   # exp1 grid T input file
-FILE2_U=l334plus_20410101_20501231_1Y_grid_U.nc   # exp1 grid U input file
-FILE2_V=l334plus_20410101_20501231_1Y_grid_V.nc   # exp1 grid V input file
-FILE2_I=l334plus_20410101_20501231_1M_icemod.nc   # exp1 ice    input file
-VAR2_T=thetao
-VAR2_S=so
-VAR2_QNET=tohfls
-VAR2_ERP=sowafldp
-VAR2_EMP=wfo
-VAR2_MLD=mldr10_1
-VAR2_U=uocetr_eff
-VAR2_V=vocetr_eff
-VAR2_Ithick=sithic
-VAR2_Ifrac=siconc
-VAR2_Isnow=sndept
-VAR2_Isal=sisali
-#
-######################### Export Variables ###############################
-#
-#===================== User PATHS =====================
-export PS_DIR PDF_DIR HTML_DIR SAXO_DIR
-export DIR_DATA DIR_CLIMATO DIR_MASK
-#===================== Model GRID =====================
-export FILE_MESH_MASK FILE_MASK_SUBDOMAIN
-#===================== DATA =====================
-export FILE_TEMP_3D   VAR_TEMP_3D
-export FILE_SAL_3D    VAR_SAL_3D
-export FILE_SST       VAR_SST
-export FILE_GEOHEAT   VAR_GEOHEAT
-export FILE_FLUX      VAR_FLUX
-export FILE_MLD       VAR_MLD
-export FILE_ICE       VAR_ICE_EXT_NH  VAR_ICE_EXT_SH VAR_ICE_area_NH VAR_ICE_area_SH
-export FILE_SNOW_ARC  VAR_SNOW_NH
-export FILE_SNOW_ANT  VAR_SNOW_SH
-#===================== EXP1 =====================
-export READ_ONLY_FIRST_RECORD
-export FILE1_T     FILE1_U  FILE1_V   FILE1_I
-export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 
-export VAR1_U      VAR1_V
-export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal
-#===================== EXP2 =====================
-export FILE2_T     FILE2_U    FILE2_V   FILE2_I
-export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 
-export VAR2_U      VAR2_V
-export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal
-#
Index: /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example1
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example1	(revision 4513)
+++ /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example1	(revision 4514)
@@ -26,7 +26,7 @@
 SAXO_DIR=/Users/sflod/SAXO_DIR
 #
-DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada      # path of data in NetCDF format
-DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data
-DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins)
+DIR_DATA=~/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format
+DIR_CLIMATO=~/idl_PLOTS/CLIMATOLOGIES    # path of climatological data
+DIR_MASK=~//idl_PLOTS/MASK  # path of mask files (ex: subbasins)
 #
 #===================== Model GRID =====================
@@ -55,20 +55,21 @@
 #===================== EXP1 =====================
 #
-FILE1_T=tr33beta_00910101_01001231_1Y_grid_T.nc   # exp1 grid T input file
-FILE1_U=tr33beta_00910101_01001231_1Y_grid_U.nc   # exp1 grid U input file
-FILE1_V=tr33beta_00910101_01001231_1Y_grid_V.nc   # exp1 grid V input file
-FILE1_I=tr33beta_00910101_01001231_1M_icemod.nc   # exp1 ice    input file
-VAR1_T=thetao
-VAR1_S=so
-VAR1_QNET=qt
-VAR1_ERP=wfcorr
-VAR1_EMP=wfo
-VAR1_MLD=mldr10_1
-VAR1_U=uocetr_eff
-VAR1_V=vocetr_eff
-VAR1_Ithick=sit
-VAR1_Ifrac=sic
-VAR1_Isnow=isnowthi
+FILE1_T=500yfwb0_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file
+FILE1_U=500yfwb0_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file
+FILE1_V=500yfwb0_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file
+FILE1_I=500yfwb0_24910101_25001231_1M_icemod.nc   # exp1 ice    input file
+VAR1_T=votemper
+VAR1_S=vosaline
+VAR1_QNET=sohflxdo
+VAR1_ERP=sowafldp
+VAR1_EMP=sowaflup
+VAR1_MLD=somxl010
+VAR1_U=vozoeftr
+VAR1_V=vomeeftr
+VAR1_Ithick=iicethic
+VAR1_Ifrac=iiceconc
+VAR1_Isnow=isnowvol
 VAR1_Isal=iicesali
+VAR1_Iage=iiceages
 #
 #===================== EXP2 =====================
@@ -90,4 +91,5 @@
 VAR2_Isnow=isnowthi
 VAR2_Isal=iicesali
+VAR2_Iage=iiceages
 #
 ######################### Export Variables ###############################
@@ -113,9 +115,9 @@
 export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 
 export VAR1_U      VAR1_V
-export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal
+export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage
 #===================== EXP2 =====================
 export FILE2_T     FILE2_U    FILE2_V   FILE2_I
 export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 
 export VAR2_U      VAR2_V
-export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal
+export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage
 #
Index: /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example2
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example2	(revision 4513)
+++ /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_plot_vardef.sh_example2	(revision 4514)
@@ -21,12 +21,12 @@
 #idl_command=idl71
 #
-PS_DIR=$( pwd )/500yfwb0_ps_plot_500y
-PDF_DIR=$( pwd )/500yfwb0_pdf_plot_500y
+PS_DIR=$( pwd )/ps_plot
+PDF_DIR=$( pwd )/pdf_plot
 HTML_DIR=$( pwd )/html_plot
 SAXO_DIR=/Users/sflod/SAXO_DIR
 #
-DIR_DATA=/Users/sflod/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format
-DIR_CLIMATO=/Users/sflod/idl_PLOTS/CLIMATOLOGIES    # path of climatological data
-DIR_MASK=/Users/sflod/idl_PLOTS/MASK  # path of mask files (ex: subbasins)
+DIR_DATA=~/idl_PLOTS/DATA_STORE/RUN_CLIMATO/lim3_ada     # path of data in NetCDF format
+DIR_CLIMATO=~/idl_PLOTS/CLIMATOLOGIES    # path of climatological data
+DIR_MASK=~//idl_PLOTS/MASK  # path of mask files (ex: subbasins)
 #
 #===================== Model GRID =====================
@@ -55,39 +55,41 @@
 #===================== EXP1 =====================
 #
-FILE1_T=tag33_20910101_21001231_1Y_grid_T.nc   # exp1 grid T input file
-FILE1_U=tag33_20910101_21001231_1Y_grid_U.nc   # exp1 grid U input file
-FILE1_V=tag33_20910101_21001231_1Y_grid_V.nc   # exp1 grid V input file
-FILE1_I=tag33_20910101_21001231_1M_icemod.nc   # exp1 ice    input file
-VAR1_T=thetao
-VAR1_S=so
-VAR1_QNET=qt
-VAR1_ERP=wfcorr
-VAR1_EMP=wfo
-VAR1_MLD=mldr10_1
-VAR1_U=uocetr_eff
-VAR1_V=vocetr_eff
-VAR1_Ithick=sit
-VAR1_Ifrac=sic
-VAR1_Isnow=isnowthi
+FILE1_T=500yfwb0_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file
+FILE1_U=500yfwb0_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file
+FILE1_V=500yfwb0_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file
+FILE1_I=500yfwb0_24910101_25001231_1M_icemod.nc   # exp1 ice    input file
+VAR1_T=votemper
+VAR1_S=vosaline
+VAR1_QNET=sohflxdo
+VAR1_ERP=sowafldp
+VAR1_EMP=sowaflup
+VAR1_MLD=somxl010
+VAR1_U=vozoeftr
+VAR1_V=vomeeftr
+VAR1_Ithick=iicethic
+VAR1_Ifrac=iiceconc
+VAR1_Isnow=isnowvol
 VAR1_Isal=iicesali
+VAR1_Iage=iiceages
 #
 #===================== EXP2 =====================
 #
-FILE2_T=tmx_2000_00910101_01001231_1Y_grid_T.nc   # exp1 grid T input file
-FILE2_U=tmx_2000_00910101_01001231_1Y_grid_U.nc   # exp1 grid U input file
-FILE2_V=tmx_2000_00910101_01001231_1Y_grid_V.nc   # exp1 grid V input file
-FILE2_I=tmx_2000_00910101_01001231_1M_icemod.nc   # exp1 ice    input file
-VAR2_T=thetao
-VAR2_S=so
-VAR2_QNET=qns+qsr
-VAR2_ERP=wfcorr
-VAR2_EMP=wfo
-VAR2_MLD=mldr10_1
-VAR2_U=uoce_eff
-VAR2_V=voce_eff
-VAR2_Ithick=sit
-VAR2_Ifrac=sic
-VAR2_Isnow=isnowthi
-VAR2_Isal=iicesali
+FILE1_T=500yfwb0_24910101_25001231_1Y_grid_T.nc   # exp1 grid T input file
+FILE1_U=500yfwb0_24910101_25001231_1Y_grid_U.nc   # exp1 grid U input file
+FILE1_V=500yfwb0_24910101_25001231_1Y_grid_V.nc   # exp1 grid V input file
+FILE1_I=500yfwb0_24910101_25001231_1M_icemod.nc   # exp1 ice    input file
+VAR1_T=votemper
+VAR1_S=vosaline
+VAR1_QNET=sohflxdo
+VAR1_ERP=sowafldp
+VAR1_EMP=sowaflup
+VAR1_MLD=somxl010
+VAR1_U=vozoeftr
+VAR1_V=vomeeftr
+VAR1_Ithick=iicethic
+VAR1_Ifrac=iiceconc
+VAR1_Isnow=isnowvol
+VAR1_Isal=iicesali
+VAR1_Iage=iiceages
 #
 ######################### Export Variables ###############################
@@ -113,9 +115,9 @@
 export VAR1_T      VAR1_S   VAR1_QNET VAR1_ERP VAR1_EMP VAR1_MLD 
 export VAR1_U      VAR1_V
-export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal
+export VAR1_Ithick VAR1_Ifrac VAR1_Isnow VAR1_Isal VAR1_Iage
 #===================== EXP2 =====================
 export FILE2_T     FILE2_U    FILE2_V   FILE2_I
 export VAR2_T      VAR2_S     VAR2_QNET VAR2_ERP VAR2_EMP VAR2_MLD 
 export VAR2_U      VAR2_V
-export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal
+export VAR2_Ithick VAR2_Ifrac VAR2_Isnow VAR2_Isal VAR2_Iage
 #
Index: /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_FRAM.pro
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_FRAM.pro	(revision 4513)
+++ /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_FRAM.pro	(revision 4514)
@@ -49,5 +49,6 @@
  ;
  ;OBSERVATIONS : mean seasonal cycle/month
-  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]
+ ; 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]
+  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]
   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]
  ;
@@ -80,9 +81,8 @@
   ICE_vol_export = reform(ICE_vol_export, 12, nyr)
   ICE_vol_export = total(ICE_vol_export,2)/nyr
-  ; ICE_vol_export = {arr:ICE_vol_export * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'}      ; annual mean
-  ICE_vol_export = {arr:ICE_vol_export * 1.e-12 * 86400 * 30 , unit : '10^3 Km^3/month'}      ; monthly mean
-
-  ;
-  ;title = 'Northern Hemisphere'+'!C'+prefix+' '+d1_d2+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)'
+  ; ICE_vol_export = {arr:ICE_vol_export * 1.e-9 * 86400 * 365 , unit : '10^3 Km^3/year'}      ; annual mean
+  ICE_vol_export = {arr:ICE_vol_export * 1.e-9 * 86400 * 30 , unit : '10^3 Km^3/month'}      ; monthly mean
+
+  ;
   title = 'Fram Strait Areal Export'+'!C'+prefix+' '+d1_d2
   jpt=12
@@ -95,10 +95,10 @@
   tot_area_expo = total(ICE_area_export.arr)
   ;
-  xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 0.851 million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2
-  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
-  xyouts, julday(7,15,1900), 0.42, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2
-  ;title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Global Annual Mean Ice Volume (Black SOLID simulation)'
+  xyouts, julday(5,15,1900), 0.49, 'Tot. Annual Export OBS = 0.851 million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2
+  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
+  xyouts, julday(7,15,1900), 0.41, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2
+  ;
   title = 'Fram Strait Volume Export'+'!C'+prefix+' '+d1_d2
-  pltt, ICE_vol_export, 't',  /REMPLI, MIN = 0., MAX = .6 , /NOERASE, XGRIDSTYLE = 1 $
+  pltt, ICE_vol_export, 't',  /REMPLI, MIN = 60., MAX = 500. , /NOERASE, XGRIDSTYLE = 1 $
        , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex
   pltt, vol_obs, 't', /REMPLI, /NOERASE, psym = 4, THICK = 4  $    ; light blue
@@ -107,8 +107,8 @@
   tot_vol_expo = total(ICE_vol_export.arr)
   tot_vol_expo_Sv = tot_vol_expo * 1.e06 * 1/86400 * 1/365 ; annual mean in Sverdrup
-  xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 2124 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2
-  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
-  xyouts, julday(9,15,1900), 0.44, 'in Sv = '+strtrim(tot_vol_expo_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2
-  xyouts, julday(7,15,1900), 0.38, 'Data from Kwok et al.(2004), 1992-1998 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2
+  xyouts, julday(5,15,1900), 490, 'Tot. Annual Export OBS = 2124 10^3 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2
+  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
+  xyouts, julday(9,15,1900), 440, 'in Sv = '+strtrim(tot_vol_expo_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2
+  xyouts, julday(7,15,1900), 360, 'Data from Kwok et al.(2004), 1992-1998 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2
 
   htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
@@ -147,5 +147,5 @@
     ICE_vol_export_2 = total(ICE_vol_export_2,2)/nyr
     ; ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-12 * 86400 * 365 , unit : '10^3 Km^3/year'}     ; annual mean
-    ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-12 * 86400 * 30 , unit : '10^3 Km^3/month'}       ; monthly mean
+    ICE_vol_export_2 = {arr:ICE_vol_export_2 * 1.e-9 * 86400 * 30 , unit : '10^3 Km^3/month'}       ; monthly mean
     ;
     ;
@@ -165,10 +165,11 @@
     tot_area_expo_2 = total(ICE_area_export_2.arr)
     ;
-    xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 0.851 million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2
-    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
-    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
-    xyouts, julday(7,15,1900), 0.42, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2
+    xyouts, julday(5,15,1900), 0.49, 'Tot. Annual Export OBS = 0.851 million Km2', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2
+    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
+    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
+    xyouts, julday(7,15,1900), 0.41, 'Data from Kwok et al.(2004), 1979-2002 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2
+    ;
     title = 'Fram Strait Volume Export'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2
-    pltt, ICE_vol_export, 't', /REMPLI,  MIN = 0., MAX = .6, /NOERASE, XGRIDSTYLE = 1 $
+    pltt, ICE_vol_export, 't', /REMPLI,  MIN = 60., MAX = 500., /NOERASE, XGRIDSTYLE = 1 $
         , small = [1, 2, 2], YTITLE = '10^3 Km^3/month', TITLE = title, DATE_FORMAT = '%M', _extra = ex
     pltt, ICE_vol_export_2, 't', /REMPLI, /NOERASE  $
@@ -181,10 +182,10 @@
     tot_vol_expo_Sv = tot_vol_expo * 1.e06 * 1/86400 * 1/365 ; annual mean in Sverdrup
     tot_vol_expo_2_Sv = tot_vol_expo_2 * 1.e06 * 1/86400 * 1/365 ; annual mean in Sverdrup
-    xyouts, julday(6,15,1900), 0.48, 'Tot. Annual Export OBS = 2124 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2
-    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
-    xyouts, julday(9,15,1900), 0.44, 'in Sv = '+strtrim(tot_vol_expo_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2
-    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
-    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
-    xyouts, julday(7,15,1900), 0.38, 'Data from Kwok et al.(2004), 1992-1998 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2
+    xyouts, julday(5,15,1900), 490, 'Tot. Annual Export OBS = 2124 10^3 Km3/year', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2
+    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
+    xyouts, julday(9,15,1900), 440, 'in Sv = '+strtrim(tot_vol_expo_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2
+    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
+    xyouts, julday(9,15,1900), 390, 'in Sv = '+strtrim(tot_vol_expo_2_Sv, 1)+' Sv', ALIGN = 0, CHARTHICK = 2, CHARSIZE=1, COLOR=2
+    xyouts, julday(7,15,1900), 360, 'Data from Kwok et al.(2004), 1992-1998 ', ALIGN = 0, CHARTHICK = 2, CHARSIZE=0.8, COLOR=2
 
     htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png  />  ' ]
Index: /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_Vel.pro
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_Vel.pro	(revision 4514)
+++ /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_ICE_Vel.pro	(revision 4514)
@@ -0,0 +1,125 @@
+function read_arr2d, filename, varname, t1, t2
+;; function that read input file and return 2d array with monthly timecounter 
+nyear = (t2-t1+1)/12
+arr2d = ncdf_lec(filename, VAR=varname)
+arr2d = arr2d[t1:t2]
+arr2d = reform(arr2d,12,nyear) ; put in 2D array
+arr2d = total(arr2d,2)/nyear ; total over 2th dimension (i.e.years)
+
+return, arr2d
+end 
+
+;; here start procedure that use function read_arr2d
+pro std_ts_ICE_Vel, masknp, s_iodir_data,  POSTSCRIPT = postscript, _extra = ex
+
+  compile_opt idl2, strictarrsubs
+  
+@common
+@std_common
+
+; get exp1 info
+  vICE_vel_1 = getenv('VAR1_Ivel')   &   prefix = getenv('V1Iv_PREF')  &   suffix = getenv('V1Iv_SUFF')
+; get exp2 info
+  vICE_vel_2 = getenv('VAR2_Ivel')   &   prefix2 = getenv('V2Iv_PREF')  &  suffix2 = getenv('V2Iv_SUFF')
+
+  cdti3 = string(cnt, format = '(i3.3)')
+  print, cdti3 + ') ' + blabla
+  filename = cdti3 + '_ts_ICE_Vel_'+prefix
+  if prefix NE prefix2 then filename = filename + '_'+prefix2
+  if KEYWORD_SET(postscript) then openps, filename+'.ps', portrait = 1
+;
+  d1_d2 = '('+strtrim(date1, 1)+' - '+strtrim(date2, 1)+')'
+;
+  iodir = std_iodir_data
+  ; ICE Velocity in NORTH Hemisphere
+  domdef, 0, jpi-1, 30, 90, /xindex
+  Velo_N = rseries_ncdf(vICE_vel_1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec,/nostruct) ;!! warning positive northward
+  
+  ;;;;;; if grossempoyenne is not present dimension of array is not ok: 54   Array[180, 54, 120] 
+  ;;;;;; with the domain dimensions [jpi/nx, jpj/ny, jpk/nz, jpt] = [180/180, 149/50, 31/31, 12]
+  
+  ICE_velo_N = grossemoyenne(Velo_N, 'xy', /integration, mask2d = masknp)
+  help, ICE_velo_N
+  print, 'max ice velo', max(ICE_velo_N)
+  
+  if jpt mod 12 ne 0 then stop
+  nyr=jpt/12.
+  ICE_velo_N = reform(ICE_velo_N, 12, nyr)
+  ICE_velo_N = total(ICE_velo_N,2)/nyr
+  ICE_velo_N = {arr:ICE_velo_N * 1.e-19 * 86400 * 365, unit : 'm/year'}  
+  help, ICE_velo_N
+  print, 'max ice velo', max(ICE_velo_N.arr)
+    
+  ;ICE Velocity in SOUTH Hemisphere
+  domdef, 0, jpi-1, -90, -30, /xindex
+  Velo_S = rseries_ncdf(vICE_vel_1, date1, date2, prefix, suffix, FIRSTONLY = 1 - allrec,/nostruct) ;!! warning positive northward
+   
+    
+  title = 'Northern Hemisphere'+'!C'+prefix+' '+d1_d2+'!C'+'Ice Velocity (Black SOLID simulation)'
+  jpt=12
+  time=julday(1,15,1900)+30*lindgen(12)
+  pltt, ICE_velo_N, 't', /REMPLI, /PORTRAIT, MIN = 0., MAX = 60000000. , XGRIDSTYLE = 1 $
+        , small = [1, 2, 1], YTITLE = varunit, TITLE = title, _extra = ex    
+  ;
+  title ='Southern Hemisphere' +'!C'+prefix+' '+d1_d2+' - '+'!C'+'Ice Velocity (Black SOLID simulation)'
+  pltt, Velo_S, 't',  /REMPLI, MIN = 0., MAX =  60000000., /NOERASE, XGRIDSTYLE = 1 $
+       , small = [1, 2, 2], YTITLE = varunit, TITLE = title, _extra = ex
+  ;
+
+  htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'.png  />  ' ]
+  if KEYWORD_SET(postscript) then closeps
+
+  if prefix NE prefix2 then BEGIN
+
+    d1_d2_2 = '('+strtrim(date1_2, 1)+' - '+strtrim(date2_2, 1)+')'
+    tsave = time
+    ; ICE Velocity in NORTH Hemisphere
+    domdef, 0, jpi-1, 30, 90, /xindex
+    Velo_N_2 = rseries_ncdf(vICE_vel_2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec,/nostruct)
+    help, Velo_N_2 
+    
+    ;;;;;; if grossempoyenne is not present dimension of array is not ok: 54   Array[180, 54, 120] 
+    ;;;;;; with the domain dimensions [jpi/nx, jpj/ny, jpk/nz, jpt] = [180/180, 149/50, 31/31, 12]
+  
+    ICE_velo_N_2 = grossemoyenne(Velo_N_2, 'xy', /integration, mask2d = masknp)
+    help, ICE_velo_N_2
+    print, 'max ice velo', max(ICE_velo_N_2)
+  
+    if jpt mod 12 ne 0 then stop
+    nyr=jpt/12.
+    ICE_velo_N_2 = reform(ICE_velo_N_2, 12, nyr)
+    ICE_velo_N_2 = total(ICE_velo_N_2,2)/nyr
+    ICE_velo_N_2 = {arr:ICE_velo_N_2 * 86400 * 365, unit : 'm/year'}  
+    help, ICE_velo_N_2
+    print, 'max ice velo', max(ICE_velo_N_2.arr)
+    
+  
+    ;ICE Velocity in SOUTH Hemisphere
+    domdef, 0, jpi-1, -90, -30, /xindex
+    Velo_S_2 = rseries_ncdf(vICE_vel_2, date1_2, date2_2, prefix2, suffix2, FIRSTONLY = 1 - allrec,/nostruct)
+;SF    help, Velo_S_2
+      
+    if KEYWORD_SET(postscript) then openps, filename+'_2.ps', portrait = 1
+
+    title = 'Northern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 +'!C'+'Ice Velocity (BLACK) '
+    pltt, ICE_velo_N , 't', MIN = 0., MAX =  60000000. , date1, date2, /REMPLI, XGRIDSTYLE = 1, window = 2  $
+          , small = [1, 2, 1], YTITLE = varunit, TITLE = title, /noerase, _extra = ex
+    pltt, ICE_velo_N_2, 't', date1, date2, /REMPLI $
+          , /ov1d, color = 250, small = [1, 2, 1], YTITLE = varunit, TITLE = title, /noerase, _extra = ex
+ 
+    title = 'Southern Hemisphere'+'!C'+prefix+' (BLACK) - '+prefix2+' (RED) '+d1_d2_2 +'!C'+'Ice Velocity (BLACK) '
+    pltt, Velo_S , 't', MIN = -.04, MAX = .04, date1, date2, /REMPLI, XGRIDSTYLE = 1 $
+          , small = [1, 2, 2], YTITLE = varunit, TITLE = title, /noerase, _extra = ex
+    pltt, Velo_S_2, 't', date1, date2, /REMPLI $
+         , /ov1d, color = 250, small = [1, 2, 2], YTITLE = varunit, TITLE = title, /noerase, _extra = ex
+
+    htmltxt = [ htmltxt, '<hr>'+blabla, '<br><img width="80%" src='+filename+'_2.png  />  ' ]
+    if KEYWORD_SET(postscript) then closeps
+
+  endif
+
+  domdef
+  
+
+  return
+end
Index: /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro
===================================================================
--- /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro	(revision 4513)
+++ /branches/2013/dev_MERGE_2013/NEMOGCM/CONFIG/ORCA2_LIM/IDL_scripts/std_ts_all.pro	(revision 4514)
@@ -78,6 +78,6 @@
   IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vol, masknp, POSTSCRIPT = postscript, _extra = ex
 ;
-;SF  cnt = cnt+1   &   blabla = 'ICE Velocity '
-;SF  IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex
+;SF    cnt = cnt+1   &   blabla = 'ICE Velocity '
+;SF    IF doplot EQ cnt OR doplot EQ 0 THEN std_ts_ICE_Vel, masknp, POSTSCRIPT = postscript, _extra = ex
 ;
   cnt = cnt+1   &   blabla = 'Sea-Ice Volume Export at Fram Strait'
