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 @ 9027

Last change on this file since 9027 was 9024, checked in by timgraham, 6 years ago
Resolved conflicts in namelists
Property svn:keywords set to `Id` Property svn:mime-type set to `text/x-fortran`
File size: 91.2 KB

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

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

Download in other formats: