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/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/SHARED – NEMO

Context Navigation

source: NEMO/branches/2020/KERNEL-03_Storkey_Coward_RK3_stage2/cfgs/SHARED/namelist_ref @ 12397

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

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

Download in other formats: