!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! NEMO/TOP :   Reference namelist 
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!!               - tracer run information                (namtrc_run)
!!               - tracer definition                     (namtrc    )
!!               - tracer data initialisation            (namtrc_dta)
!!               - tracer advection                      (namtrc_adv)
!!               - tracer lateral diffusion              (namtrc_ldf)
!!               - tracer vertical physics               (namtrc_zdf)
!!               - tracer newtonian damping              (namtrc_dmp)
!!               - dynamical tracer trends               (namtrc_trd)
!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!-----------------------------------------------------------------------
&namtrc_run      !   run information
!-----------------------------------------------------------------------
   ln_top_euler  = .false.   !  use Euler time-stepping for TOP
   ln_rsttr      = .false.   !  start from a restart file (T) or not (F)
   nn_rsttr      = 0    !  restart control = 0 initial time step is not compared to the restart file value
                             !                  = 1 do not use the value in the restart file
                             !                  = 2 calendar parameters read in the restart file
   cn_trcrst_in  = "GITLAB_CICD_00000216_restart_trc"   !  suffix of pass. sn_tracer restart name (input)
   cn_trcrst_indir = "."           !  directory from which to read input passive tracer restarts
   cn_trcrst_out = "restart_trc"   !  suffix of pass. sn_tracer restart name (output)
   cn_trcrst_outdir = "."          !  directory to which to write output passive tracer restarts
   ln_enforce_consistency = .true. ! Check namelist top for errors in tracer order (initialization + boundary conditions)
/
!-----------------------------------------------------------------------
&namtrc          !   tracers definition
!-----------------------------------------------------------------------
   jp_bgc        =  26          !  Number of passive tracers of the BGC model
   !
   ln_pisces     =  .false.     !  Run PISCES BGC model 
   ln_my_trc     =  .true.      !  Run MY_TRC BGC model
   ln_age        =  .false.     !  Run the sea water age tracer
   ln_cfc11      =  .false.     !  Run the CFC11 passive tracer
   ln_cfc12      =  .false.     !  Run the CFC12 passive tracer
   ln_sf6        =  .false.     !  Run the SF6 passive tracer
   ln_c14        =  .false.     !  Run the Radiocarbon passive tracer
   ln_benthic    =  .true.      !  Run BAMHBI benthic model
   !
   ln_trcdta     =  .true.   !  Initialisation from data input file (T) or not (F)
   ln_trcdmp     =  .false.  !  add a damping termn (T) or not (F)
   ln_trcdmp_clo =  .false.  !  damping term (T) or not (F) on closed seas
   ln_trcbc      =  .true.   !  Surface, Lateral or Open Boundaries conditions can be disabled here
   ln_trcais     =  .false.  !  Antarctic Ice Sheet nutrient supply
   !
   jp_dia3d      = 0         ! Number of 3D diagnostic variables
   jp_dia2d      = 0         ! Number of 2D diagnostic variables
   !_____________!___________!_________________________________________!____________!________________!
   !             !    name   !           title of the field            !   units    ! init from file !  sbc   !  cbc   !  obc   !  ais
!  sn_tracer(1)  = 'tracer  ', 'Tracer  Concentration                 ',   ' - '    ,   .false.
   sn_tracer(1)   = 'CFL'    ,  'variable  01'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(2)   = 'NFL'    ,  'variable  02'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(3)   = 'CEM'    ,  'variable  03'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(4)   = 'NEM'    ,  'variable  04'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(5)   = 'CDI'    ,  'variable  05'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(6)   = 'NDI'    ,  'variable  06'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(7)   = 'MIC'    ,  'variable  07'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(8)   = 'MES'    ,  'variable  08'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(9)   = 'BAC'    ,  'variable  09'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(10)  = 'DCL'    ,  'variable  10'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(11)  = 'DNL'    ,  'variable  11'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(12)  = 'DCS'    ,  'variable  12'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(13)  = 'DNS'    ,  'variable  13'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(14)  = 'NOS'    ,  'variable  14'                           ,  'mmol/m3'   ,   .true.   , .true. , .true. , .true. , .false.
   sn_tracer(15)  = 'NHS'    ,  'variable  15'                           ,  'mmol/m3'   ,   .true.   , .true. , .true. , .true. , .false.
   sn_tracer(16)  = 'SIO'    ,  'variable  16'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(17)  = 'DOX'    ,  'variable  17'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(18)  = 'DIC'    ,  'variable  18'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(19)  = 'ODU'    ,  'variable  19'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(20)  = 'POC'    ,  'variable  20'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(21)  = 'PON'    ,  'variable  21'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(22)  = 'SID'    ,  'variable  22'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(23)  = 'PHO'    ,  'variable  26'                           ,  'mmol/m3'   ,   .true.   , .true. , .true. , .true. , .false.
   sn_tracer(24)  = 'GEL'    ,  'variable  24'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(25)  = 'NOC'    ,  'variable  25'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
   sn_tracer(26)  = 'CHA'    ,  'variable  27'                           ,  'mmol/m3'   ,   .true.   , .false., .true. , .true. , .false.
/
!-----------------------------------------------------------------------
&namage          !   AGE 
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namtrc_dta      !    Initialisation from data input file
!-----------------------------------------------------------------------
!                !  file name        ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
!                !                   !  (if <0  months)  !   name    !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  ! filename      !
   sn_trcdta(1)  = 'bioinit'      ,   -12      ,  'CFL'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(1)   = 1.0
   sn_trcdta(2)  = 'bioinit'      ,   -12      ,  'NFL'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(2)   = 1.0
   sn_trcdta(3)  = 'bioinit'      ,   -12      ,  'CEM'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(3)   = 1.0
   sn_trcdta(4)  = 'bioinit'      ,   -12      ,  'NEM'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(4)   = 1.0
   sn_trcdta(5)  = 'bioinit'      ,   -12      ,  'CDI'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(5)   = 1.0
   sn_trcdta(6)  = 'bioinit'      ,   -12      ,  'NDI'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(6)   = 1.0
   sn_trcdta(7)  = 'bioinit'      ,   -12      ,  'MIC'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(7)   = 1.0
   sn_trcdta(8)  = 'bioinit'      ,   -12      ,  'MES'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(8)   = 1.0
   sn_trcdta(9)  = 'bioinit'      ,   -12      ,  'BAC'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(9)   = 1.0
   sn_trcdta(10) = 'bioinit'      ,   -12      ,  'DCL'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(10)  = 1.0
   sn_trcdta(11) = 'bioinit'      ,   -12      ,  'DNL'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(11)  = 1.0
   sn_trcdta(12) = 'bioinit'      ,   -12      ,  'DCS'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(12)  = 1.0
   sn_trcdta(13) = 'bioinit'      ,   -12      ,  'DNS'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(13)  = 1.0
   sn_trcdta(14) = 'bioinit'      ,   -12      ,  'NOS'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(14)  = 1.0
   sn_trcdta(15) = 'bioinit'      ,   -12      ,  'NHS'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(15)  = 1.0
   sn_trcdta(16) = 'bioinit'      ,   -12      ,  'SIO'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(16)  = 1.0
   sn_trcdta(17) = 'bioinit'      ,   -12      ,  'DOX'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(17)  = 1.0
   sn_trcdta(18) = 'bioinit'      ,   -12      ,  'DIC'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''   ! in gher file, use 'DICs'
   rn_trfac(18)  = 1.0
   sn_trcdta(19) = 'bioinit'      ,   -12      ,  'ODU'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(19)  = 1.0
   sn_trcdta(20) = 'bioinit'      ,   -12      ,  'POC'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(20)  = 1.0
   sn_trcdta(21) = 'bioinit'      ,   -12      ,  'PON'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(21)  = 1.0
   sn_trcdta(22) = 'bioinit'      ,   -12      ,  'SID'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(22)  = 1.0
   sn_trcdta(23) = 'bioinit'      ,   -12      ,  'PHO'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(23)  = 1.0
   sn_trcdta(24) = 'bioinit'      ,   -12      ,  'GEL'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(24)  = 1.0
   sn_trcdta(25) = 'bioinit'      ,   -12      ,  'NOC'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(25)  = 1.0
   sn_trcdta(26) = 'bioinit'      ,   -12      ,  'CHA'    ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   rn_trfac(26)  = 1.0
   cn_dir        =  '../GEO_LR/'    !  directory for the location of the data files
/
!-----------------------------------------------------------------------
&namtrc_adv      !   advection scheme for passive tracer                (default: NO selection)
!-----------------------------------------------------------------------
   ln_trcadv_fct =  .true.   !  FCT scheme
      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order 
      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order 
/
!-----------------------------------------------------------------------
&namtrc_ldf      !   lateral diffusion scheme for passive tracer        (default: NO selection)
!-----------------------------------------------------------------------
!                            !  Type of the operator:
   ln_trcldf_tra   =  .true.     !  use active tracer setting
   !                         !  Coefficient (defined with namtra_ldf coefficient)
   rn_ldf_multi    = 1.          !  multiplier of aht for TRC mixing coefficient
/
!-----------------------------------------------------------------------
&namtrc_rad      !  treatment of negative concentrations 
!-----------------------------------------------------------------------
   ln_trcrad     =  .true.  !  artificially correct negative concentrations (T) or not (F)
/
!-----------------------------------------------------------------------
&namtrc_snk      !  Sedimentation of particles
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namtrc_dcy      !  Diurnal cycle
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namtrc_opt      !  light availability in the water column (PISCES ?)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namtrc_dmp      !   passive tracer newtonian damping                   (ln_trcdmp=T)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namtrc_ice      !    Representation of sea ice growth & melt effects
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namtrc_trd      !   diagnostics on tracer trends                       ('key_trdtrc')
!                          or mixed-layer trends                        ('key_trdmld_trc')
!----------------------------------------------------------------------
   nn_trd_trc    =  5475     !  time step frequency and tracers trends
   nn_ctls_trc   =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk)
   rn_ucf_trc    =   1       !  unit conversion factor (=1 -> /seconds ; =86400. -> /day)
   ln_trdmld_trc_restart = .false. !  restart for ML diagnostics
   ln_trdmld_trc_instant = .true.  !  flag to diagnose trends of instantantaneous or mean ML T/S
   ln_trdtrc( 1) = .true.
   ln_trdtrc( 2) = .true.
   ln_trdtrc(23) = .true.
/
!----------------------------------------------------------------------
&namtrc_bc       !   data for boundary conditions
!-----------------------------------------------------------------------
   cn_dir_obc    =  '../GEO_LR/'     !  root directory for the location of SURFACE data files
   sn_trcobc(1)   = 'bdy_bio.nc'             , -12       , 'CFL'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(1)   =  1.
   sn_trcobc(2)   = 'bdy_bio.nc'             , -12       , 'NFL'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(2)   =  1.
   sn_trcobc(3)   = 'bdy_bio.nc'             , -12       , 'CEM'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(3)   =  1.
   sn_trcobc(4)   = 'bdy_bio.nc'             , -12       , 'NEM'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(4)   =  1.
   sn_trcobc(5)   = 'bdy_bio.nc'             , -12       , 'CDI'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(5)   =  1.
   sn_trcobc(6)   = 'bdy_bio.nc'             , -12       , 'NDI'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(6)   =  1.
   sn_trcobc(7)   = 'bdy_bio.nc'             , -12       , 'MIC'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(7)   =  1.
   sn_trcobc(8)   = 'bdy_bio.nc'             , -12       , 'MES'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(8)   =  1.
   sn_trcobc(9)   = 'bdy_bio.nc'             , -12       , 'BAC'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(9)   =  1.
   sn_trcobc(10)  = 'bdy_bio.nc'             , -12       , 'DCL'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(10)  =  1.
   sn_trcobc(11)  = 'bdy_bio.nc'             , -12       , 'DNL'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(11)  =  1.
   sn_trcobc(12)  = 'bdy_bio.nc'             , -12       , 'DCS'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(12)  =  1.
   sn_trcobc(13)  = 'bdy_bio.nc'             , -12       , 'DNS'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(13)  =  1.
   sn_trcobc(14)  = 'bdy_bio.nc'             , -12       , 'NOS'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(14)  =  1.
   sn_trcobc(15)  = 'bdy_bio.nc'             , -12       , 'NHS'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(15)  =  1.
   sn_trcobc(16)  = 'bdy_bio.nc'             , -12       , 'SIO'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(16)  =  1.
   sn_trcobc(17)  = 'bdy_bio.nc'             , -12       , 'DOX'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(17)  =  1.
   sn_trcobc(18)  = 'bdy_bio.nc'             , -12       , 'DIC'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(18)  =  1.
   sn_trcobc(19)  = 'bdy_bio.nc'             , -12       , 'ODU'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(19)  =  1.
   sn_trcobc(20)  = 'bdy_bio.nc'             , -12       , 'POC'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(20)  =  1.
   sn_trcobc(21)  = 'bdy_bio.nc'             , -12       , 'PON'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(21)  =  1.
   sn_trcobc(22)  = 'bdy_bio.nc'             , -12       , 'SID'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(22)  =  1.
   sn_trcobc(23)  = 'bdy_bio.nc'             , -12       , 'PHO'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(23)  =  1.
   sn_trcobc(24)  = 'bdy_bio.nc'             , -12       , 'GEL'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(24)  =  1.
   sn_trcobc(25)  = 'bdy_bio.nc'             , -12       , 'NOC'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(25)  =  1.
   sn_trcobc(26)  = 'bdy_bio.nc'             , -12       , 'CHA'    , .false.      , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trofac(26)  =  1.

   cn_dir_sbc    =  '../GEO_LR/'             !  root directory for the location of OPEN data files
   sn_trcsbc(14)    = 'deposition.nc'        ,   -12     , 'NOS'    ,  .false.     , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trsfac(14)    = 1.0
   sn_trcsbc(15)    = 'deposition.nc'        ,   -12     , 'NHS'    ,  .false.     , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trsfac(15)    = 1.0
   sn_trcsbc(23)    = 'deposition.nc'        ,   -12     , 'PHO'    ,  .false.     , .true.  , 'yearly'   , ''       , ''       , ''
   rn_trsfac(23)    = 1.0

   cn_dir_cbc    =  '../GEO_LR/rnf_bgc/'     !  root directory for the location of COASTAL data files
   sn_trccbc(1)   = 'rnf'                    , -1        , 'CFL'   ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(1)   = 1.0 
   sn_trccbc(2)   = 'rnf'                    , -1        , 'NFL'   ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(2)   = 1.0 
   sn_trccbc(3)   = 'rnf'                    , -1        , 'CEM'   ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(3)   = 1.0 
   sn_trccbc(4)   = 'rnf'                    , -1        , 'NEM'   ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(4)   = 1.0 
   sn_trccbc(5)   = 'rnf'                    , -1        , 'CDI'   ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(5)   = 1.0 
   sn_trccbc(6)   = 'rnf'                    , -1        , 'NDI'   ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(6)   = 1.0 
   sn_trccbc(7)   = 'rnf'                    , -1        , 'MIC'   ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(7)   = 1.0 
   sn_trccbc(8)   = 'rnf'                    , -1        , 'MES'   ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(8)   = 1.0 
   sn_trccbc(9)   = 'rnf'                    , -1        , 'BAC'   ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(9)   = 1.0 
   sn_trccbc(10)  = 'rnf'                    , -1        , 'DCL'    ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(10)  = 1.0 
   sn_trccbc(11)  = 'rnf'                    , -1        , 'DNL'    ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(11)  = 1.0 
   sn_trccbc(12)  = 'rnf'                    , -1        , 'DCS'    ,  .true.      , .false , 'yearly'   , ''       , ''       , ''
   rn_trcfac(12)  = 1.0 
   sn_trccbc(13)  = 'rnf'                    , -1        , 'DNS'    ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(13)  = 1.0 
   sn_trccbc(14)  = 'rnf'                    , -1        , 'NOS'    ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(14)  =  1.0 
   sn_trccbc(15)  = 'rnf'                    , -1        , 'NHS'    ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(15)  =  1.0 
   sn_trccbc(16)  = 'rnf'                    , -1        , 'SIO'    ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(16)  =  1.0 
   sn_trccbc(17)  = 'rnf'                    , -1        , 'DOX'    ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(17)  =  1.0 
   sn_trccbc(18)  = 'rnf'                    , -1        , 'DIC'    ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(18)  = 1.0 
   sn_trccbc(19)  = 'rnf'                    , -1        , 'unit'   ,  .true.     , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(19)  = 0.0
   sn_trccbc(20)  = 'rnf'                    , -1        , 'POC'    ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(20)  = 1.0 
   sn_trccbc(21)  = 'rnf'                    , -1        , 'PON'    ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(21)  = 1.0  
   sn_trccbc(22)  = 'rnf'                    , -1        , 'SID'    ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(22)  =  1.  
   sn_trccbc(23)  = 'rnf'                    , -1        , 'PHO'    ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(23)  =  1.    
   sn_trccbc(24)  = 'rnf'                    , -1        , 'GEL'   ,  .true.     , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(24)  = 1.0 
   sn_trccbc(25)  = 'rnf'                    , -1        , 'NOC'   ,  .true.     , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(25)  = 1.0 
   sn_trccbc(26)  = 'rnf'                   ,  -1       , 'CHA'    ,  .true.      , .false. , 'yearly'   , ''       , ''       , ''
   rn_trcfac(26)  = 1.0 ! 3000.0
  ln_rnf_ctl    = .true.    !  Remove runoff dilution on tracers with absent river load
   rn_sbc_time   =  1. ! 86400.   !  Time scaling factor for SBC data (seconds in a day)
   rn_cbc_time   =  1. ! 86400.   !  Time scaling factor for CBC data (seconds in a day)
/

!----------------------------------------------------------------------
&namtrc_bdy      !   Setup of tracer boundary conditions
!-----------------------------------------------------------------------
   cn_trc_dflt   = 'neumann'   ! OBC applied by default to all tracers
   cn_trc        = 'specified' ! Boundary conditions used for tracers with data files (selected in namtrc)

   nn_trcdmp_bdy = 0         !  Use damping timescales defined in nambdy of namelist
                             !  = 0 NO damping of tracers at open boudaries
                             !  = 1 Only for tracers forced with external data
                             !  = 2 Damping applied to all tracers
/
!-----------------------------------------------------------------------
&namtrc_ais      !  Representation of Antarctic Ice Sheet tracers supply
!-----------------------------------------------------------------------
   nn_ais_tr     =  1        !  tracer concentration in iceberg and ice shelf
                             !    = 0 (null concentrations)
                             !    = 1 prescribed concentrations
   rn_icbdep     =  120.     ! Mean underwater depth of iceberg (m)
/
!-------------------------------------------------------------------------!
!NAMELIST namben
!-----------------------------------------------------------------------
! General setting of the benthic module
!-----------------------------------------------------------------------
&namben
  jp_ben    = 5       ! total number of benthic tracers
  jp_benlvl = 1       ! number of benthic vertical layers
  ln_bendta = .true.  ! read input data from files . M.C: model bug, this parameter doesn't enforce anything. Ln_rsttr (above) also decides for benthic restart.
  cn_benrst_in  = "GITLAB_CICD_00000216_restart_ben"   !  suffix of pass. sn_tracer restart name (input)
  cn_benrst_out = "restart_ben"   !  suffix of pass. sn_tracer restart name (output)
  !                !    name   !           title of the field              ! initial data ! initial data !
  !                !           !                                           !  units       ! from file    !
  !                !           !                                           !              ! or not       !
  sn_bentrc(1)     = 'fCSED'   ,  'fCSED'                                  ,  'mmol/m3'   ,   .true.     , 
  sn_bentrc(2)     = 'sCSED'   ,  'sCSED'                                  ,  'mmol/m3'   ,   .true.     , 
  sn_bentrc(3)     = 'sSSED'   ,  'sSSED'                                  ,  'mmol/m3'   ,   .true.     , 
  sn_bentrc(4)     = 'fSSED'   ,  'fSSED'                                  ,  'mmol/m3'   ,   .true.     , 
  sn_bentrc(5)     = 'NCrSED'  ,  'NCrSED'                                 ,  'NONE'   ,   .true.     , 
/
!-----------------------------------------------------------------------
&namben_dta      !    Initialisation from data input file
!-----------------------------------------------------------------------
!
!                !  file name               ! frequency (hours) ! variable   ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! Sea over !
!                !                          !  (if <0  months)  !   name     !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  !  Land    !
   sn_bendta(1)  = 'bioinit'      ,        -12        ,  'fCSED'   ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   sn_bendta(2)  = 'bioinit'      ,        -12        ,  'sCSED'   ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   sn_bendta(3)  = 'bioinit'      ,        -12        ,  'sSSED'   ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   sn_bendta(4)  = 'bioinit'      ,        -12        ,  'fSSED'   ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   sn_bendta(5)  = 'bioinit'      ,        -12        ,  'NCrSED'  ,    .false.   , .true. , 'yearly'  , ''       , ''       , ''
   !
   cn_dir        = '../GEO_LR/'
   rn_benfac(1)  =   1.0      !  multiplicative factor
   rn_benfac(2)  =   1.0
   rn_benfac(3)  =   1.0
   rn_benfac(4)  =   1.0
   rn_benfac(5)  =   1.0
! rn_benfac(6)  =   1.0
/
!-----------------------------------------------------------------------
&namben_curstress      !    current bottom stress option
!-----------------------------------------------------------------------
cur_botstr     = 1    ! 0 = compute bottom stress from currents with Bamhbi formulation, logarithmic, grain size dependent
                      ! 1 = compute bottom stress as in nemo

/
!-----------------------------------------------------------------------
&namben_wbstress      !    wave bottom stress files
!-----------------------------------------------------------------------
waves_botstr     = 0  ! 0 = put bottom stress due to waves to zero
                      ! 1 = read in wave bottom stress (fill in sn_wbstress(1))
              ! 2 = compute wave bottom stress (fill in sn_wbstress(2)=period, sn_wbstress(3)=height)
              ! 3 = compute wave bottom stress and consider wave direction (fill in sn_wbstress(2)=period, sn_wbstress(3)=height, sn_wbstress(4)=vmdr)

!
!                  !  file name               ! frequency (hours) ! variable       ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! Sea over !
!                  !                          !  (if <0  months)  !   name         !   (logical)  !  (T/F) ! 'monthly' ! filename ! pairing  !  Land    !
sn_wbstress(1)     =  'WAM'                   ,         6         ,  'w_botstress' ,    .true.    , .false., 'yearly'  ,'',   ''
sn_wbstress(2)     =  'WAM'                   ,         1         ,  'VTPK'        ,    .true.    , .false., 'yearly'  ,'bslr_WAVES_bilin',   ''
sn_wbstress(3)     =  'WAM'                   ,         1         ,  'VHM0'        ,    .true.    , .false., 'yearly'  ,'bslr_WAVES_bilin',   ''
sn_wbstress(4)     =  'WAM'                   ,         1         ,  'VMDR'        ,    .true.    , .false., 'yearly'  ,'bslr_WAVES_bilin',   ''
!
cn_dir='../GEO_LR/waves/'
/
!-------------------------------------------------------------------------!
!NAMELIST namtrc_BAMHBI
!-----------------------------------------------------------------------
!   Set BAMHBI biogeochemical model parameters
!-----------------------------------------------------------------------
&namtrc_bamhbi
!
!==========
! physics
!==========
 daytosecond = 86400.
 pi=3.1415926535897932384626433832795029
 WattToPhotons = 4.56
!
!================!
! Chemical model !
!================!
! Bacterial respiration
 OCr=1.                              ! moldox/molC  Oxygen to carbon ratio of respiration process
 NCr=0.8                             ! molN/molC    Nitrogen to carbon ratio of denitRification
 ksremindox=3.                       ! mmoldox/m3     Half saturation constant for oxygen limitation in oxic mineralization,Soetaert et al., 1996
!ksdeninos=2.
 ksdeninos=0.3                       ! mmolnos/m3     Half saturation constant for nitrate limination in denitrification, Soetaert et al., 1996
 kindenidox=.5       !kindenidox=30. ! mmoldox/m3     Half saturation constant for oxygen inhibition in denitrification, Soetaert et al., 1996
 kinanoxremdox=0.0005                ! mmoldox/m3     Half saturation constant for oxygen inhibition in anoxic remineralization, fitted
 kinanoxremnos=0.0005                ! mmolnos/m3     Half saturation constant for nitrate inhibition in anoxic remineralization, fitted
 ODUCr=1.                            !   molODU/molC    mol ODU formed per mol C in anoxic mineralization
 ODU_solid=0.2                       !   - Pourcentage of formation of solid ODU

! Ammonium oxidation
 Roxnhs=3.472222222222222e-7         ! 0.03 day^-1    maximum rate of ammonium oxidation by oxygen, Soeteart et al,
 ksoxnhsdox=3.00                     ! mmoldox/m3  Half saturation constant for oxygen limitation in ammonium oxidation, Soetaert et al., 1996
!kinoxnhsdox=5.
 kinoxnhsdox=8.                      ! mmoldox/m3  Half saturation constant for oxygen inhibition in ammonium oxidation by nitrate fitted
 kinoxnhsodu=0.5                     ! inhibition of amanox by ODU .. amanox dont need nitrate !!
 ksoxnhsnos=0.3                      ! bug fix 2019-01-07
!
 ONoxnhsr=2.                         ! moldox/molnhs   mol of dox needed to oxidize one mol of NHs
 NOsNHsr=0.6                         !      molnos/molnhs mol of nos needed to oxidize one mol of nhs
!
 Roxnhsnos=5.787037037037037e-7      ! 0.05/daytosecond    !  maximum rate of ammonium oxidation by nitrate
!
! Reduced Substances oxidation
 Roxodu=1.1574074074074074e-6        ! 0.1  d-1 maximum oxidation rate of Oxygen demands units, Oguz et al. 2000 (obtained for H2S and Mn oxidation)
 ksoxodudox=1.                       ! mmoldox/m3  Half saturation constant for oxygen limitation in ODU oxidation Soetaert et al., 1996
 ksoxodunos=2
!ksoxodunos=0.3                      ! mmolNOs/m3  Half saturation constant for nitrate limitation in ODU oxidation fitted
 kinoxodudox=5.
!kinoxodudox=30.                     ! mmoldox/m3 Half saturation constant for oxygen inhibition in ODU oxidation by nitrate fitted
 OODUr=1.                            ! moldox/molODU mol of dox needed to oxidize one mol of ODU
 NODUr=0.8                           ! molnos/molodu mol of nitrate needed to oxidize one mol of ODU
 Roxodunos=5.787037037037037e-7      ! 0.05 d-1 maximum oxidation rate of Oxygen demands units by nitrate Oguz et al. 2000 (obtained for Mn oxidation)
!
!Definition of the suboxic layer
 Q10chem=2.                          ! - Q10 factor for chemical processes
!
! Phosphate precipitation by Iron3 and  Manganese3
  Fe3Mn3_density=16.1                ! density at which Fe3 and Mn3 are maximum, so phosphate precipitation occurs
  Pho_precip_fac=1.1574074074074e-8  ! 0.001 / 86400
!
!===============!
!   BACTERIA    !
!===============!
 NCrBac=0.19607843137254904          ! 1./5.1  mmolN/mmolC   Goldman ratio for Bacteria Anderson and Pondhaven, 2003
 maxhydrDOCSL=4.6296296296296294e-5  ! 4  /day  max.semi-labile DOC hydrolysis Anderson and Pondhaven, 2003
 csatdocsl=417.                      ! mmolC/m3  half-saturation for semi-labile DOC uptake Anderson and Pondhaven, 2003 instaed of 417
 maxgrowthbac=0.0001539351851851852  ! 13.3 /day  max.labile DOC or ammonium uptake Anderson and Pondhaven, 2003
 csatdocl=25.                        ! mmolC/m3  half saturation for labile DOC uptake Anderson and Pondhaven, 2003
 csatamm=0.5                         ! mmolN/m3 half saturation for ammonium uptake Anderson and Pondhaven, 2003
 mortbac=5.787037037037037e-7            ! 0.05  /day bacteria natural mortality Anderson and Pondhaven, 2003
 bactgrowtheff=0.17                  ! - bacterial gross growth efficiency on carbon  Anderson and Pondhaven, 2003
 Q10bac=2.                           ! - Q10 factor for bacteria
 Halfsaturation_Iron=100.            !  nmole/l half sauration constant for iron limitation function
 IronCsurf=10.                       ! nmole/l  Concentration of iron in surface waters
 Param1IronCurve=25000.              !  - Parameter of the curve simulating the iron concentration
 Param2IronCurve=50.                 !  -  Parameter of the curve simulating the iron c
!
!===============!
! POM
!===============!
 hydPOCmax=4.6296296296296297e-7     ! 0.04   /day detrital breakdown rate, C  Anderson and Pondhaven, 2003 0.04
 hydPONmax=6.36574074074074e-7       ! 0.055  /day detrital breakdown rate, N  Anderson and Pondhaven, 2003 0.055
 SinkingRatePOM_Constant=2.3148148148148147e-5 ! 2.0  m/day sinking rate of POM in case of not using the aggregation model
 ksatOxygenHydrolysis=2.70           !  mmolO2/m3 HAlf saturation constant for the oxygen limitation function of te hydrolysis rate instead of 0.02,0.1
!
!===============!
! DON
!===============!
 labilefraction=0.7                  !   - labile fraction phyto and non-phyto produced DOM  Anderson and Pondhaven, 2003
 labileextradocphyexcr=0.65          !   - labile fraction phyto-excreted DOC,  Anderson and Pondhaven, 2003
!
!===============!
! PHYTOPLANKTON !
!===============!
 QuantumYieldDiatoms=0.8              ! mmolC(mgChldW/m2)-1 Maximum quantun yield of diatoms Soetaert et al., 2001
 QuantumYieldFlagellates=0.6          ! mmolC(mgChldW/m2)-1 Maximum quantun yield of Flagellates Soetaert et al., 2001
 QuantumYieldEmiliana=0.6             ! mmolC(mgChldW/m2)-1 Maximum quantun yield of Emiliana Soetaert et al., 2001
!
 alphaPIFlagellates=2.491898148148148e-6  !  0.2153  molC/gchla (W/m2)-1/day   instead of 0.045
 alphaPIEmiliana=3.472222222222222e-6     !  0.3 molC/gchla  (W/m2)-1/day    instead of 0.05
 alphaPIDiatoms= 3.833333333333334e-6     ! 0.3312/daytosecond  molC/gchla (W/m2)-1/day instead of 0.09, it's 0.06 in Gregoire et al PIO2008
!
 MinNCrFlagellates=0.05               ! molN/molC  Soetaert et al., 2001
 MinNCrEmiliana=0.05                  ! molN/molC   Soetaert et al., 2001
 MinNCrDiatoms=0.05                   ! molN/molC   Soetaert et al., 2001
!
 MaxNCrFlagellates=0.2                ! molN/molC   Soetaert et al., 2001
 MaxNCrEmiliana=0.2                   ! molN/molC   Soetaert et al., 2001
 MaxNCrDiatoms=0.2                    ! molN/molC   Soetaert et al., 2001
!
 MuMaxFlagellates= 1.1574074074074073e-5    ! µmax0         ! taux de croissance maximal à 0°C du phyto (identique pour les 2 classes)
 MuMaxEmiliana=2.8935185185185186e-5    !  insted of 3 at 20 Egge and Heimdal, 1994,Sarsia, against 3 at 20 for Asknes et al 1994, 0.8 -1.9 /day Tyreell and Taylor, 1996, =mudia/1.5, Brand and Guillard, 1981
 MuMaxDiatoms= 4.050925925925926e-5    ! 3.50/daytosecond ;   4.5 instead of 3.5 mean values computed from Furnas, 1990 JPR, against 4.6 at 20 for Asknes et al 1994
!
 RespirationFlagellates=1.0416666666666665e-7         ! 0.009  /day  Lancelot et al 2002, coefficient of basal respiration
 RespirationEmiliana=1.0416666666666665e-7            ! 0.009  /day  Lancelot et al 2002, coefficient of basal respiration
 RespirationDiatoms=1.0416666666666665e-7             ! 0.009  /day  Lancelot et al 2002, coefficient of basal respiration
!
 MortalityFlagellates=3.472222222222222e-7        ! 0.03/day   0.0115 for Asknes et al. 1994
 MortalityEmiliana=3.472222222222222e-7           ! 0.03/day   0.0115 for Asknes et al. 1994
 MortalityDiatoms= 3.472222222222222e-7    ! 0.03 [ day^-1 ] , 0.0115 for Asknes et al. 1994
!
 MinGrowthRespFlagellates=0.1         ! -  respirated fraction (linked to activity) Soetaert et al., 2001
 MaxGrowthRespFlagellates=0.1         ! - max repisrated fraction
 GrowthRespEmiliana=0.1               ! -             respirated fraction (linked to activity)  Soetaert et al., 2001
 GrowthRespDiatoms=0.1                ! -             respirated fraction (linked to activity)  Soetaert et al., 2001
 GrowthRespFlagellates=0.1
 Q10Phy=2.0                           ! - temperature factor  Soetaert et al., 2001
 Q10PhyDiatoms=1.8                    ! - temperature function for diatoms instead of 1.5
!Q10PhyFlagellates=1.8                ! - temperature function for diatoms instead of 1.5
!Q10PhyEmilianna=1.8                  ! - temperature function for diatoms instead of 1.5
 MinChlNrFlagellates=1.               ! mgChla/mmolN   Soetaert et al., 2001
 MinChlNrEmiliana=1.                  ! mgChla/mmolN   Soetaert et al., 2001
 MinChlNrDiatoms=1.                   ! mgChla/mmolN   Soetaert et al., 2001
 MaxChlNrFlagellates=2.               ! mgChla/mmolN   Soetaert et al., 2001
 MaxChlNrEmiliana=2.                  ! mgChla/mmolN   Soetaert et al., 2001
 MaxChlNrDiatoms=2.                   ! mgChla/mmolN   Soetaert et al., 2001
!
 NosMaxUptakeFlagellates=5.787037037037037e-6         ! 0.50  molN/molC/d   1  maximal nitrate uptake rate
 NosMaxUptakeEmiliana=1.736111111111111e-5            ! 1.50  molN/molC/d     maximal nitrate uptake rate,
 NosMaxUptakeDiatoms=1.1574074074074073e-5            ! 1.00  molN/molC/d     maximal nitrate uptake rate
 NHsMaxUptakeFlagellates=5.787037037037037e-6         ! 0.50  molN/molC/d   1  maximal ammonium uptake rate,
 NHsMaxUptakeEmiliana=1.736111111111111e-5            ! 1.50  molN/molC/d    maximal ammonium uptake rate,
 NHsMaxUptakeDiatoms=1.1574074074074073e-5            ! 1.00  molN/molC/d     maximal ammonium uptake rate
!
 ksNOsFlagellates=3.00                ! mmolN/m3     3    half saturation constant for nitrate uptake data base
 ksNOsEmiliana=0.05                   ! mmolN/m3        half saturation constant for nitrate uptake data base
 ksNOsDiatoms=1.0                     ! mmolN/m3        half saturation constant for nitrate uptake data base 1.4 (instead of 1)
 ksNHs=1.0                            ! multiplicative factor for ksNHsFlagellates=3.0, ksNHsEmiliana=0.05, ksNHsDiatoms=1.0 ; K.Fennel uses 0.5 for all, range=0.007-1.5
 ksNHsFlagellates=3.0
 ksNHsEmiliana=0.05
 ksNHsDiatoms=1.0
 kinNHsPhy=0.5                        ! mmolNH3-N/m3   Half-satur.ct of NH3 inhibition NO3uptake phytop Soetaert et al., 2001
 leakagephy=0.02                      ! -    Phytoplankton leakage parameter  VDM et al, 2004 L&O instaed of 0.05 as in Anderson and Pondhaven, 2003
 extradocphyexcr=0.05                 ! - Extra photosynthetic carbon excretion VDM et al, 2004 L&O instaed of 0.2v6 as in Anderson and Pondhaven, 2003 instead of 0.26
 mortphydom=0.34                      ! - DOM fraction of phytoplankton mortality  instead of 0.34 Anderson and Pondhaven, 2003
!
 ksPO4Flagellates =  0.2                           ! mmolP/m3   half saturation constant for phosphate uptake by phytoplankton
 ksPO4Diatoms     =  0.1                           ! mmolP/m3   half saturation constant for phosphate uptake by phytoplankton 0.4 for green ocean code
 ksPO4Emiliana    =  0.02                          ! mmolP/m3   half saturation constant for phosphate uptake by phytoplankton 0.004 for green ocean model 0.005 for Soetaert code
 PNRedfield       =  0.0625                        ! 1.0/16.0 molP/molN   Redfield ratio
 csatpo4          =  0.03125                       ! 0.5/16.0 mmolP/m3   half saturation conatsant for ammonium uptake by bacteria computed from AW values*PNratio
 PO4MaxuptakeEmiliana = 1.0850694444444444e-6      ! 1.5/16.0 molP/molC/day    maximal PO4 uptake rate
 PO4MaxuptakeDiatoms = 7.233796296296296e-7        ! 1.0/16.0
 PO4MaxuptakeFlagellates = 3.616898148148148e-7    ! 0.5/16.0
!
 SinkingRateDiatomsmin=0.000001157    ! perturb this variable and the next ; Gregoire et al put 0.1 and 1.0 m/day ; Arthur's version has 0.5 - 2.0 m/day, CMEMS V2021 has 0.1-0.5
 SinkingRateDiatomsmax=0.000011574    !          Wang et al perturb w_phy=0.1m/day in the range 0.009-25
!                                      !   Maximum Sedimentation rate of diatoms  Tyrell and Taylor, 1996, JMS against 0.3 m/day for Asknes et al 1994
!===============!
!   SILICATE    !
!===============!
 SiMaxUptakeDiatoms=5.787037037037037e-6                  ! 0.50  molSi/molC/d
 ksSiDiatoms=3.5                                          ! mmolN/m3       0.5 half saturation constant for silicate uptake, database, between 0.1 -1 Paasche, 1980, 0.73 -1.15 data base
 SiNrDiatoms=0.833333333333333                            ! 5./6. [ molN/molSi ]  Silicate to Nitrogen ratio, database, 0.875 for Aksnes et al. 1994, data base :0.94 -1.48
 kdis_Silicious_Detritus=9.259259259259259e-7             ! 0.08  /day rate of dissolution of silicious  detritus Tusseau, 1996
 vsinkingrate_Silicious_Detritus=2.3148148148148147e-5    ! 2.0   sinking rate of Silicious_Detritus
 Q10SilicateDiss=3.3                                      ! - Q10 factor for chemical processes
!
!===============!
! ZOOPLANKTON   !
!===============!
 MaxgrazingrateMicroZoo=4.1666666666666665e-5     ! 3.6  day-1 maximum grazing rate, Strom and Morello, 1998, JPR
 MaxgrazingrateMesoZoo = 1.388888888889e-5        ! max grazing rate = 1.2/day  [ second^-1 ], it's 1.0 in Gregoire et al PIO2008
                                  ! max grazing rate of microzoo about 3 times higher than mesozoo especially, mesozoo n the BS if Calanus which is a large group
!
 Half_Saturation_MicroZoo=5.      ! mmolC/m3   half saturation constant for the grazing  Soetart et al., 2001
 Half_Saturation_MesoZoo=5.       !mmolC/m3   half saturation constant for the grazing    Soetart et al., 2001
!
 Ass_Eff_OnNitrogen=0.77          !  -   assimilation efficiency of nitrogen Anderson and Pondhaven, 2003
 Ass_Eff_OnCarbon=0.64            !   - assimilation efficiency of carbon  Anderson and Pondhaven, 2003
!
 efficiency_growth_MicroZoo=0.80  ! - Zooplankton net growth efficiency on carbon Anderson and Pondhaven, 2003
 efficiency_growth_MesoZoo=0.8    ! - MesoZooplankton net growth efficiency on carbon
!
 Messy_feeding_MicroZoo=0.23      !   - messy feeding fraction to DOM Anderson and Pondhaven, 2003
 Messy_feeding_MesoZoo=0.23       !   - messy feeding fraction to DOM Ander
!
 Capt_eff_MicroZoo_Flagellates=0. !   -    capture efficiency
 Capt_eff_MicroZoo_Emiliana=1.0   !  -    capture efficiency
 Capt_eff_MicroZoo_Diatoms=0.0    ! -    capture efficiency
 Capt_eff_MicroZoo_MicroZoo=0.    !   -    capture efficiency
 Capt_eff_MicroZoo_MesoZoo=0.     !   -    capture efficiency
 Capt_eff_MicroZoo_pom=0.         !  -    capture efficiency
 Capt_eff_MicroZoo_bac=0.7        !  -  capture efficiency
!
 Capt_eff_MesoZoo_Flagellates=0.4 !  -    capture efficiency
 Capt_eff_MesoZoo_Emiliana=0.4    !  -    capture efficiency
 Capt_eff_MesoZoo_Diatoms=1.0     !  -    capture efficiency
 Capt_eff_MesoZoo_MicroZoo=1.     !  -    capture efficiency
 Capt_eff_MesoZoo_MesoZoo=0.      !  -    capture efficiency
 Capt_eff_MesoZoo_pom=0.8         !  -    capture efficiency
 Capt_eff_MesoZoo_bac=0.          !  -  capture efficiency
!
 HalfSatMort_MicroZoo=1.             !  mmolC/m3   mortality half stauration rate
 expmortMicroZoo=2                   !  -       higher exponent of the mortality term
 NLin_Mort_MicroZoo=3.4722222222e-6  ! 0.3  day-1/mmolC/m3   mortality term in case of quadratic term (betazoo1=2) Anderson and William
 HalfSatMort_MesoZoo=1.              ! mmolC/m3   mortality half stauration rate   ! it's 5.0 in Gregoire et al PIO2008 !
 expmortMesoZoo=2                    !   -       higher exponent of the mortality term
 nlin_mort_MESozoo= 3.4722222222e-6  ! 0.3  day-1/mmolC/m3   mortality term in case of quadratic term (betazoo1=2)  Anderson and William instaed of 0.3
 Mortanoxic=2.8935185185185184e-6    ! 0.25 day-1    term of mortality due to the lack of DOX
 DOXsatmort=7.8125                   ! mmolO/m3    percentage of saturation where metabolic respiration is half the one under oxygen satyrated conditions
!
 Q10Zoo=2.                           !   -       temperature factor  Soetart et al., 2001
 NCrMicroZoo=0.18181818181818182     !  1./5.5   mmolN/mmolC   0.18 N:C molar ratio  Anderson and Pondhaven, 2003
 NCrMesoZoo=0.21413276231263384      !  1./4.67  mmolN/mmolC   0.21 N:C molar ratio
!
!
!========= PARAMETERS MISSING FROM THE MODEL, BUT THAT SHOULD BE PERTURBED SOMEDAY ========!
!r_ld=0.1           ! large detritus remineralisation rate [day^-1], range: 0.005-0.25
!tau=0.1            ! aggregation parameter [day^-1], range: 0.01-25
!theta_max=0.0535   ! max chloro-to carbon ratio [mg Chl / mg C], range= 0.005-0.15                    See Geider et al, 1997, perturb later (Catherine)
!w_LDet=1.0         ! large detritus sinking velocity [m day^-1], range=0.009-25
!
/
!====================
&namtrc_bamhbi_geider
!====================
 MinChlCrPHY=0.03                     ! MinChlNrPHY * MinNCrPHY
 MaxChlCrPHY=0.15                     ! MaxChlNrPHY * MaxNCrPHY
 MinChlNrPHY=1.                       ! mgChla/mmolN   Soetaert et al., 2001
 MaxChlNrPHY=2.                       ! mgChla/mmolN   Soetaert et al., 2001
!
 NHSNoUptakeFlagellates = 0.0
 NHSNoUptakeEmiliana = 0.0
 NHSNoUptakeDiatoms = 0.0
 PO4NoUptakeFlagellates = 0.0001
 PO4NoUptakeEmiliana = 0.0001
 PO4NoUptakeDiatoms = 0.0001
/
!
!=========================
&namtrc_bamhbi_aggregation           ! Aggregation parameters taken from Kriest, 2002, DSR
!=========================
 shear = 0.84                        ! 0.84/second        shear rate Kriest 2002, DSR
 etabio=0.62                         !       -           sinking exponent Kriest 2002, DSR
 B=0.026550925925925926              ! 2294 /m^(etabio-1)/day  sinking factor  Kriest 2002, case 9 table 2, agaisnt 206087 for single cell,DSR
 stick=0.55                          !       -           stickness  Kriest 2002, DSR
 m=0.00002                           !     m           cell size  Kriest 2002, DSR
 dzetabio=1.62                       !     -              N content exponent  Kriest 2002, table 1, ref 2a, aigains 2.28 for case 5 DSR
 C=0.4744                            !       mmolN/m^dzetabio   N content coefficient Kriest 2002 ,table 1, ref 2a, aigains 206.95 for case 5, DSR
 Maxsize=0.015                       !   m    Maximum cell size. Fitted Kriest 2002, DSR proposed a value of 1.5/100 m
 epsini=3.5                          !  -epsilon initial value Kriest 2002, DSR
 dAGG_limit=11574.074                ! maximum change per day (for numerical stability reason) -- 0.0 means no limit , default 1.e9/86400
/
!
!========================
&namtrc_bamhbi_gelatinous
!========================
 MaxgrazingrateNoctiluca  = 1.7361e-7     ! 0.020/daytosecond ! (day mmolC/m3)-1     ingestion coefficient Lancelot et al., 2002 ESCS
 MaxgrazingrateGelatinous        = 3.472E-6                ! (0.3 day mmolC/m3)-1     ingestion coefficient Lancelot et al., 2002 ESCS
!
 Ass_Eff_Noctiluca      = 0.75                     ! -   assimilation efficiency  Lancelot et al., 2002 ESCS
 Ass_Eff_Gelatinous     = 0.75                     ! -   assimilation efficiency  Lancelot et al., 2002 ESCS
!
 Capt_eff_Noctiluca_Flagellates  = 0.5             ! -    capture efficiency
 Capt_eff_Noctiluca_Emiliana  = 1.0                ! -    capture efficiency
 Capt_eff_Noctiluca_Diatoms  = 1.0                 ! -    capture efficiency
 Capt_eff_Noctiluca_microzoo = 1.0                 ! -    capture efficiency
 Capt_eff_Noctiluca_mesozoo = 0.                   ! -    capture efficiency
 Capt_eff_Noctiluca_pom = 1.                       ! -    capture efficiency
 Capt_eff_Gelatinous_Flagellates  = 0.             ! -    capture efficiency
 Capt_eff_Gelatinous_Emiliana  = 0.                ! -    capture efficiency
 Capt_eff_Gelatinous_Diatoms  = 0.                 ! -    capture efficiency
 Capt_eff_Gelatinous_microzoo = 0.                 ! -    capture efficiency
 Capt_eff_Gelatinous_mesozoo = 1.                  ! -    capture efficiency, it's 0.4 in Gregoire et al PIO2008 !
 Capt_eff_Gelatinous_pom = 0.                      ! -    capture efficiency
!
 threshold_feeding_Noctiluca = 0.833               ! mmolC/m3  feeding threshold Nakamura, 1998, JPR
 threshold_feeding_Gelatinous = 0.0                ! 0./12.  mmolC/m3  feeding threshold Lancelot et al., 2002 ESCS
!
 HalfSatMort_Noctiluca  = 0.0                      !  mmolC/m3   mortality half stauration rate
 expmortNoctiluca = 2.0                            !  -       higher exponent of the mortality term
 NLin_Mort_Noctiluca = 6.944444444444444e-7        !  0.06  day-1/mmolC/m3   mortality term in case of quadratic term (betazoo1=2)  0.04 instead of 0.17 (Lancelot et al., 2002)
 HalfSatMort_Gelatinous =0.0                       !  mmolC/m3   mortality half stauration rate
 expmortGelatinous = 2.0                           !     -       higher exponent of the mortality term
 NLin_Mort_Gelatinous = 1.0416666666666665e-7      ! 0.009  day-1/mmolC/m3   mortality term in case of quadratic term (betazoo1=2)  0.003 instead of 0.012 (Lancelot et al., 2002)
!
 efficiency_growth_Noctiluca    =  0.15            !   -    Part of the Assimilated food of zooplankton used for growth  Lancelot et al., 2002 ESCS
 efficiency_growth_Gelatinous   =  0.2             !   -    Part of the Assimilated food of zooplankton used for growth  Lancelot et al., 2002 ESCS
!
 basal_Resp_Noctiluca    = 1.1574074074074074e-9   !  0.0001  day-1  basal respiration rate
 basal_Resp_Gelatinous   =    0.00                 !   0.00   day-1  basal respiration rate
!
 NCrNoctiluca        =0.21                         !  mmolN/mmolC   N:C molar ratio Nakamura, 1998, JPR
 NCrGelatinous       =0.25                         !  mmolN/mmolC  N:C molar ratio
!
 SinkingRateNoctiluca   = 0.                       ! m/sec Sedimentation rate of zooplankton
 SinkingRateGelatinous   = 0.                      ! m/sec Sedimentation rate of zooplankton
!
 Q10Gelatinous   =3.5                              !  - Q10 fqctor for temperature dependency between 3.4-3.7 Kremer 1977
/
!
!====================
&namtrc_bamhbi_optics   ! in case SIMPLEOPTICS or FULLOPTICS are used
!====================
 kBIOABSCHL=0.03
 kBIOABSPOC=0.03
!
 a_w_short=0.0196                    ! m-1, absorption by pure sea water Smith and baker 1981
 bb_w_short=0.0015                   ! m-1, backscatering by pure sea water Smith and baker 1981
 a_chl_A_short=0.029                 ! A*CHL^B [m-1], Bricaud 95 DMITRIEV 2007
 a_chl_B_short=0.6                  ! A*CHL^B [m-1], Bricaud 95 DMITRIEV 2007
 bb_dia_short=0.0000058                 ! m2/mgPOC_DIA, backscatering by Diatoms, Vaillancourt 2004
 bb_fla_short=0.00000817                ! m2/mgPOC_FLA, backscatering by Diatoms, Vaillancourt 2004
 bb_emi_short=0.00001008               ! m2/mgPOC_EMI, backscatering by Diatoms, Vaillancourt 2004
 a_poc_short=0.5                    ! m2/g , absorption by organic matter NEUCKERMANS 2012
 bb_poc_short=0.0055                 ! m-1, backscatering by organic matter NEUCKERMANS 2012
 a_cdom_intercept_short = 0.2522 ! 0.13@10 0.02@19
 a_cdom_slope_short = -0.0122
 !
 a_w_long=0.24                    ! m-1, absorption by pure sea water Smith and baker 1981
 bb_w_long=0.0007                   ! m-1, backscatering by pure sea water Smith and baker 1981
 a_chl_A_long=0.0066                 ! A*CHL^B [m-1], DMITRIEV 2007
 a_chl_B_long=0.8                 ! A*CHL^B [m-1], Bricaud 95, DMITRIEV 2007
 bb_dia_long=0.00000441                 ! m2/mgPOC_DIA, backscatering by Diatoms, Vaillancourt 2004
 bb_fla_long=0.00000612                ! m2/mgPOC_FLA, backscatering by Diatoms, Vaillancourt 2004
 bb_emi_long=0.00000756               ! m2/mgPOC_EMI, backscatering by Diatoms, Vaillancourt 2004
 a_poc_long=0.5                    ! m2/g , absorption by organic matter NEUCKERMANS 2012
 bb_poc_long=0.005                ! m-1, backscatering by organic matter NEUCKERMANS 2012
 a_cdom_intercept_long = 0.2155
 a_cdom_slope_long = -0.0113
 ! from a, bb, to Kd !! see LEE 2005b
 mu_d = 0.8                      ! average downward cosine .. should be computed considering senith angle, a/b ratio, depth
 !
 !  Parameters for water classification - default Jerlov type I
 !  ARTCHAnge : distinction supplémentaire ente
 !  ondes visible longue (LightAbsB), LightAbsVisL,
 !  ondes visible court  (1-LightAbsB),  LightAbsVisS
 !  0.003*0.05 (abs par chl,concentration moyenne de 0.05)
 !  +0.003*0.01 (abs par det,concentration moyenne de 0.01)
 !  unit is m**-1
 LightAbsA=0.54
 LightAbsB=0.63
 LightAbsIR=4.0
 LightAbsVisL=0.23
 LightAbsVisS=0.02
 !ART les paramètres ci-dessus devraient être différents
 !dans le cas de couplage bio puisque l'absorption par la bio vient s'ajouter
 !à ces paramètres-ci qui sont l'absorption intrinsèque de l'eau
 !et qui dans le cas de la physique seule devraient sans doute représenter
 !une absorption moyenne par la bio ...
/
!
!=====================
&namtrc_bamhbi_benthic
!=====================
 pfCSED     = 0.29                  !- Fraction of carbon flux that is fast decaying
 pfSSed     = 0.5                   !- Fraction of Silice flux that is fast decaying
 fCdegrate  = 8.715277777777778e-7  !- /s Degradation rate of the fast degrated carbon  (0.0753/day)
 sCdegrate  = 3.472222222222222e-8  !- /s Degradation rate of the fast degrated carbon  (0.003/day)
 fSdisrate  = 1.6203703703703703e-7 !- /s Dissolution rate of the fast degrated silice  (0.014/day)
 sSdisrate  = 1.6203703703703705e-8 !- /s Dissolution rate of the fast degrated silice  (0.0014/day)
 psoliddepo = 0.11                  !- Fraction of precipiting ODU
 OCrdegrad  = 1                     !- molO/molC Oxygen consumption for POM degradation
 ONrnitrif  = 2                     !- molO/molN Oxygen consumption for nitrification
 NCrsedinit = 0.15                  !- molO/molN Oxygen consumption for nitrification
 Q10CDEG    = 2                     !- Q10 factor for temperature influence on C degradation
 sCburialrate = 2.5e-9  ! put zero for conservation
 sSburialrate = 2.5e-9  ! put zero for conservation
!
!
!=========RESUPENSION==========!
 criticalstress_DEP   = 0.02           ! N/m2 thresold on bottom stress for deposition
 criticalstress_ERO_S = 0.04           ! N/m2 thresold on bottom stress for erosion of slow stock, could be 0.04
 criticalstress_ERO_F = 0.02           ! N/m2 thresold on bottom stress for erosion of fast stock
 Me_sS = 1.0e-5                        ! mmol/m2s Erosion constant for slow silicious sediment
 Me_fS = 1.0e-5                        ! mmol/m2s Erosion constant for fast silicious sediment
 Me_sC = 1.0e-5                        ! mmol/m2s Erosion constant for slow carbonate sediment
 Me_fC = 1.0e-5                        ! mmol/m2s Erosion constant for fast carbonate sediment
 alphaRESUSP = 1.0                     !
/
!========== CARBONATE SYSTEM ============!
&namtrc_bamhbi_alka
 ln_reinitialize_pH=.false.
 kbulk = 1.0
 MeanpC02air = 380.0
/
!--------------------------------------------------------------------------------
&nam_pCO2          ! atmospheric pCO2 forcing files -- #define airpco2 in bamhbi.h90
!-------------------------------------------------------------------------------
!                  !  file name               ! frequency (hours) ! variable       ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! Sea over !
!                  !                          !  (if <0  months)  !   name         !   (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  !  Land    !
sn_pCO2(1)         =  'pCO2'                  ,         24        ,  'pCO2'        ,    .false.   , .false. , 'yearly'  ,    ''    ,   ''     ,   ''     !
!
cn_dir='../GEO_LR/pco2/'
/
!-----------------------------------------------------------------------
&namsto_bio             !    stochastic BGC parameters and forcings
!-----------------------------------------------------------------------
!
   ln_sto_rnf_bio  = .false.     ! multiplicative stochastic river nutriment runoff          ! L.V.: implemented in trcbc
   rn_rnfbio_std   = 0.2         ! bio runoff standard deviation
   rn_rnfbio_tcor  = 52560.0     ! bio runoff time correlation (in timesteps)
   rn_rnfbio_xycor = 72          ! bio runoff spatial correlation (radius in grid points activated > 1. ; default = 1.)
   nn_rnfbio_ord  = 1            ! order of autoregressive processes
   nn_rnfbio_flt  = 0            ! passes of Laplacian filter
   rn_rnfbio_lim  = 2.0          ! limitation factor        ! BGC runoff is provided in mmol/s --> multiplicative perturbation by both sto_rnf and sto_rnf_bio !
   !
   ln_sto_depo   = .false.       ! multiplicative stochastic atmospheric deposition runoff                     ! L.V.: implemented in trcbc
   rn_depo_std   = 0.8           ! deposition standard deviation, first try 0.2, next try 0.8
   rn_depo_tcor  = 4032.0        ! deposition time correlation (in timesteps)
   rn_depo_xycor = 36            ! deposition spatial correlation (radius in grid points activated > 1. ; default = 1.), changed from 36 to 20 (LR model)
   nn_depo_ord   = 1             ! order of autoregressive processes
   nn_depo_flt   = 0             ! passes of Laplacian filter
   rn_depo_lim   = 3.0           ! limitation factor (default = 3.0)
/
&namsto_bamhbi                         ! positive number: std/µ
   ! sto_XXXXX is a positive number: std/µ  ;  the distribution (normal, gamma, beta, loggamma) is hard-coded in stopar.F90
   ! if sto_XXXXX is zero, the parameter is not perturbed
   ! for some parameters, the special value 9999 indicates to use the same random perturbation as for another parameter
   !
   ! draw a random number in a given distribution, following Prieur et al, OM 2019
   ! - positive or negative parameters are following a gamma distribution
   ! - parameters between 0 and 1 follow a beta distribution
   ! - parameters constrained to be larger than 1: their logarithm follows a gamma distribution
   ! - non-constrained parameters follow a gaussian distribution
   ! - std dev is expected to be 5% 20% or 50% of the expected value, value is specified in the namelist
   !
   sto_MaxgrazingrateMesoZoo           = 0.0 ! 0.2
   sto_alphaPIDiatoms                  = 0.0 ! 0.3
   sto_Ass_Eff_OnCarbon                = 0.0 ! 0.2
   sto_efficiency_growth_MesoZoo       = 0.0 ! 0.1
   sto_nlin_mort_MESozoo               = 0.0 ! 0.5
   sto_HalfSatMort_MesoZoo             = 0.0 ! 0.666666
   sto_Capt_eff_Gelatinous_mesozoo     = 0.0 ! 0.2
   sto_MaxgrazingrateGelatinous        = 0.0 ! 0.2
   sto_Capt_eff_MesoZoo_Diatoms        = 0.0 ! 0.1
   sto_Ass_Eff_OnNitrogen              = 0.0 ! 0.1
   sto_MuMaxDiatoms                    = 0.0 ! 0.3
   sto_SiNrDiatoms                     = 0.0 ! 0.05
   sto_MinNCrDiatoms                   = 0.0 ! 0.1
   sto_MortalityDiatoms                = 0.0 ! 0.5
   sto_Messy_feeding_MicroZoo          = 0.0 ! 0.2
   sto_Messy_feeding_MesoZoo           = 0.0 ! 0.2
   !
   ! non-duplicate parameters from Wang, Fennel et al 2020 BGC ; see also Fennel et al 2006 doi:10.1029/2005GB002456
   sto_SinkingRateDiatomsmin           = 0.0 ! 0.5-0.75
   sto_SinkingRateDiatomsmax           = 0.0 ! 0.5-0.75
   sto_MortalityFlagellates            = 0.0 ! 9999. or 0.5   ! if 9999, use same mortality as Diatoms
   sto_MortalityEmiliana               = 0.0 ! 9999. or 0.5   ! idem  (this was the case originally in NEMO-BAMHBI)
   sto_ksNHs                           = 0.0 ! 0.5-0.75
   !
   ! 3 non-duplicate parameters from Fl. Garnier 2006
   sto_MuMaxFlagellates                = 0.0 ! 0.5
   sto_MuMaxEmiliana                   = 0.0 ! 0.3
   sto_Q10Phy                          = 0.0 ! 0.2     ! LogGamma distribution (parameter is exponential of gamma draw, so that parameter > 1) -- here sto_Q10 is std, not std/mean
   sto_Q10PhyDiatoms                   = 0.0 ! 9999.   ! if 9999, then multiply by same random coeff as Q10Phy
   sto_Q10Zoo                          = 0.0 ! 9999.   ! idem (originally, Q10Phy=2.0, Q10Diatoms=1.8, Q10Zoo=2.0)
   !
   ! non-duplicate parameters from A. Capet (Albatros project)
   sto_dzetabio                        = 0.0 ! 0.05    ! try small values such as 0.05 because dzetabio appears as exponent ; std(dzetabio)=0.081 when E(dzetabio)=1.62
   sto_etabio                          = 0.0 ! 0.05
   sto_lightAbsB                       = 0.0 ! 0.05
   sto_Q10Chem                         = 0.0 ! 9999. | 0.2
   !
   ! still missing
   ! sto_* ::  r_ld (detritus remineralisation rate), aggregation, max CHL-to-C ratio, detritus sinking rate (because it's computed, not a parameter)
   ! sto_* :: radtrans parameters
/
