1 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
---|
2 | !! NEMO/OPA Configuration namelist : used to overwrite defaults values defined in SHARED/namelist_ref |
---|
3 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
---|
4 | ! |
---|
5 | !----------------------------------------------------------------------- |
---|
6 | &namrun ! parameters of the run |
---|
7 | !----------------------------------------------------------------------- |
---|
8 | nn_no = 1 ! job number (no more used...) |
---|
9 | cn_exp = "TROP04-AMZ00" ! experience name |
---|
10 | nn_it000 = 1 ! first time step |
---|
11 | nn_itend = 3504 ! last time step (std 5475) |
---|
12 | nn_date0 = 19930101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) |
---|
13 | nn_leapy = 0 ! Leap year calendar (1) or not (0) |
---|
14 | ln_rstart = false ! start from rest (F) or from a restart file (T) |
---|
15 | nn_euler = 1 ! = 0 : start with forward time step if ln_rstart=T |
---|
16 | nn_rstctl = 0 ! restart control ==> activated only if ln_rstart=T |
---|
17 | ! ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist |
---|
18 | ! ! = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart |
---|
19 | ! ! = 2 nn_date0 read in restart ; nn_it000 : check consistancy between namelist and restart |
---|
20 | nn_istate = 0 ! output the initial state (1) or not (0) |
---|
21 | nn_stock = 99999 ! frequency of creation of a restart file (modulo referenced to 1) |
---|
22 | nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written |
---|
23 | nn_write = 99999 ! frequency of write in the output file (modulo referenced to nn_it000) |
---|
24 | |
---|
25 | / |
---|
26 | !----------------------------------------------------------------------- |
---|
27 | &namcfg ! parameters of the configuration |
---|
28 | !----------------------------------------------------------------------- |
---|
29 | cp_cfg = "trop04" ! name of the configuration |
---|
30 | jp_cfg = 004 ! resolution of the configuration |
---|
31 | jpidta = 301 ! 1st lateral dimension ( >= jpi ) |
---|
32 | jpjdta = 165 ! 2nd " " ( >= jpj ) |
---|
33 | jpkdta = 75 ! number of levels ( >= jpk ) |
---|
34 | jpiglo = 301 ! 1st dimension of global domain --> i =jpidta |
---|
35 | jpjglo = 165 ! 2nd - - --> j =jpjdta |
---|
36 | jpizoom = 1 ! left bottom (i,j) indices of the zoom |
---|
37 | jpjzoom = 1 ! in data domain indices |
---|
38 | jperio = 0 ! lateral cond. type (between 0 and 6) |
---|
39 | / |
---|
40 | !----------------------------------------------------------------------- |
---|
41 | &namzgr ! vertical coordinate |
---|
42 | !----------------------------------------------------------------------- |
---|
43 | ln_zps = .true. ! z-coordinate - partial steps |
---|
44 | ln_linssh = .true. ! linear free surface |
---|
45 | / |
---|
46 | !----------------------------------------------------------------------- |
---|
47 | &namdom ! space and time domain (bathymetry, mesh, timestep) |
---|
48 | !----------------------------------------------------------------------- |
---|
49 | nn_bathy = 1 ! compute (=0) or read (=1) the bathymetry file |
---|
50 | rn_bathy = 0. ! value of the bathymetry. if (=0) bottom flat at jpkm1 |
---|
51 | nn_closea = 0 ! remove (=0) or keep (=1) closed seas and lakes (ORCA) |
---|
52 | nn_msh = 1 ! create (=1) a mesh file or not (=0) |
---|
53 | rn_hmin = -3. ! min depth of the ocean (>0) or min number of ocean level (<0) |
---|
54 | rn_isfhmin = 1.00 ! treshold (m) to discriminate grounding ice to floating ice |
---|
55 | rn_e3zps_min= 20. ! partial step thickness is set larger than the minimum of |
---|
56 | rn_e3zps_rat= 0.1 ! rn_e3zps_min and rn_e3zps_rat*e3t, with 0<rn_e3zps_rat<1 |
---|
57 | ! |
---|
58 | rn_rdt = 1800. ! time step for the dynamics (and tracer if nn_acc=0) |
---|
59 | rn_atfp = 0.1 ! asselin time filter parameter |
---|
60 | ln_crs = .false. ! Logical switch for coarsening module |
---|
61 | jphgr_msh = 0 ! type of horizontal mesh |
---|
62 | ! = 0 curvilinear coordinate on the sphere read in coordinate.nc |
---|
63 | ! = 1 geographical mesh on the sphere with regular grid-spacing |
---|
64 | ! = 2 f-plane with regular grid-spacing |
---|
65 | ! = 3 beta-plane with regular grid-spacing |
---|
66 | ! = 4 Mercator grid with T/U point at the equator |
---|
67 | ppglam0 = 999999.0 ! longitude of first raw and column T-point (jphgr_msh = 1) |
---|
68 | ppgphi0 = 999999.0 ! latitude of first raw and column T-point (jphgr_msh = 1) |
---|
69 | ppe1_deg = 999999.0 ! zonal grid-spacing (degrees) |
---|
70 | ppe2_deg = 999999.0 ! meridional grid-spacing (degrees) |
---|
71 | ppe1_m = 999999.0 ! zonal grid-spacing (degrees) |
---|
72 | ppe2_m = 999999.0 ! meridional grid-spacing (degrees) |
---|
73 | ppsur = -3958.951371276829 ! ORCA r4, r2 and r05 coefficients |
---|
74 | ppa0 = 103.9530096000000 ! (default coefficients) |
---|
75 | ppa1 = 2.4159512690000 ! |
---|
76 | ppkth = 15.35101370000000 ! |
---|
77 | ppacr = 7.0000 ! |
---|
78 | ppdzmin = 999999. ! Minimum vertical spacing |
---|
79 | pphmax = 999999. ! Maximum depth |
---|
80 | ldbletanh = .TRUE. ! Use/do not use double tanf function for vertical coordinates |
---|
81 | ppa2 = 100.760928500000 ! Double tanh function parameters |
---|
82 | ppkth2 = 48.029893720000 ! |
---|
83 | ppacr2 = 13.000000000000 ! |
---|
84 | / |
---|
85 | !----------------------------------------------------------------------- |
---|
86 | &namcrs ! Grid coarsening for dynamics output and/or |
---|
87 | ! passive tracer coarsened online simulations |
---|
88 | !----------------------------------------------------------------------- |
---|
89 | / |
---|
90 | !----------------------------------------------------------------------- |
---|
91 | &namtsd ! data : Temperature & Salinity |
---|
92 | !----------------------------------------------------------------------- |
---|
93 | / |
---|
94 | !----------------------------------------------------------------------- |
---|
95 | &namsbc ! Surface Boundary Condition (surface module) |
---|
96 | !----------------------------------------------------------------------- |
---|
97 | nn_fsbc = 4 ! frequency of surface boundary condition computation |
---|
98 | ! (also = the frequency of sea-ice & iceberg model call) |
---|
99 | ln_dm2dc = .true. ! daily mean to diurnal cycle on short wave |
---|
100 | ln_rnf = .true. ! runoffs (T => fill namsbc_rnf) |
---|
101 | ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) |
---|
102 | nn_fwb = 0 ! FreshWater Budget: =0 unchecked |
---|
103 | ! =1 global mean of e-p-r set to zero at each time step |
---|
104 | ! =2 annual global mean of e-p-r set to zero |
---|
105 | / |
---|
106 | !----------------------------------------------------------------------- |
---|
107 | &namsbc_core ! namsbc_core CORE bulk formulae |
---|
108 | !----------------------------------------------------------------------- |
---|
109 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
---|
110 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
---|
111 | sn_wndi = 'drowned_u10_DFS5.2' , 3 , 'u10' , .true. , .false. , 'yearly' , 'weight_bicubic.nc' , '', '' |
---|
112 | sn_wndj = 'drowned_v10_DFS5.2' , 3 , 'v10' , .true. , .false. , 'yearly' , 'weight_bicubic.nc' , '', '' |
---|
113 | sn_qsr = 'drowned_radsw_DFS5.2' , 24 , 'radsw' , .false. , .false. , 'yearly' , 'weight_bilin.nc' , '', '' |
---|
114 | sn_qlw = 'drowned_radlw_DFS5.2' , 24 , 'radlw' , .false. , .false. , 'yearly' , 'weight_bilin.nc' , '', '' |
---|
115 | sn_tair = 'drowned_t2_DFS5.2' , 3 , 't2' , .true. , .false. , 'yearly' , 'weight_bilin.nc' , '', '' |
---|
116 | sn_humi = 'drowned_q2_DFS5.2' , 3 , 'q2' , .true. , .false. , 'yearly' , 'weight_bilin.nc' , '', '' |
---|
117 | sn_prec = 'drowned_precip_DFS5.2' , 24 , 'precip' , .true. , .false. , 'yearly' , 'weight_bilin.nc' , '', '' |
---|
118 | sn_snow = 'drowned_snow_DFS5.2' , 24 , 'snow' , .true. , .false. , 'yearly' , 'weight_bilin.nc' , '', '' |
---|
119 | sn_tdif = 'taudif' , 24 , 'taudif' , .true. , .false. , 'yearly' , '' , '', '' |
---|
120 | ! |
---|
121 | cn_dir = './' ! root directory for the location of the bulk files |
---|
122 | ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data |
---|
123 | rn_zqt = 2. ! Air temperature and humidity reference height (m) |
---|
124 | rn_zu = 10. ! Wind vector reference height (m) |
---|
125 | rn_pfac = 1. ! multiplicative factor for precipitation (total & snow) |
---|
126 | rn_efac = 1. ! multiplicative factor for evaporation (0. or 1.) |
---|
127 | rn_vfac = 1 ! multiplicative factor for ocean/ice velocity |
---|
128 | ! in the calculation of the wind stress (0.=absolute winds or 1.=relative winds) |
---|
129 | / |
---|
130 | !----------------------------------------------------------------------- |
---|
131 | &namtra_qsr ! penetrative solar radiation |
---|
132 | !----------------------------------------------------------------------- |
---|
133 | sn_chl ='chlorophyll', -1 , 'CHLA' , .true. , .true. , 'yearly' , '' , '' , '' |
---|
134 | |
---|
135 | cn_dir = './' ! root directory for the location of the runoff files |
---|
136 | ln_qsr_rgb = .true. ! RGB (Red-Green-Blue) light penetration |
---|
137 | nn_chldta = 0 ! RGB : Chl data (=1) or cst value (=0) |
---|
138 | / |
---|
139 | !----------------------------------------------------------------------- |
---|
140 | &namsbc_rnf ! runoffs namelist surface boundary condition |
---|
141 | !----------------------------------------------------------------------- |
---|
142 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
---|
143 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! |
---|
144 | sn_rnf = 'runoff_1m_nomask', -1 , 'sorunoff', .true. , .true. , 'yearly' , '' , '' , '' |
---|
145 | sn_cnf = 'runoff_1m_nomask', 0 , 'socoefr', .false. , .true. , 'yearly' , '' , '' , '' |
---|
146 | sn_s_rnf = 'runoffs' , 24 , 'rosaline', .true. , .true. , 'yearly' , '' , '' , '' |
---|
147 | sn_t_rnf = 'runoffs' , 24 , 'rotemper', .true. , .true. , 'yearly' , '' , '' , '' |
---|
148 | sn_dep_rnf = 'runoffs' , 0 , 'rodepth' , .false. , .true. , 'yearly' , '' , '' , '' |
---|
149 | |
---|
150 | cn_dir = './' ! root directory for the location of the runoff files |
---|
151 | ln_rnf_mouth= .true. ! specific treatment at rivers mouths |
---|
152 | rn_hrnf = 15.e0 ! depth over which enhanced vertical mixing is used (ln_rnf_mouth=T) |
---|
153 | rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T) |
---|
154 | rn_rfact = 1.e0 ! multiplicative factor for runoff |
---|
155 | ln_rnf_depth= .false. ! read in depth information for runoff |
---|
156 | ln_rnf_tem = .false. ! read in temperature information for runoff |
---|
157 | ln_rnf_sal = .false. ! read in salinity information for runoff |
---|
158 | ln_rnf_depth_ini = .false. ! compute depth at initialisation from runoff file |
---|
159 | rn_rnf_max = 5.735e-4 ! max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true ) |
---|
160 | rn_dep_max = 150. ! depth over which runoffs is spread ( ln_rnf_depth_ini = .true ) |
---|
161 | nn_rnf_depth_file = 0 ! create (=1) a runoff depth file or not (=0) |
---|
162 | |
---|
163 | / |
---|
164 | !----------------------------------------------------------------------- |
---|
165 | &namsbc_ssr ! surface boundary condition : sea surface restoring |
---|
166 | !----------------------------------------------------------------------- |
---|
167 | / |
---|
168 | !----------------------------------------------------------------------- |
---|
169 | &namsbc_alb ! albedo parameters |
---|
170 | !----------------------------------------------------------------------- |
---|
171 | / |
---|
172 | !----------------------------------------------------------------------- |
---|
173 | &namberg ! iceberg parameters |
---|
174 | !----------------------------------------------------------------------- |
---|
175 | / |
---|
176 | !----------------------------------------------------------------------- |
---|
177 | &namlbc ! lateral momentum boundary condition |
---|
178 | !----------------------------------------------------------------------- |
---|
179 | ! ! free slip ! partial slip ! no slip ! strong slip |
---|
180 | rn_shlat = 0. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat |
---|
181 | ln_vorlat = .false. ! consistency of vorticity boundary condition with analytical Eqs. |
---|
182 | / |
---|
183 | !----------------------------------------------------------------------- |
---|
184 | &namagrif ! AGRIF zoom ("key_agrif") |
---|
185 | !----------------------------------------------------------------------- |
---|
186 | nn_cln_update = 3 ! baroclinic update frequency |
---|
187 | ln_spc_dyn = .true. ! use 0 as special value for dynamics |
---|
188 | rn_sponge_tra = 200. ! coefficient for tracer sponge layer [m2/s] |
---|
189 | rn_sponge_dyn = 200. ! coefficient for dynamics sponge layer [m2/s] |
---|
190 | ln_chk_bathy = .FALSE. ! |
---|
191 | / |
---|
192 | !----------------------------------------------------------------------- |
---|
193 | &nam_tide ! tide parameters ("key_tide") |
---|
194 | !----------------------------------------------------------------------- |
---|
195 | ln_tide_pot = .true. ! use tidal potential forcing |
---|
196 | ln_tide_ramp= .false. ! |
---|
197 | rdttideramp = 0. ! |
---|
198 | clname(1) = 'DUMMY' ! name of constituent - all tidal components must be set in namelist_cfg |
---|
199 | / |
---|
200 | !----------------------------------------------------------------------- |
---|
201 | &nambdy ! cross land advection |
---|
202 | !----------------------------------------------------------------------- |
---|
203 | nb_bdy = 1 ! number of open boundary sets |
---|
204 | ln_coords_file = .true. ! =T : read bdy coordinates from file |
---|
205 | cn_coords_file = 'bdy_coordinates.nc' ! bdy coordinates files |
---|
206 | ln_mask_file = .false. ! =T : read mask from file |
---|
207 | cn_mask_file = '' ! name of mask file (if ln_mask_file=.TRUE.) |
---|
208 | cn_dyn2d = 'flather' ! |
---|
209 | nn_dyn2d_dta = 1 ! = 0, bdy data are equal to the initial state |
---|
210 | ! = 1, bdy data are read in 'bdydata .nc' files |
---|
211 | ! = 2, use tidal harmonic forcing data from files |
---|
212 | ! = 3, use external data AND tidal harmonic forcing |
---|
213 | cn_dyn3d = 'frs' ! |
---|
214 | nn_dyn3d_dta = 1 ! = 0, bdy data are equal to the initial state |
---|
215 | ! = 1, bdy data are read in 'bdydata .nc' files |
---|
216 | cn_tra = 'frs' ! |
---|
217 | nn_tra_dta = 1 ! = 0, bdy data are equal to the initial state |
---|
218 | ! = 1, bdy data are read in 'bdydata .nc' files |
---|
219 | ln_tra_dmp =.false. ! open boudaries conditions for tracers |
---|
220 | ln_dyn3d_dmp =.false. ! open boundary condition for baroclinic velocities |
---|
221 | rn_time_dmp = 1. ! Damping time scale in days |
---|
222 | rn_time_dmp_out = 1. ! Outflow damping time scale |
---|
223 | nn_rimwidth = 1 ! width of the relaxation zone |
---|
224 | ln_vol = .true. ! total volume correction (see nn_volctl parameter) |
---|
225 | nn_volctl = 0 ! = 0, the total water flux across open boundaries is zero |
---|
226 | / |
---|
227 | !----------------------------------------------------------------------- |
---|
228 | &nambdy_dta ! open boundaries - external data ("key_bdy") |
---|
229 | !----------------------------------------------------------------------- |
---|
230 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! |
---|
231 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename ! |
---|
232 | bn_ssh = 'bdyT_u2d' , 24 , 'sossheig' , .true. , .false. , 'yearly' , '' , '' , '' |
---|
233 | bn_u2d = 'bdyU_u2d' , 24 , 'vobtcrtx' , .true. , .false. , 'yearly' , '' , '' , '' |
---|
234 | bn_v2d = 'bdyV_u2d' , 24 , 'vobtcrty' , .true. , .false. , 'yearly' , '' , '' , '' |
---|
235 | bn_u3d = 'bdyU_u3d' , 24 , 'vozocrtx' , .true. , .false. , 'yearly' , '' , '' , '' |
---|
236 | bn_v3d = 'bdyV_u3d' , 24 , 'vomecrty' , .true. , .false. , 'yearly' , '' , '' , '' |
---|
237 | bn_tem = 'bdyT_tra' , 24 , 'votemper' , .true. , .false. , 'yearly' , '' , '' , '' |
---|
238 | bn_sal = 'bdyT_tra' , 24 , 'vosaline' , .true. , .false. , 'yearly' , '' , '' , '' |
---|
239 | cn_dir = './' |
---|
240 | ln_full_vel = .false. |
---|
241 | / |
---|
242 | !----------------------------------------------------------------------- |
---|
243 | &nambdy_tide ! tidal forcing at open boundaries |
---|
244 | !----------------------------------------------------------------------- |
---|
245 | filtide = './GOLFO36_bdytide_' ! file name root of tidal forcing files |
---|
246 | ln_bdytide_2ddta = .false. |
---|
247 | ln_bdytide_conj = .false. |
---|
248 | / |
---|
249 | !----------------------------------------------------------------------- |
---|
250 | &nam_tide ! tide parameters (#ifdef key_tide) |
---|
251 | !----------------------------------------------------------------------- |
---|
252 | ln_tide_pot = .true. ! use tidal potential forcing |
---|
253 | ln_tide_ramp = .false. ! |
---|
254 | rdttideramp = 0. ! |
---|
255 | clname(1) = 'M2' ! name of constituent |
---|
256 | clname(2) = 'S2' |
---|
257 | clname(3) = 'N2' |
---|
258 | clname(4) = 'K1' |
---|
259 | clname(5) = 'O1' |
---|
260 | clname(6) = 'Q1' |
---|
261 | clname(7) = 'M4' |
---|
262 | clname(8) = 'K2' |
---|
263 | clname(9) = 'P1' |
---|
264 | / |
---|
265 | !----------------------------------------------------------------------- |
---|
266 | &nambfr ! bottom friction |
---|
267 | !----------------------------------------------------------------------- |
---|
268 | nn_bfr = 1 ! type of bottom friction : = 0 : free slip, = 1 : linear friction |
---|
269 | ! = 2 : nonlinear friction |
---|
270 | rn_bfri1 = 4.e-4 ! bottom drag coefficient (linear case) |
---|
271 | rn_bfri2 = 1.e-3 ! bottom drag coefficient (non linear case). Minimum coeft if ln_loglayer=T |
---|
272 | rn_bfri2_max= 1.e-1 ! max. bottom drag coefficient (non linear case and ln_loglayer=T) |
---|
273 | rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m2/s2) |
---|
274 | rn_bfrz0 = 3.e-3 ! bottom roughness [m] if ln_loglayer=T |
---|
275 | ln_bfr2d = .false. ! horizontal variation of the bottom friction coef (read a 2D mask file ) |
---|
276 | rn_bfrien = 50. ! local multiplying factor of bfr (ln_bfr2d=T) |
---|
277 | rn_tfri1 = 4.e-4 ! top drag coefficient (linear case) |
---|
278 | rn_tfri2 = 2.5e-3 ! top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T |
---|
279 | rn_tfri2_max= 1.e-1 ! max. top drag coefficient (non linear case and ln_loglayer=T) |
---|
280 | rn_tfeb2 = 0.0 ! top turbulent kinetic energy background (m2/s2) |
---|
281 | rn_tfrz0 = 3.e-3 ! top roughness [m] if ln_loglayer=T |
---|
282 | ln_tfr2d = .false. ! horizontal variation of the top friction coef (read a 2D mask file ) |
---|
283 | rn_tfrien = 50. ! local multiplying factor of tfr (ln_tfr2d=T) |
---|
284 | |
---|
285 | ln_bfrimp = .true. ! implicit bottom friction (requires ln_zdfexp = .false. if true) |
---|
286 | ln_loglayer = .false. ! logarithmic formulation (non linear case) |
---|
287 | / |
---|
288 | !----------------------------------------------------------------------- |
---|
289 | &nambbc ! bottom temperature boundary condition (default: NO) |
---|
290 | !----------------------------------------------------------------------- |
---|
291 | ln_trabbc = .false. ! Apply a geothermal heating at the ocean bottom |
---|
292 | / |
---|
293 | !----------------------------------------------------------------------- |
---|
294 | &nambbl ! bottom boundary layer scheme |
---|
295 | !----------------------------------------------------------------------- |
---|
296 | / |
---|
297 | !----------------------------------------------------------------------- |
---|
298 | &nameos ! ocean physical parameters |
---|
299 | !----------------------------------------------------------------------- |
---|
300 | ln_teos10 = .true. ! = Use TEOS-10 equation of state |
---|
301 | / |
---|
302 | !----------------------------------------------------------------------- |
---|
303 | &namtra_adv ! advection scheme for tracer |
---|
304 | !----------------------------------------------------------------------- |
---|
305 | ln_traadv_fct = .true. ! FCT scheme |
---|
306 | nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order |
---|
307 | nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order |
---|
308 | nn_fct_zts = 0 ! > 1 , 2nd order FCT scheme with vertical sub-timestepping |
---|
309 | ! ! (number of sub-timestep = nn_fct_zts) |
---|
310 | / |
---|
311 | !----------------------------------------------------------------------- |
---|
312 | &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) |
---|
313 | !----------------------------------------------------------------------- |
---|
314 | / |
---|
315 | !---------------------------------------------------------------------------------- |
---|
316 | &namtra_ldf ! lateral diffusion scheme for tracers |
---|
317 | !---------------------------------------------------------------------------------- |
---|
318 | ! ! Operator type: |
---|
319 | ln_traldf_lap = .true. ! laplacian operator |
---|
320 | ln_traldf_blp = .false. ! bilaplacian operator |
---|
321 | ! ! Direction of action: |
---|
322 | ln_traldf_lev = .false. ! iso-level |
---|
323 | ln_traldf_hor = .false. ! horizontal (geopotential) |
---|
324 | ln_traldf_iso = .true. ! iso-neutral (standard operator) |
---|
325 | ln_traldf_triad = .false. ! iso-neutral (triad operator) |
---|
326 | ! |
---|
327 | ! ! iso-neutral options: |
---|
328 | ln_traldf_msc = .true. ! Method of Stabilizing Correction (both operators) |
---|
329 | rn_slpmax = 0.01 ! slope limit (both operators) |
---|
330 | ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) |
---|
331 | rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) |
---|
332 | ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) |
---|
333 | ! |
---|
334 | ! ! Coefficients: |
---|
335 | nn_aht_ijk_t = 20 ! space/time variation of eddy coef |
---|
336 | ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file |
---|
337 | ! ! = 0 constant |
---|
338 | ! ! = 10 F(k) =ldf_c1d |
---|
339 | ! ! = 20 F(i,j) =ldf_c2d |
---|
340 | ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation |
---|
341 | ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d |
---|
342 | ! ! = 31 F(i,j,k,t)=F(local velocity) |
---|
343 | rn_aht_0 = 300. ! lateral eddy diffusivity (lap. operator) [m2/s] |
---|
344 | rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] |
---|
345 | / |
---|
346 | !---------------------------------------------------------------------------------- |
---|
347 | &namtra_ldfeiv ! eddy induced velocity param. |
---|
348 | !---------------------------------------------------------------------------------- |
---|
349 | / |
---|
350 | !----------------------------------------------------------------------- |
---|
351 | &namtra_dmp ! tracer: T & S newtonian damping |
---|
352 | !----------------------------------------------------------------------- |
---|
353 | ln_tradmp = .false. ! add a damping termn (T) or not (F) |
---|
354 | / |
---|
355 | !----------------------------------------------------------------------- |
---|
356 | &namdyn_adv ! formulation of the momentum advection |
---|
357 | !----------------------------------------------------------------------- |
---|
358 | ln_dynadv_vec = .false. ! vector form (T) or flux form (F) |
---|
359 | nn_dynkeg = 0 ! scheme for grad(KE): =0 C2 ; =1 Hollingsworth correction |
---|
360 | ln_dynadv_cen2= .false. ! flux form - 2nd order centered scheme |
---|
361 | ln_dynadv_ubs = .true. ! flux form - 3rd order UBS scheme |
---|
362 | ln_dynzad_zts = .false. ! Use (T) sub timestepping for vertical momentum advection |
---|
363 | / |
---|
364 | !----------------------------------------------------------------------- |
---|
365 | &namdyn_vor ! option of physics/algorithm (not control by CPP keys) |
---|
366 | !----------------------------------------------------------------------- |
---|
367 | ln_dynvor_ene = .false. ! enstrophy conserving scheme |
---|
368 | ln_dynvor_ens = .false. ! energy conserving scheme |
---|
369 | ln_dynvor_mix = .false. ! mixed scheme |
---|
370 | ln_dynvor_een = .true. ! energy & enstrophy scheme |
---|
371 | nn_een_e3f = 0 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) |
---|
372 | / |
---|
373 | !----------------------------------------------------------------------- |
---|
374 | &namdyn_hpg ! Hydrostatic pressure gradient option |
---|
375 | !----------------------------------------------------------------------- |
---|
376 | / |
---|
377 | !----------------------------------------------------------------------- |
---|
378 | &namdyn_spg ! Surface pressure gradient |
---|
379 | !----------------------------------------------------------------------- |
---|
380 | ln_dynspg_ts = .true. ! Split-explicit free surface |
---|
381 | / |
---|
382 | !----------------------------------------------------------------------- |
---|
383 | &namdyn_ldf ! lateral diffusion on momentum |
---|
384 | !----------------------------------------------------------------------- |
---|
385 | ! ! Type of the operator : |
---|
386 | ! ! no diffusion: set ln_dynldf_lap=..._blp=F |
---|
387 | ln_dynldf_lap = .false. ! laplacian operator |
---|
388 | ln_dynldf_blp = .false. ! bilaplacian operator |
---|
389 | ! ! Direction of action : |
---|
390 | ln_dynldf_lev = .false. ! iso-level |
---|
391 | ln_dynldf_hor = .false. ! horizontal (geopotential) |
---|
392 | ln_dynldf_iso = .false. ! iso-neutral |
---|
393 | ! ! Coefficient |
---|
394 | nn_ahm_ijk_t = 20 ! space/time variation of eddy coef |
---|
395 | ! ! =-30 read in eddy_viscosity_3D.nc file |
---|
396 | ! ! =-20 read in eddy_viscosity_2D.nc file |
---|
397 | ! ! = 0 constant |
---|
398 | ! ! = 10 F(k)=c1d |
---|
399 | ! ! = 20 F(i,j)=F(grid spacing)=c2d |
---|
400 | ! ! = 30 F(i,j,k)=c2d*c1d |
---|
401 | ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) |
---|
402 | rn_ahm_0 = 40000. ! horizontal laplacian eddy viscosity [m2/s] |
---|
403 | rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] |
---|
404 | rn_bhm_0 = 1.e+12 ! horizontal bilaplacian eddy viscosity [m4/s] |
---|
405 | / |
---|
406 | !----------------------------------------------------------------------- |
---|
407 | &namzdf ! vertical physics |
---|
408 | !----------------------------------------------------------------------- |
---|
409 | / |
---|
410 | !----------------------------------------------------------------------- |
---|
411 | &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") |
---|
412 | !----------------------------------------------------------------------- |
---|
413 | / |
---|
414 | !----------------------------------------------------------------------- |
---|
415 | &namzdf_ddm ! double diffusive mixing parameterization ("key_zdfddm") |
---|
416 | !----------------------------------------------------------------------- |
---|
417 | / |
---|
418 | !----------------------------------------------------------------------- |
---|
419 | &namzdf_tmx ! tidal mixing parameterization ("key_zdftmx") |
---|
420 | !----------------------------------------------------------------------- |
---|
421 | / |
---|
422 | !----------------------------------------------------------------------- |
---|
423 | &nammpp ! Massively Parallel Processing ("key_mpp_mpi) |
---|
424 | !----------------------------------------------------------------------- |
---|
425 | cn_mpi_send = 'I' ! mpi send/recieve type ='S', 'B', or 'I' for standard send, |
---|
426 | ! buffer blocking send or immediate non-blocking sends, resp. |
---|
427 | nn_buffer = 0 ! size in bytes of exported buffer ('B' case), 0 no exportation |
---|
428 | ln_nnogather= .false. ! activate code to avoid mpi_allgather use at the northfold |
---|
429 | jpni = 9 ! jpni number of processors following i (set automatically if < 1) |
---|
430 | jpnj = 20 ! jpnj number of processors following j (set automatically if < 1) |
---|
431 | jpnij = 180 ! jpnij number of local domains (set automatically if < 1) |
---|
432 | / |
---|
433 | !----------------------------------------------------------------------- |
---|
434 | &namctl ! Control prints & Benchmark |
---|
435 | !----------------------------------------------------------------------- |
---|
436 | / |
---|
437 | !----------------------------------------------------------------------- |
---|
438 | &namptr ! Poleward Transport Diagnostic |
---|
439 | !----------------------------------------------------------------------- |
---|
440 | / |
---|
441 | !----------------------------------------------------------------------- |
---|
442 | &namhsb ! Heat and salt budgets |
---|
443 | !----------------------------------------------------------------------- |
---|
444 | / |
---|
445 | !----------------------------------------------------------------------- |
---|
446 | &namobs ! observation usage |
---|
447 | !----------------------------------------------------------------------- |
---|
448 | / |
---|
449 | !----------------------------------------------------------------------- |
---|
450 | &nam_asminc ! assimilation increments ('key_asminc') |
---|
451 | !----------------------------------------------------------------------- |
---|
452 | / |
---|
453 | !----------------------------------------------------------------------- |
---|
454 | &nam_dia25h ! 25h Mean Output (default F) |
---|
455 | !----------------------------------------------------------------------- |
---|
456 | / |
---|
457 | |
---|