--- MITgcm/verification/hs94.cs-32x32x5/input/data.diagnostics	2005/01/28 01:22:11	1.5
+++ MITgcm/verification/hs94.cs-32x32x5/input/data.diagnostics	2010/12/19 00:37:15	1.16
@@ -1,17 +1,62 @@
 # Diagnostic Package Choices
- &diagnostics_list
-  diag_mnc=.TRUE.,
-  frequency(1) = 10,
-   levels(1,1) = 1.,
-   fields(1,1) = 'ETAN    ','ETANSQ  ','DETADT2 ',
-   filename(1) = 'surfDiag',
-  frequency(2) = 20,
-   levels(1,2) = 1.,2.,3.,4.,5.,
-   fields(1,2) = 'UVEL    ','VVEL    ','WVEL    ','THETA   ','SALT    ',
-   filename(2) = 'dynDiag',
-  frequency(3) = 20,
-   levels(1,3) = 1.,2.,3.,4.,5.,
-   fields(1,3) = 'UVELSQ  ','VVELSQ  ','THETASQ ','UV_VEL_C',
-                 'UTHMASS ','VTHMASS ','WVELTH  ',
-   filename(3) = 'secmoms',
- &end
+#--------------------
+#  dumpAtLast (logical): always write output at the end of simulation (default=F)
+#  diag_mnc   (logical): write to NetCDF files (default=useMNC)
+#--for each output-stream:
+#  fileName(n) : prefix of the output file name (max 80c long) for outp.stream n
+#  frequency(n):< 0 : write snap-shot output every |frequency| seconds
+#               > 0 : write time-average output every frequency seconds
+#  timePhase(n)     : write at time = timePhase + multiple of |frequency|
+#    averagingFreq  : frequency (in s) for periodic averaging interval
+#    averagingPhase : phase     (in s) for periodic averaging interval
+#    repeatCycle    : number of averaging intervals in 1 cycle
+#  levels(:,n) : list of levels to write to file (Notes: declared as REAL)
+#                when this entry is missing, select all common levels of this list
+#  fields(:,n) : list of selected diagnostics fields (8.c) in outp.stream n
+#                (see "available_diagnostics.log" file for the full list of diags)
+#  missing_value(n) : missing value for real-type fields in output file "n"
+#  fileFlags(n)     : specific code (8c string) for output file "n"
+#--------------------
+ &DIAGNOSTICS_LIST
+# diag_mnc     = .FALSE.,
+#--
+  fields(1:3,1)  = 'ETAN    ','ETANSQ  ','DETADT2 ',
+   fileName(1) = 'surfDiag',
+# fileFlags(1) = 'D       ',
+  frequency(1) = 12000.,
+  fields(1:5,2)  = 'UVEL    ','VVEL    ','WVEL    ','THETA   ','SHAP_dU ',
+   fileName(2) = 'dynDiag',
+# fileFlags(2) = 'D       ',
+  frequency(2) = 2592000.,
+  fields(1:7,3)  = 'UVELSQ  ','VVELSQ  ','THETASQ ','UV_VEL_C',
+                   'UTHMASS ','VTHMASS ','WVELTH  ',
+   fileName(3) = 'secmoms',
+# fileFlags(3) = 'D       ',
+  frequency(3) = 2592000.,
+ &
+
+#--------------------
+# Parameter for Diagnostics of per level statistics:
+#--------------------
+#  diagSt_mnc (logical): write stat-diags to NetCDF files (default=diag_mnc)
+#  diagSt_regMaskFile : file containing the region-mask to read-in
+#  nSetRegMskFile   : number of region-mask sets within the region-mask file
+#  set_regMask(i)   : region-mask set-index that identifies the region "i"
+#  val_regMask(i)   : region "i" identifier value in the region mask
+#--for each output-stream:
+#  stat_fName(n) : prefix of the output file name (max 80c 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 selected diagnostics fields (8.c) in outp.stream n
+#                (see "available_diagnostics.log" file for the full list of diags)
+#--------------------
+ &DIAG_STATIS_PARMS
+#- an example just to check the agreement with MONITOR output:
+#stat_fields(1:5,1)  = 'ETAN    ','UVEL    ','VVEL    ','WVEL    ', 'THETA   ',
+# stat_fName(1) = 'dynStDiag',
+#  stat_freq(1)= -864000.,
+# stat_phase(1) = 0.,
+ &
+