1 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
---|
2 | !! NEMO/OPA : AMM12 configuration namelist used to overwrite defaults values defined in SHARED/namelist_ref |
---|
3 | !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> |
---|
4 | !----------------------------------------------------------------------- |
---|
5 | &namrun ! parameters of the run |
---|
6 | !----------------------------------------------------------------------- |
---|
7 | cn_exp = "AMM12" ! experience name |
---|
8 | nn_itend = 2880 ! last time step (std 1 day = 288) |
---|
9 | nn_date0 = 20070101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) |
---|
10 | nn_leapy = 1 ! Leap year calendar (1) or not (0) |
---|
11 | ln_rstart = .true. ! start from rest (F) or from a restart file (T) |
---|
12 | cn_ocerst_in = "amm12.restart_20070101" ! suffix of ocean restart name (input) |
---|
13 | cn_ocerst_out = "AMM12.restart.out" ! suffix of ocean restart name (input) |
---|
14 | nn_istate = 1 ! output the initial state (1) or not (0) |
---|
15 | nn_stock = 2880 ! frequency of creation of a restart file (modulo referenced to 1) |
---|
16 | nn_write = 144 ! frequency of write in the output file (modulo referenced to nit000) |
---|
17 | / |
---|
18 | !----------------------------------------------------------------------- |
---|
19 | &namzgr ! vertical coordinate |
---|
20 | !----------------------------------------------------------------------- |
---|
21 | ln_zps = .false. ! z-coordinate - partial steps (T/F) |
---|
22 | ln_sco = .true. ! s- or hybrid z-s-coordinate (T/F) |
---|
23 | / |
---|
24 | !----------------------------------------------------------------------- |
---|
25 | &namzgr_sco ! s-coordinate or hybrid z-s-coordinate |
---|
26 | !----------------------------------------------------------------------- |
---|
27 | ln_s_sh94 = .false. ! Song & Haidvogel 1994 hybrid S-sigma (T)| |
---|
28 | ln_s_sf12 = .true. ! Siddorn & Furner 2012 hybrid S-z-sigma (T)| if both are false the NEMO tanh stretching is applied |
---|
29 | ln_sigcrit = .true. ! use sigma coordinates below critical depth (T) or Z coordinates (F) for Siddorn & Furner stretch |
---|
30 | ! stretching coefficients for all functions |
---|
31 | rn_hc = 50.0 ! critical depth for transition to stretched coordinates |
---|
32 | / |
---|
33 | !----------------------------------------------------------------------- |
---|
34 | &namdom ! space and time domain (bathymetry, mesh, timestep) |
---|
35 | !----------------------------------------------------------------------- |
---|
36 | rn_rdt = 300. ! time step for the dynamics (and tracer if nn_acc=0) |
---|
37 | nn_baro = 30 ! number of barotropic time step ("key_dynspg_ts") |
---|
38 | rn_rdtmin = 300. ! minimum time step on tracers (used if nn_acc=1) |
---|
39 | rn_rdtmax = 300. ! maximum time step on tracers (used if nn_acc=1) |
---|
40 | rn_rdth = 300. ! depth variation of tracer time step (used if nn_acc=1) |
---|
41 | / |
---|
42 | !----------------------------------------------------------------------- |
---|
43 | &namtsd ! data : Temperature & Salinity |
---|
44 | !----------------------------------------------------------------------- |
---|
45 | ln_tsd_init = .false. ! Initialisation of ocean T & S with T &S input data (T) or not (F) |
---|
46 | ln_tsd_tradmp = .false. ! damping of ocean T & S toward T &S input data (T) or not (F) |
---|
47 | / |
---|
48 | !----------------------------------------------------------------------- |
---|
49 | &namsbc ! Surface Boundary Condition (surface module) |
---|
50 | !----------------------------------------------------------------------- |
---|
51 | nn_fsbc = 1 ! frequency of surface boundary condition computation |
---|
52 | ! (also = the frequency of sea-ice model call) |
---|
53 | ln_ana = .false ! analytical formulation (T => fill namsbc_ana ) |
---|
54 | ln_flx = .true. ! flux formulation (T => fill namsbc_flx ) |
---|
55 | ln_blk_core = .false. ! CORE bulk formulation (T => fill namsbc_core) |
---|
56 | nn_ice = 0 ! =0 no ice boundary condition , |
---|
57 | ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) |
---|
58 | nn_fwb = 0 ! FreshWater Budget: =0 unchecked |
---|
59 | / |
---|
60 | !----------------------------------------------------------------------- |
---|
61 | &namsbc_ana ! analytical surface boundary condition |
---|
62 | !----------------------------------------------------------------------- |
---|
63 | / |
---|
64 | !----------------------------------------------------------------------- |
---|
65 | &namsbc_flx ! surface boundary condition : flux formulation |
---|
66 | !----------------------------------------------------------------------- |
---|
67 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! |
---|
68 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! |
---|
69 | sn_utau = 'amm12_utau' , 1 , 'utau' , .false. , .false. , 'daily' , '' , '' |
---|
70 | sn_vtau = 'amm12_vtau' , 1 , 'vtau' , .false. , .false. , 'daily' , '' , '' |
---|
71 | sn_qtot = 'amm12_flx' , 3 , 'sonsfldo' , .true. , .false. , 'daily' , '' , '' |
---|
72 | sn_qsr = 'amm12_flx' , 3 , 'soshfldo' , .true. , .false. , 'daily' , '' , '' |
---|
73 | sn_emp = 'amm12_flx' , 3 , 'sowafldo' , .true. , .false. , 'daily' , '' , '' |
---|
74 | cn_dir = './fluxes/' ! root directory for the location of the flux files |
---|
75 | / |
---|
76 | !----------------------------------------------------------------------- |
---|
77 | &namsbc_clio ! namsbc_clio CLIO bulk formulae |
---|
78 | !----------------------------------------------------------------------- |
---|
79 | / |
---|
80 | !----------------------------------------------------------------------- |
---|
81 | &namsbc_core ! namsbc_core CORE bulk formulae |
---|
82 | !----------------------------------------------------------------------- |
---|
83 | / |
---|
84 | !----------------------------------------------------------------------- |
---|
85 | &namsbc_mfs ! namsbc_mfs MFS bulk formulae |
---|
86 | !----------------------------------------------------------------------- |
---|
87 | / |
---|
88 | !----------------------------------------------------------------------- |
---|
89 | &namsbc_cpl ! coupled ocean/atmosphere model ("key_coupled") |
---|
90 | !----------------------------------------------------------------------- |
---|
91 | / |
---|
92 | !----------------------------------------------------------------------- |
---|
93 | &namtra_qsr ! penetrative solar radiation |
---|
94 | !----------------------------------------------------------------------- |
---|
95 | ln_traqsr = .false. ! Light penetration (T) or not (F) |
---|
96 | nn_chldta = 0 ! RGB : Chl data (=1) or cst value (=0) |
---|
97 | / |
---|
98 | !----------------------------------------------------------------------- |
---|
99 | &namsbc_rnf ! runoffs namelist surface boundary condition |
---|
100 | !----------------------------------------------------------------------- |
---|
101 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! |
---|
102 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! |
---|
103 | sn_rnf = 'amm12_rivers' , 24 , 'rorunoff', .false. , .true. , 'yearly' , '' , '' |
---|
104 | sn_cnf = 'runoff_1m_nomask' , 0 , 'socoefr0', .false. , .true. , 'yearly' , '' , '' |
---|
105 | sn_s_rnf = 'amm12_rivers' , 24 , 'rosaline', .false. , .true. , 'yearly' , '' , '' |
---|
106 | sn_t_rnf = 'amm12_rivers' , 24 , 'rotemper', .false. , .true. , 'yearly' , '' , '' |
---|
107 | sn_dep_rnf = 'amm12_rivers' , 24 , 'rodepth' , .false. , .true. , 'yearly' , '' , '' |
---|
108 | ln_rnf_mouth = .false. ! specific treatment at rivers mouths |
---|
109 | ln_rnf_depth = .true. ! read in depth information for runoff |
---|
110 | ln_rnf_tem = .true. ! read in temperature information for runoff |
---|
111 | ln_rnf_sal = .true. ! read in salinity information for runoff |
---|
112 | / |
---|
113 | !----------------------------------------------------------------------- |
---|
114 | &namsbc_apr ! Atmospheric pressure used as ocean forcing or in bulk |
---|
115 | !----------------------------------------------------------------------- |
---|
116 | ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly'/ ! weights ! rotation ! |
---|
117 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! |
---|
118 | sn_apr = 'patm' , -1 ,'somslpre', .true. , .true. , 'yearly' , '' , '' |
---|
119 | |
---|
120 | cn_dir = './' ! root directory for the location of the bulk files |
---|
121 | rn_pref = 101000._wp ! reference atmospheric pressure [N/m2]/ |
---|
122 | ln_ref_apr = .false. ! ref. pressure: global mean Patm (T) or a constant (F) |
---|
123 | ln_apr_obc = .false. ! inverse barometer added to OBC ssh data |
---|
124 | / |
---|
125 | !----------------------------------------------------------------------- |
---|
126 | &namsbc_ssr ! surface boundary condition : sea surface restoring |
---|
127 | !----------------------------------------------------------------------- |
---|
128 | ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! |
---|
129 | ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! |
---|
130 | sn_sst = 'amm12_sstref' , 24 , 'sst' , .true. , .false. , 'daily' , .false. , '' , '' |
---|
131 | sn_sss = 'sss_data' , -1 , 'sss' , .true. , .true. , 'yearly' , .false. , '' , '' |
---|
132 | |
---|
133 | cn_dir = 'fluxes/' ! root directory for the location of the runoff files |
---|
134 | nn_sstr = 1 ! add a retroaction term in the surface heat flux (=1) or not (=0) |
---|
135 | nn_sssr = 0 ! add a damping term in the surface freshwater flux (=2) |
---|
136 | ! or to SSS only (=1) or no damping term (=0) |
---|
137 | rn_deds = -27.7 ! magnitude of the damping on salinity [mm/day] |
---|
138 | / |
---|
139 | !----------------------------------------------------------------------- |
---|
140 | &namsbc_alb ! albedo parameters |
---|
141 | !----------------------------------------------------------------------- |
---|
142 | / |
---|
143 | !----------------------------------------------------------------------- |
---|
144 | &namberg ! iceberg parameters |
---|
145 | !----------------------------------------------------------------------- |
---|
146 | / |
---|
147 | !----------------------------------------------------------------------- |
---|
148 | &namlbc ! lateral momentum boundary condition |
---|
149 | !----------------------------------------------------------------------- |
---|
150 | rn_shlat = 0 ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat |
---|
151 | ! free slip ! partial slip ! no slip ! strong slip |
---|
152 | / |
---|
153 | !----------------------------------------------------------------------- |
---|
154 | &namcla ! cross land advection |
---|
155 | !----------------------------------------------------------------------- |
---|
156 | / |
---|
157 | !----------------------------------------------------------------------- |
---|
158 | &namobc ! open boundaries parameters ("key_obc") |
---|
159 | !----------------------------------------------------------------------- |
---|
160 | / |
---|
161 | !----------------------------------------------------------------------- |
---|
162 | &namagrif ! AGRIF zoom ("key_agrif") |
---|
163 | !----------------------------------------------------------------------- |
---|
164 | / |
---|
165 | !----------------------------------------------------------------------- |
---|
166 | &nam_tide ! tide parameters (#ifdef key_tide) |
---|
167 | !----------------------------------------------------------------------- |
---|
168 | ln_tide_pot = .false. ! use tidal potential forcing |
---|
169 | clname(1) = 'Q1' ! name of constituent |
---|
170 | clname(2) = 'O1' |
---|
171 | clname(3) = 'P1' |
---|
172 | clname(4) = 'S1' |
---|
173 | clname(5) = 'K1' |
---|
174 | clname(6) = '2N2' |
---|
175 | clname(7) = 'MU2' |
---|
176 | clname(8) = 'N2' |
---|
177 | clname(9) = 'NU2' |
---|
178 | clname(10) = 'M2' |
---|
179 | clname(11) = 'L2' |
---|
180 | clname(12) = 'T2' |
---|
181 | clname(13) = 'S2' |
---|
182 | clname(14) = 'K2' |
---|
183 | clname(15) = 'M4' |
---|
184 | / |
---|
185 | !----------------------------------------------------------------------- |
---|
186 | &nambdy ! unstructured open boundaries ("key_bdy") |
---|
187 | !----------------------------------------------------------------------- |
---|
188 | / |
---|
189 | !----------------------------------------------------------------------- |
---|
190 | &nambdy_dta ! open boundaries - external data ("key_bdy") |
---|
191 | !----------------------------------------------------------------------- |
---|
192 | / |
---|
193 | !----------------------------------------------------------------------- |
---|
194 | &nambdy_tide ! tidal forcing at open boundaries |
---|
195 | !----------------------------------------------------------------------- |
---|
196 | filtide = 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files |
---|
197 | / |
---|
198 | !----------------------------------------------------------------------- |
---|
199 | &nambfr ! bottom friction |
---|
200 | !----------------------------------------------------------------------- |
---|
201 | nn_bfr = 2 ! type of bottom friction : = 0 : free slip, = 1 : linear friction |
---|
202 | ! = 2 : nonlinear friction |
---|
203 | rn_bfri2 = 2.5e-3 ! bottom drag coefficient (non linear case) |
---|
204 | rn_bfeb2 = 0.0e0 ! bottom turbulent kinetic energy background (m2/s2) |
---|
205 | ln_loglayer = .false. ! loglayer bottom friction (only effect when nn_bfr = 2) |
---|
206 | rn_bfrz0 = 0.003 ! bottom roughness (only effect when ln_loglayer = .true.) |
---|
207 | / |
---|
208 | !----------------------------------------------------------------------- |
---|
209 | &nambbc ! bottom temperature boundary condition |
---|
210 | !----------------------------------------------------------------------- |
---|
211 | ln_trabbc = .false. ! Apply a geothermal heating at the ocean bottom |
---|
212 | / |
---|
213 | !----------------------------------------------------------------------- |
---|
214 | &nambbl ! bottom boundary layer scheme |
---|
215 | !----------------------------------------------------------------------- |
---|
216 | nn_bbl_ldf = 0 ! diffusive bbl (=1) or not (=0) |
---|
217 | / |
---|
218 | !----------------------------------------------------------------------- |
---|
219 | &nameos ! ocean physical parameters |
---|
220 | !----------------------------------------------------------------------- |
---|
221 | / |
---|
222 | !----------------------------------------------------------------------- |
---|
223 | &namtra_adv ! advection scheme for tracer |
---|
224 | !----------------------------------------------------------------------- |
---|
225 | / |
---|
226 | !---------------------------------------------------------------------------------- |
---|
227 | &namtra_ldf ! lateral diffusion scheme for tracers |
---|
228 | !---------------------------------------------------------------------------------- |
---|
229 | ln_traldf_hor = .true. ! horizontal (geopotential) (needs "key_ldfslp" when ln_sco=T) |
---|
230 | ln_traldf_iso = .false. ! iso-neutral (needs "key_ldfslp") |
---|
231 | rn_aeiv_0 = 0. ! eddy induced velocity coefficient [m2/s] |
---|
232 | rn_aht_0 = 50. ! horizontal eddy diffusivity for tracers [m2/s] |
---|
233 | / |
---|
234 | !----------------------------------------------------------------------- |
---|
235 | &namtra_dmp ! tracer: T & S newtonian damping |
---|
236 | !----------------------------------------------------------------------- |
---|
237 | ln_tradmp = .false. ! add a damping termn (T) or not (F) |
---|
238 | / |
---|
239 | !----------------------------------------------------------------------- |
---|
240 | &namdyn_adv ! formulation of the momentum advection |
---|
241 | !----------------------------------------------------------------------- |
---|
242 | / |
---|
243 | !----------------------------------------------------------------------- |
---|
244 | &namdyn_vor ! option of physics/algorithm (not control by CPP keys) |
---|
245 | !----------------------------------------------------------------------- |
---|
246 | ln_dynvor_ene = .false. ! energy conserving scheme |
---|
247 | ln_dynvor_ens = .false. ! enstrophy conserving scheme |
---|
248 | / |
---|
249 | !----------------------------------------------------------------------- |
---|
250 | &namdyn_hpg ! Hydrostatic pressure gradient option |
---|
251 | !----------------------------------------------------------------------- |
---|
252 | ln_hpg_zps = .false. ! z-coordinate - partial steps (interpolation) |
---|
253 | ln_hpg_sco = .false. ! s-coordinate (standard jacobian formulation) |
---|
254 | ln_hpg_prj = .true. ! s-coordinate (Pressure Jacobian scheme) |
---|
255 | / |
---|
256 | !----------------------------------------------------------------------- |
---|
257 | !namdyn_spg ! surface pressure gradient (CPP key only) |
---|
258 | !----------------------------------------------------------------------- |
---|
259 | ! ! explicit free surface ("key_dynspg_exp") |
---|
260 | ! ! filtered free surface ("key_dynspg_flt") |
---|
261 | ! ! split-explicit free surface ("key_dynspg_ts") |
---|
262 | |
---|
263 | !----------------------------------------------------------------------- |
---|
264 | &namdyn_ldf ! lateral diffusion on momentum |
---|
265 | !----------------------------------------------------------------------- |
---|
266 | ! ! Type of the operator : |
---|
267 | ln_dynldf_bilap = .true. ! bilaplacian operator |
---|
268 | ! Coefficient |
---|
269 | rn_ahm_0_lap = 60.0 ! horizontal laplacian eddy viscosity [m2/s] |
---|
270 | rn_ahmb_0 = 0.0 ! background eddy viscosity for ldf_iso [m2/s] |
---|
271 | rn_ahm_0_blp = -1.0e+10 ! horizontal bilaplacian eddy viscosity [m4/s] |
---|
272 | / |
---|
273 | !----------------------------------------------------------------------- |
---|
274 | &namzdf ! vertical physics |
---|
275 | !----------------------------------------------------------------------- |
---|
276 | rn_avm0 = 0.1e-6 ! vertical eddy viscosity [m2/s] (background Kz if not "key_zdfcst") |
---|
277 | rn_avt0 = 0.1e-6 ! vertical eddy diffusivity [m2/s] (background Kz if not "key_zdfcst") |
---|
278 | ln_zdfevd = .false. ! enhanced vertical diffusion (evd) (T) or not (F) |
---|
279 | nn_evdm = 1 ! evd apply on tracer (=0) or on tracer and momentum (=1) |
---|
280 | / |
---|
281 | !----------------------------------------------------------------------- |
---|
282 | &namzdf_ric ! richardson number dependent vertical diffusion ("key_zdfric" ) |
---|
283 | !----------------------------------------------------------------------- |
---|
284 | / |
---|
285 | !----------------------------------------------------------------------- |
---|
286 | &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") |
---|
287 | !----------------------------------------------------------------------- |
---|
288 | / |
---|
289 | !------------------------------------------------------------------------ |
---|
290 | &namzdf_kpp ! K-Profile Parameterization dependent vertical mixing ("key_zdfkpp", and optionally: |
---|
291 | !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") |
---|
292 | / |
---|
293 | !----------------------------------------------------------------------- |
---|
294 | &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") |
---|
295 | !----------------------------------------------------------------------- |
---|
296 | rn_charn = 100000. ! Charnock constant for wb induced roughness length |
---|
297 | / |
---|
298 | !----------------------------------------------------------------------- |
---|
299 | &namzdf_ddm ! double diffusive mixing parameterization ("key_zdfddm") |
---|
300 | !----------------------------------------------------------------------- |
---|
301 | / |
---|
302 | !----------------------------------------------------------------------- |
---|
303 | &namzdf_tmx ! tidal mixing parameterization ("key_zdftmx") |
---|
304 | !----------------------------------------------------------------------- |
---|
305 | ln_tmx_itf = .FALSE. ! ITF specific parameterisation |
---|
306 | / |
---|
307 | !----------------------------------------------------------------------- |
---|
308 | &namsol ! elliptic solver / island / free surface |
---|
309 | !----------------------------------------------------------------------- |
---|
310 | / |
---|
311 | !----------------------------------------------------------------------- |
---|
312 | &nammpp ! Massively Parallel Processing ("key_mpp_mpi) |
---|
313 | !----------------------------------------------------------------------- |
---|
314 | / |
---|
315 | !----------------------------------------------------------------------- |
---|
316 | &namctl ! Control prints & Benchmark |
---|
317 | !----------------------------------------------------------------------- |
---|
318 | / |
---|
319 | !----------------------------------------------------------------------- |
---|
320 | &namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4") |
---|
321 | !----------------------------------------------------------------------- |
---|
322 | / |
---|
323 | !----------------------------------------------------------------------- |
---|
324 | &namtrd ! diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") |
---|
325 | ! ! or mixed-layer trends or barotropic vorticity ("key_trdmld" or "key_trdvor") |
---|
326 | !----------------------------------------------------------------------- |
---|
327 | / |
---|
328 | !----------------------------------------------------------------------- |
---|
329 | &namflo ! float parameters ("key_float") |
---|
330 | !----------------------------------------------------------------------- |
---|
331 | / |
---|
332 | !----------------------------------------------------------------------- |
---|
333 | &namptr ! Poleward Transport Diagnostic |
---|
334 | !----------------------------------------------------------------------- |
---|
335 | ln_diaznl = .false. ! Add zonal means and meridional stream functions |
---|
336 | ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not |
---|
337 | ! (orca configuration only, need input basins mask file named "subbasins.nc" |
---|
338 | ln_ptrcomp = .false. ! Add decomposition : overturning |
---|
339 | / |
---|
340 | !----------------------------------------------------------------------- |
---|
341 | &namhsb ! Heat and salt budgets |
---|
342 | !----------------------------------------------------------------------- |
---|
343 | / |
---|
344 | !----------------------------------------------------------------------- |
---|
345 | &nam_diaharm ! Harmonic analysis of tidal constituents ('key_diaharm') |
---|
346 | !----------------------------------------------------------------------- |
---|
347 | / |
---|
348 | !----------------------------------------------------------------------- |
---|
349 | &namdct ! transports through sections |
---|
350 | !----------------------------------------------------------------------- |
---|
351 | / |
---|
352 | !----------------------------------------------------------------------- |
---|
353 | &namobs ! observation usage switch ('key_diaobs') |
---|
354 | !----------------------------------------------------------------------- |
---|
355 | / |
---|
356 | !----------------------------------------------------------------------- |
---|
357 | &nam_asminc ! assimilation increments ('key_asminc') |
---|
358 | !----------------------------------------------------------------------- |
---|
359 | / |
---|
360 | !----------------------------------------------------------------------- |
---|
361 | &namsbc_wave ! External fields from wave model |
---|
362 | !----------------------------------------------------------------------- |
---|
363 | / |
---|
364 | !----------------------------------------------------------------------- |
---|
365 | &namdyn_nept ! Neptune effect (simplified: lateral and vertical diffusions removed) |
---|
366 | !----------------------------------------------------------------------- |
---|
367 | / |
---|