--- MITgcm/verification/aim.5l_cs/input.thSI/data.diagnostics	2005/05/22 23:19:16	1.5
+++ MITgcm/verification/aim.5l_cs/input.thSI/data.diagnostics	2005/11/25 06:07:20	1.8
@@ -13,11 +13,11 @@
  &diagnostics_list
   frequency(1) = 4500.,
    levels(1,1) = 1.,
-   fields(1,1) = 'ETAN    ','ETANSQ  ','DETADT2 ',
+   fields(1:3,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    ',
+   levels(1:5,2) = 1.,2.,3.,4.,5.,
+   fields(1:17,2) = 'UVEL    ','VVEL    ','WVEL    ','THETA   ','SALT    ',
                  'UVELSQ  ','VVELSQ  ','THETASQ ','UV_VEL_Z',
                  'UVELMASS','VVELMASS',
                  'UTHMASS ','VTHMASS ','WVELTH  ',
@@ -25,24 +25,24 @@
    filename(2) = 'dynDiag',
   frequency(3) = 4500.,
    levels(1,3) = 1.,
-   fields(1,3) = 'UFLUX   ', 'VFLUX   ',
+   fields(1:20,3) = 'UFLUX   ', 'VFLUX   ',
                  'OSR     ', 'OLR     ', 'RADSWG  ', 'RADLWG  ', 
                  'HFLUX   ', 'EVAP    ', 'PRECON  ', 'PRECLS  ', 
 #                'CLDFRC  ', 'CLDPRS  ', 'CTPCNT  ', 'CLDMAS  ',
                  'CLDFRC  ', 'CLDPRS  ', 'CLDMAS  ',
                  'KM      ', 'WINDS   ', 'TS      ', 'QS      ',
-                 'ENPREC  ', 'ALBVISDF', 'SDIAG1  ',
+                 'ENPREC  ', 'ALBVISDF', 'DTSIMPL ',
 #                'OSRCLR  ', 'OLRCLR  ', 'SWGCLR  ', 'LWGCLR  ',
    filename(3) = 'aimDiag',
   frequency(4) = 4500.,
-   levels(1,4) = 1.,2.,3.,4.,5.,
-   fields(1,4) = 'RELHUM  ', 'DIABT   ', 'DIABQ   ',
+   levels(1:5,4) = 1.,2.,3.,4.,5.,
+   fields(1:9,4) = 'RELHUM  ', 'DIABT   ', 'DIABQ   ',
                  'RADSW   ', 'RADLW   ', 'DTCONV  ', 'DTLS    ',
                  'DQCONV  ', 'DQLS    ',
 #                'SWCLR   ', 'LWCLR   ',
    filename(4) = 'aimDDT',
   frequency(5) = 4500.,
-   fields(1,5) = 'SI_Fract', 'SI_Thick', 'SI_SnowH',
+   fields(1:11,5) = 'SI_Fract', 'SI_Thick', 'SI_SnowH',
                  'SI_Tsrf ', 'SI_Tice1', 'SI_Tice2',
 #                'SI_Qice1', 'SI_Qice2', 'SIsnwAge',
                  'SIsnwPrc', 'SIalbedo',
@@ -65,9 +65,21 @@
 #-----------------
  &DIAG_STATIS_PARMS
 #- 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)= -864000.,
-#   stat_phase(1)= 0.,
+ stat_fields(1:5,1)= 'ETAN    ','UVEL    ','VVEL    ','WVEL    ', 'THETA   ',
+    stat_fname(1)= 'dynStDiag',
+     stat_freq(1)= -3600.,
+    stat_phase(1)= 0.,
+ stat_fields(1:4,4)= 'LdSnowH ', 'GrdSurfT', 'GrdTemp ', 'GrdWater',
+    stat_fname(4)= 'landStDiag',
+     stat_freq(4)= -3600.,
+    stat_phase(4)= 0.,
+ stat_fields(1:6,5)= 'SI_Fract', 'SI_Thick', 'SI_SnowH',
+                   'SI_Tsrf ', 'SI_Tice1', 'SI_Tice2',
+#                  'SI_Qice1', 'SI_Qice2', 'SIsnwAge',
+#                  'SIsnwPrc', 'SIalbedo',
+#                  'SIflx2oc', 'SIfrw2oc', 'SIsaltFx',
+    stat_fname(5)= 'thSIceStDiag',
+     stat_freq(5)= -3600.,
+    stat_phase(5)= 0.,
  &