!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! NEMO/TOP1 :  1 - tracer definition                     (namtoptrc)
!! namelists    2 - dynamical tracer trends               (namtoptrd)
!!              6 - tracer advection                      (namtopadv)
!!              7 - tracer bottom boundary                (namtopbbl)
!!              8 - tracer lateral diffusion              (namtopldf)
!!              3 - tracer vertical physics               (namtopzdf)
!!              9 - tracer newtonian damping              (namtopdmp)
!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!'''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''
&namtoptrc     !   tracers definition
!,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,
   ndttrc      =  1        !  time step frequency for passive tracers      
   nwritetrc   = 360       !  time step frequency for tracer outputs
   lrsttr      = .false.   !  start from a restart file (T) or not (F)
   nrsttr      =   0       !  restart control = 0 initial time step is not compared to the restart file value
                           !                  = 1 do not use the value in the restart file
                           !                  = 2 calendar parameters read in the restart file
   cn_trcrst_in  = "restart_trc"   !  suffix of pass. tracer restart name (input)
   cn_trcrst_out = "restart_trc"   !  suffix of pass. tracer restart name (output)
!
!              ! name  !     title of the field       !   units       ! initial data ! save   !
!              !       !                              !               ! from file    ! or not ! 
!              !       !                              !               ! or not       !        !
   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.
/
!-----------------------------------------------------------------------
&namtopadv    !   advection scheme for passive tracer 
!-----------------------------------------------------------------------
   ln_trcadv_cen2   =  .false.  !  2nd order centered scheme   
   ln_trcadv_tvd    =  .true.  !  TVD scheme   
   ln_trcadv_muscl  =  .false.   !  MUSCL scheme
   ln_trcadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries
   ln_trcadv_smolar =  .false.  !  SMOLAR scheme
   rsc              =  1.       !  tuning coefficient for Smolar. scheme 
   ncortrc          =  1        !  number of corrective phases for Smolar. scheme 
   crosster         =  .false.  !  computes Smolar crossterms (T) or not (F)
/
!-----------------------------------------------------------------------
&namtopbbl    !   bottom boundary layer scheme for passive tracer 
!-----------------------------------------------------------------------
   atrcbbl          = 1000.     !  lateral mixing coefficient in the bbl  [m2/s]
/
!-----------------------------------------------------------------------
&namtopldf    !   lateral diffusion scheme for passive tracer 
!-----------------------------------------------------------------------
   ln_trcldf_diff   =  .true.   !  performs lateral diffusion (T) or not (F)
!                               !  Type of the operator : 
   ln_trcldf_lap    =  .true.   !     laplacian operator       
   ln_trcldf_bilap  =  .false.  !     bilaplacian operator     
                                !  Direction of action  :
   ln_trcldf_level  =  .false.  !     iso-level                
   ln_trcldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T)
   ln_trcldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp")
!                               !  Coefficient
   ahtrc0      =   300.         !     horizontal eddy diffusivity for tracers [m2/s]
   ahtrb0      =     0.         !     background eddy diffusivity for ldf_iso [m2/s]
   aeivtr0     =  1000.         !     eddy induced velocity coefficient [m2/s]    (require "key_trcldf_eiv")
   trcrat      =  1.            !  ratio betweeen passive and active tracer diffusion coeff
/
!-----------------------------------------------------------------------
&namtopzdf        !   vertical physics
!-----------------------------------------------------------------------
   ln_trczdf_exp   =  .false.  !  split explicit (T) or implicit (F) time stepping
   n_trczdf_exp    =   3       !  number of sub-timestep for ln_trczdfexp=T
/
!-----------------------------------------------------------------------
&namtoprad        !  treatment of negative concentrations 
!-----------------------------------------------------------------------
   ln_trcrad   =  .true.  !  artificially correct negative concentrations (T) or not (F)
/
!-----------------------------------------------------------------------
&namtopdmp        !   passive tracer newtonian damping                 ('key_trcdmp')
!-----------------------------------------------------------------------
   ndmptr      =   20      !  type of damping in passive tracers
                           !     ='latitude', damping poleward of 'ndmp' degrees and function 
                           !                  of the distance-to-coast. Red and Med Seas as ndmptr=-1
                           !     =-1 damping only in Med and Red Seas
   ndmpftr     =    0      !  create a damping.coeff NetCDF file (=1) or not (=0)
   nmldmptr    =    1      !  type of damping: =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      =   50.     !  surface time scale for internal damping (days)
   bdmptr      =  360.     !  bottom  time scale for internal damping (days)
   hdmptr      =  800.     !  depth of transition between sdmptr and bdmptr (meters)
/
!-----------------------------------------------------------------------
&namtoptrd       !   diagnostics on tracer trends
!                          or mixed-layer trends          ('key_trdmld_trc')
!----------------------------------------------------------------------
   ntrd_trc   =  360      !  time step frequency and tracers trends
   nctls_trc  =   0       !  control surface type in mixed-layer trends (0,1 or n<jpk)
   ucf_trc    =  86400    !  unit conversion factor (=1 -> /seconds ; =86400. -> /day)
   ln_trdmld_trc_restart = .false.  !  restart for ML diagnostics
   ln_trdmld_trc_instant = .false.  !  flag to diagnose trends of instantantaneous or mean ML T/S
   luttrd(1)  =   .true.
   luttrd(2)  =   .true.
   luttrd(3)  =   .true.
   luttrd(4)  =   .true.
   luttrd(5)  =   .true.
   luttrd(6)  =   .true.
/