!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! NEMO/OCE :   Reference namelist_ref                                !!
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! NEMO/OCE  :  1 - Domain & run manager (namrun, namcfg, namdom, namtsd, namcrs, namc1d, namc1d_uvd)
!! namelists    2 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_cpl,
!!                                    namsbc_sas, namtra_qsr, namsbc_rnf,
!!                                    namisf, namsbc_apr, 
!!                                    namsbc_ssr, namsbc_wave, namberg)
!!              3 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide)
!!              4 - top/bot boundary (namdrg, namdrg_top, namdrg_bot, nambbc, nambbl)
!!              5 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_eiv, namtra_dmp)
!!              6 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf)
!!              7 - Vertical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_gls, namzdf_iwm)
!!              8 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb)
!!              9 - Obs & Assim      (namobs, nam_asminc)
!!             10 - miscellaneous    (nammpp, namctl, namsto)
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>

!!======================================================================
!!              ***  Domain & Run management namelists  ***           !!
!!                                                                    !!
!!   namrun       parameters of the run
!!   namdom       space and time domain
!!   namcfg       parameters of the configuration                       (default: user defined GYRE)
!!   namwad       Wetting and drying                                    (default: OFF)
!!   namtsd       data: temperature & salinity                          (default: OFF)
!!   namcrs       coarsened grid (for outputs and/or TOP)               (ln_crs =T)
!!   namc1d       1D configuration options                              (ln_c1d =T)
!!   namc1d_dyndmp 1D newtonian damping applied on currents             (ln_c1d =T)
!!   namc1d_uvd   1D data (currents)                                    (ln_c1d =T)
!!======================================================================
!
!-----------------------------------------------------------------------
&namrun        !   parameters of the run
!-----------------------------------------------------------------------
   nn_no       =       0   !  Assimilation cycle index
   cn_exp      =  "GITLAB_CICD"     !  experience name
   nn_it000    =       1  !  first time step
   nn_itend    = 216   ! simulate 2014-04-01 (full day)
   nn_date0    = 20140401  !  date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
   nn_stock    =      0    !  frequency of creation of a restart file (modulo referenced to 1) 
   nn_leapy    =       1   !  Leap year calendar (1) or not (0) ! Read date in namelist
   nn_istate   =       0   !  output the initial state (1) or not (0)
   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T)
      nn_rstctl   =    0   !  time info read in restart file
      cn_ocerst_in    = "GITLAB_CICD_00000216_restart"   !  suffix of ocean restart name (input)
      cn_ocerst_out   = "restart"   !  suffix of ocean restart name (output)
      nn_wxios        = 0           !  use XIOS to write restart file 0 - no, 1 - single file output, 2 - multiple file output
      ln_xios_read = .false.  !  use XIOS to read restart file (only for a single file restart)
 ln_mskland = .true.
/
!-----------------------------------------------------------------------
&namdom        !   time and space domain
!-----------------------------------------------------------------------
   ln_linssh   = .false.  !  =T  linear free surface  ==>>  model level are fixed in time
   rn_Dt       =   400.   !  time step for the dynamics and tracer
   ln_meshmask = .false.  !  =T create a mesh file
/
!-----------------------------------------------------------------------
&namcfg        !   parameters of the configuration                      (default: use namusr_def in namelist_cfg)
!-----------------------------------------------------------------------
   ln_read_cfg = .true.   !  (=T) read the domain configuration file
      cn_domcfg = "domain_cfg"  ! domain configuration filename
/
!-----------------------------------------------------------------------
&namtsd        !    Temperature & Salinity Data  (init/dmp)             (default: OFF)
!-----------------------------------------------------------------------
   !                       ! =T  read T-S fields for: 
   ln_tsd_init = .true.          !  ocean initialisation 
    
   cn_dir      = '../GEO_LR/'            !  root directory for the T-S data location 
   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 
   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 
   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   ! 
   sn_tem      = 'physinit'         ,      -12          ,'votemper' ,  .false.     , .true.,  'yearly'  ,    ''            ,    ''    ,    '' 
   sn_sal      = 'physinit'         ,      -12          ,'vosaline' ,  .false.     , .true.,  'yearly'  ,    ''            ,    ''    ,    '' 
/
!-----------------------------------------------------------------------
&namwad        !   Wetting and Drying (WaD)                             (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namcrs        !   coarsened grid (for outputs and/or TOP)              (ln_crs =T)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namc1d        !   1D configuration options                             ("key_c1d" default: PAPA station)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namc1d_dyndmp !   U & V newtonian damping                              ("key_c1d" default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namc1d_uvd    !   data: U & V currents                                 ("key_c1d" default: OFF)
!-----------------------------------------------------------------------
/
!!======================================================================
!!            ***  Surface Boundary Condition namelists  ***          !!
!!                                                                    !!
!!   namsbc          surface boundary condition manager                 (default: NO selection)
!!   namsbc_flx      flux               formulation                     (ln_flx     =T)
!!   namsbc_blk      Bulk formulae formulation                          (ln_blk     =T)
!!   namsbc_cpl      CouPLed            formulation                     ("key_oasis3" )
!!   namsbc_sas      Stand-Alone Surface module                         (SAS_SRC  only)
!!   namsbc_iif      Ice-IF: use observed ice cover                     (nn_ice = 1   )
!!   namtra_qsr      penetrative solar radiation                        (ln_traqsr  =T)
!!   namsbc_ssr      sea surface restoring term (for T and/or S)        (ln_ssr     =T)
!!   namsbc_rnf      river runoffs                                      (ln_rnf     =T)
!!   namsbc_apr      Atmospheric Pressure                               (ln_apr_dyn =T)
!!   namsbc_wave     external fields from wave model                    (ln_wave    =T)
!!   namberg         iceberg floats                                     (ln_icebergs=T)
!!======================================================================
!
!-----------------------------------------------------------------------
&namsbc        !   Surface Boundary Condition manager                   (default: NO selection)
!-----------------------------------------------------------------------
   nn_fsbc     = 1         !  frequency of SBC module call
   ln_blk      = .true.    !  Bulk formulation                          (T => fill namsbc_blk )
   ln_traqsr   = .true.    !  Light penetration in the ocean            (T => fill namtra_qsr)
   ln_rnf      = .true. 
   ln_wave     = .false.   !  Activate coupling with wave  (T => fill namsbc_wave)
/
!-----------------------------------------------------------------------
&namsbc_flx    !   surface boundary condition : flux formulation        (ln_flx =T)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namsbc_blk    !   namsbc_blk  generic Bulk formula                     (ln_blk =T)
!-----------------------------------------------------------------------
   !                    !  bulk algorithm :
   ln_NCAR      = .false.    ! "NCAR"      algorithm   (Large and Yeager 2008)
   ln_COARE_3p6 = .false.    ! "COARE 3.6" algorithm   (Edson et al. 2013)
      ln_skin_cs   = .false. !  use the cool-skin parameterization  => use at least nn_iter_algo > 10
      ln_skin_wl   = .false. !  use the warm-layer parameterization => use at least nn_iter_algo > 10
   ln_MFS       = .true.
      !
      rn_zqt      =  2.       !  Air temperature & humidity reference height (m)
   !
   ln_humi_sph = .false. !  humidity "sn_humi" is specific humidity  [kg/kg]
   ln_humi_dpt = .true.  !  humidity "sn_humi" is dew-point temperature [K]
   ln_humi_rlh = .false. !  humidity "sn_humi" is relative humidity     [%]
   ln_tair_pot = .false. !  air temperature read in "sn_tair" is already POTENTIAL TEMPERATURE, NOT ABSOLUTE (ECMWF => ln_tair_pot=.false.)
   cn_dir      = '../GEO_LR/atmo/'      !  root directory for the bulk data location
   !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ !       weights filename    ! rotation ! land/sea mask !
   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'yearly' !                           ! pairing  !    filename   !
   sn_wndi     = 'era5'                  ,         1         ,  'u10'    ,   .true.    , .false. , 'daily'  ,  'bslr_era5_bicub.nc'     ,    ''    ,     ''
   sn_wndj     = 'era5'                  ,         1         ,  'v10'    ,   .true.    , .false. , 'daily'  ,  'bslr_era5_bicub.nc'     ,    ''    ,     ''
   sn_qsr      = 'era5'              ,         1         , 'msdwswrf',   .true.    , .false. , 'daily'  ,  'bslr_era5_bilin.nc'     ,    ''    ,     ''  ! not used with MFS
   sn_qlw      = 'era5'              ,         1         , 'msdwlwrf',   .true.    , .false. , 'daily'  ,  'bslr_era5_bilin.nc'     ,    ''    ,     ''  ! not used with MFS
   sn_tair     = 'era5'                  ,         1         ,  't2m'    ,   .true.    , .false. , 'daily'  ,  'bslr_era5_bilin.nc'     ,    ''    ,     ''
   sn_humi     = 'era5'                  ,         1         ,  'd2m'    ,   .true.    , .false. , 'daily'  ,  'bslr_era5_bilin.nc'     ,    ''    ,     ''
   sn_prec     = 'era5'                  ,         1         ,  'tp'     ,   .true.    , .false. , 'daily'  ,  'bslr_era5_bilin.nc'     ,    ''    ,     ''
   sn_snow     = 'NOT USED'              ,         1         ,  'snow'   ,   .true.    , .false. , 'daily'  ,  'bslr_era5_bilin.nc'     ,    ''    ,     ''
   sn_slp      = 'era5'                  ,         1         ,  'msl'    ,   .true.    , .false. , 'daily'  ,  'bslr_era5_bicub.nc'     ,    ''    ,     ''
   sn_cc       = 'era5'                  ,         1         ,  'tcc'    ,   .true.    , .false. , 'daily'  ,  'bslr_era5_bilin.nc'     ,    ''    ,     ''
   !
   rn_pfac     = 0.27777778 ! multiplicative factor for precipitation (total & snow) : for ERA5 [m/hour -->  1000/3600 kg/s/m2];  MAR [mm/hour --> 1/3600 kg/s/m2]
   rn_msl_fac  =  0.01      ! convert Pa to hPA=millibar : 0.01 for ERA5 (given in Pa), 1.0 for MAR (given in hPa)
   rn_clc_fac  =  1.00      ! TCC data needs to be between 0.0 and 1.0 after multiplying with rn_clc_fac, which is already the case for ERA5 and MAR
   rn_windmult =  0.4       ! additional percentage of sensible and latent heat, but NOT of evaporation   
/
!-----------------------------------------------------------------------
&namsbc_cpl    !   coupled ocean/atmosphere model                       ("key_oasis3")
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namsbc_sas    !   Stand-Alone Surface module: ocean data               (SAS_SRC  only)
!-----------------------------------------------------------------------
   l_sasread   = .false.    !  =T Read in file ;  =F set all to 0. (see sbcssm)
/
!-----------------------------------------------------------------------
&namsbc_iif    !   Ice-IF : use observed ice cover                      (nn_ice = 1)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namtra_qsr    !   penetrative solar radiation                          (ln_traqsr =T)
!-----------------------------------------------------------------------
   !                       !  type of penetration                        (default: NO selection)
   ln_qsr_3bd  = .true.      !  2BD light penetration (two bands)
   rn_fracA    =  0.54        !  split between IR and Visible: % of IR
   rn_fracB    =  0.63        !  split visible between short and long: % of long wavelengths
   rn_LAbsIR   =  4.          !  IR absorption coefficient
   rn_LAbsVisL =  0.129       !  visible-long abs coef  (used to be 0.23)
   rn_LAbsVisS =  0.377       !  visible-short abs coef (used to be 0.02 or 0.09)
/
!---------------------------------------------------------------------------------
&nam_RADTRANS      ! radiative transfer module -- # define #radtrans in bamhbi.h90
!---------------------------------------------------------------------------------
! Options
!
wbedges        = 187.5,312.5,337.5,362.5,387.5,412.5,437.5,462.5,487.5,512.5,537.5,562.5,587.5,612.5,637.5,662.5,687.5,712.5,750,800,900,1000,1100,1200,1300,1400,1500,1600,1700,1800,2000,2400,3400,4000
!wbRefWLs      =     250,  325,  350,  375,  400,  425,  450,  475,  500,  525,  550,  575,  600,  625,  650,  675,  700,  725, 775,850,950,1050,1150,1250,1350,1450,1550,1650,1750,1900,2200,2900,3700
ln_solar_eqs   = .false.
!
! from Ackleson et al, 1994 (JGR)
refract_water  = 1.341
rmud_max       = 1.5
rd             = 1.5
ru             = 3.0
rmus           = 1.20481927710843373493 ! 1.0/0.83
rmuu           = 2.5                    ! 1.0/0.4
Qrs            = 3.5
rt_kmax        = 31
useMeanCosSolz = .FALSE.
sfcIrrThresh   = 1.e-4
!
darwin_part_size_p = 1.0e-15
darwin_bbmin = 0.0002
darwin_bbw = 0.5
darwin_lambda_aCDOM = 412.0         ! reference wavelength for CDOM absorption coefficient
darwin_Sdom = 0.020                  ! Elena Terzic et al (2021) estimate between  0.015 and 0.02 [/nm], or a mean value 0.017
darwin_aCDOM_fac = 0.2
darwin_rCDOM = 0.0
darwin_RPOC = 0.0
darwin_allomSpectra = .TRUE.   ! read element size (µm), in our case its present in the file
!
!SPECTRAL DATA FILES
darwin_waterAbsorbFile     = '../GEO_LR/radtrans/spectra_water.dat'
darwin_phytoAbsorbFile     = '../GEO_LR/radtrans/spectra_plankton.dat'
darwin_particleAbsorbFile  = '../GEO_LR/radtrans/spectra_particles.dat'
!
! IRRIADIANCE      !  file name               ! frequency (hours) ! variable       ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! Sea over !
!  FILES           !                          !  (if <0  months)  !   name         !   (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  !  Land    !
sn_irradiance(1)   =  'SW_spectral.g04'       ,         1      , 'PSW_SPECTRAL'     ,    .true.    , .false. , 'monthly' ,'weights_mar_nemo_bilin', '',   ''     !
sn_irradiance(2)   =  'SW_spectral.g04'       ,         1      , 'PSW_SPECTRAL_DIR' ,    .true.    , .false. , 'monthly' ,'weights_mar_nemo_bilin', '',   ''     !
cn_dir_irrad='/gpfs/projects/acad/bsmfc/MAR/v3.14-ecRAD/EUg/g04/spectral/'
!
! ABS CDOM
cdomAbsorbFile = '../GEO_LR/radtrans/cdom_sinusoidal.dat'
ndens = 15
ns = 3
/
!-----------------------------------------------------------------------
&namsbc_ssr    !   surface boundary condition : sea surface restoring   (ln_ssr =T)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namsbc_rnf    !   runoffs                                              (ln_rnf =T)
!-----------------------------------------------------------------------
   ln_rnf_depth = .true.      !  read in depth information for runoff
   ln_rnf_tem   = .true.      !  read in temperature information for runoff
   ln_rnf_sal   = .true.      !  read in salinity information for runoff
   !
   cn_dir      = '../GEO_LR/runoff/'     !  root directory for the runoff data location
   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
   !           !  file name              ! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
   !           !                         !  (if <0  months)  !   name    !   (logical) !  (T/F) ! 'monthly' !                  ! pairing  !    filename   !
   sn_rnf      = 'rnf'              ,        -1         , 'sorunoff',   .true.    , .false., 'yearly'  , ''               , ''       , ''
   sn_s_rnf    = 'rnf'              ,        -1         , 'rosaline',   .false.   , .false., 'yearly'  , ''               , ''       , ''
   sn_t_rnf    = 'rnf'              ,        -1         , 'rotemper',   .true.    , .false., 'yearly'  , ''               , ''       , ''
   sn_dep_rnf  = 'rnf'              ,        -1         , 'rodepth' ,   .false.   , .false., 'yearly'  , ''               , ''       , ''
/
!-----------------------------------------------------------------------
&namsbc_apr    !   Atmospheric pressure used as ocean forcing           (ln_apr_dyn =T)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namsbc_isf    !  Top boundary layer (ISF)                              (ln_isfcav =T : read (ln_read_cfg=T) 
!-----------------------------------------------------------------------             or set or usr_def_zgr )
/
!-----------------------------------------------------------------------
&namsbc_iscpl  !   land ice / ocean coupling option                     (ln_isfcav =T : read (ln_read_cfg=T) 
!-----------------------------------------------------------------------             or set or usr_def_zgr )
/
!-----------------------------------------------------------------------
&namsbc_wave   ! External fields from wave model                        (ln_wave=T)
!-----------------------------------------------------------------------
   ln_sdw  	= .false.   	!  get the 2D Surf Stokes Drift & Compute the 3D stokes drift
   ln_stcor	= .false.   	!  add Stokes Coriolis and tracer advection terms
   ln_cdgw 	= .false.   	!  Neutral drag coefficient read from wave model
   ln_tauoc	= .false.   	!  ocean stress is modified by wave induced stress
   ln_wave_test= .false.   	!  Test case with constant wave fields
!
   ln_charn	= .false.   	!  Charnock coefficient read from wave model (IFS only)
   ln_taw  	= .false.   	!  ocean stress is modified by wave induced stress (coupled mode)
   ln_phioc	= .false.   	!  TKE flux from wave model
   ln_bern_srfc= .false.   	!  wave induced pressure. Bernoulli head J term
   ln_breivikFV_2016 = .true. !  breivik 2016 vertical stokes profile
   ln_vortex_force = .false.   !  Vortex Force term
   ln_stshear  = .false.   	!  include stokes shear in EKE computation
!
   cn_dir  	= '../GEO_LR/waves/'  	!  root directory for the waves data location
   !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
   !       	!  file name          	! frequency (hours) ! variable  ! time interp.!  clim  ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
   !       	!                     	!  (if <0  months)  !   name	!   (logical) !  (T/F) ! 'monthly' !              	! pairing  !	filename   !
   sn_usd  	=  'WAM'	,    	1.     	, 'VSDX' 	,  .true.  , .false. , 'yearly'  ,  'bslr_WAVES_bilin.nc'          	, ''   	, ''                        	! Stokes-Drift velocity in U (x,i) direction [m/s]
   sn_vsd  	=  'WAM'	,    	1.     	, 'VSDY' 	,  .true.  , .false. , 'yearly'  ,  'bslr_WAVES_bilin.nc'          	, ''   	, ''                        	! Stokes-Drift velocity in V (y,j) direction [m/s]
   sn_hsw  	=  'WAM'	,    	1.     	, 'VHM0'     	,  .true.  , .false. , 'yearly'  ,  'bslr_WAVES_bilin.nc'          	, ''   	, ''                    	! [Significant Wave Height - m]
   sn_wmp  	=  'WAM'	,    	1.     	, 'VTPK'    	,  .true.  , .false. , 'yearly'  ,  'bslr_WAVES_bilin.nc'          	, ''   	, ''                     	! [Mean wave period - s]
   sn_tauoc	=  'WAM'	,    	1.     	, 'NWS',  .true.  , .false. , 'yearly'  ,  'bslr_WAVES_bilin.nc'          	, ''   	, ''                              	! [Normalized wave stress - no unit as it is normalized]
/
!-----------------------------------------------------------------------
&namberg       !   iceberg parameters                                   (default: OFF)
!-----------------------------------------------------------------------
/

!!======================================================================
!!               ***  Lateral boundary condition  ***                 !!
!!                                                                    !!
!!   namlbc        lateral momentum boundary condition                  (default: NO selection)
!!   namagrif      agrif nested grid   (read by child model only)       ("key_agrif")
!!   nam_tide      Tidal forcing                                        (default: OFF)
!!   nambdy        Unstructured open boundaries                         (default: OFF)
!!   nambdy_dta    Unstructured open boundaries - external data         (see  nambdy)
!!   nambdy_tide   tidal forcing at open boundaries                     (default: OFF)
!!======================================================================
!
!-----------------------------------------------------------------------
&namlbc        !   lateral momentum boundary condition                  (default: NO selection)
!-----------------------------------------------------------------------
   rn_shlat    =  2.   !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat
   ln_vorlat   = .true.   !  consistency of vorticity boundary condition with analytical Eqs.  ! COPIED FROM DIANA
/
!-----------------------------------------------------------------------
&namagrif      !  AGRIF zoom                                            ("key_agrif")
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&nam_tide      !   tide parameters                                      (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&nambdy        !  unstructured open boundaries                          (default: OFF)
!-----------------------------------------------------------------------
   ln_bdy         = .true.    !  Use unstructured open boundaries
   nb_bdy         = 1         !  number of open boundary sets
   ln_coords_file = .false.   !  =T : read bdy coordinates from file
   ln_mask_file   = .false.   !  =T : read mask from file
   cn_dyn2d    = 'twolayers'  !
   nn_dyn2d_dta   =  1        !  = 0, bdy data are equal to the initial state
      !                       !  = 1, bdy data are read in 'bdydata   .nc' files
      !                       !  = 2, use tidal harmonic forcing data from files
      !                       !  = 3, use external data AND tidal harmonic forcing
   cn_dyn3d      =  'twolayers'
   nn_dyn3d_dta  =  1
   cn_tra        =  'twolayers'
   nn_tra_dta    =  0
   ln_tra_dmp    =.false.     !  open boudaries conditions for tracers
   ln_dyn3d_dmp  =.false.     !  open boundary condition for baroclinic velocities
   nn_rimwidth   = 1        !  width of the relaxation zone
   ln_vol        = .false.    !  total volume correction (see nn_volctl parameter)
/
!-----------------------------------------------------------------------
&nambdy_index
!-----------------------------------------------------------------------
    ctypebdy = 'S'
    nbdyind  =  3      ! on the y axis
    nbdybeg  = 10      ! on the x axis: begin of the segment
    nbdyend  = 12      !                 end
/
!-----------------------------------------------------------------------
&nambdy_dta    !  open boundaries - external data                       (see nam_bdy)
!-----------------------------------------------------------------------
!              !  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      !
   bn_ssh  =    'bosph_bdyV_v1d_59levels',    -12    , 'vobtcrty' ,     .false.    , .true.  ,  'yearly' ,    ''    ,   ''     , ''            ! unused
   bn_u2d  =    'bosph_bdyV_v1d_59levels',    -12    , 'vobtcrty' ,     .false.    , .true.  ,  'yearly' ,    ''    ,   ''     , ''            ! unused
   bn_u3d  =    'bosph_bdyV_v1d_59levels',    -12    , 'vomecrty' ,     .false.    , .true.  ,  'yearly' ,    ''    ,   ''     , ''            ! unused
   bn_v2d  =    'bosph_bdyV_v1d_59levels',    -12    , 'vobtcrty' ,     .false.    , .true.  ,  'yearly' ,    ''    ,   ''     , ''
   bn_v3d  =    'bosph_bdyV_v1d_59levels',    -12    , 'vomecrty' ,     .false.    , .true.  ,  'yearly' ,    ''    ,   ''     , ''
   cn_dir  =    '../GEO_LR/'
   ln_full_vel = .false.
/
!-----------------------------------------------------------------------
&nambdy_tide   !  tidal forcing at open boundaries                      (default: OFF)
!-----------------------------------------------------------------------
/

!!======================================================================
!!                ***  Top/Bottom boundary condition  ***             !!
!!                                                                    !!
!!   namdrg        top/bottom drag coefficient                          (default: NO selection)
!!   namdrg_top    top    friction                                      (ln_OFF=F & ln_isfcav=T)
!!   namdrg_bot    bottom friction                                      (ln_OFF=F)
!!   nambbc        bottom temperature boundary condition                (default: OFF)
!!   nambbl        bottom boundary layer scheme                         (default: OFF)
!!======================================================================
!
!-----------------------------------------------------------------------
&namdrg        !   top/bottom drag coefficient                          (default: NO selection)
!-----------------------------------------------------------------------
   ln_non_lin  = .true.   !  non-linear  drag: Cd = Cd0 |U|
/
!-----------------------------------------------------------------------
&namdrg_top    !   TOP friction                                         (ln_OFF =F & ln_isfcav=T)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namdrg_bot    !   BOTTOM friction                                      (ln_OFF =F)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&nambbc        !   bottom temperature boundary condition                (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&nambbl        !   bottom boundary layer scheme                         (default: OFF)
!-----------------------------------------------------------------------
/

!!======================================================================
!!                        Tracer (T-S) namelists                      !!
!!                                                                    !!
!!   nameos        equation of state                                    (default: NO selection)
!!   namtra_adv    advection scheme                                     (default: NO selection)
!!   namtra_ldf    lateral diffusion scheme                             (default: NO selection)
!!   namtra_mle    mixed layer eddy param. (Fox-Kemper param.)          (default: OFF)
!!   namtra_eiv    eddy induced velocity param.                         (default: OFF)
!!   namtra_dmp    T & S newtonian damping                              (default: OFF)
!!======================================================================
!
!-----------------------------------------------------------------------
&nameos        !   ocean Equation Of Seawater                           (default: NO selection)
!-----------------------------------------------------------------------
   ln_eos80    = .true.         !  = Use EOS80
/
!-----------------------------------------------------------------------
&namtra_adv    !   advection scheme for tracer                          (default: NO selection)
!-----------------------------------------------------------------------
   ln_traadv_fct = .true.        !  order 2 (horizontal), 2 (vertical)
      nn_fct_h   =  2            !  =2/4, horizontal 2nd / 4th order          ! Diana uses 4th order
      nn_fct_v   =  2            !  =2/4, vertical   2nd / COMPACT 4th order  ! Diana uses 4th order
/
!-----------------------------------------------------------------------
&namtra_ldf    !   lateral diffusion scheme for tracers                 (default: NO selection)
!-----------------------------------------------------------------------
   !                       !  Operator type:
   ln_traldf_lap   = .true.   !    bilaplacian operator
   ln_traldf_iso   = .false.  !  iso-neutral (standard operator)
   ln_traldf_triad = .true.   !    iso-neutral (triad    operator)
   !		       	         !  iso-neutral options:
   ln_traldf_msc   = .false.   !  Method of Stabilizing Correction      (both operators)
   rn_slpmax       =  0.01     !  slope limit                           (both operators)
   rn_sw_triad     = 0         !  =1 switching triad ; =0 all 4 triads used (triad only)
   nn_aht_ijk_t    = 20        !  space/time variation of eddy coefficient: 20 -> replace rn_Ld with max(e1t,e2t)
      rn_Ud        = 0.001     !  lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30), A.M. : divide standard value 0.01 by 2 => 0.005
/
!-----------------------------------------------------------------------
&namtra_mle    !   mixed layer eddy parametrisation (Fox-Kemper)       (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namtra_eiv    !   eddy induced velocity param.                         (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namtra_dmp    !   tracer: T & S newtonian damping                      (default: OFF)
!-----------------------------------------------------------------------
   ln_tradmp   =  .false.   ! add a damping term (using resto.nc coef.)
      nn_zdmp  =    0         !  vertical shape =0    damping throughout the water column
      !                       !                 =1 no damping in the mixing layer (kz  criteria)
      !                       !                 =2 no damping in the mixed  layer (rho crieria)
      cn_resto = '../GEO_LR/resto.nc'   !  Name of file containing restoration coeff. field (use dmp_tools to create this)
/


!!======================================================================
!!                      ***  Dynamics namelists  ***                  !!
!!                                                                    !!
!!   nam_vvl       vertical coordinate options                          (default: z-star)
!!   namdyn_adv    formulation of the momentum advection                (default: NO selection)
!!   namdyn_vor    advection scheme                                     (default: NO selection)
!!   namdyn_hpg    hydrostatic pressure gradient                        (default: NO selection)
!!   namdyn_spg    surface pressure gradient                            (default: NO selection)
!!   namdyn_ldf    lateral diffusion scheme                             (default: NO selection)
!!   namdta_dyn    offline TOP: dynamics read in files                  (OFF_SRC only)
!!======================================================================
!
!-----------------------------------------------------------------------
&nam_vvl       !   vertical coordinate options                          (default: z-star)
!-----------------------------------------------------------------------
   ln_vvl_zstar  = .true.           !  z-star vertical coordinate
   ln_vvl_dbg    = .false.          !  debug prints    (T/F)
! L.V. 2024-01-11 --- tried z-tilde in experiment LR012
   ln_vvl_ztilde = .false.          !  z-tilde vertical coordinate: only high frequency variations
   rn_rst_e3t    = 30.0             !  ztilde to zstar restoration timescale [days]
   rn_lf_cutoff  =  5.0             !  cutoff frequency for low-pass filter  [days]
   ! add more options such as in https://doi.org/10.1029/2022MS003056
/
!-----------------------------------------------------------------------
&namdyn_adv    !   formulation of the momentum advection                (default: NO selection)
!-----------------------------------------------------------------------
   ln_dynadv_vec = .false. !  vector form - 2nd centered scheme          ! Diana uses this one instead of UBS
   ln_dynadv_ubs = .true.  !  flux form - 3rd order UBS      scheme
/
!-----------------------------------------------------------------------
&namdyn_vor    !   Vorticity / Coriolis scheme                          (default: NO selection)
!-----------------------------------------------------------------------
   ln_dynvor_een = .true. !  energy & enstrophy scheme
/
!-----------------------------------------------------------------------
&namdyn_hpg    !   Hydrostatic pressure gradient option                 (default: NO selection)
!-----------------------------------------------------------------------
   ln_hpg_zps  = .false.  !  z-coordinate - partial steps (interpolation) --> does not work with VVL
   ln_hpg_sco  = .true.
   ln_hpg_prj  = .false.  ! improved (cubic spline) version of sco, TEST ME !
/
!-----------------------------------------------------------------------
&namdyn_spg    !   surface pressure gradient                            (default: NO selection)
!-----------------------------------------------------------------------
   ln_dynspg_ts   = .true.   ! split-explicit free surface
      ln_bt_fw      = .false.    ! Forward integration of barotropic Eqs.      ! Diana uses .true.
      ln_bt_av      = .true.     ! Time filtering of barotropic variables
      ln_bt_auto    = .true.     ! Number of sub-step defined from:           ! Diana uses .false., but recommended .true. as we use
         rn_bt_cmax   =  0.8        ! =T : the Maximum Courant Number allowed
         nn_e         = 100          ! =F : the number of sub-step in rn_Dt seconds
/
!-----------------------------------------------------------------------
&namdyn_ldf    !   lateral diffusion on momentum                        (default: NO selection)
!-----------------------------------------------------------------------
   !                       !  Type of the operator :
   ln_dynldf_OFF = .false.     !  No operator (i.e. no explicit diffusion)
   nn_dynldf_typ = 0           !  =0 div-rot (default)   ;   =1 symmetric
   ln_dynldf_lap = .false.     !    laplacian operator
   ln_dynldf_blp = .true.      !  bilaplacian operator                         ! Diana uses bilaplacian momentum diffusion
   !                       !  Direction of action  :
   ln_dynldf_lev = .false.     !  iso-level
   ln_dynldf_hor = .true.      !  horizontal  (geopotential)
   ln_dynldf_iso = .false.     !  iso-neutral (lap only)
   !                       !  Coefficient
   nn_ahm_ijk_t  = 20          !  space/time variation of eddy coefficient :
      rn_Uv      = 0.01             !  lateral viscous velocity [m/s] (nn_ahm_ijk_t= 0, 10, 20, 30) 
/
!-----------------------------------------------------------------------
&namdta_dyn    !   offline ocean input files                            (OFF_SRC only)
!-----------------------------------------------------------------------
/

!!======================================================================
!!                     vertical physics namelists                     !!
!!                                                                    !!
!!    namzdf        vertical physics manager                            (default: NO selection)
!!    namzdf_ric    richardson number vertical mixing                   (ln_zdfric=T)
!!    namzdf_tke    TKE vertical mixing                                 (ln_zdftke=T)
!!    namzdf_gls    GLS vertical mixing                                 (ln_zdfgls=T)
!!    namzdf_osm    OSM vertical diffusion                              (ln_zdfosm=T)
!!    namzdf_iwm    tidal mixing parameterization                       (ln_zdfiwm=T)
!!======================================================================
!
!-----------------------------------------------------------------------
&namzdf        !   vertical physics manager                             (default: NO selection)
!-----------------------------------------------------------------------
   !                       ! adaptive-implicit vertical advection
   ln_zad_Aimp = .false.     !  Courant number dependent scheme (Shchepetkin 2015)
   !
   !                       ! type of vertical closure (required)
   ln_zdftke   = .false.     !  Turbulent Kinetic Energy closure       (T =>   fill namzdf_tke)
   ln_zdfgls   = .true.       !  GLS
   !                       ! convection
   ln_zdfevd   = .true.      !  enhanced vertical diffusion
      nn_evdm     =    1         ! apply on tracer (=0) or on tracer and momentum (=1)    !  Diana uses 0
      rn_evd      =   1.0        ! mixing coefficient [m2/s]                             !  Diana uses 100
      rn_n2_lim   =   1.e-09     ! absolute value of (instable, negative) N2 frequency below which we apply EVD, default 1.e-12
   !                       ! coefficients
   rn_avm0     =   1.0e-6     !  vertical eddy viscosity   [m2/s]       (background Kz if ln_zdfcst=F)  ! Diana uses 1.2e-5
   rn_avt0     =   1.0e-9     !  vertical eddy diffusivity [m2/s]       (background Kz if ln_zdfcst=F)  ! Diana uses 1.2e-6
/
!-----------------------------------------------------------------------
&namzdf_ric    !   richardson number dependent vertical diffusion       (ln_zdfric =T)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  (ln_zdftke =T)
!-----------------------------------------------------------------------
   rn_ediff    =   0.05    !  coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) )
   rn_emin     =   1.e-14  !  minimum value of tke [m2/s2]
   rn_emin0    =   1.e-5   !  surface minimum value of tke [m2/s2]
   rn_mxl0     =   0.01    !  surface  buoyancy lenght scale minimum value
   ln_lc       = .false.   !  Langmuir cell parameterisation (Axell 2002)
   nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to NIWs
      nn_htau  =   0       !  type of exponential decrease of tke penetration below the ML
/
!-----------------------------------------------------------------------
&namzdf_gls    !   GLS vertical diffusion                               (ln_zdfgls =T)
!-----------------------------------------------------------------------
   rn_emin       = 1.e-08  !  minimum value of e   [m2/s2], standard value is 1.e-07
   rn_epsmin     = 1.e-12  !  minimum value of eps [m2/s3]
   ln_length_lim = .true.  !  limit on the dissipation rate under stable stratification (Galperin et al., 1988)
   rn_clim_galp  = 0.267   !  galperin limit
   ln_sigpsi     = .true.  !  Activate or not Burchard 2001 mods on psi schmidt number in the wb case
   rn_crban      = 100.    !  Craig and Banner 1994 constant for wb tke flux
   rn_charn      = 70000.  !  Charnock constant for wb induced roughness length
   rn_hsro       =  0.02   !  Minimum surface roughness
   rn_hsri       =  0.03   !  Ice-ocean roughness
   rn_frac_hs    =   1.3   !  Fraction of wave height as roughness (if nn_z0_met>1)
   nn_z0_met     =     2   !  Method for surface roughness computation (0/1/2/3)
   !                       !     = 3 requires ln_wave=T
   nn_z0_ice     =   1     !  attenutaion of surface wave breaking under ice
   !                       !           = 0 no impact of ice cover
   !                       !           = 1 roughness uses rn_hsri and is weigthed by 1-TANH(10*fr_i)
   !                       !           = 2 roughness uses rn_hsri and is weighted by 1-fr_i
   !                       !           = 3 roughness uses rn_hsri and is weighted by 1-MIN(1,4*fr_i)
   nn_mxlice     =     1   !  mixing under sea ice
                           !     = 0 No scaling under sea-ice
                           !     = 1 scaling with constant Ice-ocean roughness (rn_hsri)
                           !     = 2 scaling with mean sea-ice thickness
                           !     = 3 scaling with max sea-ice thickness
   nn_bc_surf    =     1   !  surface condition (0/1=Dir/Neum)
   nn_bc_bot     =     1   !  bottom condition (0/1=Dir/Neum)
   nn_stab_func  =     2   !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB)
   nn_clos       =     1   !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
/
!-----------------------------------------------------------------------
&namzdf_osm    !   OSM vertical diffusion                               (ln_zdfosm =T)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namzdf_iwm    !    internal wave-driven mixing parameterization        (ln_zdfiwm =T)
!-----------------------------------------------------------------------
/

!!======================================================================
!!                  ***  Diagnostics namelists  ***                   !!
!!                                                                    !!
!!   namtrd       dynamics and/or tracer trends                         (default: OFF)
!!   namptr       Poleward Transport Diagnostics                        (default: OFF)
!!   namhsb       Heat and salt budgets                                 (default: OFF)
!!   namdiu       Cool skin and warm layer models                       (default: OFF)
!!   namflo       float parameters                                      ("key_float")
!!   nam_diaharm  Harmonic analysis of tidal constituents               ("key_diaharm")
!!   namdct       transports through some sections                      ("key_diadct")
!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF)
!!   nam_dia25h   25h Mean Output                                       (default: OFF)
!!   namnc4       netcdf4 chunking and compression settings             ("key_netcdf4")
!!======================================================================
!
!-----------------------------------------------------------------------
&namtrd        !   trend diagnostics                                    (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namptr        !   Poleward Transport Diagnostic                        (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namhsb        !  Heat and salt budgets                                 (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namdiu        !   Cool skin and warm layer models                      (default: OFF)
!-----------------------------------------------------------------------
   ln_diurnal      = .false.   !
   ln_diurnal_only = .false.   !
/
!-----------------------------------------------------------------------
&namflo        !   float parameters                                     ("key_float")
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&nam_diaharm   !   Harmonic analysis of tidal constituents              ("key_diaharm")
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namdct        ! transports through some sections                       ("key_diadct")
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&nam_diatmb    !  Top Middle Bottom Output                              (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&nam_dia25h    !  25h Mean Output                                       (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namnc4        !   netcdf4 chunking and compression settings            ("key_netcdf4")
!-----------------------------------------------------------------------
/
!!======================================================================
!!               ***  Observation & Assimilation  ***                 !!
!!                                                                    !!
!!   namobs       observation and model comparison                      (default: OFF)
!!   nam_asminc   assimilation increments                               ('key_asminc')
!!======================================================================
!
!-----------------------------------------------------------------------
&namobs        !  observation usage switch                              (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&nam_asminc    !   assimilation increments                              ('key_asminc')
!-----------------------------------------------------------------------
/

!!======================================================================
!!                  ***  Miscellaneous namelists  ***                 !!
!!                                                                    !!
!!   nammpp            Massively Parallel Processing                    ("key_mpp_mpi")
!!   namctl            Control prints                                   (default: OFF)
!!   namsto            Stochastic parametrization of EOS                (default: OFF)
!!======================================================================
!
!-----------------------------------------------------------------------
&nammpp        !   Massively Parallel Processing                        ("key_mpp_mpi")
!-----------------------------------------------------------------------
   jpni        =   0       !  number of processors following i (set automatically if < 1), see also ln_listonly = T
   jpnj        =   0       !  number of processors following j (set automatically if < 1), see also ln_listonly = T
   nn_hls      =   1       !  halo width (applies to both rows and columns)
   nn_comm     =   1       !  comm choice : 1=old way, 2=new way with neighbourhoods ! TEST IT !
/
!-----------------------------------------------------------------------
&namctl        !   Control prints                                       (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namsto        ! Stochastic parametrization of EOS                      (default: OFF)
!-----------------------------------------------------------------------
   ln_rst_clock  = .false.       ! generate seed from (T) system clock or (F) from hard-coded values
   ln_rststo     = .false.       ! start from mean parameter (F) or from restart file (T)
   ln_rstseed    = .true.        ! read seed of RNG from restart file
   cn_storst_in  = "restart_sto" !  suffix of stochastic parameter restart file (input)
   cn_storst_out = "restart_sto" !  suffix of stochastic parameter restart file (output)

   ln_sto_eos  = .false.   ! stochastic equation of state
   nn_sto_eos  = 1         ! number of independent random walks
   rn_eos_stdxy = 1.4      ! V.V. 1.0   ! random walk horz. standard deviation (in grid points)
   rn_eos_stdz = 0.7       ! V.V. 0.5   ! random walk vert. standard deviation (in grid points)
   rn_eos_tcor = 1440.     ! V.V. 2880. ! random walk time correlation (in timesteps)
   nn_eos_ord  = 1         ! order of autoregressive processes
   nn_eos_flt  = 0         ! passes of Laplacian filter
   rn_eos_lim  = 2.0       ! limitation factor (default = 3.0)

   ln_sto_trd = .false.          ! stochastic trend
   rn_trd_std   = 0.6            ! trend standard deviation
   rn_trd_tcor  = 5760.0         ! trend time correlation (in timesteps)
   rn_trd_xycor = 18.            ! trend spatial correlation (radius in grid points activated > 1. ; default=1.)
   nn_trd_flt  = 0               ! passes of Laplacian filter
   rn_trd_lim  = 1.0             ! limitation factor (default = 3.0)

   ln_sto_cdo  = .false.         ! stochastic ocean drag coefficient momentum Cd >>>[in NEMO3.6 sbcblk_core.F90]<<<
   ln_sto_turb = .false.         ! stochastic ocean drag coefficient turbelence Ch, Ce (activation dependency with Cd i.e. ln_sto_cdo  = .true.)
   rn_cdo_std   = 0.1            ! ocean drag coef. standard deviation
   rn_cdo_tcor  = 1728.0         ! ocean drag coef. time correlation (in timesteps)
   rn_cdo_xycor = 18.            ! ocean drag coef. spatial correlation (radius in grid points activated > 1. ; default=1.)
   nn_cdo_ord  = 1               ! order of autoregressive processes
   nn_cdo_flt  = 0               ! passes of Laplacian filter
   rn_cdo_lim  = 1.0             ! limitation factor (default = 3.0)

   ln_sto_bfr  = .false.         ! stochastic ocean bottom friction coefficient >>>[in NEMO3.6 zdfbfr.F90, in NEMO4.0 zdfdrg.F90]<<<
   rn_bfr_std   = 0.2            ! ocean bottom friction coef. standard deviation
   rn_bfr_tcor  = 17280.0        ! ocean bottom friction coef. time correlation (in timesteps)
   rn_bfr_xycor = 1.             ! ocean bottom friction coef. spatial correlation (radius in grid points activated > 1. ; default = 1.)
   nn_bfr_ord  = 1               ! order of autoregressive processes
   nn_bfr_flt  = 100             ! passes of Laplacian filter
   rn_bfr_lim  = 1.0             ! limitation factor (default = 3.0)

   ln_sto_rnf  = .false.         ! multiplicative stochastic river runoff                     ! L.V.: implemented in sbcrnf
   rn_rnf_std   = 0.2            ! ocean bottom friction coef. standard deviation
   rn_rnf_tcor  = 52560.0        ! ocean bottom friction coef. time correlation (in timesteps)
   rn_rnf_xycor = 72             ! ocean bottom friction coef. spatial correlation (radius in grid points activated > 1. ; default = 1.)
   nn_rnf_ord  = 1               ! order of autoregressive processes
   nn_rnf_flt  = 0               ! passes of Laplacian filter
   rn_rnf_lim  = 2.0             ! limitation factor (default = 3.0)

   ln_sto_wnd  = .false.         ! stochastic u,v wind forcing                                ! L.V.: implemented in sbcblk_algo_mfs
   rn_wnd_std   = 0.3            ! u,v wind forcing standard deviation
   rn_wnd_tcor  = 1728.0         ! u,v wind forcing time correlation (in timesteps)
   rn_wnd_xycor = 36.            ! u,v wind forcing spatial correlation (radius in grid points activated > 1. ; default = 1.)
   nn_wnd_ord  = 1               ! order of autoregressive processes
   nn_wnd_flt  = 0               ! passes of Laplacian filter
   rn_wnd_lim  = 1.0             ! limitation factor (default = 3.0)

   ln_sto_tair  = .false.        ! stochastic air temperature                                 ! L.V.: implemented in sbcblk_algo_mfs
   rn_tair_std   = 0.1           ! air temperature standard deviation
   rn_tair_tcor  = 8640.0        ! air temperature time correlation (in timesteps)
   rn_tair_xycor = 72.           ! air temperature spatial correlation (radius in grid points activated > 1. ; default = 1.)
   nn_tair_ord  = 1              ! order of autoregressive processes
   nn_tair_flt  = 0              ! passes of Laplacian filter
   rn_tair_lim  = 1.0            ! limitation factor (default = 3.0)

   ln_sto_slp  = .false.         ! stochastic sea level pressure >>>[in NEMO3.6 sbcapr.F90]<<<
   rn_slp_std   = 0.01           ! sea level pressure standard deviation
   rn_slp_tcor  = 2880.0         ! sea level pressure time correlation (in timesteps)
   rn_slp_xycor = 108.           ! sea level pressure spatial correlation (radius in grid points activated > 1. ; default = 1.)
   nn_slp_ord  = 1               ! order of autoregressive processes
   nn_slp_flt  = 0               ! passes of Laplacian filter
   rn_slp_lim  = 1.0             ! limitation factor (default = 3.0)

   ln_sto_wnd_eofs   = .false.   ! stochastic u,v wind forcing based on EOFs                  ! L.V.: implemented in sbcblk_algo_mfs
   nn_wnd_eofs       = 100       ! # of EOFs to use
   rn_wnd_eofs_std   = 0.003     ! u,v wind perturbation standard deviation (taking EOF magnitude into account)
   rn_wnd_eofs_tcor  = 4032.0    ! u,v wind perturbation time correlation (in timesteps) : 4032 = 7*576 = 1 week
   nn_wnd_eofs_ord   = 1         ! order of autoregressive processes
   rn_wnd_eofs_lim   = 5.0       ! limitation factor (default = 3.0): factors will be limited in +/- std*lim
   cn_wnd_eofs_file  = "../GEO/wnd_eofs"
   cn_wnd_eofs_vars  = "u10,v10"

   ln_sto_tair_eofs  = .false.   ! stochastic tair forcing based on EOFs
   nn_tair_eofs      = 100       ! # of EOFs to use
   rn_tair_eofs_std  = 0.0015    ! T2M perturbation standard deviation, 0.0015 leads to std of .87°C, non Gaussian (?)
   rn_tair_eofs_tcor = 2880.0    ! T2M perturbation time correlation (in timesteps) : 2880 = 5*576
   nn_tair_eofs_ord  = 1         ! order of autoregressive processes
   rn_tair_eofs_lim  = 5.0       ! limitation factor (default = 3.0): factors will be limited in +/- std*lim
   cn_tair_eofs_file = "../GEO/tair_eofs"
   cn_tair_eofs_vars = "t2m"

   ln_sto_tcc_eofs   = .false.   ! stochastic tcc forcing based on EOFs
   nn_tcc_eofs       = 100       ! # of EOFs to use
   rn_tcc_eofs_std   = 0.0015    ! TCC perturbation standard deviation, 0.0015 leads to std of 6%
   rn_tcc_eofs_tcor  = 1440.0    ! TCC perturbation time correlation (in timesteps) : 2.5 days
   nn_tcc_eofs_ord   = 1         ! order of autoregressive processes
   rn_tcc_eofs_lim   = 5.0       ! limitation factor (default = 3.0): factors will be limited in +/- std*lim
   cn_tcc_eofs_file  = "../GEO/tair_eofs"
   cn_tcc_eofs_vars  = "tcc"

   ln_sto_prec_eofs  = .false.   ! stochastic precipitation wind forcing based on EOFs
   nn_prec_eofs      = 100       ! # of EOFs to use
   rn_prec_eofs_std  = 0.015     ! precip perturbation standard deviation, 0.015 leads to std=4e-5
   rn_prec_eofs_tcor = 576.0     ! precip perturbation time correlation (in timesteps) : 576 = 1 day
   nn_prec_eofs_ord  = 1         ! order of autoregressive processes
   rn_prec_eofs_lim  = 5.0       ! limitation factor (default = 3.0): factors will be limited in +/- std*lim
   cn_prec_eofs_file = "../GEO/tair_eofs"
   cn_prec_eofs_vars = "tp"
/


