!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ! 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) ! ctrcnm tracer name ! ctrcun tracer unit ! lutini logical to read initial value from file or not ! &nattrc nwritetrc = 360, lrsttr = .false. nrsttr = 0, ! ijulian = 010101, ctrcnm(1) = 'DET', ctrcnl(1) = 'Detritus', ctrcun(1) = 'mmole-N/m3', lutini(1) = .true., ctrcnm(2) = 'ZOO', ctrcnl(2) = 'Zooplankton concentration', ctrcun(2) = 'mmole-N/m3', lutini(2) = .true., ctrcnm(3) = 'PHY', ctrcnl(3) = 'Phytoplankton concentration', ctrcun(3) = 'mmole-N/m3', lutini(3) = .true., ctrcnm(4) = 'NO3', ctrcnl(4) = 'Nitrate concentration', ctrcun(4) = 'mmole-N/m3', lutini(4) = .true., ctrcnm(5) = 'NH4', ctrcnl(5) = 'Ammonium concentration', ctrcun(5) = 'mmole-N/m3', lutini(5) = .true., ctrcnm(6) = 'DOM', ctrcnl(6) = 'Dissolved organic matter', ctrcun(6) = 'mmole-N/m3', lutini(6) = .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., / !----------------------------------------------------------------------- ! natadd 2/3 D diagnostics (#ifdef key_trc_diaadd) !----------------------------------------------------------------------- ! nwriteadd : time step frequency for additional arrays outputs ! ctrc3d short title (max 8 characters) for 3d output arrays ! ctrc3l long title (max 80 characters) for 3d output arrays ! ctrc3u units ! ctrc2d short title (max 8 characters) for 2d output arrays ! ctrc2l long title (max 80 characters) for 2d output arrays ! ctrc2u units ! &natadd nwriteadd = 360, ctrc3d(1) = 'FNO3PHY', ctrc3l(1) = 'FNO3PHY', ctrc3u(1) = '-', ctrc3d(2) = 'FNH4PHY', ctrc3l(2) = 'FNH4PHY', ctrc3u(2) = '-', ctrc3d(3) = 'FNH4NO3', ctrc3l(3) = 'FNH4NO3', ctrc3u(3) = '-', ctrc2d(1) = 'TNO3PHY' ctrc2l(1) = 'TNO3PHY' ctrc2u(1) = '-' ctrc2d(2) = 'TNH4PHY' ctrc2l(2) = 'TNH4PHY' ctrc2u(2) = '-' ctrc2d(3) = 'TPHYDOM' ctrc2l(3) = 'TPHYDOM' ctrc2u(3) = '-' ctrc2d(4) = 'TPHYNH4' ctrc2l(4) = 'TPHYNH4' ctrc2u(4) = '-' ctrc2d(5) = 'TPHYZOO' ctrc2l(5) = 'TPHYZOO' ctrc2u(5) = '-' ctrc2d(6) = 'TPHYDET' ctrc2l(6) = 'TPHYDET' ctrc2u(6) = '-' ctrc2d(7) = 'TDETZOO' ctrc2l(7) = 'TDETZOO' ctrc2u(7) = '-' ctrc2d(8) = 'TDETSED' ctrc2l(8) = 'TDETSED' ctrc2u(8) = '-' ctrc2d(9) = 'TZOODET' ctrc2l(9) = 'TZOODET' ctrc2u(9) = '-' ctrc2d(10) = 'TZOOBOD' ctrc2l(10) = 'TZOOBOD' ctrc2u(10) = '-' ctrc2d(11) = 'TZOONH4' ctrc2l(11) = 'TZOONH4' ctrc2u(11) = '-' ctrc2d(12) = 'TZOODOM' ctrc2l(12) = 'TZOODOM' ctrc2u(12) = '-' ctrc2d(13) = 'TNH4NO3' ctrc2l(13) = 'TNH4NO3' ctrc2u(13) = '-' ctrc2d(14) = 'TDOMNH4' ctrc2l(14) = 'TDOMNH4' ctrc2u(14) = '-' ctrc2d(15) = 'TDETNH4' ctrc2l(15) = 'TDETNH4' ctrc2u(15) = '-' ctrc2d(16) = 'TPHYTOT' ctrc2l(16) = 'TPHYTOT' ctrc2u(16) = '-' ctrc2d(17) = 'TZOOTOT' ctrc2l(17) = 'TZOOTOT' ctrc2u(17) = '-' ctrc2d(18) = 'TDETDOM' ctrc2l(18) = 'TDETDOM' ctrc2u(18) = '-' ctrc2d(19) = 'SEDPOC' ctrc2l(19) = 'SEDPOC' ctrc2u(19) = '-' / !----------------------------------------------------------------------- ! 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 = 0, 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. / !----------------------------------------------------------------------- ! namtrcldf lateral diffusion scheme for tracer (option not control by CPP keys) !----------------------------------------------------------------------- ! 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