# Diagnostic Package Choices
#-----------------
# for each output-stream:
#  filename(n) : prefix of the output file name (only 8.c long) for outp.stream n
#  frequency(n):< 0 : write snap-shot output every multiple of |frequency| (iter)
#               > 0 : write time-average output every multiple of frequency (iter)
#  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 diagnostics fields (8.c) (see "available_diagnostics" file 
#                 for the list of all available diag. in this particular config)
#--------------------------------------------------------------------
#
 &diagnostics_list
#
   dumpatlast = .TRUE.,
#
#---
  frequency(1) = 2635200.0,
   fields(1:4,1) = 'ETAN    ','SIarea  ','SIheff  ','SIhsnow ',
   filename(1) = 'diags/state_2d_set1',
#---
  frequency(3) = 2635200.0,
   fields(1:2,3) = 'THETA   ','SALT    ',
   filename(3) = 'diags/state_3d_set1',
#---
  frequency(5) = 2635200.0,
   fields(1:3,5) = 'DETADT2 ','MXLDEPTH','PHIBOT  ',
   filename(5) = 'diags/other_2d_set1',
#---
  frequency(6) = 2635200.0,
   fields(1:4,6) = 'SIatmQnt','SIatmFW ','sIceLoad','oceSPDep',
   filename(6) = 'diags/other_2d_set2',
#---
  frequency(7) = 2635200.0,
   fields(1:4,7) = 'oceTAUX ','oceTAUY ','oceQnet ','oceFWflx',
   filename(7) = 'diags/other_2d_set3',
#---
  frequency(9) = 2635200.0,
   fields(1:2,9) = 'DRHODR  ','RHOAnoma',
   filename(9) = 'diags/other_3d_set1',
#---
  frequency(11) = 2635200.0,
   fields(1:3,11) = 'UVELMASS','VVELMASS','WVELMASS',
   filename(11) = 'diags/trsp_3d_set1',
#---
  frequency(12) = 2635200.0,
   fields(1:3,12) = 'ADVx_TH ','ADVy_TH ','ADVr_TH ',
   filename(12) = 'diags/trsp_3d_set2',
#---
  frequency(13) = 2635200.0,
   fields(1:3,13) = 'ADVx_SLT','ADVy_SLT','ADVr_SLT',
   filename(13) = 'diags/trsp_3d_set3',
#---
  frequency(14) = 2635200.0,
   fields(1:2,14) = 'DFxE_TH ','DFyE_TH ',
   filename(14) = 'diags/trsp_3d_set4',
#---
  frequency(15) = 2635200.0,
   fields(1:2,15) = 'DFrE_TH ','DFrI_TH ',
   filename(15) = 'diags/trsp_3d_set5',
#---
  frequency(16) = 2635200.0,
   fields(1:2,16) = 'DFxE_SLT','DFyE_SLT',
   filename(16) = 'diags/trsp_3d_set6',
#---
  frequency(17) = 2635200.0,
   fields(1:2,17) = 'DFrE_SLT','DFrI_SLT',
   filename(17) = 'diags/trsp_3d_set7',
#---
  frequency(18) = 2635200.0,
   fields(1:2,18) = 'GM_PsiX ','GM_PsiY ',
   filename(18) = 'diags/trsp_3d_set8',
#---
#---
#---
#---
  frequency(20) = -604800.0,
   fields(1:3,20) = 'ETAN    ','SIheff  ','SIhsnow ',
   filename(20) = 'diags/budg2d_snap_set1',
  timePhase(20) = 0.,
  fileFlags(20) = 'D       ',
#---
  frequency(21) = -604800.0,
   fields(1:2,21) = 'THETA   ','SALT    ',
   filename(21) = 'diags/budg2d_snap_set2',
  timePhase(21) = 0.,
  fileFlags(21) = 'DI      ',
#---
  frequency(22) = -604800.0,
   fields(1:2,22) = 'sIceLoad','PHIBOT  ',
   filename(22) = 'diags/budg2d_snap_set3',
  timePhase(22) = 0.,
  fileFlags(22) = 'D       ',
#---
  frequency(24) = 604800.0,
   fields(1:4,24) = 'oceFWflx','SIatmFW ',
# SDIAG5 will become SItflux
                 'TFLUX   ','SDIAG5  ',
   filename(24) = 'diags/budg2d_zflux_set1',
  fileFlags(24) = 'D       ',
#---
  frequency(25) = 604800.0,
# Need SIdHbFLO, along with SIheff, to match oceSflux since 
# flooding takes SIsal0 from ocean, whereas the other terms 
# in dHEFF also take SIsal0 from ocean!!!BUG???
# This point needs clarification, and is likely a BUG : 
# I think I should count FLO positively now.
# Also it is oceSflux and not SFLUX that matches dSALT 
# because SFLUX lacks +oceSPflx.
#                'SFLUX   ','oceSflux','SIdHbFLO','oceSPflx',
# after bug fixes
   fields(1:6,25) = 'SFLUX   ','oceSPflx',
# to check effect of LWDOWN*emissi with regard to evap
#                 'EXFpreci','EXFevap ',
# when SSS relax without global mean removal
                 'SRELAX  ','SDIAG1  ',
#when removing global mean fresh water and heat
                 'SDIAG2  ','SDIAG3  ',
   filename(25) = 'diags/budg2d_zflux_set2',
  fileFlags(25) = 'D       ',
#---
  frequency(27) = 604800.0,
   fields(1:2,27) = 'UVELMASS','VVELMASS',
   filename(27) = 'diags/budg2d_hflux_set1',
  fileFlags(27) = 'DI      ',
#---
  frequency(28) = 604800.0,
   fields(1:4,28) = 'ADVx_TH ','ADVy_TH ','DFxE_TH ','DFyE_TH ',
   filename(28) = 'diags/budg2d_hflux_set2',
  fileFlags(28) = 'DI      ',
#---
  frequency(29) = 604800.0,
   fields(1:4,29) = 'ADVx_SLT','ADVy_SLT','DFxE_SLT','DFyE_SLT',
   filename(29) = 'diags/budg2d_hflux_set3',
  fileFlags(29) = 'DI      ',
#---
  frequency(30) = 604800.0,
   fields(1:4,30) ='ADVxHEFF','ADVyHEFF','DFxEHEFF','DFyEHEFF',
   filename(30) = 'diags/budg2d_hflux_set4',
  fileFlags(30) = 'D       ',
#---
  frequency(31) = 604800.0,
   fields(1:4,31) ='ADVxSNOW','ADVySNOW','DFxESNOW','DFyESNOW',
   filename(31) = 'diags/budg2d_hflux_set5',
  fileFlags(31) = 'D       ',
#---
  frequency(33) = 604800.0,
   fields(1:4,33) ='EXFpreci','EXFevap ','EXFroff ','EXFempmr',
   filename(33) = 'diags/exf_zflux_set1',
#---
  frequency(34) = 604800.0,
   fields(1:5,34) ='EXFswdn ','EXFlwdn ','EXFswnet','EXFlwnet','EXFqnet ',
   filename(34) = 'diags/exf_zflux_set2',
#---
#---
#---
#---
#
 /
#
#
# 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
# diagSt_regMaskFile='basin_masks_eccollc_90x50.bin',
# nSetRegMskFile=1,
# set_regMask(1)= 1,  1,  1,  1,  1,  1,  1,  1,  1,
#                 1,  1,  1,  1,  1,  1,  1,  1
# val_regMask(1)= 1., 2., 3., 4., 5., 6., 7., 8., 9.,
#                10.,11.,12.,13.,14.,15.,16.,17.
##---
# stat_fields(1,1)= 'ETAN    ','ETANSQ  ','DETADT2 ',
#                   'UVEL    ','VVEL    ','WVEL    ',
#                   'THETA   ','SALT    ',
#    stat_fname(1)= 'dynStDiag',
#     stat_freq(1)= 2635200.,
# stat_region(1,1)=  1, 2, 3, 4, 5, 6, 7, 8, 9,
#                   10,11,12,13,14,15,16,17
##---
# stat_fields(1,2)= 'oceTAUX ','oceTAUY ',
#                   'surForcT','surForcS','TFLUX   ','SFLUX   ',
#                   'oceQnet ','oceSflux','oceFWflx',
#    stat_fname(2)= 'surfStDiag',
#     stat_freq(2)= 2635200.,
# stat_region(1,2)=  1, 2, 3, 4, 5, 6, 7, 8, 9,
#                   10,11,12,13,14,15,16,17
 /