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/2016/dev_r6522_SIMPLIF_3/NEMOGCM/CONFIG/SHARED – NEMO

Context Navigation

source: branches/2016/dev_r6522_SIMPLIF_3/NEMOGCM/CONFIG/SHARED/namelist_ref @ 6862

Last change on this file since 6862 was 6862, checked in by lovato, 8 years ago
#1729 - trunk: removed key_bdy from the code and set usage of ln_bdy. Tested with SETTE.
Property svn:keywords set to `Id`
File size: 92.6 KB

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

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

Download in other formats:

Original Format