source: trunk/CONFIG/GYRE_LOBSTER/EXP00/namelist_top @ 1090

!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
! OPA MODEL general namelist for passive tracers
! -------------
!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!
!-----------------------------------------------------------------------
!       nattrc  
!-----------------------------------------------------------------------
!
!       NATTRC   
!       nwritetrc time step frequency for tracer outputs
!       lrsttr    boolean term for tracer model restart (true or false)
!       nrsttr    control of the time step for tracer model restart (0, 1 or 2)
!       tracer type defined by :
!                 *  short name 
!                 *  long_name
!                 *  units
!                 *  logical to read initial value from file or not
!                 *  logical to save value  
! 
&nattrc
   nwritetrc   = 360 
   lrsttr      = .false.
   nrsttr      = 0
   tracer(1)   = 'DET' , 'Detritus                   ',  'mmole-N/m3' ,  .false. ,  .true.
   tracer(2)   = 'ZOO' , 'Zooplankton concentration  ',  'mmole-N/m3' ,  .false. ,  .true.
   tracer(3)   = 'PHY' , 'Phytoplankton concentration',  'mmole-N/m3' ,  .false. ,  .true.
   tracer(4)   = 'NO3' , 'Nitrate concentration      ',  'mmole-N/m3' ,  .false. ,  .true.
   tracer(5)   = 'NH4' , 'Ammonium concentration     ',  'mmole-N/m3' ,  .false. ,  .true.
   tracer(6)   = 'DOM' , 'Dissolved organic matter   ',  'mmole-N/m3' ,  .false. ,  .true.
/
!-----------------------------------------------------------------------
!       natrtd   dynamical tracers trends (#ifdef key_trc_diatrd) 
!-----------------------------------------------------------------------
!   nwritetrd : time step frequency for dynamical trends outputs
!   luttrd    : logical to keep large diagnostics with trends or not
!               one value per tracer
&natrtd
   nwritetrd  = 360
   luttrd(1)  = .false.
   luttrd(2)  = .false.
   luttrd(3)  = .false.
   luttrd(4)  = .false.
   luttrd(5)  = .false.
   luttrd(6)  = .false.
/
!-----------------------------------------------------------------------
!       natdia  additional 2D/3D  (#ifdef key_trc_diaadd)
!-----------------------------------------------------------------------
!  nwritedia : time step frequency for additional arrays outputs
!  2D/3D diagnostic type defined by :
!                 *  short name 
!                 *  long_name
!                 *  units
! 
&natdia
   nwritedia   = 360
   diag2d(1)   = 'TNO3PHY' , 'TNO3PHY',  '-'
   diag2d(2)   = 'TNH4PHY' , 'TNH4PHY',  '-'
   diag2d(3)   = 'TPHYDOM' , 'TPHYDOM',  '-'
   diag2d(4)   = 'TPHYNH4' , 'TPHYNH4',  '-'
   diag2d(5)   = 'TPHYZOO' , 'TPHYZOO',  '-'
   diag2d(6)   = 'TPHYDET' , 'TPHYDET',  '-'
   diag2d(7)   = 'TDETZOO' , 'TDETZOO',  '-'
   diag2d(8)   = 'TDETSED' , 'TDETSED',  '-'
   diag2d(9)   = 'TZOODET' , 'TZOODET',  '-'
   diag2d(10)  = 'TZOOBOD' , 'TZOOBOD',  '-'
   diag2d(11)  = 'TZOONH4' , 'TZOONH4',  '-'
   diag2d(12)  = 'TZOODOM' , 'TZOODOM',  '-'
   diag2d(13)  = 'TNH4NO3' , 'TNH4NO3',  '-'
   diag2d(14)  = 'TDOMNH4' , 'TDOMNH4',  '-'
   diag2d(15)  = 'TDETNH4' , 'TDETNH4',  '-'
   diag2d(16)  = 'TPHYTOT' , 'TPHYTOT',  '-'
   diag2d(17)  = 'TZOOTOT' , 'TZOOTOT',  '-'
   diag2d(18)  = 'TDETDOM' , 'TDETDOM',  '-'
   diag2d(19)  = 'SEDPOC ' , 'SEDPOC ',  '-'
   diag3d(1)   = 'FNO3PHY' , 'FNO3PHY',  '-'
   diag3d(2)   = 'FNH4PHY' , 'FNH4PHY',  '-'
   diag3d(3)   = 'FNH4NO3' , 'FNH4NO3',  '-'
/
!-----------------------------------------------------------------------
!       natnum  numerical schemes 
!-----------------------------------------------------------------------
!  ndttrc time step frequency for passive tracers
!  lhdf logical if true computes horizontal diffusion
!  rsc tuning coefficient for Smolar. 
!  rtrn truncation value for Smolar.
!  ncor number of corrective phases for Smolar. 
!  crosster logical if true computes Smolar crossterms
&natnum
   ndttrc   = 1
   rsc      = 1.
   rtrn     = 1.e-15
   ncortrc  = 1
   crosster = .false.
/
!-----------------------------------------------------------------------
!       namtrcadv   advection scheme for tracer (option not control by CPP keys)
!-----------------------------------------------------------------------
!  ln_trcadv_cen2     2nd order centered scheme    (default F)
!  ln_trcadv_tvd      TVD scheme                   (default F)
!  ln_trcadv_muscl    MUSCL scheme                 (default F)
!  ln_trcadv_muscl2   MUSCL2 scheme                (default F)
!  ln_trcadv_smolar   SMOLAR scheme                (default T)
&namtrcadv
   ln_trcadv_cen2   =  .false.
   ln_trcadv_tvd    =  .true.
   ln_trcadv_muscl  =  .false.
   ln_trcadv_muscl2 =  .false.
   ln_trcadv_smolar =  .false.
/
!
!-----------------------------------------------------------------------
!       namtrcbbl   bottom boundary layer scheme
!-----------------------------------------------------------------------
!  atrcbbl   lateral tracer coeff. for bottom boundary layer scheme(m2/s)
&namtrcbbl
   atrcbbl = 1000.
/
!-----------------------------------------------------------------------
!       namtrcldf   lateral diffusion scheme for tracer (option not control by CPP keys)
!-----------------------------------------------------------------------
!  Flag to performs lateral diffusion or not :
!     ln_trcldf_diff
!  Type of the operator : 
!     ln_trcldf_lap    laplacian operator          (default T)
!     ln_trcldf_bilap  bilaplacian operator        (default F)
!  Direction of action  : 
!     ln_trcldf_level  iso-level                   (default F)
!     ln_trcldf_hor    horizontal (geopotential)   (default F)^**
!     ln_trcldf_iso    iso-neutral                 (default T)^*
! ^* require key_ldfslp to compute the direction of the lateral diffusion
! ^** require key_ldfslp in s-coordinate
!  ahtrb0     background eddy diffusivity for isopycnal diffusion (m2/s)
!  trcrat     ratio betweeen passive and active tracer diffusion coeff
!  ahtrc0     horizontal eddy diffus. for passive tracer
!  aeivtr0    eddy induced veloc. coef. for passive tracer
&namtrcldf
   ln_trcldf_diff   =  .true.
   ln_trcldf_lap    =  .true.
   ln_trcldf_bilap  =  .false.
   ln_trcldf_level  =  .false.
   ln_trcldf_hor    =  .false.
   ln_trcldf_iso    =  .true.
   ahtrb0           =  0. 
   trcrat           =  1.
   ahtrc0           =  300.
   aeivtr0          =  1000.
/
!-----------------------------------------------------------------------
!       namtrczdf   vertical physics
!-----------------------------------------------------------------------
!  ln_zdfexp   vertical physics: (=T)  time splitting (T)     (Default=F)
!                               (=F)  euler backward (F)
!  n_zdfexp   number of sub-timestep for time splitting scheme
&namtrczdf
   ln_trczdf_exp =  .false.
   n_trczdf_exp =      3
/
!-----------------------------------------------------------------------
!       namtrcdmp   tracer newtonian damping ('key_trcdmp')
!-----------------------------------------------------------------------
!  ndmptr    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)
!  ndmpftr   =1 create a damping.coeff NetCDF file (the 3D damping array)
!  nmldmptr  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 )
!  sdmptr    surface time scale for internal damping (days)
!  bdmptr    bottom time scale for internal damping (days)
!  hdmptr    depth of transition between sdmp and bdmp (meters)
&namtrcdmp
   ndmptr   =  20
   ndmpftr  =   0
   nmldmptr =   0
   sdmptr   =  50.
   bdmptr   = 360.
   hdmptr   = 800.
/