! LIMITS NAMELIST
!
! CONTENTS:
!
! NDAY1    : First day of integration (i)
! NDAY2    : Last day of integration (i)
! IDATE    : Model date in YYYYMMDD (i)
! IDATE0   : Initial experiment date in YYYYMMDD (i)
! RUNID    : Experiment name (a)
! EXPCNF   : Experiment configuration (a)
! RUNTYP   : CESM run type (a)
! GRFILE   : Name of file containing grid specification (a)
! ICFILE   : Name of file containing initial conditions, that is either
!            a valid restart file or 'inicon.nc' if climatological based
!            initial conditions are desired (a)
! BACLIN   : Baroclinic time step (sec) (f)
! BATROP   : Barotropic time step (sec) (f)
! MDV2HI   : Laplacian diffusion velocity for momentum dissipation (cm/s) (f)
! MDV2LO   : Laplacian diffusion velocity for momentum dissipation (cm/s) (f)
! MDV4HI   : Biharmonic diffusion velocity for momentum dissipation (cm/s) (f)
! MDV4LO   : Biharmonic diffusion velocity for momentum dissipation (cm/s) (f)
! MDC2HI   : Laplacian diffusivity for momentum dissipation (cm**2/s) (f)
! MDC2LO   : Laplacian diffusivity for momentum dissipation (cm**2/s) (f)
! VSC2HI   : Parameter in deformation-dependent Laplacian viscosity (f)
! VSC2LO   : Parameter in deformation-dependent Laplacian viscosity (f)
! VSC4HI   : Parameter in deformation-dependent Biharmonic viscosity (f)
! VSC4LO   : Parameter in deformation-dependent Biharmonic viscosity (f)
! CBAR     : rms flow speed for linear bottom friction law (cm/s) (f)
! CB       : Nondiemnsional coefficient of quadratic bottom friction (f)
! CWBDTS   : Coastal wave breaking damping resiprocal time scale (1/s) (f)
! CWBDLS   : Coastal wave breaking damping length scale (m) (f)
! MOMMTH   : Momentum equation discretization method. Valid methods:
!            'enscon' (Sadourny (1975) enstrophy conserving), 'enecon'
!            (Sadourny (1975) energy conserving), 'enedis' (Sadourny
!            (1975) energy conserving with some dissipation) (a)
! EITMTH   : Eddy-induced transport parameterization method. Valid
!            methods: 'intdif', 'gm' (a)
! EDRITP   : Type of Richardson number used in eddy diffusivity
!            computation. Valid types: 'shear', 'large scale' (a)
! BMCMTH   : Baroclinic mass flux correction method. Valid methods:
!            'uc' (upstream column), 'dluc' (depth limited upstream
!            column) (a)
! RMPMTH   : Method of applying eddy-induced transport in the remap
!            transport algorithm. Valid methods: 'eitvel', 'eitflx' (a)
! EDWMTH   : Method to estimate eddy diffusivity weight as a function of
!            the ration of Rossby radius of deformation to the
!            horizontal grid spacing. Valid methods: 'smooth', 'step' (a)
! MLRTTP   : Type of mixed layer restratification time scale. Valid
!            types: 'variable', 'constant', 'limited' (a)
! EDSPRS   : Apply eddy mixing suppression away from steering level (l)
! EGC      : Parameter c in Eden and Greatbatch (2008) parameterization (f)
! EGGAM    : Parameter gamma in E. & G. (2008) param. (f)
! EGLSMN   : Minimum eddy length scale in  E. & G. (2008) param. (cm) (f)
! EGMNDF   : Minimum diffusivity in E. & G. (2008) param. (cm**2/s) (f)
! EGMXDF   : Maximum diffusivity in E. & G. (2008) param. (cm**2/s) (f)
! EGIDFQ   : Factor relating the isopycnal diffusivity to the layer
!            interface diffusivity in the Eden and Greatbatch (2008)
!            parameterization. egidfq=difint/difiso () (f)
! RI0      : Critical gradient richardson number for shear driven
!            vertical mixing () (f)
! RM0      : Efficiency factor of wind TKE generation in the Oberhuber
!            (1993) TKE closure () (f)
! RM5      : Efficiency factor of TKE generation by momentum
!            entrainment in the Oberhuber (1993) TKE closure () (f)
! CE       : Efficiency factor for the restratification by mixed layer
!            eddies (Fox-Kemper et al., 2008) () (f)
! BDMTYP   : Type of background diapycnal mixing. If bdmtyp=1 the
!            background diffusivity is a constant divided by the
!            Brunt-Vaisala frequency, if bdmtyp=2 the background
!            diffusivity is constant () (i)
! BDMC1    : Background diapycnal diffusivity times buoyancy frequency
!            frequency (cm**2/s**2) (f)
! BDMC2    : Background diapycnal diffusivity (cm**2/s) (f)
! TDFILE   : Name of file containing tidal wave energy dissipation
!            divided by by bottom buoyancy frequency (a)
! TKEPF    : Fraction of surface TKE that penetrates beneath mixed layer
!            () (f)
! NIWGF    : Global factor applied to the energy input by near-intertial
!            motions () (f)
! NIWBF    : Fraction of near-inertial energy dissipated in the boundary
!            layer () (f)
! NIWLF    : Fraction of near-inertial energy dissipated locally beneath
!            the boundary layer () (f)
! SWAMTH   : Shortwave radiation absorption method. Valid methods:
!            'top-layer', 'jerlov', 'chlorophyll' (a)
! JWTYPE   : Number indicating the Jerlov (1968) water type (i)
! CHLOPT   : Chlorophyll concentration option. Valid options:
!            'climatology' (a)
! CCFILE   : Name of file containing chlorophyll concentration climatology (a)
! TRXDAY   : e-folding time scale (days) for SST relax., if 0 no relax. (f)
! SRXDAY   : e-folding time scale (days) for SSS relax., if 0 no relax. (f)
! TRXDPT   : Maximum mixed layer depth for e-folding SST relaxation (m) (f)
! SRXDPT   : Maximum mixed layer depth for e-folding SSS relaxation (m) (f)
! TRXLIM   : Max. absolute value of SST difference in relaxation (degC) (f)
! SRXLIM   : Max. absolute value of SSS difference in relaxation (psu) (f)
! APTFLX   : Apply diagnosed heat flux flag (l)
! APSFLX   : Apply diagnosed freshwater flux flag (l)
! DITFLX   : Diagnose heat flux flag (l)
! DISFLX   : Diagnose freshwater flux flag (l)
! SRXBAL   : Balance the SSS relaxation (l)
! SCFILE   : Name of file containing SSS climatology used for relaxation (a)
! SMTFRC   : Smooth CESM forcing (l)
! SPRFAC   : Send precipitation/runoff factor to CESM coupler (l)
! ATM_PATH : Path to forcing fields in case of EXPCNF 'ben02clim' or
!            'ben02syn' (a)
! RSTFRQ   : Restart frequency in days (30=1month,365=1year) (i)
! RSTFMT   : Format of restart file (valid arguments are 0 for classic,
!            1 for 64-bit offset and 2 for netcdf4/hdf5 format) (i)
! RSTCMP   : Compression flag for restart file (i)
! IOTYPE   : 0 = netcdf, 1 = pnetcdf
&LIMITS
  NDAY1    = 0
  NDAY2    = 0
  IDATE    = 20000401
  IDATE0   = 20000401
  RUNID    = 'unset'
  EXPCNF   = 'cesm'
  RUNTYP   = 'hybrid'
  GRFILE   = '/work/shared/noresm/inputdata/ocn/blom/grid/grid_tnx1v4_20170622.nc'
  ICFILE   = '/work/shared/noresm/inputdata/ocn/blom/inicon/inicon_tnx1v4_20170622.nc'
  BACLIN   = 3600.
  BATROP   = 60.
  MDV2HI   = 2.
  MDV2LO   = .4
  MDV4HI   = 0.
  MDV4LO   = 0.
  MDC2HI   = 5000.e4
  MDC2LO   = 300.e4
  VSC2HI   = .5
  VSC2LO   = .5
  VSC4HI   = 0.
  VSC4LO   = 0.
  CBAR     = 5.
  CB       = .002
  CWBDTS   = .75e-4
  CWBDLS   = 25.
  MOMMTH   = 'enscon'
  EITMTH   = 'gm'
  EDRITP   = 'large scale'
  BMCMTH   = 'uc'
  RMPMTH   = 'eitvel'
  EDWMTH   = 'smooth'
  MLRTTP   = 'constant'
  EDSPRS   = .true.
  EGC      = 0.85
  EGGAM    = 200.
  EGLSMN   = 4000.e2
  EGMNDF   = 100.e4
  EGMXDF   = 1500.e4
  EGIDFQ   = 1.
  RI0      = 1.2
  RM0      = 1.2
  RM5      = 0.
  CE       = .5
  BDMTYP   = 2
  BDMC1    = 5.e-4
  BDMC2    = .1
  TDFILE   = '/work/shared/noresm/inputdata/ocn/blom/bndcon/tidal_dissipation_tnx1v4_20170605.nc'
  TKEPF    = .006
  NIWGF    = .4
  NIWBF    = .35
  NIWLF    = .5
  SWAMTH   = 'jerlov'
  JWTYPE   = 3
  CHLOPT   = 'climatology'
  CCFILE   = '/work/shared/noresm/inputdata/ocn/blom/bndcon/chlorophyll_concentration_tnx1v4_20170608.nc'
  TRXDAY   = 0.
  SRXDAY   = 0.
  TRXDPT   = 1.
  SRXDPT   = 1.
  TRXLIM   = 1.5
  SRXLIM   = .5
  APTFLX   = .false.
  APSFLX   = .false.
  DITFLX   = .false.
  DISFLX   = .false.
  SRXBAL   = .false.
  SCFILE   = '/work/shared/noresm/inputdata/ocn/blom/bndcon/sss_clim_core_tnx1v4_20170604.nc'
  SMTFRC   = .false.
  SPRFAC   = .false.
  ATM_PATH = 'unset'
  RSTFRQ   = 1
  RSTFMT   = 1
  RSTCMP   = 0
  IOTYPE   = 1
/
 
! NAMELIST FOR CHANNEL WIDTH MODIFICATIONS
!
! CONTENTS:
!
! CWMTAG   : Array of geographical names of channels to be modified (a)
! CWMEDG   : Array of C grid cell edges to be modified. Valid options:
!            'u' or 'v' (a)
! CWMI     : Array of grid cell i-indices (i)
! CWMJ     : Array of grid cell j-indices (i)
! CWMWTH   : Array of modified grid cell widths (m) (f)
! 
&CWMOD
  CWMTAG = 'Gibraltar','Gibraltar'
  CWMEDG =         'u',        'u'
  CWMI   =         105,        106
  CWMJ   =         273,        273
  CWMWTH =       30.e3,      30.e3
/

! NAMELIST FOR MERIDIONAL OVERTURNING AND FLUX DIAGNOSTICS
!
! CONTENTS:
!
! MER_ORFILE : Name of file containing ocean region specification (a)
! MER_MIFILE : Name of file containing zonal section specification for
!              meridional transport computation (a)
! MER_REGNAM : Array of region names for meridional overturning and flux
!              diagnostics (a)
! MER_REGFLG : Array of mask flags in ocean regions file to be included
!              for each region (i)
! MER_MINLAT : Minimum latitude to be considered for each region (f)
! MER_MAXLAT : Maximum latitude to be considered for each region (f)
! 
&MERDIA
  MER_ORFILE = '/work/shared/noresm/inputdata/ocn/blom/grid/ocean_regions_tnx1v4_20190729.nc'
  MER_MIFILE = '/work/shared/noresm/inputdata/ocn/blom/grid/mertra_index_tnx1v4_20190615.dat'
  MER_REGNAM = 'atlantic_arctic_ocean', 'atlantic_arctic_extended_ocean', 'indian_pacific_ocean', 'global_ocean'
  MER_REGFLG(1,:) = 2, 4
  MER_REGFLG(2,:) = 2, 4, 6, 7, 8, 9
  MER_REGFLG(3,:) = 3, 5 
  MER_REGFLG(4,:) = 0
  MER_MINLAT = -34., -34., -34., -90.
  MER_MAXLAT =  90.,  90.,  90.,  90.
/

! NAMELIST FOR SECTION TRANSPORT DIAGNOSTICS
!
! CONTENTS:
!
! SEC_SIFILE : Name of file containing section specification for section
!              transport computation (a)
! 
&SECDIA
  SEC_SIFILE = '/work/shared/noresm/inputdata/ocn/blom/grid/section_index_tnx1v4_20190611.dat'
/

! IO-NAMELIST FOR DIAGNOSTIC OUTPUT
!
! Description:
!   BLOM supports multiple output groups for its diagnostic output.
!   Each output group is represented by one column in the namlist and may
!   have its own output format, averaging period, and file frequency.
!   The maximum number of output groups is currently limited to 10 but
!   can be changed easily in mod_dia.F.
!
!   The output precision can be choosen on a per-variable basis.
!
!   Multiple time-slices can be written to the same output file
!   provided that no variable is written in packed data format
!   (i.e. as int2 with scale factor and offset).
!
!   Compression of the output (i.e. storage of only wet points)
!   and the file format can be choosen on a per-file basis.
!
!   All time periods are specified in number of days for positive
!   integer values and fraction of day for negative integer values.
!   The length of the actual calendar month is used if 30 is written.
!   The length of the actual calendar year is used if 365 is written.
!   A variable is not written when 0 is specified.
!
! Namelist acronyms:
!   GLB_     - global parameters i.e. valid for entire output group
!   H2D_     - 2d fields with the horizontal dimensions
!   LYR_     - 3d fields with sigma layers as vertical coordinate
!   LVL_     - 3d fields with levitus leves as vertical coordinate
!   MSC_     - miscellanous, non-gridded fields
!
! Global parameters:
!   FNAMETAG - tag used in file name (c10) 
!   AVEPERIO - average period in days (i) 
!   FILEFREQ - how often to start a new file in days (i) 
!   COMPFLAG - switch for compressed/uncompressed output (i) 
!   NCFORMAT - netcdf format (valid arguments are 0 for classic,
!              1 for 64-bit offset and 2 for netcdf4/hdf5 format)
!
! Arguments for diagnostic variables:
!   0        - variable is not written
!   2        - variable is written as int2 with scale factor and offset
!   4        - variable is written as real4
!   8        - variable is written as real8
!
! Output variables:
!   ABSWND   - absolute wind speed [m s-1]
!   ALB      - surface albedo []
!   BTMSTR   - Barotropic mass streamfunction [kg s-1]
!   BRNFLX   - brine flux [kg m-2 s-1]
!   BRNPD    - brine plume depth [m]
!   DFL      - non-solar heat flux derivative [W m-2 K-1]
!   EVA      - evaporation [kg m-2 s-1]
!   FICE     - ice concentration [%]
!   FMLTFZ   - fresh water flux due to melting/freezing [kg m-2 s-1]
!   HICE     - ice thickness [m]
!   HMLTFZ   - heat flux due to melting/freezing [W m-2]
!   HSNW     - snow depth [m]
!   IAGE     - ice age [d]
!   IDKEDT   - mixed layer inertial kinetic energy tendency [kg s-3]
!   LIP      - liquid precipitation [kg m-2 s-1]
!   MAXMLD   - maximum mixed layer depth [m]
!   MLD      - mixed layer depth [m]
!   MLDU     - mixed layer depth at u-point [m]
!   MLDV     - mixed layer depth at v-point [m]
!   MLTS     - mixed layer thickness using "sigma-t" criterion [m]
!   MLTSMN   - minimum mixed layer thickness using "sigma-t" criterion [m]
!   MLTSMX   - maximum mixed layer thickness using "sigma-t" criterion [m]
!   MLTSSQ   - mixed layer thickness squared using "sigma-t" criterion [m2]
!   MTKEUS   - mixed layer TKE tendency related to friction velocity [kg s-3]
!   MTKENI   - mixed layer TKE tendency related to near inertial mot. [kg s-3]
!   MTKEBF   - mixed layer TKE tendency related to buoyancy forcing [kg s-3]
!   MTKERS   - mixed layer TKE tendency related to eddy restrat. [kg s-3]
!   MTKEPE   - mixed layer TKE tendency related to pot. energy change [kg s-3]
!   MTKEKE   - mixed layer TKE tendency related to kin. energy change [kg s-3]
!   MTY      - wind stress y-component [N m-2]
!   MXLU     - mixed layer velocity x-component [m s-1]
!   MXLV     - mixed layer velocity y-component [m s-1]
!   NSF      - non-solar heat flux [W m-2]
!   PBOT     - bottom pressure [Pa]
!   PSRF     - surface pressure [Pa]
!   RFIFLX   - frozen runoff [kg m-2 s-1]
!   RNFFLX   - liquid runoff [kg m-2 s-1]
!   SALFLX   - salt flux received by ocean [kg m-2 s-1]
!   SALRLX   - restoring salt flux received by ocean [kg m-2 s-1]
!   SBOT     - bottom salinity [g kg-1]
!   SEALV    - sea level [m]
!   SLVSQ    - sea level squared [m2]
!   SFL      - salt flux [kg m-2 s-1]
!   SOP      - solid precipitation [kg m-2 s-1]
!   SIGMX    - mixed layer density [kg m-3]
!   SSS      - ocean surface salinity [g kg-1]
!   SSSSQ    - ocean surface salinity squared [g2 kg-2]
!   SST      - ocean surface temperature [degC]
!   SSTSQ    - ocean surface temperature squared [degC2]
!   SURFLX   - heat flux received by ocean [W m-2]
!   SURRLX   - restoring heat flux received by ocean [W m-2]
!   SWA      - short-wave heat flux [W m-2]
!   T20D     - 20C isoterm depth [m]
!   TAUX     - momentum flux received by ocean x-component [N m-2]
!   TAUY     - momentum flux received by ocean y-component [N m-2]
!   TBOT     - bottom temperature [degC]
!   TICE     - ice temperature [degC]
!   TSRF     - surface temperature [degC]
!   UB       - barotropic velocity x-component [m s-1]
!   UICE     - ice velocity x-component [m s-1]
!   USTAR    - friction velocity [m s-1]
!   USTAR3   - friction velocity cubed [m3 s-3]
!   VB       - barotropic velocity y-component [m s-1]
!   VICE     - ice velocity y-component [m s-1]
!   ZTX      - wind stress x-component [N m-2]
!   BFSQ     - buoyancy frequency squared [s-1]
!   DIFDIA   - diapycnal diffusivity [log10(m2 s-1)]
!   DIFINT   - layer interface diffusivity [log10(m2 s-1)]
!   DIFISO   - isopycnal diffusivity [log10(m2 s-1)]
!   DP       - layer pressure thickness [Pa]
!   DZ       - layer thickness [m]
!   SALN     - salinity [g kg-1]
!   TEMP     - temperature [degC]
!   TRC      - tracer []
!   UFLX     - mass flux in x-direction [kg s-1]
!   UTFLX    - heat flux in x-direction [W]
!   USFLX    - salt flux in x-direction [kg s-1]
!   UMFLTD   - mass flux due to thickness diffusion in x-direction [kg s-1]
!   UTFLTD   - heat flux due to thickness diffusion in x-direction [W]
!   UTFLLD   - heat flux due to lateral diffusion in x-direction [W]
!   USFLTD   - salt flux due to thickness diffusion in x-direction [kg s-1]
!   USFLLD   - salt flux due to lateral diffusion in x-direction [kg s-1]
!   UVEL     - velocity x-component [m s-1]
!   VFLX     - mass flux in y-direction [kg s-1]
!   VTFLX    - heat flux in y-direction [W]
!   VSFLX    - salt flux in y-direction [kg s-1]
!   VMFLTD   - mass flux due to thickness diffusion in y-direction [kg s-1]
!   VTFLTD   - heat flux due to thickness diffusion in y-direction [W]
!   VTFLLD   - heat flux due to lateral diffusion in y-direction [W]
!   VSFLTD   - salt flux due to thickness diffusion in y-direction [kg s-1]
!   VSFLLD   - salt flux due to lateral diffusion in y-direction [kg s-1]
!   VVEL     - velocity x-component [m s-1]
!   WFLX     - vertical mass flux [kg s-1]
!   WFLX2    - vertical mass flux squared [kg2 s-2]
!   PV       - potential vorticity [m-1 s-1]
!   TKE      - turbulent kinetic energy [m2 s-2]
!   GLS_PSI  - generic length scale [m2 s-3]
!   IDLAGE   - ideal age [year]
!   MMFLXL   - meridional overturning circ. (MOC) on isopycnic layers [kg s-1]
!   MMFLXD   - MOC on z-levels [kg s-1]
!   MMFTDL   - MOC due to thickness diffusion on isopycnic layers [kg s-1]
!   MMFTDD   - MOC due to thickness diffusion on z-levels [kg s-1]
!   MHFLX    - meridional heat flux [W]
!   MHFTD    - meridional heat flux due to thickness diffusion [W]
!   MHFLD    - meridional heat flux due to lateral diffusion [W]
!   MSFLX    - meridional salt flux [kg s-1]
!   MSFTD    - meridional salt flux due to thickness diffusion [kg s-1]
!   MSFLD    - meridional salt flux due to lateral diffusion [kg s-1]
!   VOLTR    - section transports [kg s-1]
!   MASSGS   - global sum of mass [kg]
!   VOLGS    - global sum of volume [m3]
!   SALNGA   - global average temperature [degC]
!   TEMPGA   - global average temperature [degC]
!   SSSGA    - global average sea surface salinity [g kg-1]
!   SSTGA    - global average sea surface temperature [degC]
!
&DIAPHY
  GLB_FNAMETAG = 'hd','hm','hy'
  GLB_AVEPERIO = 1,  30, 365
  GLB_FILEFREQ = 30,  30, 365
  GLB_COMPFLAG = 0,   0,   0
  GLB_NCFORMAT = 1,   1,   1
  H2D_ABSWND   = 0,   4,   0
  H2D_ALB      = 0,   0,   0
  H2D_BTMSTR   = 0,   4,   0
  H2D_BRNFLX   = 0,   4,   0
  H2D_BRNPD    = 0,   4,   0
  H2D_DFL      = 0,   0,   0
  H2D_EVA      = 0,   4,   0
  H2D_FICE     = 0,   4,   0
  H2D_FMLTFZ   = 0,   4,   0
  H2D_HICE     = 0,   0,   0
  H2D_HMLTFZ   = 0,   4,   0
  H2D_HSNW     = 0,   0,   0
  H2D_IAGE     = 0,   0,   0
  H2D_IDKEDT   = 0,   4,   0
  H2D_LIP      = 0,   4,   0
  H2D_MAXMLD   = 4,   4,   0
  H2D_MLD      = 0,   4,   0
  H2D_MLDU     = 0,   0,   0
  H2D_MLDV     = 0,   0,   0
  H2D_MLTS     = 4,   4,   0
  H2D_MLTSMN   = 0,   4,   0
  H2D_MLTSMX   = 0,   4,   0
  H2D_MLTSSQ   = 0,   4,   0
  H2D_MTKEUS   = 0,   4,   0
  H2D_MTKENI   = 0,   4,   0
  H2D_MTKEBF   = 0,   4,   0
  H2D_MTKERS   = 0,   4,   0
  H2D_MTKEPE   = 0,   4,   0
  H2D_MTKEKE   = 0,   4,   0
  H2D_MTY      = 0,   4,   0
  H2D_MXLU     = 0,   4,   0
  H2D_MXLV     = 0,   4,   0
  H2D_NSF      = 0,   4,   0
  H2D_PBOT     = 0,   4,   0
  H2D_PSRF     = 0,   4,   0
  H2D_RFIFLX   = 0,   4,   0
  H2D_RNFFLX   = 0,   4,   0
  H2D_SALFLX   = 0,   4,   0
  H2D_SALRLX   = 0,   4,   0
  H2D_SBOT     = 0,   4,   0
  H2D_SEALV    = 4,   4,   0
  H2D_SLVSQ    = 0,   4,   0
  H2D_SFL      = 0,   4,   0
  H2D_SOP      = 0,   4,   0
  H2D_SIGMX    = 0,   4,   0
  H2D_SSS      = 4,   4,   0
  H2D_SSSSQ    = 4,   4,   0
  H2D_SST      = 4,   4,   0
  H2D_SSTSQ    = 4,   4,   0
  H2D_SURFLX   = 0,   4,   0
  H2D_SURRLX   = 0,   4,   0
  H2D_SWA      = 0,   4,   0
  H2D_T20D     = 4,   4,   0
  H2D_TAUX     = 0,   4,   0
  H2D_TAUY     = 0,   4,   0
  H2D_TBOT     = 0,   4,   0
  H2D_TICE     = 0,   0,   0
  H2D_TSRF     = 0,   0,   0
  H2D_UB       = 0,   4,   0
  H2D_UICE     = 0,   0,   0
  H2D_USTAR    = 0,   4,   0
  H2D_USTAR3   = 0,   4,   0
  H2D_VB       = 0,   4,   0
  H2D_VICE     = 0,   0,   0
  H2D_ZTX      = 0,   4,   0
  LYR_BFSQ     = 0,   4,   0
  LYR_DIFDIA   = 0,   4,   0
  LYR_DIFINT   = 0,   4,   0
  LYR_DIFISO   = 0,   4,   0
  LYR_DP       = 0,   4,   0
  LYR_DZ       = 0,   4,   0
  LYR_SALN     = 0,   4,   0
  LYR_TEMP     = 0,   4,   0
  LYR_TRC      = 0,   0,   0
  LYR_UFLX     = 0,   4,   0
  LYR_UTFLX    = 0,   4,   0
  LYR_USFLX    = 0,   4,   0
  LYR_UMFLTD   = 0,   0,   4
  LYR_UTFLTD   = 0,   0,   4
  LYR_UTFLLD   = 0,   0,   4
  LYR_USFLTD   = 0,   0,   4
  LYR_USFLLD   = 0,   0,   4
  LYR_UVEL     = 0,   4,   0
  LYR_VFLX     = 0,   4,   0
  LYR_VTFLX    = 0,   4,   0
  LYR_VSFLX    = 0,   4,   0
  LYR_VMFLTD   = 0,   4,   0
  LYR_VTFLTD   = 0,   4,   0
  LYR_VTFLLD   = 0,   4,   0
  LYR_VSFLTD   = 0,   4,   0
  LYR_VSFLLD   = 0,   4,   0
  LYR_VVEL     = 0,   4,   0
  LYR_WFLX     = 0,   4,   0
  LYR_WFLX2    = 0,   4,   0
  LYR_PV       = 0,   4,   0
  LYR_TKE      = 0,   4,   0
  LYR_GLS_PSI  = 0,   4,   0
  LYR_IDLAGE   = 0,   4,   0
  LVL_BFSQ     = 0,   4,   0
  LVL_DIFDIA   = 0,   4,   0
  LVL_DIFINT   = 0,   4,   0
  LVL_DIFISO   = 0,   4,   0
  LVL_DZ       = 0,   4,   0
  LVL_SALN     = 0,   4,   0
  LVL_TEMP     = 0,   4,   0
  LVL_TRC      = 0,   0,   0
  LVL_UFLX     = 0,   4,   0
  LVL_UTFLX    = 0,   4,   0
  LVL_USFLX    = 0,   4,   0
  LVL_UMFLTD   = 0,   0,   4
  LVL_UTFLTD   = 0,   0,   4
  LVL_UTFLLD   = 0,   0,   4
  LVL_USFLTD   = 0,   0,   4
  LVL_USFLLD   = 0,   0,   4
  LVL_UVEL     = 0,   4,   0
  LVL_VFLX     = 0,   4,   0
  LVL_VTFLX    = 0,   4,   0
  LVL_VSFLX    = 0,   4,   0
  LVL_VMFLTD   = 0,   0,   4
  LVL_VTFLTD   = 0,   0,   4
  LVL_VTFLLD   = 0,   0,   4
  LVL_VSFLTD   = 0,   0,   4
  LVL_VSFLLD   = 0,   0,   4
  LVL_VVEL     = 0,   4,   0
  LVL_WFLX     = 0,   4,   0
  LVL_WFLX2    = 0,   4,   0
  LVL_PV       = 0,   4,   0
  LVL_TKE      = 0,   4,   0
  LVL_GLS_PSI  = 0,   4,   0
  LVL_IDLAGE   = 0,   4,   0
  MSC_MMFLXL   = 0,   4,   0
  MSC_MMFLXD   = 0,   4,   0
  MSC_MMFTDL   = 0,   4,   0
  MSC_MMFTDD   = 0,   4,   0
  MSC_MHFLX    = 0,   4,   0
  MSC_MHFTD    = 0,   4,   0
  MSC_MHFLD    = 0,   4,   0
  MSC_MSFLX    = 0,   4,   0
  MSC_MSFTD    = 0,   4,   0
  MSC_MSFLD    = 0,   4,   0
  MSC_VOLTR    = 0,   4,   0
  MSC_MASSGS   = 0,   4,   0
  MSC_VOLGS    = 0,   4,   0
  MSC_SALNGA   = 0,   4,   0
  MSC_TEMPGA   = 0,   4,   0
  MSC_SSSGA    = 0,   4,   0
  MSC_SSTGA    = 0,   4,   0
/

! NAMELIST FOR iHAMOCC OPTIONS
!
! CONTENTS:
!
! ATM_CO2   : Atmospheric CO2 concentration [ppmv]
! FEDEPFILE : File name (incl. full path) for iron (dust) deposition data
! DO_RIVINPT: Logical switch to activate riverine input
! RIVINFILE : File name (incl. full path) for riverine input data
! DO_NDEP   : Logical switch to activate N-deposition
! NDEPFILE  : File name (incl. full path) for atmopheric N-deposition data
! INIXXX    : Initial condition file for iHAMOCC, where XXX=DIC, ALK, PO4,
!             OXY, NO3, SIL, D13C, and D14C
&BGCNML
  ATM_CO2    = 367.0
  FEDEPFILE  = '/work/shared/noresm/inputdata/ocn/blom/bndcon/dustdep_mhw2006_tnx1v4_20171107.nc'
  DO_RIVINPT = .true.
  RIVINFILE  = '/work/shared/noresm/inputdata/ocn/blom/bndcon/river_nutrients_GNEWS2000_tnx1v4_20170820.nc' 
  DO_NDEP    = .true.
  NDEPFILE   = '/work/shared/noresm/inputdata/ocn/blom/bndcon/ndep_2000_CMIP6_tnx1v4_20200826.nc'
  INIDIC     = '/work/shared/noresm/inputdata/ocn/blom/inicon/glodapv2_Ct_preind_OMIPinit_20171107.nc'
  INIALK     = '/work/shared/noresm/inputdata/ocn/blom/inicon/glodapv2_At_OMIPinit_20171107.nc'
  INIPO4     = '/work/shared/noresm/inputdata/ocn/blom/inicon/woa13_phosphate_OMIPinit_20171107.nc'
  INIOXY     = '/work/shared/noresm/inputdata/ocn/blom/inicon/woa13_oxygen_OMIPinit_20171107.nc'
  ININO3     = '/work/shared/noresm/inputdata/ocn/blom/inicon/woa13_nitrate_OMIPinit_20171107.nc'
  INISIL     = '/work/shared/noresm/inputdata/ocn/blom/inicon/woa13_silicate_OMIPinit_20171107.nc'
  INID13C    = ''
  INID14C    = ''
/

! IO-NAMELIST FOR DIAGNOSTIC iHAMOCC OUTPUT
!
! Namelist acronyms:
!   GLB_     - global parameters i.e. valid for entire output group
!   SRF_     - surface variables (includes some non-surface 2d fields)
!   FLX_     - interior ocean particle fuxes at different depth (2d)
!   INT_     - vertically integrated fields (2d)
!   LYR_     - 3d fields with sigma layers as vertical coordinate
!   LVL_     - 3d fields with levitus levels as vertical coordinate
!   SDM_     - 3d sediment fields
!   BUR_     - 2d fields of sediment burial
!
! Global parameters:
!   FNAMETAG - tag used in file name (c10) 
!   AVEPERIO - average period in days (i) 
!   FILEFREQ - how often to start a new file in days (i) 
!   COMPFLAG - switch for compressed/uncompressed output (i) 
!   NCFORMAT - netcdf format (valid arguments are 0 for classic,
!              1 for 64-bit offset and 2 for netcdf4/hdf5 format)
!   INVENTORY- how often to write an inventory of tracers to ocean 
!              output file
!  
! Arguments for diagnostic variables:
!   0        - variable is not written
!   2        - variable is written as int2 with scale factor and offset
!   4        - variable is written as real4
!   8        - variable is written as real8
!
! Output variables:
!
! Main BGC tracer fields (SFR, LYR or LVL)
!   DIC            - Dissolved carbon (dissic) [mol C m-3]
!   ALKALI         - Alkalinity (talk) [eq m-3]
!   OXYGEN         - Oxygen (o2) [mol O2 m-3]
!   ANO3           - Nitrate (no3) [mol N m-3]
!   PHOSPH         - Phosphorus (po4) [mol P m-3]
!   IRON           - Dissolved iron (dfe) [mol Fe m-3]
!   SILICA         - Silicate (si) [mol Si m-3]
!   PHYTO          - Phytoplankton (phyc) [mol C m-3]
!   NATDIC         - Natural dissolved carbon (natdissic) [mol C m-3]
!   NATALKALI      - Natural alkalinity (nattalk) [eq m-3]
!
! Other 3d tracer or diagnostic variables (LYR or LVL)
!   DP             - Layer thickness (pddpo) [m]
!   PREFO2         - preformed oxygen (p_o2) [mol O2 m-3]
!   O2SAT          - Saturated oxygen (satoxy) [mol O2 m-3]
!   PREFPO4        - preformed phosphate (p_po4) [mol PO4 m-3]
!   PREFALK        - preformed alkalinity (p_talk) [eq m-3]
!   PREFDIC        - preformed DIC (p_dic) [mol C m-3]
!   DICSAT         - saturated DIC (dic_sat) [mol C m-3]
!   POC            - Detrius (detoc) [mol C m-3]
!   CALC           - CaCO3 shells (calc) [mol C m-3]
!   OPAL           - Opal shells (opal) [mol Si m-3]
!   DOC            - Dissolved organic carbon (dissoc) [mol C m-3]
!   GRAZER         - Zooplankton (zooc) [mol C m-3]
!   PHOSY          - Primary production (pp) [mol C m-3 s-1]
!   CO3            - Carbonate ions (co3) [mol C m-3]
!   N2O            - Nitrous oxide concentration [mol N2O m-3]
!   PH             - pH (ph) [-log10([h+])]
!   OMEGAA         - Aragonite saturation state (omegaa) [1]
!   OMEGAC         - Calcite saturation state (omegac) [1]
!   NATDIC         - Natural dissolved carbon (natdissic) [mol C m-3]
!   NATALKALI      - Natural alkalinity (nattalk) [eq m-3]
!   NATCO3         - Natural carbonate ion concentration (natco3) [mol C m-3]
!   NATCALC        - Natural CaCO3 shells (natcalc) [mol C m-3]
!   NATPH          - Natural pH (natph) [-log10([h+])]
!   NATOMEGAA      - Natural aragonite saturation state (natomegaa) [1]
!   NATOMEGAC      - Natural calcite saturation state (natomegac) [1]
!   DIC13          - Dissolved C13 (dissic13) [mol C m-3]
!   DIC14          - Dissolved C14 (dissic14) [mol C m-3]
!   D13C           - delta 13C of DIC [1]
!   D14C           - delta 14C of DIC [1]
!   BIGD14C        - Delta 14C of DIC [1]
!   POC13          - Particulate organic carbon 13 (detoc13) [mol C m-3]
!   DOC13          - Dissolved organic carbon 13 (dissoc13) [mol C m-3]
!   CALC13         - Particulate inorganic carbon 13 (calc13) [mol C m-3]
!   PHYTO13        - 13C of phytoplankton biomass (phyc13) [mol C m-3]
!   GRAZER13       - 13C of zootoplankton biomass (zooc13) [mol C m-3]
!   CFC11          - CFC11 concentration [mol CFC11 m-3]
!   CFC12          - CFC12 concentration [mol CFC12 m-3]
!   SF6            - SF6 concentration [mol SF6 m-3]
!   NOS            - Number of sinking particles (aggregate scheme, nos) [cm-3]
!   WPHY           - Mass sinking velocity (aggregate scheme) [m d-1]
!   WNOS           - Number sinking velocity (aggregate scheme) [m d-1]
!   EPS            - Epsilon exponent (aggregate scheme) [1]
!   ASIZE          - Average particle size (aggregate scheme)
!
! 2d-fields (SFR)
!   EXPORT         - Export production (epc100) [mol C m-2 s-1]
!   EXPOCA         - Ca export production (epcalc100) [mol Ca m-2 s-1]
!   EXPOSI         - Si export production (epsi100) [mol Si m-2 s-1]
!   PCO2           - Surface PCO2 (spco2) [uatm]
!   KWCO2          - kwco2 x solubility
!   CO2FXD         - Downward CO2 flux (co2fxd) [kg C m-2 s-1]
!   CO2FXU         - Upward CO2 flux (co2fxu) [kg C m-2 s-1]
!   NIFLUX         - Nitrogen flux (fgn2) [mol N2 m-2 s-1]
!   OXFLUX         - Oxygen flux (fgo2) [mol O2 m-2 s-1]
!   N2OFX          - Nitrous oxide flux [mol N2O m-2 s-1]
!   DMSFLUX        - DMS flux (dmsflux) [mol DMS m-2 s-1]
!   DMS            - surface DMS concentration (dms) [mol DMS m-3]
!   DMSPROD        - DMS production (dmsprod) [mol DMS m-2 s-1]
!   DMS_BAC        - DMS bacterial consuption (dms_bac) [mol DMS m-2 s-1]
!   DMS_UV         - DMS decomposition by UV (dms_uv) [mol DMS m-2 s-1]
!   ATMCO2         - Atmospheric CO2 (atmco2) [ppm]
!   ATMO2          - Atmospheric O2 (atmo2) [ppm]
!   ATMN2          - Atmospheric N2 (atmn2) [ppm]
!   NATPCO2        - Natural surface PCO2 (spco2) [uatm]
!   NATCO2FX       - Natural net CO2 flux (natco2fx) [kg C m-2 s-1]
!   CO213FXD       - Downward 13CO2 flux (co213fxd) [kg C m-2 s-1]
!   CO213FXU       - Upward 13CO2 flux (co213fxu) [kg C m-2 s-1]
!   CO214FXD       - Downward 14CO2 flux (co214fxd) [kg C m-2 s-1]
!   CO214FXU       - Upward 14CO2 flux (co214fxu) [kg C m-2 s-1]
!   CFC11          - CFC11 flux [mol CFC11 m-2 s-1]
!   CFC12          - CFC12 flux [mol CFC12 m-2 s-1]
!   SF6            - SF6 flux [mol SF6 m-2 s-1]
!
! Vertically integrated fields (INT)
!   PHOSY          - Vertically integrated primary production [mol C m-3 s-1]
!   NFIX           - Vertically integrated nitrogen fixation
!   DNIT           - Vertically integrated denitrification
!
! Particle fluxes (FLX, where ****=0100,0500,1000,2000,4000, or _BOT)
!   CARFLX****     - POC flux at **** metres depth [mol C m-2 s-1]
!   BSIFLX****     - Biogenic silica flux at **** metres depth [mol Si m-2 s-1]
!   CALFLX****     - Calcium carbonate flux at **** metres depth [mol C m-2 s-1]
!
! Sediment fields (SDM)
!   POWAIC         - (powdic) [mol C m-3]
!   POWAAL         - (powalk) [eq m-3]
!   POWAPH         - (powpho) [eq m-3]
!   POWAOX         - (powox) [mol O2 m-3]
!   POWN2          - (pown2) [mol N2 m-3]
!   POWNO3         - (powno3)[mol N m-3]
!   POWASI         - (powsi) [mol Si m-3]
!   SSSO12         - (ssso12) [mol m-3]
!   SSSSIL         - (ssssil) [mol Si m-3]
!   SSSC12         - (sssc12) [mol C m-3]
!   SSSTER         - (ssster) [mol m-3]
!
! Burial fields (BUR)
!   SSSO12         - Solid sediment organic carbon [mol P m-2]
!   SSSSIL         - Solid sediment silica [mol Si m-2]
!   SSSC12         - Solid sediment inorganic carbon [mol C m-2]
!   SSSTER         - 
!
&DIABGC
  GLB_FNAMETAG  = 'hbgcd','hbgcm','hbgcy'
  GLB_AVEPERIO  = 1,30,365
  GLB_FILEFREQ  = 30,30,365
  GLB_COMPFLAG  = 0, 0, 0
  GLB_NCFORMAT  = 1, 1, 1
  GLB_INVENTORY = 0, 1, 0
  SRF_PHOSPH    = 0, 2, 2
  SRF_OXYGEN    = 0, 2, 2
  SRF_IRON      = 0, 2, 2
  SRF_ANO3      = 0, 2, 2
  SRF_ALKALI    = 4, 2, 2
  SRF_SILICA    = 0, 2, 2
  SRF_DIC       = 4, 2, 2
  SRF_PHYTO     = 4, 2, 2
  SRF_EXPORT    = 0, 2, 2
  SRF_EXPOSI    = 0, 2, 2
  SRF_EXPOCA    = 0, 2, 2
  SRF_KWCO2     = 0, 2, 2
  SRF_PCO2      = 0, 2, 2
  SRF_CO2FXD    = 4, 2, 2
  SRF_CO2FXU    = 4, 2, 2
  SRF_OXFLUX    = 0, 2, 2
  SRF_NIFLUX    = 0, 2, 2
  SRF_N2OFX     = 0, 0, 2
  SRF_DMSFLUX   = 0, 2, 2
  SRF_DMS       = 0, 2, 2
  SRF_DMSPROD   = 0, 2, 2
  SRF_DMS_BAC   = 0, 2, 2
  SRF_DMS_UV    = 0, 2, 2
  SRF_ATMCO2    = 0, 2, 2
  SRF_ATMO2     = 0, 2, 2
  SRF_ATMN2     = 0, 2, 2
  SRF_NATDIC    = 0, 2, 2
  SRF_NATALKALI = 0, 2, 2
  SRF_NATPCO2   = 0, 2, 2
  SRF_NATCO2FX  = 0, 2, 2
  SRF_CO213FXD  = 0, 2, 2
  SRF_CO213FXU  = 0, 2, 2
  SRF_CO214FXD  = 0, 2, 2
  SRF_CO214FXU  = 0, 2, 2
  SRF_CFC11     = 0, 2, 2
  SRF_CFC12     = 0, 2, 2
  SRF_SF6       = 0, 2, 2
  INT_PHOSY     = 4, 2, 2
  INT_NFIX      = 0, 2, 2
  INT_DNIT      = 0, 2, 2
  FLX_CAR0100   = 0, 2, 2
  FLX_CAR0500   = 0, 2, 2
  FLX_CAR1000   = 0, 2, 2
  FLX_CAR2000   = 0, 2, 2
  FLX_CAR4000   = 0, 2, 2
  FLX_CAR_BOT   = 0, 2, 2
  FLX_BSI0100   = 0, 2, 2
  FLX_BSI0500   = 0, 2, 2
  FLX_BSI1000   = 0, 2, 2
  FLX_BSI2000   = 0, 2, 2
  FLX_BSI4000   = 0, 2, 2
  FLX_BSI_BOT   = 0, 2, 2
  FLX_CAL0100   = 0, 2, 2
  FLX_CAL0500   = 0, 2, 2
  FLX_CAL1000   = 0, 2, 2
  FLX_CAL2000   = 0, 2, 2
  FLX_CAL4000   = 0, 2, 2
  FLX_CAL_BOT   = 0, 2, 2
  LYR_PHYTO     = 0, 0, 2
  LYR_GRAZER    = 0, 0, 2
  LYR_DOC       = 0, 0, 2
  LYR_PHOSY     = 0, 0, 2
  LYR_PHOSPH    = 0, 0, 2
  LYR_OXYGEN    = 0, 0, 4
  LYR_IRON      = 0, 0, 2
  LYR_ANO3      = 0, 0, 2
  LYR_ALKALI    = 0, 0, 2
  LYR_SILICA    = 0, 0, 2
  LYR_DIC       = 0, 0, 2
  LYR_POC       = 0, 0, 2
  LYR_CALC      = 0, 0, 2
  LYR_OPAL      = 0, 0, 2
  LYR_CO3       = 0, 0, 2
  LYR_N2O       = 0, 0, 0
  LYR_PH        = 0, 0, 2
  LYR_OMEGAC    = 0, 0, 2
  LYR_OMEGAA    = 0, 0, 2
  LYR_PREFO2    = 0, 0, 4
  LYR_O2SAT     = 0, 0, 4
  LYR_PREFPO4   = 0, 0, 2
  LYR_PREFALK   = 0, 0, 2
  LYR_PREFDIC   = 0, 0, 2
  LYR_DICSAT    = 0, 0, 2
  LYR_NATDIC    = 0, 0, 2
  LYR_NATALKALI = 0, 0, 2
  LYR_NATCO3    = 0, 0, 2
  LYR_NATCALC   = 0, 0, 2
  LYR_NATPH     = 0, 0, 2
  LYR_NATOMEGAC = 0, 0, 2
  LYR_NATOMEGAA = 0, 0, 2
  LYR_DIC13     = 0, 0, 2
  LYR_DIC14     = 0, 0, 2
  LYR_D13C      = 0, 0, 2
  LYR_D14C      = 0, 0, 2
  LYR_BIGD14C   = 0, 0, 2
  LYR_POC13     = 0, 0, 2
  LYR_DOC13     = 0, 0, 2
  LYR_CALC13    = 0, 0, 2
  LYR_PHYTO13   = 0, 0, 2
  LYR_GRAZER13  = 0, 0, 2
  LYR_CFC11     = 0, 0, 2
  LYR_CFC12     = 0, 0, 2
  LYR_SF6       = 0, 0, 2
  LYR_NOS       = 0, 0, 2
  LYR_WPHY      = 0, 0, 2
  LYR_WNOS      = 0, 0, 2
  LYR_EPS       = 0, 0, 0
  LYR_ASIZE     = 0, 0, 0
  LYR_DP        = 0, 2, 2
  LVL_PHYTO     = 0, 2, 2
  LVL_GRAZER    = 0, 2, 2
  LVL_DOC       = 0, 2, 2
  LVL_PHOSY     = 0, 2, 2
  LVL_PHOSPH    = 0, 2, 2
  LVL_OXYGEN    = 0, 4, 4
  LVL_IRON      = 0, 2, 2
  LVL_ANO3      = 0, 2, 2
  LVL_ALKALI    = 0, 2, 2
  LVL_SILICA    = 0, 2, 2
  LVL_DIC       = 0, 2, 2
  LVL_POC       = 0, 2, 2
  LVL_CALC      = 0, 2, 2
  LVL_OPAL      = 0, 2, 2
  LVL_CO3       = 0, 2, 2
  LVL_N2O       = 0, 0, 2
  LVL_PH        = 0, 2, 2
  LVL_OMEGAC    = 0, 2, 2
  LVL_OMEGAA    = 0, 2, 2
  LVL_PREFO2    = 0, 4, 4
  LVL_O2SAT     = 0, 4, 4
  LVL_PREFPO4   = 0, 2, 2
  LVL_PREFALK   = 0, 2, 2
  LVL_PREFDIC   = 0, 2, 2
  LVL_DICSAT    = 0, 2, 2
  LVL_NATDIC    = 0, 2, 2
  LVL_NATALKALI = 0, 2, 2
  LVL_NATCO3    = 0, 2, 2
  LVL_NATCALC   = 0, 2, 2
  LVL_NATPH     = 0, 2, 2
  LVL_NATOMEGAC = 0, 2, 2
  LVL_NATOMEGAA = 0, 2, 2
  LVL_DIC13     = 0, 2, 2
  LVL_DIC14     = 0, 2, 2
  LVL_D13C      = 0, 2, 2
  LVL_POC13     = 0, 2, 2
  LVL_DOC13     = 0, 2, 2
  LVL_CALC13    = 0, 2, 2
  LVL_PHYTO13   = 0, 2, 2
  LVL_GRAZER13  = 0, 2, 2
  LVL_CFC11     = 0, 2, 2
  LVL_CFC12     = 0, 2, 2
  LVL_SF6       = 0, 2, 2
  LVL_NOS       = 0, 2, 2
  LVL_WPHY      = 0, 2, 2
  LVL_WNOS      = 0, 2, 2
  LVL_EPS       = 0, 0, 0
  LVL_ASIZE     = 0, 0, 0
  SDM_POWAIC    = 0, 0, 2
  SDM_POWAAL    = 0, 0, 2
  SDM_POWAPH    = 0, 0, 2
  SDM_POWAOX    = 0, 0, 2
  SDM_POWN2     = 0, 0, 2
  SDM_POWNO3    = 0, 0, 2
  SDM_POWASI    = 0, 0, 2
  SDM_SSSO12    = 0, 0, 2
  SDM_SSSSIL    = 0, 0, 2
  SDM_SSSC12    = 0, 0, 2
  SDM_SSSTER    = 0, 0, 2
  BUR_SSSO12    = 0, 0, 2
  BUR_SSSSIL    = 0, 0, 2
  BUR_SSSC12    = 0, 0, 2
  BUR_SSSTER    = 0, 0, 2
/
