# Diagnostic Package Choices
#--------------------
#  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"
#--------------------
# This example dumps EXF diagnostics as snapshot after 10 time-steps
# Note: EXF air-sea fluxes over Sea-Ice are wrong
 &DIAGNOSTICS_LIST
   dumpAtLast  = .TRUE.,
#--
  fields(1:11,1) = 'EXFtaux ','EXFtauy ','EXFqnet ','EXFempmr',
                   'EXFhl   ','EXFhs   ','EXFswnet','EXFlwnet',
                   'EXFuwind','EXFvwind','EXFatemp',
#  fileName(1) = 'exfDiag',
  frequency(1) = 86400.,

#- with pkg/thsice:
# fields(1:11,2) = 'SI_Fract','SI_Thick','THETA   ','SI_Tsrf ',
##                 'SI_Tsrf ','SI_Tice1','SI_Tice2',
##                 'SI_Qice1','SI_Qice2',
##                 'SIsnwPrc','SIalbedo','SIsnwAge',
##                 'oceQnet ','oceQsw  ','oceFWflx','oceSflux',
#                  'SIflx2oc','SIfrw2oc','SIsaltFx',
#                  'SIflxAtm','SIfrwAtm',
#                  'EXFqnet ','EXFempmr',
#- with pkg/seaice:
  fields(1:10,2) = 'SIarea  ','SIheff  ','THETA   ','SItices ',
#                  'SIuice  ','SIvice  ','SIhsnow ',
#                  'oceQnet ','oceQsw  ','oceFWflx','oceSflux',
                   'SIqnet  ','SIqsw   ','SIempmr ','oceSflux',
                   'SIatmQnt','SIatmFW ',
   fileName(2) = 'iceDiag',
  frequency(2) =  86400.,
  missing_value(2) = -999.,

  fields(1:10,3) = 'SIarea  ','SIheff  ','THETA   ','SItices ',
                   'SIqnet  ','SIqsw   ','SIempmr ','oceSflux',
                   'SIatmQnt','SIatmFW ',
   fileName(3) = 'snapshot',
  frequency(3) = -86400.,
  timePhase(3) =  -3600.,
  missing_value(3) = -999.,

 &

#--------------------
# 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
#- with pkg/thsice:
#stat_fields(1:12,1) = 'SI_Fract','SI_Thick','SI_SnowH','THETA   ',
#stat_fields(1:11,1) = 'SI_Fract','SI_Thick','THETA   ',
#                      'SI_Tsrf ','SI_Tice1','SI_Tice2',
#                      'SIflx2oc','SIfrw2oc','SIsaltFx',
#                      'SIflxAtm','SIfrwAtm',
#- with pkg/seaice:
#stat_fields(1:11,1) = 'SIarea  ','SIheff  ','SIhsnow ',
 stat_fields(1:10,1) = 'SIarea  ','SIheff  ',
                       'THETA   ','SItices ',
#                      'oceQnet ','oceQsw  ','oceFWflx','oceSflux',
                       'SIqnet  ','SIqsw   ','SIempmr ','oceSflux',
                       'SIatmQnt','SIatmFW ',
  stat_fName(1) = 'iceStDiag',
   stat_freq(1) = 43200.,
  stat_phase(1) = 3600.,
 &