--- MITgcm/verification/aim.5l_cs/input.thSI/data.diagnostics	2005/05/14 20:56:46	1.3
+++ MITgcm/verification/aim.5l_cs/input.thSI/data.diagnostics	2006/02/21 01:52:15	1.12
@@ -11,14 +11,15 @@
 #                 for the list of all available diag. in this particular config)
 #-----------------
  &diagnostics_list
-  frequency(1) = 4500,,
+  frequency(1) = 4500.,
    levels(1,1) = 1.,
    fields(1,1) = 'ETAN    ','ETANSQ  ','DETADT2 ',
    filename(1) = 'surfDiag',
   frequency(2) = 4500.,
    levels(1,2) = 1.,2.,3.,4.,5.,
    fields(1,2) = 'UVEL    ','VVEL    ','WVEL    ','THETA   ','SALT    ',
-                 'UVELSQ  ','VVELSQ  ','THETASQ ','UV_VEL_Z',
+                 'UVELSQ  ','VVELSQ  ','WVELSQ  ','THETASQ ',
+                 'UV_VEL_Z','PHIHYD  ','PHIHYDSQ',
                  'UVELMASS','VVELMASS',
                  'UTHMASS ','VTHMASS ','WVELTH  ',
                  'USLTMASS','VSLTMASS','WVELSLT ',
@@ -26,13 +27,13 @@
   frequency(3) = 4500.,
    levels(1,3) = 1.,
    fields(1,3) = 'UFLUX   ', 'VFLUX   ',
-                 'OSR     ', 'OLR     ', 'RADSWG  ', 'RADLWG  ', 
+                 'TSR     ', 'OLR     ', 'RADSWG  ', 'RADLWG  ', 
                  'HFLUX   ', 'EVAP    ', 'PRECON  ', 'PRECLS  ', 
-#                'CLDFRC  ', 'CLDPRS  ', 'CTPCNT  ', 'CLDMAS  ',
                  'CLDFRC  ', 'CLDPRS  ', 'CLDMAS  ',
-                 'KM      ', 'WINDS   ', 'TS      ', 'QS      ',
-                 'ENPREC  ', 'ALBVISDF', 'SDIAG1  ',
-#                'OSRCLR  ', 'OLRCLR  ', 'SWGCLR  ', 'LWGCLR  ',
+                 'DRAG    ', 'WINDS   ', 'TS      ', 'QS      ',
+                 'ENPREC  ', 'ALBVISDF', 'DTSIMPL ',
+#                'DWNLWG  ',
+#                'TSRCLR  ', 'OLRCLR  ', 'SWGCLR  ', 'LWGCLR  ',
    filename(3) = 'aimDiag',
   frequency(4) = 4500.,
    levels(1,4) = 1.,2.,3.,4.,5.,
@@ -44,10 +45,49 @@
   frequency(5) = 4500.,
    fields(1,5) = 'SI_Fract', 'SI_Thick', 'SI_SnowH',
                  'SI_Tsrf ', 'SI_Tice1', 'SI_Tice2',
-#                'SI_Qice1', 'SI_Qice2', 'SIsnwAge',
+                 'SI_Qice1', 'SI_Qice2', 'SIsnwAge',
                  'SIsnwPrc', 'SIalbedo',
                  'SIflx2oc', 'SIfrw2oc', 'SIsaltFx',
 #                'SIflxAtm', 'SIfrwAtm',
 #                'SItOcMxL', 'SIsOcMxL',
    filename(5) = 'thSIceDiag',
- &end
+ &
+
+# Parameter for Diagnostics of per level statistics:
+#-----------------
+# for each output-stream:
+#  stat_fname(n) : prefix of the output file name (only 8.c long) for outp.stream n
+#  stat_freq(n):< 0 : write snap-shot output every |stat_freq| seconds
+#               > 0 : write time-average output every stat_freq seconds
+#  stat_phase(n)    : write at time = stat_phase + multiple of |stat_freq|
+#  stat_region(:,n) : list of "regions" (default: 1 region only=global)
+#  stat_fields(:,n) : list of diagnostics fields (8.c) (see "available_diagnostics.log"
+#                 file for the list of all available diag. in this particular config)
+#-----------------
+ &DIAG_STATIS_PARMS
+#- regional mask: 3 lat. band: 1 : y <= -24 ; 2 : -24<y<24 ; 3 : 24 <= y
+ diagSt_regMaskFile='regMask_lat24.bin',
+ nSetRegMskFile=1,
+ set_regMask(1)= 1,  1,  1,
+ val_regMask(1)= 1., 2., 3.,
+#- an example just to check the agreement with MONITOR output:
+ stat_fields(1,1)= 'ETAN    ','UVEL    ','VVEL    ','WVEL    ', 'THETA   ',
+    stat_fname(1)= 'dynStDiag',
+     stat_freq(1)= -3600.,
+    stat_phase(1)= 0.,
+ stat_fields(1,4)= 'LdSnowH ', 'GrdSurfT', 'GrdTemp ', 'GrdWater',
+ stat_region(1,4)= 1, 2, 3,
+    stat_fname(4)= 'landStDiag',
+     stat_freq(4)= -3600.,
+    stat_phase(4)= 0.,
+ stat_fields(1,5)= 'SI_Fract', 'SI_Thick', 'SI_SnowH',
+                   'SI_Tsrf ', 'SI_Tice1', 'SI_Tice2',
+#                  'SI_Qice1', 'SI_Qice2',
+                   'SIsnwPrc', 'SIalbedo', 'SIsnwAge',
+                   'SIflx2oc', 'SIfrw2oc', 'SIsaltFx',
+ stat_region(1,5)= 1, 3, 0,
+    stat_fname(5)= 'thSIceStDiag',
+     stat_freq(5)= -3600.,
+    stat_phase(5)= 0.,
+ &
+