Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! NEMO/OCE Configuration namelist : overwrite default values defined in SHARED/namelist_ref
!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
!! SPITZ 1/12 configuration !!
!!======================================================================
!! *** Domain & Run management namelists *** !!
!! !!
!! namrun parameters of the run
!! namdom space and time domain
!! namcfg parameters of the configuration (default: user defined GYRE)
!! namwad Wetting and drying (default: OFF)
!! namtsd data: temperature & salinity (default: OFF)
!! namcrs coarsened grid (for outputs and/or TOP) (ln_crs =T)
!! namc1d 1D configuration options (ln_c1d =T)
!! namc1d_dyndmp 1D newtonian damping applied on currents (ln_c1d =T)
!! namc1d_uvd 1D data (currents) (ln_c1d =T)
!!======================================================================
!
!-----------------------------------------------------------------------
&namrun ! parameters of the run
!-----------------------------------------------------------------------
cn_exp = "SPITZ" ! experience name
nn_it000 = 1 ! first time step
nn_itend = 100 ! last time step (std 5475)
/
!-----------------------------------------------------------------------
&namdom ! time and space domain
!-----------------------------------------------------------------------
rn_Dt = 720. ! time step for the dynamics and tracer
!
/
!-----------------------------------------------------------------------
&namcfg ! parameters of the configuration (default: use namusr_def in namelist_cfg)
!-----------------------------------------------------------------------
ln_read_cfg = .true. ! (=T) read the domain configuration file
! ! (=F) user defined configuration (F => create/check namusr_def)
cn_domcfg = "domain_cfg" ! domain configuration filename
/
!-----------------------------------------------------------------------
&namtile ! parameters of the tiling
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namtsd ! Temperature & Salinity Data (init/dmp) (default: OFF)
!-----------------------------------------------------------------------
! ! =T read T-S fields for:
ln_tsd_init = .true. ! ocean initialisation
cn_dir = './' ! root directory for the T-S data location
!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
sn_tem = 'T_SPITZ12' , 24. ,'votemper', .false. , .false., 'daily' , '' , '' , ''
sn_sal = 'S_SPITZ12' , 24. ,'vosaline', .false. , .false., 'daily' , '' , '' , ''
/
!!======================================================================
!! *** Surface Boundary Condition namelists *** !!
!! !!
!! namsbc surface boundary condition manager (default: NO selection)
!! namsbc_flx flux formulation (ln_flx =T)
!! namsbc_blk Bulk formulae formulation (ln_blk =T)
!! namsbc_cpl CouPLed formulation ("key_oasis3" )
!! namsbc_sas Stand-Alone Surface module (SAS_SRC only)
!! namsbc_iif Ice-IF: use observed ice cover (nn_ice = 1 )
!! namtra_qsr penetrative solar radiation (ln_traqsr =T)
!! namsbc_ssr sea surface restoring term (for T and/or S) (ln_ssr =T)
!! namsbc_rnf river runoffs (ln_rnf =T)
!! namsbc_apr Atmospheric Pressure (ln_apr_dyn =T)
!! namsbc_isf ice shelf melting/freezing (ln_isfcav =T : read (ln_read_cfg=T) or set or usr_def_zgr )
!! namsbc_iscpl coupling option between land ice model and ocean (ln_isfcav =T)
!! namsbc_wave external fields from wave model (ln_wave =T)
!! namberg iceberg floats (ln_icebergs=T)
!!======================================================================
!
!-----------------------------------------------------------------------
&namsbc ! Surface Boundary Condition manager (default: NO selection)
!-----------------------------------------------------------------------
nn_fsbc = 1 ! frequency of SBC module call
! ! (control sea-ice & iceberg model call)
! Type of air-sea fluxes
ln_usr = .false. ! user defined formulation (T => check usrdef_sbc)
ln_flx = .false. ! flux formulation (T => fill namsbc_flx )
ln_blk = .true. ! Bulk formulation (T => fill namsbc_blk )
! Sea-ice :
nn_ice = 2 ! SI3
ln_ice_embd = .false. ! =T embedded sea-ice (pressure + mass and salt exchanges)
! ! =F levitating ice (no pressure, mass and salt exchanges)
! Misc. options of sbc :
ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr )
/
!-----------------------------------------------------------------------
&namsbc_blk ! namsbc_blk generic Bulk formula (ln_blk =T)
!-----------------------------------------------------------------------
! ! bulk algorithm :
ln_NCAR = .true. ! "NCAR" algorithm (Large and Yeager 2008)
!
ln_Cx_ice_LU12 = .false. ! air-ice drags = F(ice concentration) (Lupkes et al. 2012)
ln_Cx_ice_LG15 = .true. ! air-ice drags = F(ice concentration) (Lupkes et al. 2015)
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
!
cn_dir = './' ! root directory for the bulk data location
!___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
sn_wndi = 'MARv3.6-9km-Svalbard-2hourly_spitz' , 2. , 'u10' , .true. , .false. , 'yearly' , 'weights_bicub', 'Uwnd' , ''
sn_wndj = 'MARv3.6-9km-Svalbard-2hourly_spitz' , 2. , 'v10' , .true. , .false. , 'yearly' , 'weights_bicub', 'Vwnd' , ''
sn_qsr = 'MARv3.6-9km-Svalbard-2hourly_spitz' , 2. , 'ssrd' , .true. , .false. , 'yearly' , 'weights_bilin', '' , ''
sn_qlw = 'MARv3.6-9km-Svalbard-2hourly_spitz' , 2. , 'strd' , .true. , .false. , 'yearly' , 'weights_bilin', '' , ''
sn_tair = 'MARv3.6-9km-Svalbard-2hourly_spitz' , 2. , 't10' , .true. , .false. , 'yearly' , 'weights_bilin', '' , ''
sn_humi = 'MARv3.6-9km-Svalbard-2hourly_spitz' , 2. , 'humi' , .true. , .false. , 'yearly' , 'weights_bilin', '' , ''
sn_prec = 'MARv3.6-9km-Svalbard-2hourly_spitz' , 2. , 'precip' , .true. , .false. , 'yearly' , 'weights_bilin', '' , ''
sn_snow = 'MARv3.6-9km-Svalbard-2hourly_spitz' , 2. , 'snow' , .true. , .false. , 'yearly' , 'weights_bilin', '' , ''
sn_slp = 'MARv3.6-9km-Svalbard-2hourly_spitz' , 2. , 'slp' , .true. , .false. , 'yearly' , 'weights_bilin', '' , ''
/
!-----------------------------------------------------------------------
&namsbc_sas ! Stand-Alone Surface module: ocean data (SAS_SRC only)
!-----------------------------------------------------------------------
l_sasread = .false. ! =T Read in file ; =F set all to 0. (see sbcssm)
/
!-----------------------------------------------------------------------
&namtra_qsr ! penetrative solar radiation (ln_traqsr =T)
!-----------------------------------------------------------------------
! ! type of penetration (default: NO selection)
ln_qsr_rgb = .true. ! RGB light penetration (Red-Green-Blue)
/
!!======================================================================
!! *** Lateral boundary condition *** !!
!! !!
!! namlbc lateral momentum boundary condition (default: NO selection)
!! namagrif agrif nested grid (read by child model only) ("key_agrif")
!! nam_tide Tidal forcing (default: OFF)
!! nambdy Unstructured open boundaries (default: OFF)
!! nambdy_dta Unstructured open boundaries - external data (see nambdy)
!! nambdy_tide tidal forcing at open boundaries (default: OFF)
!!======================================================================
!
!-----------------------------------------------------------------------
&namlbc ! lateral momentum boundary condition (default: NO selection)
!-----------------------------------------------------------------------
rn_shlat = 2. ! no slip
/
!-----------------------------------------------------------------------
&nam_tide ! tide parameters (default: OFF)
!-----------------------------------------------------------------------
ln_tide = .true. ! Activate tides
ln_tide_pot = .false. ! use tidal potential forcing
sn_tide_cnames(1) = 'M2' ! name of constituent
sn_tide_cnames(2) = 'S2'
sn_tide_cnames(3) = 'N2'
sn_tide_cnames(4) = 'K2'
sn_tide_cnames(5) = 'K1'
sn_tide_cnames(6) = 'O1'
sn_tide_cnames(7) = 'P1'
sn_tide_cnames(8) = 'Q1'
/
!-----------------------------------------------------------------------
&nambdy ! unstructured open boundaries (default: OFF)
!-----------------------------------------------------------------------
ln_bdy = .true. ! Use unstructured open boundaries
nb_bdy = 1 ! number of open boundary sets
ln_coords_file = .true. ! =T : read bdy coordinates from file
cn_coords_file = 'coordinates_bdyr1_grid12.nc' ! bdy coordinates files
cn_dyn2d = 'flather' !
nn_dyn2d_dta = 3 ! = 0, bdy data are equal to the initial state
! ! = 1, bdy data are read in 'bdydata .nc' files
! ! = 2, use tidal harmonic forcing data from files
! ! = 3, use external data AND tidal harmonic forcing
cn_dyn3d = 'frs' !
nn_dyn3d_dta = 1 ! = 0, bdy data are equal to the initial state
! ! = 1, bdy data are read in 'bdydata .nc' files
cn_tra = 'frs' !
nn_tra_dta = 1 ! = 0, bdy data are equal to the initial state
! ! = 1, bdy data are read in 'bdydata .nc' files
cn_ice = 'frs' !
nn_ice_dta = 1 ! = 0, bdy data are equal to the initial state
! ! = 1, bdy data are read in 'bdydata .nc' files
nn_rimwidth = 1 ! width of the relaxation zone
ln_vol = .false. ! total volume correction (see nn_volctl parameter)
/
!-----------------------------------------------------------------------
&nambdy_dta ! open boundaries - external data (see nam_bdy)
!-----------------------------------------------------------------------
ln_zinterp = .false. ! T if a vertical interpolation is required. Variables gdep[tuv] and e3[tuv] must exist in the file
! ! automatically defined to T if the number of vertical levels in bdy dta /= jpk
ln_full_vel = .false. ! T if [uv]3d are "full" velocities and not only its baroclinic components
! ! in this case, baroclinic and barotropic velocities will be recomputed -> [uv]2d not needed
!
cn_dir = './'
! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !
! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename !
bn_ssh = 'bdyT_u2d_grid12' , 72. , 'sossheig' , .true. , .false. , 'yearly' , '' , '' , ''
bn_u2d = 'bdyU_u2d_grid12' , 72. , 'vobtcrtx' , .true. , .false. , 'yearly' , '' , '' , ''
bn_v2d = 'bdyV_u2d_grid12' , 72. , 'vobtcrty' , .true. , .false. , 'yearly' , '' , '' , ''
bn_u3d = 'bdyU_u3d_grid12' , 72. , 'vozocrtx' , .true. , .false. , 'yearly' , '' , '' , ''
bn_v3d = 'bdyV_u3d_grid12' , 72. , 'vomecrty' , .true. , .false. , 'yearly' , '' , '' , ''
bn_tem = 'bdyT_tem_grid12' , 72. , 'votemper' , .true. , .false. , 'yearly' , '' , '' , ''
bn_sal = 'bdyT_sal_grid12' , 72. , 'vosaline' , .true. , .false. , 'yearly' , '' , '' , ''
bn_a_i = 'bdyT_ice_grid12' , 72. , 'ileadfra' , .true. , .false. , 'yearly' , '' , '' , ''
bn_h_i = 'bdyT_ice_grid12' , 72. , 'iicethic' , .true. , .false. , 'yearly' , '' , '' , ''
bn_h_s = 'bdyT_ice_grid12' , 72. , 'isnowthi' , .true. , .false. , 'yearly' , '' , '' , ''
/
!-----------------------------------------------------------------------
&nambdy_tide ! tidal forcing at open boundaries (default: OFF)
!-----------------------------------------------------------------------
filtide = 'bdytide_grid12_' ! file name root of tidal forcing files
/
!!======================================================================
!! *** Top/Bottom boundary condition *** !!
!! !!
!! namdrg top/bottom drag coefficient (default: NO selection)
!! namdrg_top top friction (ln_drg_OFF=F & ln_isfcav=T)
!! namdrg_bot bottom friction (ln_drg_OFF=F)
!! nambbc bottom temperature boundary condition (default: OFF)
!! nambbl bottom boundary layer scheme (default: OFF)
!!======================================================================
!
!-----------------------------------------------------------------------
&namdrg ! top/bottom drag coefficient (default: NO selection)
!-----------------------------------------------------------------------
ln_loglayer = .true. ! logarithmic drag: Cd = vkarmn/log(z/z0) |U|
ln_drgimp = .true. ! implicit top/bottom friction flag
ln_drgice_imp = .true. ! implicit ice-ocean drag
/
!-----------------------------------------------------------------------
&namdrg_bot ! BOTTOM friction (ln_drg_OFF =F)
!-----------------------------------------------------------------------
rn_Cd0 = 2.5e-3 ! drag coefficient [-]
rn_Cdmax = 0.1 ! drag value maximum [-] (logarithmic drag)
rn_ke0 = 0. ! background kinetic energy [m2/s2] (non-linear cases)
rn_z0 = 3.e-3 ! roughness [m] (ln_loglayer=T)
/
!-----------------------------------------------------------------------
&nambbl ! bottom boundary layer scheme (default: OFF)
!-----------------------------------------------------------------------
ln_trabbl = .true. ! Bottom Boundary Layer parameterisation flag
nn_bbl_ldf = 1 ! diffusive bbl (=1) or not (=0)
nn_bbl_adv = 0 ! advective bbl (=1/2) or not (=0)
/
!!======================================================================
!! Tracer (T-S) namelists !!
!! !!
!! nameos equation of state (default: NO selection)
!! namtra_adv advection scheme (default: NO selection)
!! namtra_ldf lateral diffusion scheme (default: NO selection)
!! namtra_mle mixed layer eddy param. (Fox-Kemper param.) (default: OFF)
!! namtra_eiv eddy induced velocity param. (default: OFF)
!! namtra_dmp T & S newtonian damping (default: OFF)
!!======================================================================
!
!-----------------------------------------------------------------------
&nameos ! ocean Equation Of Seawater (default: NO selection)
!-----------------------------------------------------------------------
ln_eos80 = .true. ! = Use EOS80
ln_teos10 = .false. ! = Use TEOS-10 equation of state
/
!-----------------------------------------------------------------------
&namtra_adv ! advection scheme for tracer (default: NO selection)
!-----------------------------------------------------------------------
ln_traadv_fct = .true. ! FCT scheme
nn_fct_h = 4 ! =2/4, horizontal 2nd / 4th order
nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order
/
!-----------------------------------------------------------------------
&namtra_ldf ! lateral diffusion scheme for tracers (default: NO selection)
!-----------------------------------------------------------------------
ln_traldf_blp = .true. ! bilaplacian operator
ln_traldf_iso = .true. ! iso-neutral (Standard operator)
ln_traldf_msc = .true. ! Method of Stabilizing Correction (both operators)
! ! Coefficients:
nn_aht_ijk_t = 31 ! space/time variation of eddy coefficient:
! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing)
/
!!======================================================================
!! *** Dynamics namelists *** !!
!! !!
!! nam_vvl vertical coordinate options (default: z-star)
!! namdyn_adv formulation of the momentum advection (default: NO selection)
!! namdyn_vor advection scheme (default: NO selection)
!! namdyn_hpg hydrostatic pressure gradient (default: NO selection)
!! namdyn_spg surface pressure gradient (default: NO selection)
!! namdyn_ldf lateral diffusion scheme (default: NO selection)
!! namdta_dyn offline TOP: dynamics read in files (OFF_SRC only)
!!======================================================================
!
!-----------------------------------------------------------------------
&nam_vvl ! vertical coordinate options (default: z-star)
!-----------------------------------------------------------------------
ln_vvl_zstar = .true. ! z-star vertical coordinate
/
!-----------------------------------------------------------------------
&namdyn_adv ! formulation of the momentum advection (default: NO selection)
!-----------------------------------------------------------------------
ln_dynadv_ubs = .true. ! flux form - 3rd order UBS scheme
/
!-----------------------------------------------------------------------
&namdyn_vor ! Vorticity / Coriolis scheme (default: NO selection)
!-----------------------------------------------------------------------
ln_dynvor_eeT = .true. ! energy conserving scheme (een using e3t)
/
!-----------------------------------------------------------------------
&namdyn_hpg ! Hydrostatic pressure gradient option (default: NO selection)
!-----------------------------------------------------------------------
ln_hpg_sco = .true. ! s-coordinate (standard jacobian formulation)
/
!-----------------------------------------------------------------------
&namdyn_spg ! surface pressure gradient (default: NO selection)
!-----------------------------------------------------------------------
ln_dynspg_ts = .true. ! split-explicit free surface
ln_bt_fw = .false. ! Forward integration of barotropic Eqs.
/
!-----------------------------------------------------------------------
&namdyn_ldf ! lateral diffusion on momentum (default: NO selection)
!-----------------------------------------------------------------------
ln_dynldf_OFF = .true. ! No operator (i.e. no explicit diffusion)
/
!!======================================================================
!! vertical physics namelists !!
!! !!
!! namzdf vertical physics manager (default: NO selection)
!! namzdf_ric richardson number vertical mixing (ln_zdfric=T)
!! namzdf_tke TKE vertical mixing (ln_zdftke=T)
!! namzdf_gls GLS vertical mixing (ln_zdfgls=T)
!! namzdf_osm OSM vertical diffusion (ln_zdfosm=T)
!! namzdf_iwm tidal mixing parameterization (ln_zdfiwm=T)
!!======================================================================
!
!-----------------------------------------------------------------------
&namzdf ! vertical physics manager (default: NO selection)
!-----------------------------------------------------------------------
! ! adaptive-implicit vertical advection
ln_zad_Aimp = .true. ! Courant number dependent scheme (Shchepetkin 2015)
! ! type of vertical closure (required)
ln_zdftke = .true. ! Turbulent Kinetic Energy closure (T => fill namzdf_tke)
ln_zdfgls = .false. ! Generic Length Scale closure (T => fill namzdf_gls)
! ! convection
ln_zdfevd = .true. ! enhanced vertical diffusion
!
ln_zdfddm = .true. ! double diffusive mixing
!
! ! Coefficients
rn_avm0 = 1.2e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F)
rn_avt0 = 2.0e-6 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F)
nn_avb = 0 ! profile for background avt & avm (=1) or not (=0)
nn_havtb = 1 ! horizontal shape for avtb (=1) or not (=0)
/
!-----------------------------------------------------------------------
&namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion (ln_zdftke =T)
!-----------------------------------------------------------------------
ln_mxl0 = .true. ! surface mixing length scale = F(wind stress) (T) or not (F)
nn_mxlice = 0 ! type of scaling under sea-ice
! = 0 no scaling under sea-ice
! = 1 scaling with constant sea-ice thickness
! = 2 scaling with mean sea-ice thickness ( only with SI3 sea-ice model )
! = 3 scaling with maximum sea-ice thickness
nn_eice = 0 ! attenutaion of langmuir & surface wave breaking under ice
! ! = 0 no impact of ice cover on langmuir & surface wave breaking
! ! = 1 weigthed by 1-TANH(10*fr_i)
! ! = 2 weighted by 1-fr_i
! ! = 3 weighted by 1-MIN(1,4*fr_i)
/
!!======================================================================
!! *** Diagnostics namelists *** !!
!! !!
!! namtrd dynamics and/or tracer trends (default: OFF)
!! namhsb Heat and salt budgets (default: OFF)
!! namdiu Cool skin and warm layer models (default: OFF)
!! namdiu Cool skin and warm layer models (default: OFF)
!! namflo float parameters (default: OFF)
!! nam_diadct transports through some sections (default: OFF)
!! nam_dia25h 25h Mean Output (default: OFF)
!! namnc4 netcdf4 chunking and compression settings ("key_netcdf4")
!!======================================================================
!
!-----------------------------------------------------------------------
&namhsb ! Heat and salt budgets (default: OFF)
!-----------------------------------------------------------------------
ln_diahsb = .false. ! check the heat and salt budgets (T) or not (F)
/
!-----------------------------------------------------------------------
&namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4")
!-----------------------------------------------------------------------
nn_nchunks_k = 75 ! number of chunks in k-dimension
! ! setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which
! ! is optimal for postprocessing which works exclusively with horizontal slabs
/
!!======================================================================
!! *** Observation & Assimilation *** !!
!! !!
!! namobs observation and model comparison (default: OFF)
!! nam_asminc assimilation increments ('key_asminc')
!!======================================================================
!
!!======================================================================
!! *** Miscellaneous namelists *** !!
!! !!
!! nammpp Massively Parallel Processing
!! namctl Control prints (default: OFF)
!! namsto Stochastic parametrization of EOS (default: OFF)
!!======================================================================
!
!-----------------------------------------------------------------------
&nammpp ! Massively Parallel Processing
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namctl ! Control prints (default: OFF)
!-----------------------------------------------------------------------
/
!-----------------------------------------------------------------------
&namsto ! Stochastic parametrization of EOS (default: OFF)
!-----------------------------------------------------------------------
/