1 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
---|
2 | !! NEMO/TOP1 : 1 - tracer definition (namtrc ) |
---|
3 | !! 2 - tracer data initialisation (namtrc_dta) |
---|
4 | !! 3 - tracer advection (namtrc_adv) |
---|
5 | !! 4 - tracer lateral diffusion (namtrc_ldf) |
---|
6 | !! 5 - tracer vertical physics (namtrc_zdf) |
---|
7 | !! 6 - tracer newtonian damping (namtrc_dmp) |
---|
8 | !! 7 - dynamical tracer trends (namtrc_trd) |
---|
9 | !! 8 - tracer output diagonstics (namtrc_dia) |
---|
10 | !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
---|
11 | !''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''' |
---|
12 | &namtrc ! tracers definition |
---|
13 | !,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,, |
---|
14 | nn_dttrc = 1 ! time step frequency for passive sn_tracers |
---|
15 | ln_top_euler = .true. ! Euler timestepping |
---|
16 | ! |
---|
17 | ! ! name ! title of the field ! initial data ! initial data ! save ! |
---|
18 | ! ! ! ! units ! from file ! or not ! |
---|
19 | ! ! ! ! ! or not ! ! |
---|
20 | sn_tracer(1) = 'DUMMY' , 'Dummy tracer' , 'dummy units', .false. , .true. |
---|
21 | / |
---|
22 | !----------------------------------------------------------------------- |
---|
23 | &namtrc_dta ! Initialisation from data input file |
---|
24 | !----------------------------------------------------------------------- |
---|
25 | ! |
---|
26 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! |
---|
27 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! |
---|
28 | sn_trcdta(1) = 'data_1m_O2_nomask' , -1 , 'O2' , .false. , .true. , 'yearly' , '' , '' |
---|
29 | sn_trcdta(2) = 'data_1m_PO4_nomask' , -1 , 'PO4' , .false. , .true. , 'yearly' , '' , '' |
---|
30 | sn_trcdta(3) = 'data_1m_NO3_nomask' , -1 , 'NO3' , .false. , .true. , 'yearly' , '' , '' |
---|
31 | sn_trcdta(6) = 'data_1m_Si_nomask' , -1 , 'Si' , .false. , .true. , 'yearly' , '' , '' |
---|
32 | sn_trcdta(8) = 'data_1m_Fer_nomask' , -12 , 'Fer' , .false. , .true. , 'yearly' , '' , '' |
---|
33 | sn_trcdta(45) = 'data_1m_DOC_nomask' , -12 , 'DOC' , .false. , .true. , 'yearly' , '' , '' |
---|
34 | sn_trcdta(56) = 'data_1m_DIC_nomask' , -12 , 'DIC' , .false. , .true. , 'yearly' , '' , '' |
---|
35 | sn_trcdta(57) = 'data_1m_Alkalini_nomask', -12 , 'Alkalini', .false. , .true. , 'yearly' , '' , '' |
---|
36 | ! |
---|
37 | cn_dir = './' ! root directory for the location of the data files |
---|
38 | rn_trfac(1) = 22.4 ! conversion factor from ml/l to mmol m3 (1 mol O2 = 22.4 l) |
---|
39 | rn_trfac(2) = 1.0 ! - - - - |
---|
40 | rn_trfac(3) = 1.0 ! - - - - |
---|
41 | rn_trfac(6) = 1.0 ! - - - - |
---|
42 | rn_trfac(8) = 1.e9 ! conversion factor from mol/l to umol/m3 |
---|
43 | rn_trfac(45) = 12.e6 ! conversion factor from mol/l to mg/m3 |
---|
44 | rn_trfac(56) = 12.0 ! conversion factor from umol/l to mg/m3 |
---|
45 | rn_trfac(57) = 1.0 ! - - - - |
---|
46 | / |
---|
47 | !----------------------------------------------------------------------- |
---|
48 | &namtrc_bc |
---|
49 | !----------------------------------------------------------------------- |
---|
50 | ! ! file name ! frequency (hr) ! variable ! time interp. ! clim !'yearly' or ! weights ! rotation |
---|
51 | ! ! ! (if <0 mon) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing |
---|
52 | sn_trcsbc(3) = 'ndeposition.orca' , -1 , 'ndep' , .false. , .true. , 'yearly' , '' , '' |
---|
53 | sn_trcsbc(8) = 'dust.orca' , -1 , 'dust' , .false. , .true. , 'yearly' , '' , '' |
---|
54 | |
---|
55 | rn_trsfac(3) = 1.0 |
---|
56 | rn_trsfac(8) = 5.8866e+05 ! conversion factor from kg Fe/m2/mon to umol/m2/d = 1.e9/(55.85)*12/365 |
---|
57 | ! |
---|
58 | sn_trccbc(2) = 'river' , -12 , 'RIVERPHOS' , .false. , .true. , 'yearly' , '' , '' |
---|
59 | sn_trccbc(3) = 'river' , -12 , 'RIVERNIT' , .false. , .true. , 'yearly' , '' , '' |
---|
60 | sn_trccbc(6) = 'river' , -12 , 'RIVERSIL' , .false. , .true. , 'yearly' , '' , '' |
---|
61 | sn_trccbc(8) = 'river' , -12 , 'RIVERFe80' , .false. , .true. , 'yearly' , '' , '' |
---|
62 | sn_trccbc(45) = 'river' , -12 , 'RIVERDOC' , .false. , .true. , 'yearly' , '' , '' |
---|
63 | sn_trccbc(56) = 'river' , -12 , 'RIVERDIC' , .false. , .true. , 'yearly' , '' , '' |
---|
64 | 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 |
---|
65 | 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 |
---|
66 | 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 |
---|
67 | 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 |
---|
68 | rn_trcfac(45) = 2.7397e+12 ! 1.e15/365 ! conversion factor from Tg/y to mg/d |
---|
69 | rn_trcfac(56) = 2.7397e+12 ! 1.e15/365 ! conversion factor from Tg/y to mg/d |
---|
70 | cn_dir = './' ! root directory for the location of the boundary condition files |
---|
71 | / |
---|
72 | !----------------------------------------------------------------------- |
---|
73 | &namtrc_adv ! advection scheme for passive tracer |
---|
74 | !----------------------------------------------------------------------- |
---|
75 | ln_trcadv_cen2 = .false. ! 2nd order centered scheme |
---|
76 | ln_trcadv_tvd = .false. ! TVD scheme |
---|
77 | ln_trcadv_muscl = .true. ! MUSCL scheme |
---|
78 | ln_trcadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries |
---|
79 | ln_trcadv_ubs = .false. ! UBS scheme |
---|
80 | ln_trcadv_qck = .false. ! QUICKEST scheme |
---|
81 | / |
---|
82 | !----------------------------------------------------------------------- |
---|
83 | &namtrc_ldf ! lateral diffusion scheme for passive tracer |
---|
84 | !----------------------------------------------------------------------- |
---|
85 | ln_trcldf_diff = .true. ! performs lateral diffusion (T) or not (F) |
---|
86 | ! ! Type of the operator : |
---|
87 | ln_trcldf_lap = .true. ! laplacian operator |
---|
88 | ln_trcldf_bilap = .false. ! bilaplacian operator |
---|
89 | ! Direction of action : |
---|
90 | ln_trcldf_level = .false. ! iso-level |
---|
91 | ln_trcldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" when ln_sco=T) |
---|
92 | ln_trcldf_iso = .true. ! iso-neutral (require "key_ldfslp") |
---|
93 | ! ! Coefficient |
---|
94 | rn_ahtrc_0 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] |
---|
95 | rn_ahtrb_0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] |
---|
96 | / |
---|
97 | !----------------------------------------------------------------------- |
---|
98 | &namtrc_zdf ! vertical physics |
---|
99 | !----------------------------------------------------------------------- |
---|
100 | ln_trczdf_exp = .false. ! split explicit (T) or implicit (F) time stepping |
---|
101 | nn_trczdf_exp = 3 ! number of sub-timestep for ln_trczdfexp=T |
---|
102 | / |
---|
103 | !----------------------------------------------------------------------- |
---|
104 | &namtrc_rad ! treatment of negative concentrations |
---|
105 | !----------------------------------------------------------------------- |
---|
106 | ln_trcrad = .true. ! artificially correct negative concentrations (T) or not (F) |
---|
107 | / |
---|
108 | !----------------------------------------------------------------------- |
---|
109 | &namtrc_dmp ! passive tracer newtonian damping |
---|
110 | !----------------------------------------------------------------------- |
---|
111 | ln_trcdmp = .false. ! add a damping term (T) or not (F) |
---|
112 | nn_hdmp_tr = -1 ! horizontal shape =-1, damping in Med and Red Seas only |
---|
113 | ! =XX, damping poleward of XX degrees (XX>0) |
---|
114 | ! + F(distance-to-coast) + Red and Med Seas |
---|
115 | nn_zdmp_tr = 1 ! vertical shape =0 damping throughout the water column |
---|
116 | ! =1 no damping in the mixing layer (kz criteria) |
---|
117 | ! =2 no damping in the mixed layer (rho crieria) |
---|
118 | rn_surf_tr = 50. ! surface time scale of damping [days] |
---|
119 | rn_bot_tr = 360. ! bottom time scale of damping [days] |
---|
120 | rn_dep_tr = 800. ! depth of transition between rn_surf and rn_bot [meters] |
---|
121 | nn_file_tr = 0 ! create a damping.coeff NetCDF file (=1) or not (=0) |
---|
122 | / |
---|
123 | !----------------------------------------------------------------------- |
---|
124 | &namtrc_trd ! diagnostics on tracer trends ('key_trdtrc') |
---|
125 | ! or mixed-layer trends ('key_trdmld_trc') |
---|
126 | !---------------------------------------------------------------------- |
---|
127 | nn_trd_trc = 360 ! time step frequency and tracers trends |
---|
128 | nn_ctls_trc = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) |
---|
129 | rn_ucf_trc = 86400 ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) |
---|
130 | ln_trdmld_trc_restart = .false. ! restart for ML diagnostics |
---|
131 | ln_trdmld_trc_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S |
---|
132 | ln_trdtrc(1) = .true. |
---|
133 | ln_trdtrc(2) = .true. |
---|
134 | ln_trdtrc(23) = .true. |
---|
135 | / |
---|
136 | !----------------------------------------------------------------------- |
---|
137 | &namtrc_dia ! parameters for passive tracer additional diagnostics |
---|
138 | !---------------------------------------------------------------------- |
---|
139 | ln_diatrc = .false. ! save additional diag. (T) or not (F) |
---|
140 | ln_diabio = .false. |
---|
141 | nn_writedia = 5475 ! time step frequency for diagnostics |
---|
142 | nn_writebio = 5475 |
---|
143 | / |
---|