# 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:24,1) = 'ETAN    ','SIarea  ','SIheff ','SIhsnow ',
#stuff that is not quite state variables (and may not be quite
#synchroneous) but are added here to reduce number of files
                 'DETADT2 ','PHIBOT  ','sIceLoad',
                 'MXLDEPTH',
                 'SIatmQnt','SIatmFW ','oceQnet ','oceFWflx',
                 'oceTAUX ','oceTAUY ',
                 'ADVxHEFF','ADVyHEFF','DFxEHEFF','DFyEHEFF',
                 'ADVxSNOW','ADVySNOW','DFxESNOW','DFyESNOW',
                 'SIuice  ','SIvice  ',
   filename(1) = 'diags/state_2d_set1',
#---
  frequency(2) = 2635200.0,
   fields(1:17,2) = 'THETA   ','SALT    ', 'TRAC01    ',
                   'TRAC02  ', 'TRAC03 ', 'TRAC04 ', 'TRAC05 ',
                   'TRAC06  ', 'TRAC07 ', 'TRAC08 ', 'TRAC09 ',
                   'TRAC10  ', 'TRAC11 ', 'TRAC12 ', 'TRAC13 ',
                   'TRAC14  ', 'TRAC15 ', 
#stuff that is not quite state variables (and may not be quite
#synchronous) but are added here to reduce file number
#                 'DRHODR  ',
#                'KPPviscA','KPPdiffT','KPPdiffTS',
#                'GGL90TKE','GGL90Lmx','GGL90Prl',
#                'GGL90ArU','GGL90ArV','GGL90Kr ',
   filename(2) = 'diags/state_3d_set1',
#---
  frequency(3) = 2635200.0,
   fields(1:3,3) = 'UVELMASS','VVELMASS','WVELMASS',
#                 'GM_PsiX ','GM_PsiY ',
#uvel is approximately uvelmass divided by initial hfac, so we dont output the former
#                'UVEL    ','VVEL    ','WVEL     ',
#full 3D temperature fluxes : 'DFxE_TH ','DFyE_TH ','DFrE_TH ','DFrI_TH ','ADVx_TH ','ADVy_TH ','ADVr_TH ',
#but for present computations I only need the vertically integrated horizontal components (see trsp_3d_set2)
   filename(3) = 'diags/trsp_3d_set1',
#---
#  frequency(4) = 2635200.0,
#   fields(1:8,4) = 'DFxE_TH ','DFyE_TH ','ADVx_TH ','ADVy_TH ',
#                 'DFxE_SLT','DFyE_SLT','ADVx_SLT','ADVy_SLT',
#   filename(4) = 'diags/trsp_3d_set2',
#vertically integrate fields we only use to compute vertically integr.
#meridional transports (also omit vertical transports, both to save space)
#---
#  frequency(5) = -2635200.0,
#   fields(1:6,5) = 'ETAN    ','SIheff  ','SIhsnow ',
#                 'SIarea  ','sIceLoad','PHIBOT  ',
#   filename(5) = 'diags/budg2d_snap_set1',
#   timePhase(5)= 0.,
#   fileFlags(5) = 'D       ',
#---
#  frequency(6) = -2635200.0,
#   fields(1:2,6) = 'THETA   ','SALT    ',
#   filename(6) = 'diags/budg2d_snap_set2',
#   timePhase(6)= 0.,
#   fileFlags(6) = 'DI      ',
#---
  frequency(4) = 2635200.0,
   fields(1:6,4) = 'oceFWflx','SIatmFW ','TFLUX   ','SItflux ',
                   'SFLUX   ','oceQsw  ',
   filename(4) = 'diags/budg2d_zflux_set1',
   fileFlags(4) = 'D       ',
#---
  frequency(5) = 2635200.0,
   fields(1:21,5) ='EXFhs   ', 'EXFhl   ', 'EXFlwnet', 'EXFswnet', 'EXFlwdn ',
                  'EXFswdn ','EXFqnet ','EXFtaux ','EXFtauy ','EXFuwind',
                  'EXFvwind','EXFwspee','EXFatemp','EXFaqh ','EXFevap ',
                  'EXFpreci','EXFsnow ','EXFempmr','EXFpress','EXFroff ',
                  'EXFroft ',
   filename(5) = 'diags/EXF_forcing_set1',
   fileFlags(5) = 'D       ',
#---
  frequency(6) = 2635200.0,
   fields(1:2,6) = 'RHOAnoma','PHIHYD  ',
   filename(6) = 'diags/state_3d_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)= 3153600.,
# 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)= 3153600.,
# stat_region(1,2)=  1, 2, 3, 4, 5, 6, 7, 8, 9,
#                   10,11,12,13,14,15,16,17
 /