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 branches/2017/dev_merge_2017/NEMOGCM/CONFIG/SHARED – NEMO

Context Navigation

source: branches/2017/dev_merge_2017/NEMOGCM/CONFIG/SHARED/namelist_ref @ 9490

Last change on this file since 9490 was 9490, checked in by gm, 6 years ago
#2075 - dev_merge_2017: scale-aware setting of lateral viscous and diffusive coefficient
Property svn:keywords set to `Id` Property svn:mime-type set to `text/x-fortran`
File size: 101.2 KB

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

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

Download in other formats: