!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! OPA namelists :   1 - run manager      (namrun, nam_ctl, nam_mpp)
!!                   2 - Domain           (nam_zgr, nam_zgr_sco, namdom)
!!                   3 - Surface boundary (namtau, namflx, namqsr, namfwb, namalb, namcpl)
!!                   4 - lateral boundary (namlbc, namobc, namcla)
!!                   5 - bottom  boundary (nambfr, nambbc, nambbl)
!!                   6 - Tracer           (nam_traadv, nam_traldf, namtdp)
!!                   7 - dynamics         (nam_dynvor, nam_dynhpg, namcfg, nam_dynldf)
!!                   8 - Verical physics  (namzdf, namnpc, namgap, namznl, namspr, namflo)
!!                   9 - diagnostics      (namtrd, namgap, namznl, namptr, namspr, namflo)
!!                   9 - miscellaneous    (solver)
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!  CAUTION: some scripts does not support CAPITALs for logical definition 
!  *******       use .true. or .false. and NOT .TRUE. or .FALSE.

!!======================================================================
!!                   ***  Run management namelists  ***
!!======================================================================
!!   namrun       parameters of the run
!!   nam_ctl      Control prints & Benchmark
!!   nam_mpp      Massively Parallel Processing
!!======================================================================

!-----------------------------------------------------------------------
&namrun   !   parameters of the run
!-----------------------------------------------------------------------
   no         =       0    !  job number
   cexper     =  "ORCA2"   !  experience name for vairmer format
   ln_rstart  = .false.    !  boolean term for restart (true or false)
   nrstdt     =       0    !  restart control = 0 restart, do not control nit000 in the restart file.
                           !                  = 1 restart, control nit000 in the restart file. Do not
                           !                      use the date in the restart file (use ndate0 in namelist)
                           !                  = 2 restart, control nit000 in the restart file, use the date
                           !                      in the restart file. ndate0 in the namelist is ignored.
   nit000     =       1    !  number of the first time step
   nitend     =    5475    !  number of the first time step
   ndate0     =  010101    !  initial calendar date aammjj
   nleapy     =       0    !  Leap year calendar (0/1)
   ninist     =       0    !  initial state output flag (0/1)
   nstock     =    5475    !  number of the first time step
   nwrite     =    5475    !  frequency of OUTPUT file
   nrunoff    =       2    !  = 0 no, 1 runoff, 2 runoff+river mouth ups adv
/
!-----------------------------------------------------------------------
&nam_ctl      !   Control prints & Benchmark
!-----------------------------------------------------------------------
   ln_ctl =  .false.       !  trends control print (expensive!)
   nprint =       0        !  level of print (0 no print)
   nictls =       0        !  start i indice to make the control SUM (very usefull to compare mono-
   nictle =       0        !  end   i indice to make the control SUM (-versus multi processor runs)
   njctls =       0        !  start j indice to make the control SUM (very usefull to compare mono-
   njctle =       0        !  end   j indice to make the control SUM (-versus multi processor runs)
   isplt  =       1        !  number of processors following i
   jsplt  =       1        !  number of processors following j
   nbench =       0        !  Bench parameter (0/1): CAUTION it must be zero except for bench
                           !     for which we don't care about physical meaning of the results
   nbit_cmp =     0        !  bit comparison mode parameter (0/1): enables bit comparison between 
                           !             single and multiple processor runs
/
!-----------------------------------------------------------------------
&nam_mpp      !   Massively Parallel Processing
!-----------------------------------------------------------------------
   c_mpi_send =  'S'       !  mpi send/recieve type   ='S', 'B', or 'I' for standard send,
                           !  buffer blocking send or immediate non-blocking sends, resp.
/

!!======================================================================
!!                      ***  Domain namelists  ***
!!======================================================================
!!   nam_zgr       vertical coordinate
!!   nam_zgr_sco   s-coordinate or hybrid z-s-coordinate
!!   namdom        space and time domain (bathymetry, mesh, timestep)
!!======================================================================
!-----------------------------------------------------------------------
&nam_zgr       !   vertical coordinate
!-----------------------------------------------------------------------
   ln_zco   =  .false.     !  z-coordinate - full    steps   (T/F)      ("key_zco" may also be defined)
   ln_zps   =  .true.      !  z-coordinate - partial steps   (T/F)
   ln_sco   =  .false.     !  s- or hybrid z-s-coordinate    (T/F)
/
!-----------------------------------------------------------------------
&nam_zgr_sco   !   s-coordinate or hybrid z-s-coordinate
!-----------------------------------------------------------------------
   sbot_min =  300.        !  minimum depth of s-bottom surface (>0) (m)
   sbot_max = 5250.        !  maximum depth of s-bottom surface (= ocean depth) (>0) (m)
   theta    =    6.0       !  surface control parameter (0<=theta<=20)
   thetb    =    0.75      !  bottom control parameter  (0<=thetb<= 1)
   r_max    =    0.15      !  maximum cut-off r-value allowed (0<r_max<1)
/
!-----------------------------------------------------------------------
&namdom   !   space and time domain (bathymetry, mesh, timestep)
!-----------------------------------------------------------------------
   ntopo     =     1       !  = 0/1 ,compute/read the bathymetry file (mbathy, nb of T-ocean levels)
   e3zps_min =     5.      !  the thickness of the partial step is set larger than the minimum
   e3zps_rat =     0.1     !  of e3zps_min and e3zps_rat * e3t   (N.B. 0<e3zps_rat<1)
   nmsh      =     0       !  =1 create a mesh file (coordinates, scale factors, masks)
   nacc      =     0       !  the acceleration of convergence method   = 0, no acceleration, rdt = rdttra
                           !                                           = 1, acceleration used, rdt < rdttra(k)
   atfp      =     0.1     !  asselin time filter parameter
   rdt       =  5760.      !  time step for the dynamics (and tracer if nacc=0)
   rdtmin    =  5760.      !  minimum time step on tracers
   rdtmax    =  5760.      !  maximum time step on tracers
   rdth      =   800.      !  depth variation of tracer time step
   rdtbt     =    90.      !  barotropic time step (for the time splitting algorithm)
   nfice     =     5       !  frequency of ice model call
   nfbulk    =     5       !  frequency of bulk formulea call (not used if ice used)
   nclosea   =     0       !  = 0 no closed sea in the model domain
                           !  = 1 closed sea (Black Sea, Caspian Sea, Great US Lakes...) 
/

!!======================================================================
!!            ***  Surface Boundary Condition namelists  ***
!!======================================================================
!!   namtau      surface wind stress
!!   namflx      surface fluxes
!!   namqsr      penetrative solar radiation
!!   namfwb      freshwater budget correction
!!   namalb      albedo parameters
!!   namcpl      coupled ocean/atmosphere model (#ifdef "key_coupled")
!!======================================================================
!-----------------------------------------------------------------------
&namtau   !   surface wind stress
!-----------------------------------------------------------------------
   ntau000 =      0        !  gently increase the stress over the first ntau_rst time-steps
   tau0x   =      0.e0     !  uniform value used as default surface heat flux
   tau0y   =      0.e0     !  uniform value used as default solar radiation flux
/
!-----------------------------------------------------------------------
&namflx   !   surface fluxes
!-----------------------------------------------------------------------
   q0      =      0.e0     !  uniform value used as default surface heat flux
   qsr0    =      0.e0     !  uniform value used as default solar radiation flux
   emp0    =      0.e0     !  uniform value used as default surface freswater budget (E-P)
   dqdt0   =      -40.     !  feedback coefficient for SST damping (W/m2/K)
   deds0   =        0.     !  feedback coefficient for SSS damping (mm/day)
/
!-----------------------------------------------------------------------
&namqsr   !   penetrative solar radiation
!-----------------------------------------------------------------------
   ln_traqsr = .true.      !  penetrative solar radiation (T) or not (F)     (Default=T)
   rabs     =   0.58       !  fraction of qsr associated with xsi1
   xsi1     =   0.35       !  first depth of extinction
   xsi2     =   23.0       !  second depth of extinction
/
!-----------------------------------------------------------------------
&namfwb   !   freshwater budget correction
!-----------------------------------------------------------------------
   ln_fwb    = .true.      !  flag for freshwater budget correction (0 annual mean)
/
!-----------------------------------------------------------------------
&namalb   !   albedo parameters
!-----------------------------------------------------------------------
   cgren    =      0.06    !  correction of the snow or ice albedo to take into account
   albice   =      0.5     !  albedo of melting ice in the arctic and antarctic
   alphd    =      0.80    !  coefficients for linear interpolation used to
   alphc    =      0.65    !  compute albedo between two extremes values 
   alphdi   =      0.72    !  (Pyane, 1972)
/
!-----------------------------------------------------------------------
&namcpl   !   coupled ocean/atmosphere model                            ("key_coupled")
!-----------------------------------------------------------------------
/

!!======================================================================
!!               ***  Lateral boundary condition  ***
!!======================================================================
!!   namlbc        lateral momentum boundary condition
!!   namcro        cross land advection
!!   namobc        open boundaries parameters                           ("key_obc")
!!======================================================================

!-----------------------------------------------------------------------
&namlbc   !   lateral momentum boundary condition
!-----------------------------------------------------------------------
   shlat  =      2.        !      shlat = 0 : free slip
                           !  0 < shlat < 2 : partial slip
                           !      shlat = 2 : no slip
                           !  2 < shlat     : strong slip
/
!-----------------------------------------------------------------------
&namcro   !   cross land advection
!-----------------------------------------------------------------------
!  n_cla   advection between 2 ocean pts separates by land 
   n_cla   = 1
/
!-----------------------------------------------------------------------
&namobc   !   open boundaries parameters                            ("key_obc")
!-----------------------------------------------------------------------
    nobc_dta    =    0     !  = 0 the obc data are equal to the initial state
                           !  = 1 the obc data are read in 'obc   .dta' files
    rdpein      =    1.    !  ???
    rdpwin      =    1.    !  ???
    rdpnin      =   30.    !  ???
    rdpsin      =    1.    !  ???
    rdpeob      = 1500.    !  time relaxation (days) for the east  open boundary
    rdpwob      =   15.    !    "        "           for the west  open boundary
    rdpnob      =  150.    !    "        "           for the north open boundary
    rdpsob      =   15.    !    "        "           for the south open boundary
    zbsic1      =  140.e+6 !  barotropic stream function on first  isolated coastline
    zbsic2      =    1.e+6 !    "                   "    on second isolated coastline
    zbsic3      =    0.    !    "                   "    on thrid  isolated coastline
    ln_obc_clim = .true.   !  climatological obc data files (default T)
    ln_vol_cst  = .true.   !  total volume conserved
/

!!======================================================================
!!                 ***  Bottom boundary condition  ***
!!======================================================================
!!   nambfr        bottom friction
!!   nambbc        bottom temperature boundary condition                ("key_trabbc")
!!======================================================================

!-----------------------------------------------------------------------
&nambfr   !   bottom friction
!-----------------------------------------------------------------------
   nbotfr =       1        !  type of bottom friction :   = 0 : no slip        ,  = 2 : nonlinear friction
                           !                              = 1 : linear friction,  = 3 : free slip
   bfri1  =   4.e-4        !  bottom drag coefficient (linear case)
   bfri2  =   1.e-3        !  bottom drag coefficient (non linear case)
   bfeb2  =  2.5e-3        !  bottom turbulent kinetic energy background  (m^2/s^2)
/
!-----------------------------------------------------------------------
&nambbc   !   bottom temperature boundary condition
!-----------------------------------------------------------------------
   ngeo_flux =  2             !  = 0 no geothermal heat flux
                              !  = 1 constant geothermal heat flux
                              !  = 2 variable geothermal heat flux (read in geothermal_heating.nc in mW/m2) 
   ngeo_flux_const = 86.4e-3  !  Constant value of geothermal heat flux (W/m2)
/

!!======================================================================
!!                        Tracer (T & S ) namelists
!!======================================================================
!!   nameos        equation of state
!!   nam_traadv    advection scheme
!!   nam_traldf    lateral diffusion scheme
!!   nambbl        bottom boundary layer scheme                         ("key_trabbl_dif","key_trabbl_adv")
!!   namdmp        tracer newtonian damping                             ("key_tradmp")
!!======================================================================

!-----------------------------------------------------------------------
&nameos   !   ocean physical parameters
!-----------------------------------------------------------------------
   neos   =      0         !  type of equation of state and Brunt-Vaisala frequency
                           !     = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) )
                           !     = 1, linear: rho(T)   = rau0 * ( 1.028 - ralpha * T )
                           !     = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T )
   ralpha =  2.e-4         !  thermal expension coefficient (neos= 1 or 2)
   rbeta  =  0.001         !  saline  expension coefficient (neos= 2)
/
!-----------------------------------------------------------------------
&nam_traadv   !   advection scheme for tracer 
!-----------------------------------------------------------------------
   ln_traadv_cen2   =  .true.    !  2nd order centered scheme    (default T)
   ln_traadv_tvd    =  .false.   !  TVD scheme                   (default F)
   ln_traadv_muscl  =  .false.   !  MUSCL scheme                 (default F)
   ln_traadv_muscl2 =  .false.   !  MUSCL2 scheme                (default F)
   ln_traadv_ubs    =  .false.   !  UBS scheme                   (default F)
/
!-----------------------------------------------------------------------
&nam_traldf   !   lateral diffusion scheme for tracer 
!-----------------------------------------------------------------------
!                                !  Type of the operator : 
   ln_traldf_lap    =  .true.    !     laplacian operator          (default T)
   ln_traldf_bilap  =  .false.   !     bilaplacian operator        (default F)
                                 !  Direction of action  :
   ln_traldf_level  =  .false.   !     iso-level                   (default F)
   ln_traldf_hor    =  .false.   !     horizontal (geopotential)   (default F)   (require "key_ldfslp" when ln_sco=T)
   ln_traldf_iso    =  .true.    !     iso-neutral                 (default T)   (require "key_ldfslp")
!                                !  Coefficient
   aht0    =  2000.              !     horizontal eddy diffusivity for tracers (m2/s)
   ahtb0   =     0.              !     background eddy diffusivity for isopycnal diffusion (m2/s)
   aeiv0   =  2000.              !     eddy induced velocity coefficient (m2/s)  (require "key_traldf_eiv")
/
!-----------------------------------------------------------------------
&nam_trabbl   !   bottom boundary layer scheme                          ("key_trabbl_dif","key_trabbl_adv")
!-----------------------------------------------------------------------
!  "key_trabbl	                  ! diffusive bbl 
!  "key_trabbl_adv               ! advective bbl 
   atrbbl = 10000.               !  lateral tracer coeff. for bottom boundary layer scheme (m2/s)
/
!-----------------------------------------------------------------------
&namdmp   tracer newtonian damping ('key_tradmp')
!-----------------------------------------------------------------------
   ndmp   =   -1                 !  type of damping in temperature and salinity 
                                 !     ='latitude', damping poleward of 'ndmp' degrees and function 
                                 !                  of the distance-to-coast. Red and Med Seas as ndmp=-1
                                 !     =-1 damping only in Med and Red Seas
   ndmpf  =    1                 !  =1 create a damping.coeff NetCDF file (the 3D damping array)
   nmldmp =    1                 !  type of damping in the mixed layer
                                 !     =0 damping throughout the water column
                                 !     =1 no damping in the mixed layer defined by avt >5cm2/s )
                                 !     =2 no damping in the mixed layer defined rho<rho(surf)+.01 )
   sdmp   =  50.                 !  surface time scale for internal damping (days)
   bdmp   = 360.                 !  bottom  time scale for internal damping (days)
   hdmp   = 800.                 !  depth of transition between sdmp and bdmp (meters)
/

!!======================================================================
!!                      ***  Dynamics namelists  ***
!!======================================================================
!!   nam_dynvor    advection scheme
!!   nam_dynhpg    hydrostatic pressure gradient
!!   namflg        hydrostatic pressure gradient time stepping
!!   nam_dynldf    lateral diffusion scheme
!!======================================================================

!-----------------------------------------------------------------------
&nam_dynvor   !   option of physics/algorithm 
!-----------------------------------------------------------------------
   ln_dynvor_ene = .false.    !  enstrophy conserving scheme (default T)
   ln_dynvor_ens = .true.     !  energy conserving scheme    (default F)
   ln_dynvor_mix = .false.    !  mixed scheme                (default F)
   ln_dynvor_een = .false.    !  energy & enstrophy scheme   (default F)
/
!-----------------------------------------------------------------------
&nam_dynhpg   !   Hydrostatic pressure gradient option
!-----------------------------------------------------------------------
   ln_hpg_zco = .false.       !  z-coordinate - full steps                   
   ln_hpg_zps = .true.        !  z-coordinate - partial steps (interpolation)
   ln_hpg_sco = .false.       !  s-coordinate (standard jacobian formulation)
   ln_hpg_hel = .false.       !  s-coordinate (helsinki modification)
   ln_hpg_wdj = .false.       !  s-coordinate (weighted density jacobian)
   ln_hpg_djc = .false.       !  s-coordinate (Density Jacobian with Cubic polynomial)
   ln_hpg_rot = .false.       !  s-coordinate (ROTated axes scheme)
   gamm       = 0.e0          !  weighting coefficient (wdj scheme)
/
!-----------------------------------------------------------------------
&namflg   !   algorithm flags (algorithm not control by CPP keys)
!-----------------------------------------------------------------------
   ln_dynhpg_imp   =  .false. !  hydrostatic pressure gradient: semi-implicit time scheme  (T)
                              !                                 centered      time scheme  (F)
   nn_dynhpg_rst   =  0       !  add dynhpg implicit variables in restart ot not (1/0)
/
!-----------------------------------------------------------------------
!nam_dynspg   !   surface pressure gradient   (CPP key only)
!-----------------------------------------------------------------------
!                                !  explicit free surface               ("key_dynspg_exp")
!                                !  filtered free surface               ("key_dynspg_flt")
!                                !  split-explicit free surface         ("key_dynspg_ts")
!                                !  rigid-lid                           ("key_dynspg_rl")

!-----------------------------------------------------------------------
&nam_dynldf   !   lateral diffusion on momentum
!-----------------------------------------------------------------------
                                 !  Type of the operator : 
   ln_dynldf_lap    =  .true.    !     laplacian operator          (default T)
   ln_dynldf_bilap  =  .false.   !     bilaplacian operator        (default F)
                                 !  Direction of action  : 
   ln_dynldf_level  =  .false.   !     iso-level                   (default F)
   ln_dynldf_hor    =  .true.    !     horizontal (geopotential)   (default F)   (require "key_ldfslp" in s-coord.)
   ln_dynldf_iso    =  .false.   !     iso-neutral                 (default T)   (require "key_ldfslp")
                                 !  Coefficient
   ahm0    = 40000.              !     horizontal eddy viscosity for the dynamics (m2/s)
   ahmb0   =     0.              !     background eddy viscosity for isopycnal diffusion (m2/s)
/

!!======================================================================
!!             Tracers & Dynamics vertical physics namelists
!!======================================================================
!!       namzdf        vertical physics
!!       namnpc        non penetrative convection                       
!!       nambbl        bottom boundary layer scheme                     ("key_trabbl_dif", "key_trabbl_adv")
!!       namric        richardson number dependent vertical mixing      ("key_zdfric"      )
!!       namtke        TKE dependent vertical mixing                    ("key_zdftke"      )
!!       namkpp        KPP dependent vertical mixing                    ("key_zdfkpp"      )
!!       namddm        double diffusive mixing parameterization         ("key_zdfddm"      )
!!======================================================================

!-----------------------------------------------------------------------
&namzdf   !   vertical physics
!-----------------------------------------------------------------------
   ln_zdfnpc = .false.        !  Non-Penetrative Convection          (default T)
   avm0      = 1.2e-4         !  vertical eddy viscosity   for the dynamic (m2/s)   (background Kz if not "key_zdfcst")
   avt0      = 1.2e-5         !  vertical eddy diffusivity for tracers     (m2/s)   (background Kz if not "key_zdfcst")
   ln_zdfevd = .true.         !  enhanced vertical diffusion         (default T)
   avevd     =   100.         !  vertical coefficient for enhanced diffusion scheme (m2/s)
   n_evdm    =     0          !  enhanced mixing apply on tracer (=0) or on tracer and momentum (=1)
   ln_zdfexp =  .false.       !  =T/F  split explicit / implicit      (Default=F)
   n_zdfexp  =     3          !  number of sub-timestep for ln_zdfexp=T
/
!-----------------------------------------------------------------------
&namnpc   !   non penetrative convection
!-----------------------------------------------------------------------
   nnpc1  =      1            !  non penetrative convective scheme computation frequency
   nnpc2  =    365            !  non penetrative convective scheme print frequency
/
!-----------------------------------------------------------------------
&namric   richardson number dependent vertical diffusion                ("key_zdfric" )
!-----------------------------------------------------------------------
   avmri = 100.e-4            !  maximum value of the vertical viscosity
   alp   =      5.            !  coefficient of the parameterization
   nric  =       2            !  coefficient of the parameterization
/
!-----------------------------------------------------------------------
&namtke   !   turbulent eddy kinetic dependent vertical diffusion       ("key_zdftke" )
!-----------------------------------------------------------------------
   ln_rstke = .false.         !  flag to restart with tke from a run without tke (default F)
   ediff =       0.1          !  coef. to compute vertical eddy coef. (avt=ediff*mxl*sqrt(e) )
   ediss =       0.7          !  coef. of the Kolmogoroff dissipation
   ebb   =      3.75          !  coef. of the surface input of tke
   efave =        1.          !  coef. to applied to the tke diffusion ( avtke=efave*avm )
   emin  =     1.e-6          !  minimum value of tke (m^2/s^2)
   emin0 =     1.e-4          !  surface minimum value of tke (m^2/s^2)
   nitke =        50          !  number of restart iterative loops
!! ri_c  =         0.22222222 !  critic richardson number (default = 2/9)
   nmxl  =         2          !  mixing length    = 0 bounded by the distance to surface and bottom
                              !                   = 1 bounded by the local vertical scale factor
                              !                   = 2 first vertical derivative of mixing length bounded by 1
   npdl  =         1          !  prandtl number   = 0 no vertical prandtl number (avt=avm)
                              !                   = 1 prandtl number function of richarson number (avt=pdl*avm)
                              !                   = 2 same as = 1 but a shapiro filter is applied on pdl
   nave  =         1          !  horizontal averaged (=1) or not (=0) on avt and amv   (default =1)
   navb  =         0          !  constant (=0) or profile (=1) background on avt
/
!-----------------------------------------------------------------------
&namkpp   !   K-Profile Parameterization dependent vertical mixing      ("key_zdfkpp", and optionnally:
!                                                                        "key_kppcustom" or "key_kpplktb")
!-----------------------------------------------------------------------
   ln_kpprimix  = .true.      !  shear instability mixing  (default T)
   difmiw       =  1.0e-04    !  constant internal wave viscosity (m2/s)
   difsiw       =  0.1e-04    !  constant internal wave diffusivity (m2/s)
   Riinfty      =  0.8        !  local Richardson Number limit for shear instability
   difri        =  0.0050     !  maximum shear mixing at Rig = 0    (m2/s)
   bvsqcon      = -0.01e-07   !  Brunt-Vaisala squared (1/s**2) for maximum convection
   difcon       =  1.         !  maximum mixing in interior convection (m2/s)
   navb         =  0          !  horizontal averaged (=1) or not (=0) on avt and amv   (default =1)
   nave         =  1          !  constant (=0) or profile (=1) background on avt
/
!-----------------------------------------------------------------------
&namddm   !   double diffusive mixing parameterization                  ("key_zdfddm")
!-----------------------------------------------------------------------
      avts  = 1.e-4           !  maximum avs (vertical mixing on salinity)
      hsbfr = 1.6             !  heat/salt buoyancy flux ratio
/

!!======================================================================
!!                  ***  Miscelaneous namelists  ***
!!======================================================================
!!   namsol      elliptic solver / island / free surface 
!!======================================================================

!-----------------------------------------------------------------------
&namsol   !   elliptic solver / island / free surface 
!-----------------------------------------------------------------------
   nsolv     =      1         !  type of elliptic solver  =1 preconditioned conjugate gradient (pcg)
                              !                           =2 successive-over-relaxation (sor)
                              !                           =3 FETI (fet)                               ("key_feti")
                              !                           =4 sor with extra outer halo
   nsol_arp  =      0         !  absolute/relative (0/1) precision convergence test
   nmin      =    300         !  minimum of iterations for the SOR solver
   nmax      =    800         !  maximum of iterations for the SOR solver
   nmod      =     10         !  frequency of test for the SOR solver
   eps       =  1.E-6         !  absolute precision of the solver
   resmax    = 1.E-10         !  absolute precision for the SOR solver
   sor       =   1.92         !  optimal coefficient for SOR solver
   epsisl    = 1.e-10         !  absolute precision on stream function solver
   nmisl     =   4000         !  maximum pcg iterations for island                                    ("key_islands")
   rnu       =     1.         !  strength of the additional force used in filtered free surface
/

!!======================================================================
!!                       ***  Diagnostics namelists  ***
!!======================================================================
!!   namtrd       dynamics and/or tracer trends                         ("key_trddyn","key_trdtra","key_trdmld")
!!   namgap       level mean model-data gap                             ("key_diagap")
!!   namptr       Poleward Transport Diagnostics                        ("key_diaspr")
!!   namspr       surface pressure diagnosed in rigid-lid               ("key_diaspr")
!!   namflo       float parameters                                      ("key_float")
!!======================================================================

!-----------------------------------------------------------------------
&namtrd    !   diagnostics on dynamics and/or tracer trends             ("key_trddyn" and/or "key_trdtra")
!                          or mixed-layer trends                        ('key_trdmld')
!                          or barotropic vorticity                      ("key_trdvor")
!-----------------------------------------------------------------------
   ntrd  = 365                   !  time step frequency dynamics and tracers trends
   nctls =   0                   !  control surface type in mixed-layer trends (0,1 or n<jpk)
   ln_trdmld_restart = .false.   !  restart for ML diagnostics
   ucf   =  1.                   !  unit conversion factor (=1 -> /seconds | =86400. -> /day)
   ln_trdmld_instant = .false.   !  flag to diagnose trends of instantantaneous or mean ML T/S
/
!-----------------------------------------------------------------------
&namgap    !   level mean model-data gap                                ('key_diagap')
!-----------------------------------------------------------------------
   ngap =  15                    !  time-step frequency of model-data gap computation
   nprg =  10                    !  time-step frequency of gap print in model output
/
!-----------------------------------------------------------------------
&namptr   !   Poleward Transport Diagnostic
!-----------------------------------------------------------------------
   ln_diaptr = .false.           !  logical flag for Poleward transport computation
   ln_subbas = .false.           !  logical flag for Atlantic/Pacific/Indian basins computation
                                 !  need input basins mask file named "subbasins.nc"
   nf_ptr    =  15               !  Frequency of computation
/
!-----------------------------------------------------------------------
&namspr  !   surface pressure diagnostic
!-----------------------------------------------------------------------
   nmaxp   =   1000              !  maximum of iterations for the solver
   epsp    =  1.e-3              !  absolute precision of the solver
   niterp  =    400              !  number of iteration done by the solver
/
!-----------------------------------------------------------------------
&namflo    !   float parameters                                         ("key_float")
!-----------------------------------------------------------------------
    ln_rstflo = .false.          !  boolean term for float restart (true or false)
    nwritefl  =      75          !  frequency of float output file 
    nstockfl  =    5475          !  frequency of float restart file 
    ln_argo   = .false.          !  Argo type floats (stay at the surface each 10 days)
    ln_flork4 = .false.          !  = T trajectories computed with a 4th order Runge-Kutta
                                 !  = F  (default)   computed with Blanke' scheme
/