!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> !! NEMO/TOP1 : 1 - tracer definition (namtrc ) !! 2 - tracer data initialisation (namtrc_dta) !! 3 - tracer advection (namtrc_adv) !! 4 - tracer lateral diffusion (namtrc_ldf) !! 5 - tracer vertical physics (namtrc_zdf) !! 6 - tracer newtonian damping (namtrc_dmp) !! 7 - dynamical tracer trends (namtrc_trd) !! 8 - tracer output diagonstics (namtrc_dia) !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' &namtrc ! tracers definition !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, nn_dttrc = 1 ! time step frequency for passive sn_tracers ln_top_euler = .true. ! Euler timestepping ! ! ! name ! title of the field ! initial data ! initial data ! save ! ! ! ! ! units ! from file ! or not ! ! ! ! ! ! or not ! ! sn_tracer(1) = 'DUMMY' , 'Dummy tracer' , 'dummy units', .false. , .true. / !----------------------------------------------------------------------- &namtrc_dta ! Initialisation from data input file !----------------------------------------------------------------------- ! ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! sn_trcdta(1) = 'data_1m_O2_nomask' , -1 , 'O2' , .false. , .true. , 'yearly' , '' , '' sn_trcdta(2) = 'data_1m_PO4_nomask' , -1 , 'PO4' , .false. , .true. , 'yearly' , '' , '' sn_trcdta(3) = 'data_1m_NO3_nomask' , -1 , 'NO3' , .false. , .true. , 'yearly' , '' , '' sn_trcdta(6) = 'data_1m_Si_nomask' , -1 , 'Si' , .false. , .true. , 'yearly' , '' , '' sn_trcdta(8) = 'data_1m_Fer_nomask' , -12 , 'Fer' , .false. , .true. , 'yearly' , '' , '' sn_trcdta(45) = 'data_1m_DOC_nomask' , -12 , 'DOC' , .false. , .true. , 'yearly' , '' , '' sn_trcdta(56) = 'data_1m_DIC_nomask' , -12 , 'DIC' , .false. , .true. , 'yearly' , '' , '' sn_trcdta(57) = 'data_1m_Alkalini_nomask', -12 , 'Alkalini', .false. , .true. , 'yearly' , '' , '' ! cn_dir = './' ! root directory for the location of the data files rn_trfac(1) = 22.4 ! conversion factor from ml/l to mmol m3 (1 mol O2 = 22.4 l) rn_trfac(2) = 1.0 ! - - - - rn_trfac(3) = 1.0 ! - - - - rn_trfac(6) = 1.0 ! - - - - rn_trfac(8) = 1.e9 ! conversion factor from mol/l to umol/m3 rn_trfac(45) = 12.e6 ! conversion factor from mol/l to mg/m3 rn_trfac(56) = 12.0 ! conversion factor from umol/l to mg/m3 rn_trfac(57) = 1.0 ! - - - - / !----------------------------------------------------------------------- &namtrc_bc !----------------------------------------------------------------------- ! ! file name ! frequency (hr) ! variable ! time interp. ! clim !'yearly' or ! weights ! rotation ! ! ! (if <0 mon) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing sn_trcsbc(3) = 'ndeposition.orca' , -1 , 'ndep' , .false. , .true. , 'yearly' , '' , '' sn_trcsbc(8) = 'dust.orca' , -1 , 'dust' , .false. , .true. , 'yearly' , '' , '' rn_trsfac(3) = 1.0 rn_trsfac(8) = 5.8866e+05 ! conversion factor from kg Fe/m2/mon to umol/m2/d = 1.e9/(55.85)*12/365 ! sn_trccbc(2) = 'river' , -12 , 'RIVERPHOS' , .false. , .true. , 'yearly' , '' , '' sn_trccbc(3) = 'river' , -12 , 'RIVERNIT' , .false. , .true. , 'yearly' , '' , '' sn_trccbc(6) = 'river' , -12 , 'RIVERSIL' , .false. , .true. , 'yearly' , '' , '' sn_trccbc(8) = 'river' , -12 , 'RIVERFe80' , .false. , .true. , 'yearly' , '' , '' sn_trccbc(45) = 'river' , -12 , 'RIVERDOC' , .false. , .true. , 'yearly' , '' , '' sn_trccbc(56) = 'river' , -12 , 'RIVERDIC' , .false. , .true. , 'yearly' , '' , '' rn_trcfac(2) = 8.8464e+10 ! 1 mol P = 30.97 g; 1.e12*1.e3/30.97/365 conversion factor from Tg/y to mmol/d rn_trcfac(3) = 1.9556e+11 ! 1 mol N = 14.01 g; 1.e12*1.e3/14.01/365 conversion factor from Tg/y to mmol/d rn_trcfac(6) = 9.7534e+10 ! 1 mol Si = 28.09 g; 1.e12*1.e3/28.09/365 conversion factor from Tg/y to mmol/d rn_trcfac(8) = 4.9055e+13 ! 1 mol Fe = 55.85 g; 1.e12*1.e6/55.85/365 conversion factor from Tg/y to umol/d rn_trcfac(45) = 2.7397e+12 ! 1.e15/365 ! conversion factor from Tg/y to mg/d rn_trcfac(56) = 2.7397e+12 ! 1.e15/365 ! conversion factor from Tg/y to mg/d cn_dir = './' ! root directory for the location of the boundary condition files / !----------------------------------------------------------------------- &namtrc_adv ! advection scheme for passive tracer !----------------------------------------------------------------------- ln_trcadv_cen2 = .false. ! 2nd order centered scheme ln_trcadv_tvd = .false. ! TVD scheme ln_trcadv_muscl = .true. ! MUSCL scheme ln_trcadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries ln_trcadv_ubs = .false. ! UBS scheme ln_trcadv_qck = .false. ! QUICKEST scheme / !----------------------------------------------------------------------- &namtrc_ldf ! 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 rn_ahtrc_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] rn_ahtrb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] / !----------------------------------------------------------------------- &namtrc_zdf ! vertical physics !----------------------------------------------------------------------- ln_trczdf_exp = .false. ! split explicit (T) or implicit (F) time stepping nn_trczdf_exp = 3 ! number of sub-timestep for ln_trczdfexp=T / !----------------------------------------------------------------------- &namtrc_rad ! treatment of negative concentrations !----------------------------------------------------------------------- ln_trcrad = .true. ! artificially correct negative concentrations (T) or not (F) / !----------------------------------------------------------------------- &namtrc_dmp ! passive tracer newtonian damping !----------------------------------------------------------------------- ln_trcdmp = .false. ! add a damping term (T) or not (F) nn_hdmp_tr = -1 ! horizontal shape =-1, damping in Med and Red Seas only ! =XX, damping poleward of XX degrees (XX>0) ! + F(distance-to-coast) + Red and Med Seas nn_zdmp_tr = 1 ! vertical shape =0 damping throughout the water column ! =1 no damping in the mixing layer (kz criteria) ! =2 no damping in the mixed layer (rho crieria) rn_surf_tr = 50. ! surface time scale of damping [days] rn_bot_tr = 360. ! bottom time scale of damping [days] rn_dep_tr = 800. ! depth of transition between rn_surf and rn_bot [meters] nn_file_tr = 0 ! create a damping.coeff NetCDF file (=1) or not (=0) / !----------------------------------------------------------------------- &namtrc_trd ! diagnostics on tracer trends ('key_trdtrc') ! or mixed-layer trends ('key_trdmld_trc') !---------------------------------------------------------------------- nn_trd_trc = 360 ! time step frequency and tracers trends nn_ctls_trc = 0 ! control surface type in mixed-layer trends (0,1 or n /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 ln_trdtrc(1) = .true. ln_trdtrc(2) = .true. ln_trdtrc(23) = .true. / !----------------------------------------------------------------------- &namtrc_dia ! parameters for passive tracer additional diagnostics !---------------------------------------------------------------------- ln_diatrc = .false. ! save additional diag. (T) or not (F) ln_diabio = .false. nn_writedia = 5475 ! time step frequency for diagnostics nn_writebio = 5475 /