New URL for NEMO forge! http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.

namelist_ref in NEMO/branches/UKMO/NEMO_4.0.1_GO8_package/cfgs/SHARED – NEMO

Context Navigation

source: NEMO/branches/UKMO/NEMO_4.0.1_GO8_package/cfgs/SHARED/namelist_ref @ 12088

Last change on this file since 12088 was 12088, checked in by deazer, 4 years ago
Updating GO8 Package branch to bring in required BDY bug fixes frouse with CO8 The mirror branch is already updated to have this change, where we merge in the mirror to the package branch
Property svn:mime-type set to `text/x-fortran`
File size: 104.8 KB

Line
1	!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
2	!! NEMO/OCE : Reference namelist_ref !!
3	!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
4	!! NEMO/OCE : 1 - Domain & run manager (namrun, namcfg, namdom, namtsd, namcrs, namc1d, namc1d_uvd)
5	!! namelists 2 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_cpl,
6	!! namsbc_sas, namtra_qsr, namsbc_rnf,
7	!! namsbc_isf, namsbc_iscpl, namsbc_apr,
8	!! namsbc_ssr, namsbc_wave, namberg)
9	!! 3 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide)
10	!! 4 - top/bot boundary (namdrg, namdrg_top, namdrg_bot, nambbc, nambbl)
11	!! 5 - Tracer (nameos, namtra_adv, namtra_ldf, namtra_eiv, namtra_dmp)
12	!! 6 - dynamics (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf)
13	!! 7 - Vertical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_gls, namzdf_iwm)
14	!! 8 - diagnostics (namnc4, namtrd, namspr, namflo, namhsb)
15	!! 9 - Obs & Assim (namobs, nam_asminc)
16	!! 10 - miscellaneous (nammpp, namctl, namsto)
17	!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
18
19	!!======================================================================
20	!! * Domain & Run management namelists * !!
21	!! !!
22	!! namrun parameters of the run
23	!! namdom space and time domain
24	!! namcfg parameters of the configuration (default: user defined GYRE)
25	!! namwad Wetting and drying (default: OFF)
26	!! namtsd data: temperature & salinity (default: OFF)
27	!! namcrs coarsened grid (for outputs and/or TOP) (ln_crs =T)
28	!! namc1d 1D configuration options ("key_c1d")
29	!! namc1d_dyndmp 1D newtonian damping applied on currents ("key_c1d")
30	!! namc1d_uvd 1D data (currents) ("key_c1d")
31	!!======================================================================
32	!
33	!-----------------------------------------------------------------------
34	&namrun ! parameters of the run
35	!-----------------------------------------------------------------------
36	nn_no = 0 ! Assimilation cycle index
37	cn_exp = "ORCA2" ! experience name
38	nn_it000 = 1 ! first time step
39	nn_itend = 5840 ! last time step (std 5840)
40	nn_date0 = 010101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1)
41	nn_time0 = 0 ! initial time of day in hhmm
42	nn_leapy = 0 ! Leap year calendar (1) or not (0)
43	ln_rstart = .false. ! start from rest (F) or from a restart file (T)
44	nn_euler = 1 ! = 0 : start with forward time step if ln_rstart=T
45	nn_rstctl = 0 ! restart control ==> activated only if ln_rstart=T
46	! ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist
47	! ! = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart
48	! ! = 2 nn_date0 read in restart ; nn_it000 : check consistancy between namelist and restart
49	cn_ocerst_in = "restart" ! suffix of ocean restart name (input)
50	cn_ocerst_indir = "." ! directory from which to read input ocean restarts
51	cn_ocerst_out = "restart" ! suffix of ocean restart name (output)
52	cn_ocerst_outdir = "." ! directory in which to write output ocean restarts
53	ln_iscpl = .false. ! cavity evolution forcing or coupling to ice sheet model
54	nn_istate = 0 ! output the initial state (1) or not (0)
55	ln_rst_list = .false. ! output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F)
56	nn_stock = 0 ! used only if ln_rst_list = F: output restart freqeuncy (modulo referenced to 1)
57	! ! = 0 force to write restart files only at the end of the run
58	! ! = -1 do not do any restart
59	nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written
60	nn_write = 0 ! used only if key_iomput is not defined: output frequency (modulo referenced to nn_it000)
61	! ! = 0 force to write output files only at the end of the run
62	! ! = -1 do not do any output file
63	ln_mskland = .false. ! mask land points in NetCDF outputs
64	ln_cfmeta = .false. ! output additional data to netCDF files required for compliance with the CF metadata standard
65	ln_clobber = .true. ! clobber (overwrite) an existing file
66	nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (works only with iom_nf90 routines)
67	ln_xios_read = .FALSE. ! use XIOS to read restart file (only for a single file restart)
68	nn_wxios = 0 ! use XIOS to write restart file 0 - no, 1 - single file output, 2 - multiple file output
69	/
70	!-----------------------------------------------------------------------
71	&namdom ! time and space domain
72	!-----------------------------------------------------------------------
73	ln_linssh = .false. ! =T linear free surface ==>> model level are fixed in time
74	rn_isfhmin = 1.00 ! treshold [m] to discriminate grounding ice from floating ice
75	!
76	rn_rdt = 5400. ! time step for the dynamics and tracer
77	rn_atfp = 0.1 ! asselin time filter parameter
78	!
79	ln_crs = .false. ! Logical switch for coarsening module (T => fill namcrs)
80	!
81	ln_meshmask = .false. ! =T create a mesh file
82	/
83	!-----------------------------------------------------------------------
84	&namcfg ! parameters of the configuration (default: use namusr_def in namelist_cfg)
85	!-----------------------------------------------------------------------
86	ln_read_cfg = .false. ! (=T) read the domain configuration file
87	! ! (=F) user defined configuration (F => create/check namusr_def)
88	cn_domcfg = "domain_cfg" ! domain configuration filename
89	!
90	ln_closea = .false. ! T => keep closed seas (defined by closea_mask field) in the
91	! ! domain and apply special treatment of freshwater fluxes.
92	! ! F => suppress closed seas (defined by closea_mask field)
93	! ! from the bathymetry at runtime.
94	! ! If closea_mask field doesn't exist in the domain_cfg file
95	! ! then this logical does nothing.
96	ln_write_cfg = .false. ! (=T) create the domain configuration file
97	cn_domcfg_out = "domain_cfg_out" ! newly created domain configuration filename
98	!
99	ln_use_jattr = .false. ! use (T) the file attribute: open_ocean_jstart, if present
100	! ! in netcdf input files, as the start j-row for reading
101	/
102	!-----------------------------------------------------------------------
103	&namtsd ! Temperature & Salinity Data (init/dmp) (default: OFF)
104	!-----------------------------------------------------------------------
105	! ! =T read T-S fields for:
106	ln_tsd_init = .false. ! ocean initialisation
107	ln_tsd_dmp = .false. ! T-S restoring (see namtra_dmp)
108
109	cn_dir = './' ! root directory for the T-S data location
110	!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
111	! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
112	! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
113	sn_tem = 'data_1m_potential_temperature_nomask', -1. , 'votemper', .true. , .true. , 'yearly' , '' , '' , ''
114	sn_sal = 'data_1m_salinity_nomask' , -1. , 'vosaline', .true. , .true. , 'yearly' , '' , '' , ''
115	/
116	!-----------------------------------------------------------------------
117	&namwad ! Wetting and Drying (WaD) (default: OFF)
118	!-----------------------------------------------------------------------
119	ln_wd_il = .false. ! T/F activation of iterative limiter
120	ln_wd_dl = .false. ! T/F activation of directional limiter
121	ln_wd_dl_bc = .false. ! T/F Directional limiteer Baroclinic option
122	ln_wd_dl_rmp = .false. ! T/F Turn on directional limiter ramp
123	rn_wdmin0 = 0.30 ! depth at which WaD starts
124	rn_wdmin1 = 0.2 ! Minimum wet depth on dried cells
125	rn_wdmin2 = 0.0001 ! Tolerance of min wet depth on dried cells
126	rn_wdld = 2.5 ! Land elevation below which WaD is allowed
127	nn_wdit = 20 ! Max iterations for WaD limiter
128	rn_wd_sbcdep = 5.0 ! Depth at which to taper sbc fluxes
129	rn_wd_sbcfra = 0.999 ! Fraction of SBC fluxes at taper depth (Must be <1)
130	/
131	!-----------------------------------------------------------------------
132	&namcrs ! coarsened grid (for outputs and/or TOP) (ln_crs =T)
133	!-----------------------------------------------------------------------
134	nn_factx = 3 ! Reduction factor of x-direction
135	nn_facty = 3 ! Reduction factor of y-direction
136	nn_binref = 0 ! Bin centering preference: NORTH or EQUAT
137	! ! 0, coarse grid is binned with preferential treatment of the north fold
138	! ! 1, coarse grid is binned with centering at the equator
139	! ! Symmetry with nn_facty being odd-numbered. Asymmetry with even-numbered nn_facty.
140	ln_msh_crs = .false. ! =T create a mesh & mask file
141	nn_crs_kz = 0 ! 0, MEAN of volume boxes
142	! ! 1, MAX of boxes
143	! ! 2, MIN of boxes
144	ln_crs_wn = .true. ! wn coarsened (T) or computed using horizontal divergence ( F )
145	/
146	!-----------------------------------------------------------------------
147	&namc1d ! 1D configuration options ("key_c1d" default: PAPA station)
148	!-----------------------------------------------------------------------
149	rn_lat1d = 50 ! Column latitude
150	rn_lon1d = -145 ! Column longitude
151	ln_c1d_locpt = .true. ! Localization of 1D config in a grid (T) or independant point (F)
152	/
153	!-----------------------------------------------------------------------
154	&namc1d_dyndmp ! U & V newtonian damping ("key_c1d" default: OFF)
155	!-----------------------------------------------------------------------
156	ln_dyndmp = .false. ! add a damping term (T) or not (F)
157	/
158	!-----------------------------------------------------------------------
159	&namc1d_uvd ! data: U & V currents ("key_c1d" default: OFF)
160	!-----------------------------------------------------------------------
161	! ! =T read U-V fields for:
162	ln_uvd_init = .false. ! ocean initialisation
163	ln_uvd_dyndmp = .false. ! U-V restoring
164
165	cn_dir = './' ! root directory for the U-V data location
166	!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
167	! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
168	! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
169	sn_ucur = 'ucurrent' , -1. ,'u_current', .false. , .true. , 'monthly' , '' , 'Ume' , ''
170	sn_vcur = 'vcurrent' , -1. ,'v_current', .false. , .true. , 'monthly' , '' , 'Vme' , ''
171	/
172
173	!!======================================================================
174	!! * Surface Boundary Condition namelists * !!
175	!! !!
176	!! namsbc surface boundary condition manager (default: NO selection)
177	!! namsbc_flx flux formulation (ln_flx =T)
178	!! namsbc_blk Bulk formulae formulation (ln_blk =T)
179	!! namsbc_cpl CouPLed formulation ("key_oasis3" )
180	!! namsbc_sas Stand-Alone Surface module (SAS_SRC only)
181	!! namsbc_iif Ice-IF: use observed ice cover (nn_ice = 1 )
182	!! namtra_qsr penetrative solar radiation (ln_traqsr =T)
183	!! namsbc_ssr sea surface restoring term (for T and/or S) (ln_ssr =T)
184	!! namsbc_rnf river runoffs (ln_rnf =T)
185	!! namsbc_apr Atmospheric Pressure (ln_apr_dyn =T)
186	!! namsbc_isf ice shelf melting/freezing (ln_isfcav =T : read (ln_read_cfg=T) or set or usr_def_zgr )
187	!! namsbc_iscpl coupling option between land ice model and ocean (ln_isfcav =T)
188	!! namsbc_wave external fields from wave model (ln_wave =T)
189	!! namberg iceberg floats (ln_icebergs=T)
190	!!======================================================================
191	!
192	!-----------------------------------------------------------------------
193	&namsbc ! Surface Boundary Condition manager (default: NO selection)
194	!-----------------------------------------------------------------------
195	nn_fsbc = 2 ! frequency of SBC module call
196	! ! (control sea-ice & iceberg model call)
197	! Type of air-sea fluxes
198	ln_usr = .false. ! user defined formulation (T => check usrdef_sbc)
199	ln_flx = .false. ! flux formulation (T => fill namsbc_flx )
200	ln_blk = .false. ! Bulk formulation (T => fill namsbc_blk )
201	! ! Type of coupling (Ocean/Ice/Atmosphere) :
202	ln_cpl = .false. ! atmosphere coupled formulation ( requires key_oasis3 )
203	ln_mixcpl = .false. ! forced-coupled mixed formulation ( requires key_oasis3 )
204	nn_components = 0 ! configuration of the opa-sas OASIS coupling
205	! ! =0 no opa-sas OASIS coupling: default single executable config.
206	! ! =1 opa-sas OASIS coupling: multi executable config., OPA component
207	! ! =2 opa-sas OASIS coupling: multi executable config., SAS component
208	! Sea-ice :
209	nn_ice = 0 ! =0 no ice boundary condition
210	! ! =1 use observed ice-cover ( => fill namsbc_iif )
211	! ! =2 or 3 automatically for SI3 or CICE ("key_si3" or "key_cice")
212	! ! except in AGRIF zoom where it has to be specified
213	ln_ice_embd = .false. ! =T embedded sea-ice (pressure + mass and salt exchanges)
214	! ! =F levitating ice (no pressure, mass and salt exchanges)
215	! Misc. options of sbc :
216	ln_traqsr = .false. ! Light penetration in the ocean (T => fill namtra_qsr)
217	ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave
218	ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr)
219	nn_fwb = 0 ! FreshWater Budget: =0 unchecked
220	! ! =1 global mean of e-p-r set to zero at each time step
221	! ! =2 annual global mean of e-p-r set to zero
222	ln_rnf = .false. ! runoffs (T => fill namsbc_rnf)
223	ln_apr_dyn = .false. ! Patm gradient added in ocean & ice Eqs. (T => fill namsbc_apr )
224	ln_isf = .false. ! ice shelf (T => fill namsbc_isf & namsbc_iscpl)
225	ln_wave = .false. ! Activate coupling with wave (T => fill namsbc_wave)
226	ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave)
227	ln_sdw = .false. ! Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave)
228	nn_sdrift = 0 ! Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift
229	! ! = 0 Breivik 2015 parameterization: v_z=v_0[exp(2kz)/(1-8k*z)]
230	! ! = 1 Phillips: v_z=v_o[exp(2kz)-betasqrt(-2kpiz)erfc(sqrt(-2kz))]
231	! ! = 2 Phillips as (1) but using the wave frequency from a wave model
232	ln_tauwoc = .false. ! Activate ocean stress modified by external wave induced stress (T => ln_wave=.true. & fill namsbc_wave)
233	ln_tauw = .false. ! Activate ocean stress components from wave model
234	ln_stcor = .false. ! Activate Stokes Coriolis term (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave)
235	nn_lsm = 0 ! =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) ,
236	! =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field)
237	/
238	!-----------------------------------------------------------------------
239	&namsbc_flx ! surface boundary condition : flux formulation (ln_flx =T)
240	!-----------------------------------------------------------------------
241	cn_dir = './' ! root directory for the fluxes data location
242	!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
243	! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
244	! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
245	sn_utau = 'utau' , 24. , 'utau' , .false. , .false., 'yearly' , '' , '' , ''
246	sn_vtau = 'vtau' , 24. , 'vtau' , .false. , .false., 'yearly' , '' , '' , ''
247	sn_qtot = 'qtot' , 24. , 'qtot' , .false. , .false., 'yearly' , '' , '' , ''
248	sn_qsr = 'qsr' , 24. , 'qsr' , .false. , .false., 'yearly' , '' , '' , ''
249	sn_emp = 'emp' , 24. , 'emp' , .false. , .false., 'yearly' , '' , '' , ''
250	/
251	!-----------------------------------------------------------------------
252	&namsbc_blk ! namsbc_blk generic Bulk formula (ln_blk =T)
253	!-----------------------------------------------------------------------
254	! ! bulk algorithm :
255	ln_NCAR = .false. ! "NCAR" algorithm (Large and Yeager 2008)
256	ln_COARE_3p0 = .false. ! "COARE 3.0" algorithm (Fairall et al. 2003)
257	ln_COARE_3p5 = .false. ! "COARE 3.5" algorithm (Edson et al. 2013)
258	ln_ECMWF = .false. ! "ECMWF" algorithm (IFS cycle 31)
259	!
260	rn_zqt = 10. ! Air temperature & humidity reference height (m)
261	rn_zu = 10. ! Wind vector reference height (m)
262	ln_Cd_L12 = .false. ! air-ice drags = F(ice concentration) (Lupkes et al. 2012)
263	ln_Cd_L15 = .false. ! air-ice drags = F(ice concentration) (Lupkes et al. 2015)
264	ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data
265	rn_pfac = 1. ! multiplicative factor for precipitation (total & snow)
266	rn_efac = 1. ! multiplicative factor for evaporation (0. or 1.)
267	rn_vfac = 0. ! multiplicative factor for ocean & ice velocity used to
268	! ! calculate the wind stress (0.=absolute or 1.=relative winds)
269
270	cn_dir = './' ! root directory for the bulk data location
271	!___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________!
272	! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
273	! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
274	sn_wndi = 'u_10.15JUNE2009_fill' , 6. , 'U_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bicubic_noc.nc' , 'Uwnd' , ''
275	sn_wndj = 'v_10.15JUNE2009_fill' , 6. , 'V_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bicubic_noc.nc' , 'Vwnd' , ''
276	sn_qsr = 'ncar_rad.15JUNE2009_fill' , 24. , 'SWDN_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''
277	sn_qlw = 'ncar_rad.15JUNE2009_fill' , 24. , 'LWDN_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''
278	sn_tair = 't_10.15JUNE2009_fill' , 6. , 'T_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''
279	sn_humi = 'q_10.15JUNE2009_fill' , 6. , 'Q_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''
280	sn_prec = 'ncar_precip.15JUNE2009_fill', -1. , 'PRC_MOD1', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''
281	sn_snow = 'ncar_precip.15JUNE2009_fill', -1. , 'SNOW' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''
282	sn_slp = 'slp.15JUNE2009_fill' , 6. , 'SLP' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''
283	sn_tdif = 'taudif_core' , 24 , 'taudif' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''
284	/
285	!-----------------------------------------------------------------------
286	&namsbc_cpl ! coupled ocean/atmosphere model ("key_oasis3")
287	!-----------------------------------------------------------------------
288	nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentially sending/receiving data
289	ln_usecplmask = .false. ! use a coupling mask file to merge data received from several models
290	! ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel)
291	nn_cats_cpl = 5 ! Number of sea ice categories over which coupling is to be carried out (if not 1)
292	!_____________!__________________________!____________!_____________!______________________!________!
293	! ! description ! multiple ! vector ! vector ! vector !
294	! ! ! categories ! reference ! orientation ! grids !
295	!* send *
296	sn_snd_temp = 'weighted oce and ice' , 'no' , '' , '' , ''
297	sn_snd_alb = 'weighted ice' , 'no' , '' , '' , ''
298	sn_snd_thick = 'none' , 'no' , '' , '' , ''
299	sn_snd_crt = 'none' , 'no' , 'spherical' , 'eastward-northward' , 'T'
300	sn_snd_co2 = 'coupled' , 'no' , '' , '' , ''
301	sn_snd_crtw = 'none' , 'no' , '' , '' , 'U,V'
302	sn_snd_ifrac = 'none' , 'no' , '' , '' , ''
303	sn_snd_wlev = 'coupled' , 'no' , '' , '' , ''
304	sn_snd_cond = 'weighted ice' , 'no' , '' , '' , ''
305	sn_snd_thick1 = 'ice and snow' , 'no' , '' , '' , ''
306	sn_snd_mpnd = 'weighted ice' , 'no' , '' , '' , ''
307	sn_snd_sstfrz = 'coupled' , 'no' , '' , '' , ''
308	sn_snd_ttilyr = 'weighted ice' , 'no' , '' , '' , ''
309	!* receive *
310	sn_rcv_w10m = 'none' , 'no' , '' , '' , ''
311	sn_rcv_taumod = 'coupled' , 'no' , '' , '' , ''
312	sn_rcv_tau = 'oce only' , 'no' , 'cartesian' , 'eastward-northward' , 'U,V'
313	sn_rcv_dqnsdt = 'coupled' , 'no' , '' , '' , ''
314	sn_rcv_qsr = 'oce and ice' , 'no' , '' , '' , ''
315	sn_rcv_qns = 'oce and ice' , 'no' , '' , '' , ''
316	sn_rcv_emp = 'conservative' , 'no' , '' , '' , ''
317	sn_rcv_rnf = 'coupled' , 'no' , '' , '' , ''
318	sn_rcv_cal = 'coupled' , 'no' , '' , '' , ''
319	sn_rcv_co2 = 'coupled' , 'no' , '' , '' , ''
320	sn_rcv_hsig = 'none' , 'no' , '' , '' , ''
321	sn_rcv_iceflx = 'none' , 'no' , '' , '' , ''
322	sn_rcv_mslp = 'none' , 'no' , '' , '' , ''
323	sn_rcv_phioc = 'none' , 'no' , '' , '' , ''
324	sn_rcv_sdrfx = 'none' , 'no' , '' , '' , ''
325	sn_rcv_sdrfy = 'none' , 'no' , '' , '' , ''
326	sn_rcv_wper = 'none' , 'no' , '' , '' , ''
327	sn_rcv_wnum = 'none' , 'no' , '' , '' , ''
328	sn_rcv_wfreq = 'none' , 'no' , '' , '' , ''
329	sn_rcv_wdrag = 'none' , 'no' , '' , '' , ''
330	sn_rcv_ts_ice = 'none' , 'no' , '' , '' , ''
331	sn_rcv_isf = 'none' , 'no' , '' , '' , ''
332	sn_rcv_icb = 'none' , 'no' , '' , '' , ''
333	sn_rcv_tauwoc = 'none' , 'no' , '' , '' , ''
334	sn_rcv_tauw = 'none' , 'no' , '' , '' , ''
335	sn_rcv_wdrag = 'none' , 'no' , '' , '' , ''
336	/
337	!-----------------------------------------------------------------------
338	&namsbc_sas ! Stand-Alone Surface module: ocean data (SAS_SRC only)
339	!-----------------------------------------------------------------------
340	l_sasread = .true. ! =T Read in file ; =F set all to 0. (see sbcssm)
341	ln_3d_uve = .false. ! specify whether we are supplying a 3D u,v and e3 field
342	ln_read_frq = .false. ! specify whether we must read frq or not
343
344	cn_dir = './' ! root directory for the ocean data location
345	!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
346	! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
347	! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
348	sn_usp = 'sas_grid_U' , 120. , 'uos' , .true. , .true. , 'yearly' , '' , '' , ''
349	sn_vsp = 'sas_grid_V' , 120. , 'vos' , .true. , .true. , 'yearly' , '' , '' , ''
350	sn_tem = 'sas_grid_T' , 120. , 'sosstsst', .true. , .true. , 'yearly' , '' , '' , ''
351	sn_sal = 'sas_grid_T' , 120. , 'sosaline', .true. , .true. , 'yearly' , '' , '' , ''
352	sn_ssh = 'sas_grid_T' , 120. , 'sossheig', .true. , .true. , 'yearly' , '' , '' , ''
353	sn_e3t = 'sas_grid_T' , 120. , 'e3t_m' , .true. , .true. , 'yearly' , '' , '' , ''
354	sn_frq = 'sas_grid_T' , 120. , 'frq_m' , .true. , .true. , 'yearly' , '' , '' , ''
355	/
356	!-----------------------------------------------------------------------
357	&namsbc_iif ! Ice-IF : use observed ice cover (nn_ice = 1)
358	!-----------------------------------------------------------------------
359	cn_dir = './' ! root directory for the ice cover data location
360	!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
361	! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
362	! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
363	sn_ice ='ice_cover_clim.nc' , -12. ,'ice_cover', .true. , .true. , 'yearly' , '' , '' , ''
364	/
365	!-----------------------------------------------------------------------
366	&namtra_qsr ! penetrative solar radiation (ln_traqsr =T)
367	!-----------------------------------------------------------------------
368	! ! type of penetration (default: NO selection)
369	ln_qsr_rgb = .false. ! RGB light penetration (Red-Green-Blue)
370	ln_qsr_2bd = .false. ! 2BD light penetration (two bands)
371	ln_qsr_bio = .false. ! bio-model light penetration
372	! ! RGB & 2BD choices:
373	rn_abs = 0.58 ! RGB & 2BD: fraction absorbed in the very near surface
374	rn_si0 = 0.35 ! RGB & 2BD: shortess depth of extinction
375	nn_chldta = 0 ! RGB : Chl data (=1) or cst value (=0)
376	rn_si1 = 23.0 ! 2BD : longest depth of extinction
377
378	cn_dir = './' ! root directory for the chlorophyl data location
379	!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
380	! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
381	! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
382	sn_chl ='chlorophyll' , -1. , 'CHLA' , .true. , .true. , 'yearly' , '' , '' , ''
383	/
384	!-----------------------------------------------------------------------
385	&namsbc_ssr ! surface boundary condition : sea surface restoring (ln_ssr =T)
386	!-----------------------------------------------------------------------
387	nn_sstr = 0 ! add a retroaction term to the surface heat flux (=1) or not (=0)
388	rn_dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K]
389	nn_sssr = 0 ! add a damping term to the surface freshwater flux (=2)
390	! ! or to SSS only (=1) or no damping term (=0)
391	rn_deds = -166.67 ! magnitude of the damping on salinity [mm/day]
392	ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2)
393	rn_sssr_bnd = 4.e0 ! ABS(Max/Min) value of the damping erp term [mm/day]
394
395	cn_dir = './' ! root directory for the SST/SSS data location
396	!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
397	! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
398	! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
399	sn_sst = 'sst_data' , 24. , 'sst' , .false. , .false., 'yearly' , '' , '' , ''
400	sn_sss = 'sss_data' , -1. , 'sss' , .true. , .true. , 'yearly' , '' , '' , ''
401	/
402	!-----------------------------------------------------------------------
403	&namsbc_rnf ! runoffs (ln_rnf =T)
404	!-----------------------------------------------------------------------
405	ln_rnf_mouth = .false. ! specific treatment at rivers mouths
406	rn_hrnf = 15.e0 ! depth over which enhanced vertical mixing is used (ln_rnf_mouth=T)
407	rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T)
408	rn_rfact = 1.e0 ! multiplicative factor for runoff
409	ln_rnf_depth = .false. ! read in depth information for runoff
410	ln_rnf_tem = .false. ! read in temperature information for runoff
411	ln_rnf_sal = .false. ! read in salinity information for runoff
412	ln_rnf_depth_ini = .false. ! compute depth at initialisation from runoff file
413	rn_rnf_max = 5.735e-4 ! max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true )
414	rn_dep_max = 150. ! depth over which runoffs is spread ( ln_rnf_depth_ini = .true )
415	nn_rnf_depth_file = 0 ! create (=1) a runoff depth file or not (=0)
416
417	cn_dir = './' ! root directory for the runoff data location
418	!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
419	! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
420	! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
421	sn_rnf = 'runoff_core_monthly' , -1. , 'sorunoff', .true. , .true. , 'yearly' , '' , '' , ''
422	sn_cnf = 'runoff_core_monthly' , 0. , 'socoefr0', .false. , .true. , 'yearly' , '' , '' , ''
423	sn_s_rnf = 'runoffs' , 24. , 'rosaline', .true. , .true. , 'yearly' , '' , '' , ''
424	sn_t_rnf = 'runoffs' , 24. , 'rotemper', .true. , .true. , 'yearly' , '' , '' , ''
425	sn_dep_rnf = 'runoffs' , 0. , 'rodepth' , .false. , .true. , 'yearly' , '' , '' , ''
426	/
427	!-----------------------------------------------------------------------
428	&namsbc_apr ! Atmospheric pressure used as ocean forcing (ln_apr_dyn =T)
429	!-----------------------------------------------------------------------
430	rn_pref = 101000. ! reference atmospheric pressure [N/m2]/
431	ln_ref_apr = .false. ! ref. pressure: global mean Patm (T) or a constant (F)
432	ln_apr_obc = .false. ! inverse barometer added to OBC ssh data
433
434	cn_dir = './' ! root directory for the Patm data location
435	!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
436	! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
437	! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
438	sn_apr = 'patm' , -1. ,'somslpre' , .true. , .true. , 'yearly' , '' , '' , ''
439	/
440	!-----------------------------------------------------------------------
441	&namsbc_isf ! Top boundary layer (ISF) (ln_isfcav =T : read (ln_read_cfg=T)
442	!----------------------------------------------------------------------- or set or usr_def_zgr )
443	! ! type of top boundary layer
444	nn_isf = 1 ! ice shelf melting/freezing
445	! 1 = presence of ISF ; 2 = bg03 parametrisation
446	! 3 = rnf file for ISF ; 4 = ISF specified freshwater flux
447	! options 1 and 4 need ln_isfcav = .true. (domzgr)
448	! ! nn_isf = 1 or 2 cases:
449	rn_gammat0 = 1.e-4 ! gammat coefficient used in blk formula
450	rn_gammas0 = 1.e-4 ! gammas coefficient used in blk formula
451	! ! nn_isf = 1 or 4 cases:
452	rn_hisf_tbl = 30. ! thickness of the top boundary layer (Losh et al. 2008)
453	! ! 0 => thickness of the tbl = thickness of the first wet cell
454	! ! nn_isf = 1 case
455	nn_isfblk = 1 ! 1 ISOMIP like: 2 equations formulation (Hunter et al., 2006)
456	! ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015)
457	nn_gammablk = 1 ! 0 = cst Gammat (= gammat/s)
458	! ! 1 = velocity dependend Gamma (u* * gammat/s) (Jenkins et al. 2010)
459	! ! 2 = velocity and stability dependent Gamma (Holland et al. 1999)
460
461	!___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________!
462	! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !
463	! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !
464	!* nn_isf = 4 case
465	sn_fwfisf = 'rnfisf' , -12. ,'sowflisf' , .false. , .true. , 'yearly' , '' , '' , ''
466	!* nn_isf = 3 case
467	sn_rnfisf = 'rnfisf' , -12. ,'sofwfisf' , .false. , .true. , 'yearly' , '' , '' , ''
468	!* nn_isf = 2 and 3 cases
469	sn_depmax_isf ='rnfisf' , -12. ,'sozisfmax', .false. , .true. , 'yearly' , '' , '' , ''
470	sn_depmin_isf ='rnfisf' , -12. ,'sozisfmin', .false. , .true. , 'yearly' , '' , '' , ''
471	!* nn_isf = 2 case
472	sn_Leff_isf = 'rnfisf' , -12. ,'Leff' , .false. , .true. , 'yearly' , '' , '' , ''
473	/
474	!-----------------------------------------------------------------------
475	&namsbc_iscpl ! land ice / ocean coupling option (ln_isfcav =T : read (ln_read_cfg=T)
476	!----------------------------------------------------------------------- or set or usr_def_zgr )
477	nn_drown = 10 ! number of iteration of the extrapolation loop (fill the new wet cells)
478	ln_hsb = .false. ! activate conservation module (conservation exact after a time of rn_fiscpl)
479	nn_fiscpl = 43800 ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency)
480	/
481	!-----------------------------------------------------------------------
482	&namsbc_wave ! External fields from wave model (ln_wave=T)
483	!-----------------------------------------------------------------------
484	cn_dir = './' ! root directory for the waves data location
485	!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
486	! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
487	! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
488	sn_cdg = 'sdw_ecwaves_orca2' , 6. , 'drag_coeff' , .true. , .true. , 'yearly' , '' , '' , ''
489	sn_usd = 'sdw_ecwaves_orca2' , 6. , 'u_sd2d' , .true. , .true. , 'yearly' , '' , '' , ''
490	sn_vsd = 'sdw_ecwaves_orca2' , 6. , 'v_sd2d' , .true. , .true. , 'yearly' , '' , '' , ''
491	sn_hsw = 'sdw_ecwaves_orca2' , 6. , 'hs' , .true. , .true. , 'yearly' , '' , '' , ''
492	sn_wmp = 'sdw_ecwaves_orca2' , 6. , 'wmp' , .true. , .true. , 'yearly' , '' , '' , ''
493	sn_wfr = 'sdw_ecwaves_orca2' , 6. , 'wfr' , .true. , .true. , 'yearly' , '' , '' , ''
494	sn_wnum = 'sdw_ecwaves_orca2' , 6. , 'wave_num' , .true. , .true. , 'yearly' , '' , '' , ''
495	sn_tauwoc = 'sdw_ecwaves_orca2' , 6. , 'wave_stress', .true. , .true. , 'yearly' , '' , '' , ''
496	sn_tauwx = 'sdw_ecwaves_orca2' , 6. , 'wave_stress', .true. , .true. , 'yearly' , '' , '' , ''
497	sn_tauwy = 'sdw_ecwaves_orca2' , 6. , 'wave_stress', .true. , .true. , 'yearly' , '' , '' , ''
498	/
499	!-----------------------------------------------------------------------
500	&namberg ! iceberg parameters (default: OFF)
501	!-----------------------------------------------------------------------
502	ln_icebergs = .false. ! activate iceberg floats (force =F with "key_agrif")
503	!
504	! ! diagnostics:
505	ln_bergdia = .true. ! Calculate budgets
506	nn_verbose_level = 0 ! Turn on more verbose output if level > 0
507	nn_verbose_write = 15 ! Timesteps between verbose messages
508	nn_sample_rate = 1 ! Timesteps between sampling for trajectory storage
509	!
510	! ! iceberg setting:
511	! ! Initial mass required for an iceberg of each class
512	rn_initial_mass = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11
513	! ! Proportion of calving mass to apportion to each class
514	rn_distribution = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02
515	! ! Ratio between effective and real iceberg mass (non-dim)
516	! ! i.e. number of icebergs represented at a point
517	rn_mass_scaling = 2000., 200., 50., 20., 10., 5., 2., 1., 1., 1.
518	! thickness of newly calved bergs (m)
519	rn_initial_thickness = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250.
520	!
521	rn_rho_bergs = 850. ! Density of icebergs
522	rn_LoW_ratio = 1.5 ! Initial ratio L/W for newly calved icebergs
523	ln_operator_splitting = .true. ! Use first order operator splitting for thermodynamics
524	rn_bits_erosion_fraction = 0. ! Fraction of erosion melt flux to divert to bergy bits
525	rn_sicn_shift = 0. ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1)
526	ln_passive_mode = .false. ! iceberg - ocean decoupling
527	nn_test_icebergs = 10 ! Create test icebergs of this class (-1 = no)
528	! ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2)
529	rn_test_box = 108.0, 116.0, -66.0, -58.0
530	ln_use_calving = .false. ! Use calving data even when nn_test_icebergs > 0
531	rn_speed_limit = 0. ! CFL speed limit for a berg
532
533	cn_dir = './' ! root directory for the calving data location
534	!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
535	! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
536	! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
537	sn_icb = 'calving' , -1. ,'calvingmask', .true. , .true. , 'yearly' , '' , '' , ''
538	/
539
540	!!======================================================================
541	!! * Lateral boundary condition * !!
542	!! !!
543	!! namlbc lateral momentum boundary condition (default: NO selection)
544	!! namagrif agrif nested grid (read by child model only) ("key_agrif")
545	!! nam_tide Tidal forcing (default: OFF)
546	!! nambdy Unstructured open boundaries (default: OFF)
547	!! nambdy_dta Unstructured open boundaries - external data (see nambdy)
548	!! nambdy_tide tidal forcing at open boundaries (default: OFF)
549	!!======================================================================
550	!
551	!-----------------------------------------------------------------------
552	&namlbc ! lateral momentum boundary condition (default: NO selection)
553	!-----------------------------------------------------------------------
554	! ! free slip ! partial slip ! no slip ! strong slip
555	rn_shlat = -9999. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat
556	ln_vorlat = .false. ! consistency of vorticity boundary condition with analytical Eqs.
557	/
558	!-----------------------------------------------------------------------
559	&namagrif ! AGRIF zoom ("key_agrif")
560	!-----------------------------------------------------------------------
561	ln_spc_dyn = .true. ! use 0 as special value for dynamics
562	rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [m2/s]
563	rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [m2/s]
564	ln_chk_bathy = .false. ! =T check the parent bathymetry
565	/
566	!-----------------------------------------------------------------------
567	&nam_tide ! tide parameters (default: OFF)
568	!-----------------------------------------------------------------------
569	ln_tide = .false. ! Activate tides
570	ln_tide_pot = .true. ! use tidal potential forcing
571	ln_scal_load = .false. ! Use scalar approximation for
572	rn_scal_load = 0.094 ! load potential
573	ln_read_load = .false. ! Or read load potential from file
574	cn_tide_load = 'tide_LOAD_grid_T.nc' ! filename for load potential
575	!
576	ln_tide_ramp = .false. ! Use linear ramp for tides at startup
577	rdttideramp = 0. ! ramp duration in days
578	clname(1) = 'DUMMY' ! name of constituent - all tidal components must be set in namelist_cfg
579	/
580	!-----------------------------------------------------------------------
581	&nambdy ! unstructured open boundaries (default: OFF)
582	!-----------------------------------------------------------------------
583	ln_bdy = .false. ! Use unstructured open boundaries
584	nb_bdy = 0 ! number of open boundary sets
585	ln_coords_file = .true. ! =T : read bdy coordinates from file
586	cn_coords_file = 'coordinates.bdy.nc' ! bdy coordinates files
587	ln_mask_file = .false. ! =T : read mask from file
588	cn_mask_file = '' ! name of mask file (if ln_mask_file=.TRUE.)
589	cn_dyn2d = 'none' !
590	nn_dyn2d_dta = 0 ! = 0, bdy data are equal to the initial state
591	! ! = 1, bdy data are read in 'bdydata .nc' files
592	! ! = 2, use tidal harmonic forcing data from files
593	! ! = 3, use external data AND tidal harmonic forcing
594	cn_dyn3d = 'none' !
595	nn_dyn3d_dta = 0 ! = 0, bdy data are equal to the initial state
596	! ! = 1, bdy data are read in 'bdydata .nc' files
597	cn_tra = 'none' !
598	nn_tra_dta = 0 ! = 0, bdy data are equal to the initial state
599	! ! = 1, bdy data are read in 'bdydata .nc' files
600	cn_ice = 'none' !
601	nn_ice_dta = 0 ! = 0, bdy data are equal to the initial state
602	! ! = 1, bdy data are read in 'bdydata .nc' files
603	!
604	ln_tra_dmp =.false. ! open boudaries conditions for tracers
605	ln_dyn3d_dmp =.false. ! open boundary condition for baroclinic velocities
606	rn_time_dmp = 1. ! Damping time scale in days
607	rn_time_dmp_out = 1. ! Outflow damping time scale
608	nn_rimwidth = 10 ! width of the relaxation zone
609	ln_vol = .false. ! total volume correction (see nn_volctl parameter)
610	nn_volctl = 1 ! = 0, the total water flux across open boundaries is zero
611	/
612	!-----------------------------------------------------------------------
613	&nambdy_dta ! open boundaries - external data (see nam_bdy)
614	!-----------------------------------------------------------------------
615	ln_zinterp = .false. ! T if a vertical interpolation is required. Variables gdep[tuv] and e3[tuv] must exist in the file
616	! ! automatically defined to T if the number of vertical levels in bdy dta /= jpk
617	ln_full_vel = .false. ! T if [uv]3d are "full" velocities and not only its baroclinic components
618	! ! in this case, baroclinic and barotropic velocities will be recomputed -> [uv]2d not needed
619	!
620	cn_dir = 'bdydta/' ! root directory for the BDY data location
621	!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
622	! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
623	! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
624	bn_ssh = 'amm12_bdyT_u2d' , 24. , 'sossheig', .true. , .false., 'daily' , '' , '' , ''
625	bn_u2d = 'amm12_bdyU_u2d' , 24. , 'vobtcrtx', .true. , .false., 'daily' , '' , '' , ''
626	bn_v2d = 'amm12_bdyV_u2d' , 24. , 'vobtcrty', .true. , .false., 'daily' , '' , '' , ''
627	bn_u3d = 'amm12_bdyU_u3d' , 24. , 'vozocrtx', .true. , .false., 'daily' , '' , '' , ''
628	bn_v3d = 'amm12_bdyV_u3d' , 24. , 'vomecrty', .true. , .false., 'daily' , '' , '' , ''
629	bn_tem = 'amm12_bdyT_tra' , 24. , 'votemper', .true. , .false., 'daily' , '' , '' , ''
630	bn_sal = 'amm12_bdyT_tra' , 24. , 'vosaline', .true. , .false., 'daily' , '' , '' , ''
631	!* for si3
632	bn_a_i = 'amm12_bdyT_ice' , 24. , 'siconc' , .true. , .false., 'daily' , '' , '' , ''
633	bn_h_i = 'amm12_bdyT_ice' , 24. , 'sithic' , .true. , .false., 'daily' , '' , '' , ''
634	bn_h_s = 'amm12_bdyT_ice' , 24. , 'snthic' , .true. , .false., 'daily' , '' , '' , ''
635	bn_t_i = 'NOT USED' , 24. , 'sitemp' , .true. , .false., 'daily' , '' , '' , ''
636	bn_t_s = 'NOT USED' , 24. , 'sntemp' , .true. , .false., 'daily' , '' , '' , ''
637	bn_tsu = 'NOT USED' , 24. , 'sittop' , .true. , .false., 'daily' , '' , '' , ''
638	bn_s_i = 'NOT USED' , 24. , 'sisalt' , .true. , .false., 'daily' , '' , '' , ''
639	! melt ponds (be careful, bn_aip is the pond concentration (not fraction), so it differs from rn_iceapnd)
640	bn_aip = 'NOT USED' , 24. , 'siapnd' , .true. , .false., 'daily' , '' , '' , ''
641	bn_hip = 'NOT USED' , 24. , 'sihpnd' , .true. , .false., 'daily' , '' , '' , ''
642	! if bn_t_i etc are "not used", then define arbitrary temperatures and salinity and ponds
643	rn_ice_tem = 270. ! arbitrary temperature of incoming sea ice
644	rn_ice_sal = 10. ! -- salinity --
645	rn_ice_age = 30. ! -- age --
646	rn_ice_apnd = 0.2 ! -- pond fraction = a_ip/a_i --
647	rn_ice_hpnd = 0.05 ! -- pond depth --
648	/
649	!-----------------------------------------------------------------------
650	&nambdy_tide ! tidal forcing at open boundaries (default: OFF)
651	!-----------------------------------------------------------------------
652	filtide = 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files
653	ln_bdytide_2ddta = .false. !
654	ln_bdytide_conj = .false. !
655	/
656
657	!!======================================================================
658	!! * Top/Bottom boundary condition * !!
659	!! !!
660	!! namdrg top/bottom drag coefficient (default: NO selection)
661	!! namdrg_top top friction (ln_OFF=F & ln_isfcav=T)
662	!! namdrg_bot bottom friction (ln_OFF=F)
663	!! nambbc bottom temperature boundary condition (default: OFF)
664	!! nambbl bottom boundary layer scheme (default: OFF)
665	!!======================================================================
666	!
667	!-----------------------------------------------------------------------
668	&namdrg ! top/bottom drag coefficient (default: NO selection)
669	!-----------------------------------------------------------------------
670	ln_OFF = .false. ! free-slip : Cd = 0 (F => fill namdrg_bot
671	ln_lin = .false. ! linear drag: Cd = Cd0 Uc0 & namdrg_top)
672	ln_non_lin = .false. ! non-linear drag: Cd = Cd0 \|U\|
673	ln_loglayer = .false. ! logarithmic drag: Cd = vkarmn/log(z/z0) \|U\|
674	!
675	ln_drgimp = .true. ! implicit top/bottom friction flag
676	/
677	!-----------------------------------------------------------------------
678	&namdrg_top ! TOP friction (ln_OFF =F & ln_isfcav=T)
679	!-----------------------------------------------------------------------
680	rn_Cd0 = 1.e-3 ! drag coefficient [-]
681	rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0)
682	rn_Cdmax = 0.1 ! drag value maximum [-] (logarithmic drag)
683	rn_ke0 = 2.5e-3 ! background kinetic energy [m2/s2] (non-linear cases)
684	rn_z0 = 3.0e-3 ! roughness [m] (ln_loglayer=T)
685	ln_boost = .false. ! =T regional boost of Cd0 ; =F constant
686	rn_boost = 50. ! local boost factor [-]
687	/
688	!-----------------------------------------------------------------------
689	&namdrg_bot ! BOTTOM friction (ln_OFF =F)
690	!-----------------------------------------------------------------------
691	rn_Cd0 = 1.e-3 ! drag coefficient [-]
692	rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0)
693	rn_Cdmax = 0.1 ! drag value maximum [-] (logarithmic drag)
694	rn_ke0 = 2.5e-3 ! background kinetic energy [m2/s2] (non-linear cases)
695	rn_z0 = 3.e-3 ! roughness [m] (ln_loglayer=T)
696	ln_boost = .false. ! =T regional boost of Cd0 ; =F constant
697	rn_boost = 50. ! local boost factor [-]
698	/
699	!-----------------------------------------------------------------------
700	&nambbc ! bottom temperature boundary condition (default: OFF)
701	!-----------------------------------------------------------------------
702	ln_trabbc = .false. ! Apply a geothermal heating at the ocean bottom
703	nn_geoflx = 2 ! geothermal heat flux: = 1 constant flux
704	! ! = 2 read variable flux [mW/m2]
705	rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [mW/m2]
706
707	cn_dir = './' ! root directory for the geothermal data location
708	!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
709	! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
710	! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
711	sn_qgh ='geothermal_heating.nc' , -12. , 'heatflow', .false. , .true. , 'yearly' , '' , '' , ''
712	/
713	!-----------------------------------------------------------------------
714	&nambbl ! bottom boundary layer scheme (default: OFF)
715	!-----------------------------------------------------------------------
716	ln_trabbl = .false. ! Bottom Boundary Layer parameterisation flag
717	nn_bbl_ldf = 1 ! diffusive bbl (=1) or not (=0)
718	nn_bbl_adv = 0 ! advective bbl (=1/2) or not (=0)
719	rn_ahtbbl = 1000. ! lateral mixing coefficient in the bbl [m2/s]
720	rn_gambbl = 10. ! advective bbl coefficient [s]
721	/
722
723	!!======================================================================
724	!! Tracer (T-S) namelists !!
725	!! !!
726	!! nameos equation of state (default: NO selection)
727	!! namtra_adv advection scheme (default: NO selection)
728	!! namtra_ldf lateral diffusion scheme (default: NO selection)
729	!! namtra_mle mixed layer eddy param. (Fox-Kemper param.) (default: OFF)
730	!! namtra_eiv eddy induced velocity param. (default: OFF)
731	!! namtra_dmp T & S newtonian damping (default: OFF)
732	!!======================================================================
733	!
734	!-----------------------------------------------------------------------
735	&nameos ! ocean Equation Of Seawater (default: NO selection)
736	!-----------------------------------------------------------------------
737	ln_teos10 = .false. ! = Use TEOS-10
738	ln_eos80 = .false. ! = Use EOS80
739	ln_seos = .false. ! = Use S-EOS (simplified Eq.)
740	!
741	! ! S-EOS coefficients (ln_seos=T):
742	! ! rd(T,S,Z)rau0 = -a0(1+.5lambdadT+muZ+nudS)dT+b0dS
743	rn_a0 = 1.6550e-1 ! thermal expension coefficient
744	rn_b0 = 7.6554e-1 ! saline expension coefficient
745	rn_lambda1 = 5.9520e-2 ! cabbeling coeff in T^2 (=0 for linear eos)
746	rn_lambda2 = 7.4914e-4 ! cabbeling coeff in S^2 (=0 for linear eos)
747	rn_mu1 = 1.4970e-4 ! thermobaric coeff. in T (=0 for linear eos)
748	rn_mu2 = 1.1090e-5 ! thermobaric coeff. in S (=0 for linear eos)
749	rn_nu = 2.4341e-3 ! cabbeling coeff in T*S (=0 for linear eos)
750	/
751	!-----------------------------------------------------------------------
752	&namtra_adv ! advection scheme for tracer (default: NO selection)
753	!-----------------------------------------------------------------------
754	ln_traadv_OFF = .false. ! No tracer advection
755	ln_traadv_cen = .false. ! 2nd order centered scheme
756	nn_cen_h = 4 ! =2/4, horizontal 2nd order CEN / 4th order CEN
757	nn_cen_v = 4 ! =2/4, vertical 2nd order CEN / 4th order COMPACT
758	ln_traadv_fct = .false. ! FCT scheme
759	nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order
760	nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order
761	ln_traadv_mus = .false. ! MUSCL scheme
762	ln_mus_ups = .false. ! use upstream scheme near river mouths
763	ln_traadv_ubs = .false. ! UBS scheme
764	nn_ubs_v = 2 ! =2 , vertical 2nd order FCT / COMPACT 4th order
765	ln_traadv_qck = .false. ! QUICKEST scheme
766	/
767	!-----------------------------------------------------------------------
768	&namtra_ldf ! lateral diffusion scheme for tracers (default: NO selection)
769	!-----------------------------------------------------------------------
770	! ! Operator type:
771	ln_traldf_OFF = .false. ! No explicit diffusion
772	ln_traldf_lap = .false. ! laplacian operator
773	ln_traldf_blp = .false. ! bilaplacian operator
774	!
775	! ! Direction of action:
776	ln_traldf_lev = .false. ! iso-level
777	ln_traldf_hor = .false. ! horizontal (geopotential)
778	ln_traldf_iso = .false. ! iso-neutral (standard operator)
779	ln_traldf_triad = .false. ! iso-neutral (triad operator)
780	!
781	! ! iso-neutral options:
782	ln_traldf_msc = .false. ! Method of Stabilizing Correction (both operators)
783	rn_slpmax = 0.01 ! slope limit (both operators)
784	ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only)
785	rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only)
786	ln_botmix_triad = .false. ! lateral mixing on bottom (triad only)
787	!
788	! ! Coefficients:
789	nn_aht_ijk_t = 0 ! space/time variation of eddy coefficient:
790	! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file
791	! ! = 0 constant
792	! ! = 10 F(k) =ldf_c1d
793	! ! = 20 F(i,j) =ldf_c2d
794	! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation
795	! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d
796	! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing)
797	! ! time invariant coefficients: aht0 = 1/2 Ud*Ld (lap case)
798	! ! or = 1/12 Ud*Ld^3 (blp case)
799	rn_Ud = 0.01 ! lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30)
800	rn_Ld = 200.e+3 ! lateral diffusive length [m] (nn_aht_ijk_t= 0, 10)
801	/
802	!-----------------------------------------------------------------------
803	&namtra_mle ! mixed layer eddy parametrisation (Fox-Kemper) (default: OFF)
804	!-----------------------------------------------------------------------
805	ln_mle = .false. ! (T) use the Mixed Layer Eddy (MLE) parameterisation
806	rn_ce = 0.06 ! magnitude of the MLE (typical value: 0.06 to 0.08)
807	nn_mle = 1 ! MLE type: =0 standard Fox-Kemper ; =1 new formulation
808	rn_lf = 5.e+3 ! typical scale of mixed layer front (meters) (case rn_mle=0)
809	rn_time = 172800. ! time scale for mixing momentum across the mixed layer (seconds) (case rn_mle=0)
810	rn_lat = 20. ! reference latitude (degrees) of MLE coef. (case rn_mle=1)
811	nn_mld_uv = 0 ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max)
812	nn_conv = 0 ! =1 no MLE in case of convection ; =0 always MLE
813	rn_rho_c_mle = 0.01 ! delta rho criterion used to calculate MLD for FK
814	/
815	!-----------------------------------------------------------------------
816	&namtra_eiv ! eddy induced velocity param. (default: OFF)
817	!-----------------------------------------------------------------------
818	ln_ldfeiv = .false. ! use eddy induced velocity parameterization
819	!
820	! ! Coefficients:
821	nn_aei_ijk_t = 0 ! space/time variation of eddy coefficient:
822	! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file
823	! ! = 0 constant
824	! ! = 10 F(k) =ldf_c1d
825	! ! = 20 F(i,j) =ldf_c2d
826	! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation
827	! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d
828	! ! time invariant coefficients: aei0 = 1/2 Ue*Le
829	rn_Ue = 0.02 ! lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30)
830	rn_Le = 200.e+3 ! lateral diffusive length [m] (nn_aht_ijk_t= 0, 10)
831	!
832	ln_ldfeiv_dia =.false. ! diagnose eiv stream function and velocities
833	/
834	!-----------------------------------------------------------------------
835	&namtra_dmp ! tracer: T & S newtonian damping (default: OFF)
836	!-----------------------------------------------------------------------
837	ln_tradmp = .false. ! add a damping term (using resto.nc coef.)
838	nn_zdmp = 0 ! vertical shape =0 damping throughout the water column
839	! ! =1 no damping in the mixing layer (kz criteria)
840	! ! =2 no damping in the mixed layer (rho crieria)
841	cn_resto = 'resto.nc' ! Name of file containing restoration coeff. field (use dmp_tools to create this)
842	/
843
844	!!======================================================================
845	!! * Dynamics namelists * !!
846	!! !!
847	!! nam_vvl vertical coordinate options (default: z-star)
848	!! namdyn_adv formulation of the momentum advection (default: NO selection)
849	!! namdyn_vor advection scheme (default: NO selection)
850	!! namdyn_hpg hydrostatic pressure gradient (default: NO selection)
851	!! namdyn_spg surface pressure gradient (default: NO selection)
852	!! namdyn_ldf lateral diffusion scheme (default: NO selection)
853	!! namdta_dyn offline TOP: dynamics read in files (OFF_SRC only)
854	!!======================================================================
855	!
856	!-----------------------------------------------------------------------
857	&nam_vvl ! vertical coordinate options (default: z-star)
858	!-----------------------------------------------------------------------
859	ln_vvl_zstar = .true. ! z-star vertical coordinate
860	ln_vvl_ztilde = .false. ! z-tilde vertical coordinate: only high frequency variations
861	ln_vvl_layer = .false. ! full layer vertical coordinate
862	ln_vvl_ztilde_as_zstar = .false. ! ztilde vertical coordinate emulating zstar
863	ln_vvl_zstar_at_eqtor = .false. ! ztilde near the equator
864	rn_ahe3 = 0.0 ! thickness diffusion coefficient
865	rn_rst_e3t = 30.0 ! ztilde to zstar restoration timescale [days]
866	rn_lf_cutoff = 5.0 ! cutoff frequency for low-pass filter [days]
867	rn_zdef_max = 0.9 ! maximum fractional e3t deformation
868	ln_vvl_dbg = .true. ! debug prints (T/F)
869	/
870	!-----------------------------------------------------------------------
871	&namdyn_adv ! formulation of the momentum advection (default: NO selection)
872	!-----------------------------------------------------------------------
873	ln_dynadv_OFF = .false. ! linear dynamics (no momentum advection)
874	ln_dynadv_vec = .false. ! vector form - 2nd centered scheme
875	nn_dynkeg = 0 ! grad(KE) scheme: =0 C2 ; =1 Hollingsworth correction
876	ln_dynadv_cen2 = .false. ! flux form - 2nd order centered scheme
877	ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme
878	/
879	!-----------------------------------------------------------------------
880	&namdyn_vor ! Vorticity / Coriolis scheme (default: NO selection)
881	!-----------------------------------------------------------------------
882	ln_dynvor_ene = .false. ! energy conserving scheme
883	ln_dynvor_ens = .false. ! enstrophy conserving scheme
884	ln_dynvor_mix = .false. ! mixed scheme
885	ln_dynvor_enT = .false. ! energy conserving scheme (T-point)
886	ln_dynvor_eeT = .false. ! energy conserving scheme (een using e3t)
887	ln_dynvor_een = .false. ! energy & enstrophy scheme
888	nn_een_e3f = 0 ! =0 e3f = mi(mj(e3t))/4
889	! ! =1 e3f = mi(mj(e3t))/mi(mj( tmask))
890	ln_dynvor_msk = .false. ! vorticity multiplied by fmask (=T) ==>>> PLEASE DO NOT ACTIVATE
891	! ! (f-point vorticity schemes only)
892	/
893	!-----------------------------------------------------------------------
894	&namdyn_hpg ! Hydrostatic pressure gradient option (default: NO selection)
895	!-----------------------------------------------------------------------
896	ln_hpg_zco = .false. ! z-coordinate - full steps
897	ln_hpg_zps = .false. ! z-coordinate - partial steps (interpolation)
898	ln_hpg_sco = .false. ! s-coordinate (standard jacobian formulation)
899	ln_hpg_isf = .false. ! s-coordinate (sco ) adapted to isf
900	ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial)
901	ln_hpg_prj = .false. ! s-coordinate (Pressure Jacobian scheme)
902	/
903	!-----------------------------------------------------------------------
904	&namdyn_spg ! surface pressure gradient (default: NO selection)
905	!-----------------------------------------------------------------------
906	ln_dynspg_exp = .false. ! explicit free surface
907	ln_dynspg_ts = .false. ! split-explicit free surface
908	ln_bt_fw = .true. ! Forward integration of barotropic Eqs.
909	ln_bt_av = .true. ! Time filtering of barotropic variables
910	nn_bt_flt = 1 ! Time filter choice = 0 None
911	! ! = 1 Boxcar over nn_baro sub-steps
912	! ! = 2 Boxcar over 2*nn_baro " "
913	ln_bt_auto = .true. ! Number of sub-step defined from:
914	rn_bt_cmax = 0.8 ! =T : the Maximum Courant Number allowed
915	nn_baro = 30 ! =F : the number of sub-step in rn_rdt seconds
916	rn_bt_alpha = 0. ! Temporal diffusion parameter (if ln_bt_av=F)
917	/
918	!-----------------------------------------------------------------------
919	&namdyn_ldf ! lateral diffusion on momentum (default: NO selection)
920	!-----------------------------------------------------------------------
921	! ! Type of the operator :
922	ln_dynldf_OFF = .false. ! No operator (i.e. no explicit diffusion)
923	ln_dynldf_lap = .false. ! laplacian operator
924	ln_dynldf_blp = .false. ! bilaplacian operator
925	! ! Direction of action :
926	ln_dynldf_lev = .false. ! iso-level
927	ln_dynldf_hor = .false. ! horizontal (geopotential)
928	ln_dynldf_iso = .false. ! iso-neutral (lap only)
929	! ! Coefficient
930	nn_ahm_ijk_t = 0 ! space/time variation of eddy coefficient :
931	! ! =-30 read in eddy_viscosity_3D.nc file
932	! ! =-20 read in eddy_viscosity_2D.nc file
933	! ! = 0 constant
934	! ! = 10 F(k)=c1d
935	! ! = 20 F(i,j)=F(grid spacing)=c2d
936	! ! = 30 F(i,j,k)=c2d*c1d
937	! ! = 31 F(i,j,k)=F(grid spacing and local velocity)
938	! ! = 32 F(i,j,k)=F(local gridscale and deformation rate)
939	! ! time invariant coefficients : ahm = 1/2 Uv*Lv (lap case)
940	! ! or = 1/12 Uv*Lv^3 (blp case)
941	rn_Uv = 0.1 ! lateral viscous velocity [m/s] (nn_ahm_ijk_t= 0, 10, 20, 30)
942	rn_Lv = 10.e+3 ! lateral viscous length [m] (nn_ahm_ijk_t= 0, 10)
943	! ! Smagorinsky settings (nn_ahm_ijk_t= 32) :
944	rn_csmc = 3.5 ! Smagorinsky constant of proportionality
945	rn_minfac = 1.0 ! multiplier of theorectical lower limit
946	rn_maxfac = 1.0 ! multiplier of theorectical upper limit
947	! ! iso-neutral laplacian operator (ln_dynldf_iso=T) :
948	rn_ahm_b = 0.0 ! background eddy viscosity [m2/s]
949	/
950	!-----------------------------------------------------------------------
951	&namdta_dyn ! offline ocean input files (OFF_SRC only)
952	!-----------------------------------------------------------------------
953	ln_dynrnf = .false. ! runoffs option enabled (T) or not (F)
954	ln_dynrnf_depth = .false. ! runoffs is spread in vertical (T) or not (F)
955	! fwbcorr = 3.786e-06 ! annual global mean of empmr for ssh correction
956
957	cn_dir = './' ! root directory for the ocean data location
958	!___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________!
959	! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask !
960	! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename !
961	sn_tem = 'dyna_grid_T' , 120. , 'votemper' , .true. , .true. , 'yearly' , '' , '' , ''
962	sn_sal = 'dyna_grid_T' , 120. , 'vosaline' , .true. , .true. , 'yearly' , '' , '' , ''
963	sn_mld = 'dyna_grid_T' , 120. , 'somixhgt' , .true. , .true. , 'yearly' , '' , '' , ''
964	sn_emp = 'dyna_grid_T' , 120. , 'sowaflup' , .true. , .true. , 'yearly' , '' , '' , ''
965	sn_fmf = 'dyna_grid_T' , 120. , 'iowaflup' , .true. , .true. , 'yearly' , '' , '' , ''
966	sn_ice = 'dyna_grid_T' , 120. , 'soicecov' , .true. , .true. , 'yearly' , '' , '' , ''
967	sn_qsr = 'dyna_grid_T' , 120. , 'soshfldo' , .true. , .true. , 'yearly' , '' , '' , ''
968	sn_wnd = 'dyna_grid_T' , 120. , 'sowindsp' , .true. , .true. , 'yearly' , '' , '' , ''
969	sn_uwd = 'dyna_grid_U' , 120. , 'uocetr_eff', .true. , .true. , 'yearly' , '' , '' , ''
970	sn_vwd = 'dyna_grid_V' , 120. , 'vocetr_eff', .true. , .true. , 'yearly' , '' , '' , ''
971	sn_wwd = 'dyna_grid_W' , 120. , 'wocetr_eff', .true. , .true. , 'yearly' , '' , '' , ''
972	sn_avt = 'dyna_grid_W' , 120. , 'voddmavs' , .true. , .true. , 'yearly' , '' , '' , ''
973	sn_ubl = 'dyna_grid_U' , 120. , 'sobblcox' , .true. , .true. , 'yearly' , '' , '' , ''
974	sn_vbl = 'dyna_grid_V' , 120. , 'sobblcoy' , .true. , .true. , 'yearly' , '' , '' , ''
975	/
976
977	!!======================================================================
978	!! vertical physics namelists !!
979	!! !!
980	!! namzdf vertical physics manager (default: NO selection)
981	!! namzdf_ric richardson number vertical mixing (ln_zdfric=T)
982	!! namzdf_tke TKE vertical mixing (ln_zdftke=T)
983	!! namzdf_gls GLS vertical mixing (ln_zdfgls=T)
984	!! namzdf_osm OSM vertical diffusion (ln_zdfosm=T)
985	!! namzdf_iwm tidal mixing parameterization (ln_zdfiwm=T)
986	!!======================================================================
987	!
988	!-----------------------------------------------------------------------
989	&namzdf ! vertical physics manager (default: NO selection)
990	!-----------------------------------------------------------------------
991	! ! adaptive-implicit vertical advection
992	ln_zad_Aimp = .false. ! Courant number dependent scheme (Shchepetkin 2015)
993	!
994	! ! type of vertical closure (required)
995	ln_zdfcst = .false. ! constant mixing
996	ln_zdfric = .false. ! local Richardson dependent formulation (T => fill namzdf_ric)
997	ln_zdftke = .false. ! Turbulent Kinetic Energy closure (T => fill namzdf_tke)
998	ln_zdfgls = .false. ! Generic Length Scale closure (T => fill namzdf_gls)
999	ln_zdfosm = .false. ! OSMOSIS BL closure (T => fill namzdf_osm)
1000	!
1001	! ! convection
1002	ln_zdfevd = .false. ! enhanced vertical diffusion
1003	nn_evdm = 0 ! apply on tracer (=0) or on tracer and momentum (=1)
1004	rn_evd = 100. ! mixing coefficient [m2/s]
1005	ln_zdfnpc = .false. ! Non-Penetrative Convective algorithm
1006	nn_npc = 1 ! frequency of application of npc
1007	nn_npcp = 365 ! npc control print frequency
1008	!
1009	ln_zdfddm = .false. ! double diffusive mixing
1010	rn_avts = 1.e-4 ! maximum avs (vertical mixing on salinity)
1011	rn_hsbfr = 1.6 ! heat/salt buoyancy flux ratio
1012	!
1013	! ! gravity wave-driven vertical mixing
1014	ln_zdfiwm = .false. ! internal wave-induced mixing (T => fill namzdf_iwm)
1015	ln_zdfswm = .false. ! surface wave-induced mixing (T => ln_wave=ln_sdw=T )
1016	!
1017	! ! coefficients
1018	rn_avm0 = 1.2e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F)
1019	rn_avt0 = 1.2e-5 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F)
1020	nn_avb = 0 ! profile for background avt & avm (=1) or not (=0)
1021	nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0)
1022	/
1023	!-----------------------------------------------------------------------
1024	&namzdf_ric ! richardson number dependent vertical diffusion (ln_zdfric =T)
1025	!-----------------------------------------------------------------------
1026	rn_avmri = 100.e-4 ! maximum value of the vertical viscosity
1027	rn_alp = 5. ! coefficient of the parameterization
1028	nn_ric = 2 ! coefficient of the parameterization
1029	ln_mldw = .false. ! enhanced mixing in the Ekman layer
1030	rn_ekmfc = 0.7 ! Factor in the Ekman depth Equation
1031	rn_mldmin = 1.0 ! minimum allowable mixed-layer depth estimate (m)
1032	rn_mldmax = 1000.0 ! maximum allowable mixed-layer depth estimate (m)
1033	rn_wtmix = 10.0 ! vertical eddy viscosity coeff [m2/s] in the mixed-layer
1034	rn_wvmix = 10.0 ! vertical eddy diffusion coeff [m2/s] in the mixed-layer
1035	/
1036	!-----------------------------------------------------------------------
1037	&namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion (ln_zdftke =T)
1038	!-----------------------------------------------------------------------
1039	rn_ediff = 0.1 ! coef. for vertical eddy coef. (avt=rn_ediffmxlsqrt(e) )
1040	rn_ediss = 0.7 ! coef. of the Kolmogoroff dissipation
1041	rn_ebb = 67.83 ! coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T)
1042	rn_emin = 1.e-6 ! minimum value of tke [m2/s2]
1043	rn_emin0 = 1.e-4 ! surface minimum value of tke [m2/s2]
1044	rn_bshear = 1.e-20 ! background shear (>0) currently a numerical threshold (do not change it)
1045	nn_pdl = 1 ! Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm)
1046	nn_mxl = 2 ! mixing length: = 0 bounded by the distance to surface and bottom
1047	! ! = 1 bounded by the local vertical scale factor
1048	! ! = 2 first vertical derivative of mixing length bounded by 1
1049	! ! = 3 as =2 with distinct dissipative an mixing length scale
1050	ln_mxl0 = .true. ! surface mixing length scale = F(wind stress) (T) or not (F)
1051	rn_mxl0 = 0.04 ! surface buoyancy lenght scale minimum value
1052	ln_drg = .false. ! top/bottom friction added as boundary condition of TKE
1053	ln_lc = .true. ! Langmuir cell parameterisation (Axell 2002)
1054	rn_lc = 0.15 ! coef. associated to Langmuir cells
1055	nn_etau = 1 ! penetration of tke below the mixed layer (ML) due to NIWs
1056	! = 0 none ; = 1 add a tke source below the ML
1057	! = 2 add a tke source just at the base of the ML
1058	! = 3 as = 1 applied on HF part of the stress (ln_cpl=T)
1059	rn_efr = 0.05 ! fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2)
1060	nn_htau = 1 ! type of exponential decrease of tke penetration below the ML
1061	! = 0 constant 10 m length scale
1062	! = 1 0.5m at the equator to 30m poleward of 40 degrees
1063	rn_eice = 4 ! below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4
1064	/
1065	!-----------------------------------------------------------------------
1066	&namzdf_gls ! GLS vertical diffusion (ln_zdfgls =T)
1067	!-----------------------------------------------------------------------
1068	rn_emin = 1.e-7 ! minimum value of e [m2/s2]
1069	rn_epsmin = 1.e-12 ! minimum value of eps [m2/s3]
1070	ln_length_lim = .true. ! limit on the dissipation rate under stable stratification (Galperin et al., 1988)
1071	rn_clim_galp = 0.267 ! galperin limit
1072	ln_sigpsi = .true. ! Activate or not Burchard 2001 mods on psi schmidt number in the wb case
1073	rn_crban = 100. ! Craig and Banner 1994 constant for wb tke flux
1074	rn_charn = 70000. ! Charnock constant for wb induced roughness length
1075	rn_hsro = 0.02 ! Minimum surface roughness
1076	rn_frac_hs = 1.3 ! Fraction of wave height as roughness (if nn_z0_met>1)
1077	nn_z0_met = 2 ! Method for surface roughness computation (0/1/2/3)
1078	! ! =3 requires ln_wave=T
1079	nn_bc_surf = 1 ! surface condition (0/1=Dir/Neum)
1080	nn_bc_bot = 1 ! bottom condition (0/1=Dir/Neum)
1081	nn_stab_func = 2 ! stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB)
1082	nn_clos = 1 ! predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
1083	/
1084	!-----------------------------------------------------------------------
1085	&namzdf_osm ! OSM vertical diffusion (ln_zdfosm =T)
1086	!-----------------------------------------------------------------------
1087	ln_use_osm_la = .false. ! Use namelist rn_osm_la
1088	rn_osm_la = 0.3 ! Turbulent Langmuir number
1089	rn_osm_dstokes = 5. ! Depth scale of Stokes drift (m)
1090	nn_ave = 0 ! choice of horizontal averaging on avt, avmu, avmv
1091	ln_dia_osm = .true. ! output OSMOSIS-OBL variables
1092	rn_osm_hbl0 = 10. ! initial hbl value
1093	ln_kpprimix = .true. ! Use KPP-style Ri# mixing below BL
1094	rn_riinfty = 0.7 ! Highest local Ri_g permitting shear instability
1095	rn_difri = 0.005 ! max Ri# diffusivity at Ri_g = 0 (m^2/s)
1096	ln_convmix = .true. ! Use convective instability mixing below BL
1097	rn_difconv = 1. ! diffusivity when unstable below BL (m2/s)
1098	nn_osm_wave = 0 ! Method used to calculate Stokes drift
1099	! ! = 2: Use ECMWF wave fields
1100	! ! = 1: Pierson Moskowitz wave spectrum
1101	! ! = 0: Constant La# = 0.3
1102	/
1103	!-----------------------------------------------------------------------
1104	&namzdf_iwm ! internal wave-driven mixing parameterization (ln_zdfiwm =T)
1105	!-----------------------------------------------------------------------
1106	nn_zpyc = 1 ! pycnocline-intensified dissipation scales as N (=1) or N^2 (=2)
1107	ln_mevar = .true. ! variable (T) or constant (F) mixing efficiency
1108	ln_tsdiff = .true. ! account for differential T/S mixing (T) or not (F)
1109	/
1110
1111	!!======================================================================
1112	!! * Diagnostics namelists * !!
1113	!! !!
1114	!! namtrd dynamics and/or tracer trends (default: OFF)
1115	!! namptr Poleward Transport Diagnostics (default: OFF)
1116	!! namhsb Heat and salt budgets (default: OFF)
1117	!! namdiu Cool skin and warm layer models (default: OFF)
1118	!! namdiu Cool skin and warm layer models (default: OFF)
1119	!! namflo float parameters (default: OFF)
1120	!! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)
1121	!! nam_diadct transports through some sections (default: OFF)
1122	!! nam_diatmb Top Middle Bottom Output (default: OFF)
1123	!! nam_dia25h 25h Mean Output (default: OFF)
1124	!! namnc4 netcdf4 chunking and compression settings ("key_netcdf4")
1125	!!======================================================================
1126	!
1127	!-----------------------------------------------------------------------
1128	&namtrd ! trend diagnostics (default: OFF)
1129	!-----------------------------------------------------------------------
1130	ln_glo_trd = .false. ! (T) global domain averaged diag for T, T^2, KE, and PE
1131	ln_dyn_trd = .false. ! (T) 3D momentum trend output
1132	ln_dyn_mxl = .false. ! (T) 2D momentum trends averaged over the mixed layer (not coded yet)
1133	ln_vor_trd = .false. ! (T) 2D barotropic vorticity trends (not coded yet)
1134	ln_KE_trd = .false. ! (T) 3D Kinetic Energy trends
1135	ln_PE_trd = .false. ! (T) 3D Potential Energy trends
1136	ln_tra_trd = .false. ! (T) 3D tracer trend output
1137	ln_tra_mxl = .false. ! (T) 2D tracer trends averaged over the mixed layer (not coded yet)
1138	nn_trd = 365 ! print frequency (ln_glo_trd=T) (unit=time step)
1139	/
1140	!!gm nn_ctls = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk)
1141	!!gm rn_ucf = 1. ! unit conversion factor (=1 -> /seconds ; =86400. -> /day)
1142	!!gm cn_trdrst_in = "restart_mld" ! suffix of ocean restart name (input)
1143	!!gm cn_trdrst_out = "restart_mld" ! suffix of ocean restart name (output)
1144	!!gm ln_trdmld_restart = .false. ! restart for ML diagnostics
1145	!!gm ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S
1146	!!gm
1147	!-----------------------------------------------------------------------
1148	&namptr ! Poleward Transport Diagnostic (default: OFF)
1149	!-----------------------------------------------------------------------
1150	ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F)
1151	ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not
1152	/
1153	!-----------------------------------------------------------------------
1154	&namhsb ! Heat and salt budgets (default: OFF)
1155	!-----------------------------------------------------------------------
1156	ln_diahsb = .false. ! output the heat and salt budgets (T) or not (F)
1157	/
1158	!-----------------------------------------------------------------------
1159	&namdiu ! Cool skin and warm layer models (default: OFF)
1160	!-----------------------------------------------------------------------
1161	ln_diurnal = .false. !
1162	ln_diurnal_only = .false. !
1163	/
1164	!-----------------------------------------------------------------------
1165	&namflo ! float parameters (default: OFF)
1166	!-----------------------------------------------------------------------
1167	ln_floats = .false. ! activate floats or not
1168	jpnfl = 1 ! total number of floats during the run
1169	jpnnewflo = 0 ! number of floats for the restart
1170	ln_rstflo = .false. ! float restart (T) or not (F)
1171	nn_writefl = 75 ! frequency of writing in float output file
1172	nn_stockfl = 5475 ! frequency of creation of the float restart file
1173	ln_argo = .false. ! Argo type floats (stay at the surface each 10 days)
1174	ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T)
1175	! ! or computed with Blanke' scheme (F)
1176	ln_ariane = .true. ! Input with Ariane tool convention(T)
1177	ln_flo_ascii= .true. ! Output with Ariane tool netcdf convention(F) or ascii file (T)
1178	/
1179	!-----------------------------------------------------------------------
1180	&nam_diaharm ! Harmonic analysis of tidal constituents (default: OFF)
1181	!-----------------------------------------------------------------------
1182	ln_diaharm = .false. ! Choose tidal harmonic output or not
1183	nit000_han = 1 ! First time step used for harmonic analysis
1184	nitend_han = 75 ! Last time step used for harmonic analysis
1185	nstep_han = 15 ! Time step frequency for harmonic analysis
1186	tname(1) = 'M2' ! Name of tidal constituents
1187	tname(2) = 'K1' ! ---
1188	/
1189	!-----------------------------------------------------------------------
1190	&nam_diadct ! transports through some sections (default: OFF)
1191	!-----------------------------------------------------------------------
1192	ln_diadct = .false. ! Calculate transport thru sections or not
1193	nn_dct = 15 ! time step frequency for transports computing
1194	nn_dctwri = 15 ! time step frequency for transports writing
1195	nn_secdebug = 112 ! 0 : no section to debug
1196	! ! -1 : debug all section
1197	! ! 0 < n : debug section number n
1198	/
1199	!-----------------------------------------------------------------------
1200	&nam_diatmb ! Top Middle Bottom Output (default: OFF)
1201	!-----------------------------------------------------------------------
1202	ln_diatmb = .false. ! Choose Top Middle and Bottom output or not
1203	/
1204	!-----------------------------------------------------------------------
1205	&nam_dia25h ! 25h Mean Output (default: OFF)
1206	!-----------------------------------------------------------------------
1207	ln_dia25h = .false. ! Choose 25h mean output or not
1208	/
1209	!-----------------------------------------------------------------------
1210	&namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4")
1211	!-----------------------------------------------------------------------
1212	nn_nchunks_i = 4 ! number of chunks in i-dimension
1213	nn_nchunks_j = 4 ! number of chunks in j-dimension
1214	nn_nchunks_k = 31 ! number of chunks in k-dimension
1215	! ! setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which
1216	! ! is optimal for postprocessing which works exclusively with horizontal slabs
1217	ln_nc4zip = .true. ! (T) use netcdf4 chunking and compression
1218	! ! (F) ignore chunking information and produce netcdf3-compatible files
1219	/
1220
1221	!!======================================================================
1222	!! * Observation & Assimilation * !!
1223	!! !!
1224	!! namobs observation and model comparison (default: OFF)
1225	!! nam_asminc assimilation increments ('key_asminc')
1226	!!======================================================================
1227	!
1228	!-----------------------------------------------------------------------
1229	&namobs ! observation usage switch (default: OFF)
1230	!-----------------------------------------------------------------------
1231	ln_diaobs = .false. ! Logical switch for the observation operator
1232	!
1233	ln_t3d = .false. ! Logical switch for T profile observations
1234	ln_s3d = .false. ! Logical switch for S profile observations
1235	ln_sla = .false. ! Logical switch for SLA observations
1236	ln_sst = .false. ! Logical switch for SST observations
1237	ln_sss = .false. ! Logical swithc for SSS observations
1238	ln_sic = .false. ! Logical switch for Sea Ice observations
1239	ln_vel3d = .false. ! Logical switch for velocity observations
1240	ln_altbias = .false. ! Logical switch for altimeter bias correction
1241	ln_sstbias = .false. ! Logical switch for SST bias correction
1242	ln_nea = .false. ! Logical switch for rejection of observations near land
1243	ln_grid_global = .true. ! Logical switch for global distribution of observations
1244	ln_grid_search_lookup = .false. ! Logical switch for obs grid search w/lookup table
1245	ln_ignmis = .true. ! Logical switch for ignoring missing files
1246	ln_s_at_t = .false. ! Logical switch for computing model S at T obs if not there
1247	ln_sstnight = .false. ! Logical switch for calculating night-time average for SST obs
1248	ln_bound_reject = .false. ! Logical to remove obs near boundaries in LAMs.
1249	ln_sla_fp_indegs = .true. ! Logical for SLA: T=> averaging footprint is in degrees, F=> in metres
1250	ln_sst_fp_indegs = .true. ! Logical for SST: T=> averaging footprint is in degrees, F=> in metres
1251	ln_sss_fp_indegs = .true. ! Logical for SSS: T=> averaging footprint is in degrees, F=> in metres
1252	ln_sic_fp_indegs = .true. ! Logical for SIC: T=> averaging footprint is in degrees, F=> in metres
1253	! All of the files variables below are arrays. Use namelist_cfg to add more files
1254	cn_profbfiles = 'profiles_01.nc' ! Profile feedback input observation file names
1255	cn_slafbfiles = 'sla_01.nc' ! SLA feedback input observation file names
1256	cn_sstfbfiles = 'sst_01.nc' ! SST feedback input observation file names
1257	cn_sssfbfiles = 'sss_01.nc' ! SSS feedback input observation file names
1258	cn_sicfbfiles = 'sic_01.nc' ! SIC feedback input observation file names
1259	cn_velfbfiles = 'vel_01.nc' ! Velocity feedback input observation file names
1260	cn_altbiasfile = 'altbias.nc' ! Altimeter bias input file name
1261	cn_sstbiasfiles = 'sstbias.nc' ! SST bias input file name
1262	cn_gridsearchfile ='gridsearch.nc' ! Grid search file name
1263	rn_gridsearchres = 0.5 ! Grid search resolution
1264	rn_mdtcorr = 1.61 ! MDT correction
1265	rn_mdtcutoff = 65.0 ! MDT cutoff for computed correction
1266	rn_dobsini = 00010101.000000 ! Initial date in window YYYYMMDD.HHMMSS
1267	rn_dobsend = 00010102.000000 ! Final date in window YYYYMMDD.HHMMSS
1268	rn_sla_avglamscl = 0. ! E/W diameter of SLA observation footprint (metres/degrees)
1269	rn_sla_avgphiscl = 0. ! N/S diameter of SLA observation footprint (metres/degrees)
1270	rn_sst_avglamscl = 0. ! E/W diameter of SST observation footprint (metres/degrees)
1271	rn_sst_avgphiscl = 0. ! N/S diameter of SST observation footprint (metres/degrees)
1272	rn_sss_avglamscl = 0. ! E/W diameter of SSS observation footprint (metres/degrees)
1273	rn_sss_avgphiscl = 0. ! N/S diameter of SSS observation footprint (metres/degrees)
1274	rn_sic_avglamscl = 0. ! E/W diameter of SIC observation footprint (metres/degrees)
1275	rn_sic_avgphiscl = 0. ! N/S diameter of SIC observation footprint (metres/degrees)
1276	nn_1dint = 0 ! Type of vertical interpolation method
1277	nn_2dint = 0 ! Default horizontal interpolation method
1278	nn_2dint_sla = 0 ! Horizontal interpolation method for SLA
1279	nn_2dint_sst = 0 ! Horizontal interpolation method for SST
1280	nn_2dint_sss = 0 ! Horizontal interpolation method for SSS
1281	nn_2dint_sic = 0 ! Horizontal interpolation method for SIC
1282	nn_msshc = 0 ! MSSH correction scheme
1283	nn_profdavtypes = -1 ! Profile daily average types - array
1284	/
1285	!-----------------------------------------------------------------------
1286	&nam_asminc ! assimilation increments ('key_asminc')
1287	!-----------------------------------------------------------------------
1288	ln_bkgwri = .false. ! Logical switch for writing out background state
1289	ln_trainc = .false. ! Logical switch for applying tracer increments
1290	ln_dyninc = .false. ! Logical switch for applying velocity increments
1291	ln_sshinc = .false. ! Logical switch for applying SSH increments
1292	ln_asmdin = .false. ! Logical switch for Direct Initialization (DI)
1293	ln_asmiau = .false. ! Logical switch for Incremental Analysis Updating (IAU)
1294	nitbkg = 0 ! Timestep of background in [0,nitend-nit000-1]
1295	nitdin = 0 ! Timestep of background for DI in [0,nitend-nit000-1]
1296	nitiaustr = 1 ! Timestep of start of IAU interval in [0,nitend-nit000-1]
1297	nitiaufin = 15 ! Timestep of end of IAU interval in [0,nitend-nit000-1]
1298	niaufn = 0 ! Type of IAU weighting function
1299	ln_salfix = .false. ! Logical switch for ensuring that the sa > salfixmin
1300	salfixmin = -9999 ! Minimum salinity after applying the increments
1301	nn_divdmp = 0 ! Number of iterations of divergence damping operator
1302	/
1303
1304	!!======================================================================
1305	!! * Miscellaneous namelists * !!
1306	!! !!
1307	!! nammpp Massively Parallel Processing ("key_mpp_mpi")
1308	!! namctl Control prints (default: OFF)
1309	!! namsto Stochastic parametrization of EOS (default: OFF)
1310	!!======================================================================
1311	!
1312	!-----------------------------------------------------------------------
1313	&nammpp ! Massively Parallel Processing ("key_mpp_mpi")
1314	!-----------------------------------------------------------------------
1315	ln_listonly = .false. ! do nothing else than listing the best domain decompositions (with land domains suppression)
1316	! ! if T: the largest number of cores tested is defined by max(mppsize, jpni*jpnj)
1317	ln_nnogather = .true. ! activate code to avoid mpi_allgather use at the northfold
1318	jpni = 0 ! number of processors following i (set automatically if < 1), see also ln_listonly = T
1319	jpnj = 0 ! number of processors following j (set automatically if < 1), see also ln_listonly = T
1320	/
1321	!-----------------------------------------------------------------------
1322	&namctl ! Control prints (default: OFF)
1323	!-----------------------------------------------------------------------
1324	ln_ctl = .FALSE. ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T
1325	sn_cfctl%l_config = .TRUE. ! IF .true. then control which reports are written with the following
1326	sn_cfctl%l_runstat = .FALSE. ! switches and which areas produce reports with the proc integer settings.
1327	sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure
1328	sn_cfctl%l_oceout = .FALSE. ! that all areas report.
1329	sn_cfctl%l_layout = .FALSE. !
1330	sn_cfctl%l_mppout = .FALSE. !
1331	sn_cfctl%l_mpptop = .FALSE. !
1332	sn_cfctl%procmin = 0 ! Minimum area number for reporting [default:0]
1333	sn_cfctl%procmax = 1000000 ! Maximum area number for reporting [default:1000000]
1334	sn_cfctl%procincr = 1 ! Increment for optional subsetting of areas [default:1]
1335	sn_cfctl%ptimincr = 1 ! Timestep increment for writing time step progress info
1336	nn_print = 0 ! level of print (0 no extra print)
1337	nn_ictls = 0 ! start i indice of control sum (use to compare mono versus
1338	nn_ictle = 0 ! end i indice of control sum multi processor runs
1339	nn_jctls = 0 ! start j indice of control over a subdomain)
1340	nn_jctle = 0 ! end j indice of control
1341	nn_isplt = 1 ! number of processors in i-direction
1342	nn_jsplt = 1 ! number of processors in j-direction
1343	ln_timing = .false. ! timing by routine write out in timing.output file
1344	ln_diacfl = .false. ! CFL diagnostics write out in cfl_diagnostics.ascii
1345	/
1346	!-----------------------------------------------------------------------
1347	&namsto ! Stochastic parametrization of EOS (default: OFF)
1348	!-----------------------------------------------------------------------
1349	ln_sto_eos = .false. ! stochastic equation of state
1350	nn_sto_eos = 1 ! number of independent random walks
1351	rn_eos_stdxy = 1.4 ! random walk horz. standard deviation (in grid points)
1352	rn_eos_stdz = 0.7 ! random walk vert. standard deviation (in grid points)
1353	rn_eos_tcor = 1440. ! random walk time correlation (in timesteps)
1354	nn_eos_ord = 1 ! order of autoregressive processes
1355	nn_eos_flt = 0 ! passes of Laplacian filter
1356	rn_eos_lim = 2.0 ! limitation factor (default = 3.0)
1357	ln_rststo = .false. ! start from mean parameter (F) or from restart file (T)
1358	ln_rstseed = .true. ! read seed of RNG from restart file
1359	cn_storst_in = "restart_sto" ! suffix of stochastic parameter restart file (input)
1360	cn_storst_out = "restart_sto" ! suffix of stochastic parameter restart file (output)
1361	/

Note: See TracBrowser for help on using the repository browser.

Download in other formats: