# 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"
#--------------------
 &DIAGNOSTICS_LIST
# diag_mnc     = .FALSE.,
# dumpAtLast   = .TRUE.,
#--
  fields(1:12,1) = 'ETAN    ','ETANSQ  ','DETADT2 ','PHIBOT  ','PHIBOTSQ',
                   'oceTAUX ','oceTAUY ',
                   'oceQnet ','oceSflux','oceFWflx',
#                  'surForcT','surForcS','TFLUX   ','SFLUX   ','oceFreez',
                   'TRELAX  ','SRELAX  ',
#                  'GM_VisbK',
   fileName(1) = 'surfDiag',
  frequency(1) = 1555200000.,
  fields(1:15,2) = 'UVEL    ','VVEL    ','WVEL    ','PHIHYD  ',
                   'VVELMASS','UVELMASS','WVELSQ  ',
                   'THETA   ','UTHMASS ','VTHMASS ','WTHMASS ',
                   'SALT    ','USLTMASS','VSLTMASS','WSLTMASS',
# do not specify levels => all levels are selected
#  fileName(2) = 'dynDiag',
  frequency(2) = 1555200000.,

  fields(1:8,3)  = 'DRHODR  ','RHOAnoma','CONVADJ ',
#                  'GM_PsiX ','GM_PsiY ',
                   'GM_Kwx  ','GM_Kwy  ','GM_Kwz  ',
                   'GM_Kux  ','GM_Kvy  ',
#                  'GM_Kuz  ','GM_Kvz  ',
#- disable this output list by commenting out the file name
#  fileName(3) = 'oceDiag',
  frequency(3) = 1555200000.,

  fields(1:7,4)  = 'ADVx_TH ','ADVy_TH ','ADVr_TH ',
                   'DIFx_TH ','DIFy_TH ','DFrE_TH ',
                   'DFrI_TH ',
#                  'ADVx_SLT',
#  fileName(4) = 'flxDiag',
  frequency(4) = 1728000.,

  fields(1:5,5)  = 'SIheff  ','SIarea  ','SIhsnow ',
                   'SIuice  ','SIvice  ',
  fileName(5) = 'seaiceDiag',
  frequency(5) = 1555200000.,
 averagingFreq(5) = 2592000.,
   repeatCycle(5) = 12,
#---
  frequency(6) = 86400.0,
  fields(1:12,6) ='UVELMASS','VVELMASS',
                 'ADVx_TH ','ADVy_TH ','DFxE_TH ','DFyE_TH ',
                 'ADVx_SLT','ADVy_SLT','DFxE_SLT','DFyE_SLT',
                 'AB_gT   ','AB_gS   ',
#  filename(6) = 'budg2d_hflux_set2',
   fileFlags(6) = 'DI      ',
#---
  frequency(7) = -86400.0,
   fields(1:3,7) = 'ETAN    ','SIheff  ','SIhsnow ',
#  filename(7) = 'budg2d_snap_set1',
   timePhase(7)= 0.,
   fileFlags(7) = 'D       ',
#---
  frequency(8) = -86400.0,
   fields(1:2,8) = 'THETA   ','SALT    ',
#  filename(8) = 'budg2d_snap_set2',
   timePhase(8)= 0.,
   fileFlags(8) = 'DI      ',
#---
  frequency(9) = 86400.0,
   fields(1:16,9) = 'oceFWflx','SIatmFW ','TFLUX   ','SItflux ','SFLUX   ',
                  'SRELAX  ','TRELAX  ','oceSflux','oceQnet ','SIatmQnt','oceQsw  ',
                  'WTHMASS ','WSLTMASS','SIaaflux','SIsnPrcp','SIacSubl',
   levels(1,9) = 1.,
#  filename(9) = 'budg2d_zflux_set1',
   fileFlags(9) = 'D       ',
#---
  frequency(10) = 86400.0,
   fields(1:8,10) ='ADVxHEFF','ADVyHEFF','DFxEHEFF','DFyEHEFF',
                 'ADVxSNOW','ADVySNOW','DFxESNOW','DFyESNOW',
#  filename(10) = 'budg2d_hflux_set1',
   fileFlags(10) = 'D       ',
 &

#--------------------
# 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
#- 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:3) = 1,  1,  1,
 val_regMask(1:3) = 1., 2., 3.,
#---
 stat_fields(1:8,1)  = 'ETAN    ','UVEL    ','VVEL    ','WVEL    ',
                       'THETA   ','SALT    ','CONVADJ ','DETADT2 ',
  stat_fName(1) = 'dynStDiag',
   stat_freq(1) = 864000.,
 stat_fields(1:5,5)  = 'SIheff  ','SIarea  ','SIhsnow ',
                       'SIuice  ','SIvice  ',
 stat_region(1:3,5)  = 1, 3, 0,
  stat_fName(5) = 'seaiceStDiag',
   stat_freq(5) = 864000.,
 &