!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! NEMO/OPA  :  1 - run manager      (namrun)
!! namelists    2 - miscellaneous    (namctl,nammpp)
!!              3 - Domain           (namzgr, namzgr_sco, namdom)
!!              6 - Tracer           (nameos, namcla, namqsr)
!!              7 - Inputs dynamics  (namdyna)
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!  CAUTION: some scripts does not support CAPITALs for logical use .true./.false., not .TRUE./.FALSE.

!!======================================================================
!!                   ***  Run management namelists  ***
!!======================================================================
!!   namrun            parameters of the run
!!======================================================================

!-----------------------------------------------------------------------
&namrun        !   parameters of the run
!-----------------------------------------------------------------------
   no          =       0   !  job number
   cexper      =  "PISCES"  !  experience name 
   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T)
   nrstdt      =       0   !  restart control = 0 nit000 is not compared to the restart file value
                           !                  = 1 use ndate0 in namelist (not the value in the restart file)
                           !                  = 2 calendar parameters read in the restart file
   nit000      =       1   !  first time step
   nitend      =     1460  !  last  time step
   ndate0      =  010101   !  initial calendar date yymmdd (used if nrstdt=1)
   nleapy      =       0   !  Leap year calendar (1) or not (0)
   ninist      =       0   !  output the initial state (1) or not (0)
   nstock      =     1460  !  frequency of creation of a restart file (modulo referenced to 1)
   nwrite      =     1460  !  frequency of write in the output file   (modulo referenced to nit000)
   ln_mskland  = .false.   !  mask land points in NetCDF outputs (costly: + ~15%)
/
!-----------------------------------------------------------------------
&namctl       !   Control prints & Benchmark
!-----------------------------------------------------------------------
   ln_ctl     = .false.    !  trends control print (expensive!)
   nprint     =    0       !  level of print (0 no extra print)
   nictls     =    1       !  start i indice of control sum (use to compare mono versus
   nictle     =    182       !  end   i indice of control sum        multi processor runs
   njctls     =    1       !  start j indice of control               over a subdomain)
   njctle     =    149       !  end   j indice of control 
   isplt      =    1       !  number of processors in i-direction
   jsplt      =    1       !  number of processors in j-direction
   nbench     =    0       !  Bench mode (1/0): CAUTION use zero except for bench
/
!-----------------------------------------------------------------------
&nammpp      !   Massively Parallel Processing                         ("key_mpp_mpi)
!-----------------------------------------------------------------------
   c_mpi_send =  'S'       !  mpi send/recieve type   ='S', 'B', or 'I' for standard send,
                           !  buffer blocking send or immediate non-blocking sends, resp.
/
!-----------------------------------------------------------------------
&namzgr       !   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)
/
!-----------------------------------------------------------------------
&namzgr_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)
!-----------------------------------------------------------------------
   e3zps_min   =    25.    !  the thickness of the partial step is set larger than the minimum
   e3zps_rat   =    0.2    !  of e3zps_min and e3zps_rat * e3t   (N.B. 0<e3zps_rat<1)
   nmsh        =    1      !  create (=1) a mesh file (coordinates, scale factors, masks) or not (=0)
   nacc        =    0      !  =1 acceleration of convergence method used, rdt < rdttra(k)
                           !  =0, no acceleration, rdt = rdttra
   atfp        =    0.1    !  asselin time filter parameter
   rdt         = 21600.    !  time step for the dynamics (and tracer if nacc=0)
   rdtmin      = 21600.    !  minimum time step on tracers (used if nacc=1)
   rdtmax      = 21600.    !  maximum time step on tracers (used if nacc=1)
   rdth        =  800.     !  depth variation of tracer time step  (used if nacc=1)
/
!-----------------------------------------------------------------------
&namtraldf    !   lateral diffusion scheme for tracer
!-----------------------------------------------------------------------
!                               !  Type of the operator :
   ln_traldf_lap    =  .true.   !     laplacian operator
   ln_traldf_bilap  =  .false.  !     bilaplacian operator
                                !  Direction of action  :
   ln_traldf_level  =  .false.  !     iso-level
   ln_traldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T)
   ln_traldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp")
   ln_traldf_grif   =  .false.  !     griffies skew flux formulation    (require "key_ldfslp")
!                               !  Coefficient
   aht0        =  2000.         !     horizontal eddy diffusivity for tracers [m2/s]
   ahtb0       =     0.         !     background eddy diffusivity for ldf_iso [m2/s]
   aeiv0       =  2000.         !     eddy induced velocity coefficient [m2/s]    (require "key_traldf_eiv")
/
!-----------------------------------------------------------------------
&namcla        !   cross land advection
!-----------------------------------------------------------------------
   n_cla       =    0      !  advection between 2 ocean pts separates by land
/
!-----------------------------------------------------------------------
&namqsr        !   penetrative solar radiation
!-----------------------------------------------------------------------
   ln_qsr_bio  = .false.   !  bio-model light penetration
   rn_abs      =   0.58    !  RGB & 2 bands: fraction of light (rn_si1)
   rn_si0      =   0.35    !  RGB & 2 bands: shortess depth of extinction
   rn_si2      =   62.0    !  3 bands: longest depth of extinction (for blue waveband & 0.01 mg/m2 Chl)
/  
!-----------------------------------------------------------------------
&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)
/
!-----------------------------------------------------------------------
&namdyn        !   offline parameters 
!-----------------------------------------------------------------------
    ndtadyn   =  73        ! number of period in the file for one year    
    ndtatot   =  73        ! total number of period in the file    
    nsptint   =  1         ! indicator for time interpolation    
    nficdyn   =  2         ! number of file to read    
    lperdyn   = .true.     ! periodicity of the unique file (T)
!                          ! F  (default)   computed with Blanke' scheme 
    cfile_grid_T = 'NEMOV3_5d_21210101_21211231_grid_T.nc' ! name of grid_T file
    cfile_grid_U = 'NEMOV3_5d_21210101_21211231_grid_U.nc' ! name of grid_U file
    cfile_grid_V = 'NEMOV3_5d_21210101_21211231_grid_V.nc' ! name of grid_V file
    cfile_grid_W = 'NEMOV3_5d_21210101_21211231_grid_W.nc' ! name of grid_W file
/