!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> !! NEMO/OPA : 1 - run manager (namrun) !! namelists 2 - Domain (nam_zgr, nam_zgr_sco, namdom) !! 3 - Surface boundary (namsbc, namsbc_ana, namsbc_flx, namsbc_clio, namsbc_core !! namsbc_cpl, namqsr, namsbc_rnf, namsbc_ssr, namalb) !! 4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, namtide) !! 5 - bottom boundary (nambfr, nambbc, nambbl) !! 6 - Tracer (nameos, nam_traadv, nam_traldf, namtdp) !! 7 - dynamics (nam_dynadv, nam_dynvor, nam_dynhpg, namflg, nam_dynspg, nam_dynldf) !! 8 - Verical physics (namzdf, namnpc, namric, namtke, namkpp, namddm) !! 9 - diagnostics (namtrd, namgap, namspr, namflo, namptr) !! 9 - miscellaneous (namsol, nam_mpp, nam_mpp_dyndist, namctl) !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ! CAUTION: some scripts does not support CAPITALs for logical use .true./.false., not .TRUE./.FALSE. !!====================================================================== !! *** Run management namelists *** !!====================================================================== !! namrun parameters of the run !!====================================================================== !----------------------------------------------------------------------- &namrun ! parameters of the run !----------------------------------------------------------------------- no = 0 ! job number cexper = "GYRE" ! experience name ln_rstart = .false. ! start from rest (F) or from a restart file (T) nrstdt = 0 ! restart control = 0 nit000 is not compared to the restart file value ! = 1 use ndate0 in namelist (not the value in the restart file) ! = 2 calendar parameters read in the restart file nit000 = 1 ! first time step nitend = 4320 ! last time step ndate0 = 010101 ! initial calendar date yymmdd (used if nrstdt=1) nleapy = 30 ! Leap year calendar (1) or not (0) ninist = 0 ! output the initial state (1) or not (0) nstock = 4320 ! frequency of creation of a restart file nwrite = 60 ! frequency of write in the output file ln_dimgnnn = .false. ! DIMG file format: 1 file for all processors (F) or by processor (T) / !!====================================================================== !! *** Domain namelists *** !!====================================================================== !! nam_zgr vertical coordinate !! nam_zgr_sco s-coordinate or hybrid z-s-coordinate !! namdom space and time domain (bathymetry, mesh, timestep) !!====================================================================== !----------------------------------------------------------------------- &nam_zgr ! vertical coordinate !----------------------------------------------------------------------- ln_zco = .true. ! z-coordinate - full steps (T/F) ("key_zco" may also be defined) ln_zps = .false. ! z-coordinate - partial steps (T/F) ln_sco = .false. ! s- or hybrid z-s-coordinate (T/F) / !----------------------------------------------------------------------- &nam_zgr_sco ! s-coordinate or hybrid z-s-coordinate !----------------------------------------------------------------------- sbot_min = 300. ! minimum depth of s-bottom surface (>0) (m) sbot_max = 5250. ! maximum depth of s-bottom surface (= ocean depth) (>0) (m) theta = 6.0 ! surface control parameter (0<=theta<=20) thetb = 0.75 ! bottom control parameter (0<=thetb<= 1) r_max = 0.15 ! maximum cut-off r-value allowed (0 fill namsbc_ana ) ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) ln_blk_clio = .false. ! CLIO bulk formulation (T => fill namsbc_clio) ln_blk_core = .false. ! CORE bulk formulation (T => fill namsbc_core) ln_cpl = .false. ! Coupled formulation (T => fill namsbc_cpl ) nn_ice = 0 ! =0 no ice boundary condition , ! =1 use observed ice-cover , ! =2 ice-model used ("key_lim3" or "key_lim2) nn_ico_cpl = 0 ! ice-ocean coupling : =0 each nn_fsbc ! =1 stresses recomputed each ocean time step ("key_lim3" only) ! =2 combination of 0 and 1 cases ("key_lim3" only) ln_dm2dc = .false. ! daily mean to diurnal cycle short wave (qsr) ln_rnf = .false. ! runoffs (T => fill namsbc_rnf) ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) nn_fwb = 0 ! FreshWater Budget: =0 unchecked , ! =1 annual global mean of e-p-r set to zero , ! =2 global mean of e-p-r set to zero at each nn_fsbc time step / !----------------------------------------------------------------------- &namsbc_ana ! analytical surface boundary condition !----------------------------------------------------------------------- nn_tau000 = 100 ! gently increase the stress over the first ntau_rst time-steps rn_utau0 = 0.1e0 ! uniform value for the i-stress rn_vtau0 = 0.e0 ! uniform value for the j-stress rn_q0 = 0.e0 ! uniform value for the total heat flux rn_qsr0 = 0.e0 ! uniform value for the solar radiation rn_emp0 = 0.e0 ! uniform value for the freswater budget (E-P) / !----------------------------------------------------------------------- &namsbc_flx ! surface boundary condition : flux formulation !----------------------------------------------------------------------- ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! sn_utau = 'utau' , 24. , 'utau' , .false. , .false. , 'yearly' sn_vtau = 'vtau' , 24. , 'vtau' , .false. , .false. , 'yearly' sn_qtot = 'qtot' , 24. , 'qtot' , .false. , .false. , 'yearly' sn_qsr = 'qsr' , 24. , 'qsr' , .false. , .false. , 'yearly' sn_emp = 'emp' , 24. , 'emp' , .false. , .false. , 'yearly' ! cn_dir = './' ! root directory for the location of the flux files / !----------------------------------------------------------------------- &namsbc_clio ! namsbc_clio CLIO bulk formulea !----------------------------------------------------------------------- ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! sn_utau = 'taux_1m' , -1. , 'sozotaux' , .false. , .true. , 'yearly' sn_vtau = 'tauy_1m' , -1. , 'sometauy' , .false. , .true. , 'yearly' sn_wndm = 'flx' , -1. , 'socliowi' , .false. , .true. , 'yearly' sn_tair = 'flx' , -1. , 'socliot2' , .false. , .true. , 'yearly' sn_humi = 'flx' , -1. , 'socliohu' , .false. , .true. , 'yearly' sn_ccov = 'flx' , -1. , 'socliocl' , .false. , .true. , 'yearly' sn_prec = 'flx' , -1. , 'socliopl' , .false. , .true. , 'yearly' ! cn_dir = './' ! root directory for the location of the bulk files are / !----------------------------------------------------------------------- &namsbc_core ! namsbc_core CORE bulk formulea !----------------------------------------------------------------------- ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! sn_wndi = 'u10' , 24. , 'U_10_MOD' , .false. , .true. , 'yearly' sn_wndj = 'v10' , 24. , 'V_10_MOD' , .false. , .true. , 'yearly' sn_qsr = 'rad' , 24. , 'SWDN_MOD' , .false. , .true. , 'yearly' sn_qlw = 'rad' , 24. , 'LWDN_MOD' , .false. , .true. , 'yearly' sn_tair = 't10' , 24. , 'T_10_MOD' , .false. , .true. , 'yearly' sn_humi = 'q10' , 24. , 'Q_10_MOD' , .false. , .true. , 'yearly' sn_prec = 'precip' , -1. , 'RAIN' , .false. , .true. , 'yearly' sn_snow = 'precip' , -1. , 'SNOW' , .false. , .true. , 'yearly' ! cn_dir = './' ! root directory for the location of the bulk files ln_2m = .false. ! air temperature and humidity referenced at 2m (T) instead 10m (F) alpha_precip= 1. ! multiplicative factor for precipitation (total & snow) / !----------------------------------------------------------------------- &namsbc_cpl ! coupled ocean/atmosphere model ("key_coupled") !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &namqsr ! penetrative solar radiation !----------------------------------------------------------------------- ln_traqsr = .true. ! penetrative solar radiation (T) or not (F) rabs = 0.58 ! fraction of qsr associated with xsi1 xsi1 = 0.35 ! first depth of extinction xsi2 = 23.0 ! second depth of extinction / !----------------------------------------------------------------------- &namsbc_rnf ! runoffs namelist surface boundary condition !----------------------------------------------------------------------- ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! sn_rnf = 'runoff_1m_nomask' , -1. , 'sorunoff' , .true. , .true. , 'yearly' sn_cnf = 'runoff_1m_nomask' , 0. , 'socoefr' , .false. , .true. , 'yearly' ! cn_dir = './' ! root directory for the location of the runoff files ln_rnf_emp = .false. ! runoffs included into precipitation field (T) or into a file (F) ln_rnf_mouth = .false. ! specific treatment at rivers mouths rn_hrnf = 0.e0 ! depth over which enhanced vertical mixing is used rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] / !----------------------------------------------------------------------- &namsbc_ssr ! surface boundary condition : sea surface restoring !----------------------------------------------------------------------- ! ! file name ! frequency (hours) ! variable ! time interpol. ! clim ! 'yearly' or ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! sn_sst = 'sst_data' , 24. , 'sst' , .false. , .false. , 'yearly' sn_sss = 'sss_data' , -1. , 'sss' , .true. , .false. , 'yearly' ! cn_dir = './' ! root directory for the location of the runoff files nn_sstr = 0 ! add a retroaction term in the surface heat flux (=1) or not (=0) nn_sssr = 0 ! add a damping term in the surface freshwater flux (=1) or not (=0) dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K] deds = -27.7 ! magnitude of the damping on salinity [mm/day/psu] / !----------------------------------------------------------------------- &namalb ! albedo parameters !----------------------------------------------------------------------- cgren = 0.06 ! correction of the snow or ice albedo to take into account the albice = 0.5 ! albedo of melting ice in the arctic and antarctic alphd = 0.80 ! coefficients for linear interpolation used to alphc = 0.65 ! compute albedo between two extremes values alphdi = 0.72 ! (Pyane, 1972) / !!====================================================================== !! *** Lateral boundary condition *** !!====================================================================== !! namlbc lateral momentum boundary condition !! namcla cross land advection !! namobc open boundaries parameters ("key_obc") !! namagrif agrif nested grid ( read by child model only ) ("key_agrif") !! nambdy Unstructured open boundaries ("key_bdy") !! namtide Tidal forcing at open boundaries ("key_bdy_tides") !!====================================================================== !----------------------------------------------------------------------- &namlbc ! lateral momentum boundary condition !----------------------------------------------------------------------- shlat = 0. ! shlat = 0 : free slip ! 0 < shlat < 2 : partial slip ! shlat = 2 : no slip ! 2 < shlat : strong slip / !----------------------------------------------------------------------- &namcla ! cross land advection !----------------------------------------------------------------------- n_cla = 0 ! advection between 2 ocean pts separates by land / !----------------------------------------------------------------------- &namobc ! open boundaries parameters ("key_obc") !----------------------------------------------------------------------- nobc_dta = 1 ! = 0 the obc data are equal to the initial state ! = 1 the obc data are read in 'obc.dta' files cffile = 'annual' ! set to annual if obc datafile hold 1 year of data ! set to monthly if obc datafile hold 1 month of data rdpein = 1. ! ??? rdpwin = 1. ! ??? rdpnin = 1. ! ??? rdpsin = 1. ! ??? rdpeob = 3000. ! time relaxation (days) for the east open boundary rdpwob = 15. ! " " " west " rdpnob = 3000. ! " " " north " rdpsob = 15. ! " " " south " zbsic1 = 140.e+6 ! barotropic stream function on first isolated coastline zbsic2 = 1.e+6 ! " " second " zbsic3 = 0. ! " " thrid " ln_obc_clim= .false. ! climatological obc data files (T) or not (F) ln_vol_cst = .true. ! impose the total volume conservation (T) or not (F) / !----------------------------------------------------------------------- &namagrif ! ("key_agrif") !----------------------------------------------------------------------- nbclineupdate = 3 ! baroclinic update frequency ln_spc_dyn = .true. ! use 0 as special value for dynamics visc_tra = 2880. ! viscosity coeeficient for tracers sponge layer visc_dyn = 2880. ! viscosity coeeficient for dynamics sponge layer / !----------------------------------------------------------------------- &nambdy ! unstructured open boundaries parameters ("key_bdy") !----------------------------------------------------------------------- filbdy_mask = '' ! name of mask file (if ln_bdy_mask=.TRUE.) filbdy_data_T = 'bdydata_grid_T.nc' ! name of data file (T-points) filbdy_data_U = 'bdydata_grid_U.nc' ! name of data file (U-points) filbdy_data_V = 'bdydata_grid_V.nc' ! name of data file (V-points) ln_bdy_clim = .false. ! contain 1 (T) or 12 (F) time dumps and be cyclic ln_bdy_vol = .true. ! total volume correction (see volbdy parameter) ln_bdy_mask = .false. ! boundary mask from filbdy_mask (T) or boundaries are on edges of domain (F) ln_bdy_tides = .true. ! Apply tidal harmonic forcing with Flather condition ln_bdy_dyn_fla = .true. ! Apply Flather condition to velocities ln_bdy_tra_frs = .false. ! Apply FRS condition to temperature and salinity ln_bdy_dyn_frs = .false. ! Apply FRS condition to velocities nbdy_dta = 1 ! = 0, bdy data are equal to the initial state ! = 1, bdy data are read in 'bdydata .nc' files nb_rimwidth = 9 ! width of the relaxation zone volbdy = 0 ! = 0, the total water flux across open boundaries is zero ! = 1, the total volume of the system is conserved / !----------------------------------------------------------------------- &namtide ! tidal forcing at unstructured boundaries !----------------------------------------------------------------------- filtide = 'bdytide_' ! file name root of tidal forcing files tide_cpt = 'M2','S1' ! names of tidal components used tide_speed = 28.984106, 15.000001 ! phase speeds of tidal components (deg/hour) ln_tide_date = .false. ! adjust tidal harmonics for start date of run / !!====================================================================== !! *** Bottom boundary condition *** !!====================================================================== !! nambfr bottom friction !! nambbc bottom temperature boundary condition ("key_trabbc") !! nambbl bottom boundary layer scheme ("key_trabbl_dif","key_trabbl_adv") !!====================================================================== !----------------------------------------------------------------------- &nambfr ! bottom friction !----------------------------------------------------------------------- nbotfr = 2 ! type of bottom friction : = 0 : no slip, = 2 : nonlinear friction ! = 3 : free slip, = 1 : linear friction bfri1 = 4.e-4 ! bottom drag coefficient (linear case) bfri2 = 1.e-3 ! bottom drag coefficient (non linear case) bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m^2/s^2) / !----------------------------------------------------------------------- &nambbc ! bottom temperature boundary condition !----------------------------------------------------------------------- ngeo_flux = 0 ! geothermal heat flux = 0 no flux considered ! = 1 constant flux ! = 2 variable flux (read in geothermal_heating.nc in mW/m2) ngeo_flux_const = 86.4e-3 ! Constant value of geothermal heat flux [W/m2] / !----------------------------------------------------------------------- &nambbl ! bottom boundary layer scheme !----------------------------------------------------------------------- ! ! diffusive bbl ("key_trabbl") ! ! advective bbl ("key_trabbl_adv") atrbbl = 10000. ! lateral mixing coefficient in the bbl [m2/s] / !!====================================================================== !! Tracer (T & S ) namelists !!====================================================================== !! nameos equation of state !! nam_traadv advection scheme !! nam_traldf lateral diffusion scheme !! namtdp tracer newtonian damping ("key_tradmp") !!====================================================================== !----------------------------------------------------------------------- &nameos ! ocean physical parameters !----------------------------------------------------------------------- neos = 2 ! type of equation of state and Brunt-Vaisala frequency ! = 0, UNESCO (formulation of Jackett and McDougall (1994) and of McDougall (1987) ) ! = 1, linear: rho(T) = rau0 * ( 1.028 - ralpha * T ) ! = 2, linear: rho(T,S) = rau0 * ( rbeta * S - ralpha * T ) ralpha = 2.e-4 ! thermal expension coefficient (neos= 1 or 2) rbeta = 7.7e-4 ! saline expension coefficient (neos= 2) / !----------------------------------------------------------------------- &nam_traadv ! advection scheme for tracer !----------------------------------------------------------------------- ln_traadv_cen2 = .false. ! 2nd order centered scheme ln_traadv_tvd = .true. ! TVD scheme ln_traadv_muscl = .false. ! MUSCL scheme ln_traadv_muscl2 = .false. ! MUSCL2 scheme + cen2 at boundaries ln_traadv_ubs = .false. ! UBS scheme / !----------------------------------------------------------------------- &nam_traldf ! lateral diffusion scheme for tracer !----------------------------------------------------------------------- ! ! Type of the operator : ln_traldf_lap = .true. ! laplacian operator ln_traldf_bilap = .false. ! bilaplacian operator ! Direction of action : ln_traldf_level = .false. ! iso-level ln_traldf_hor = .false. ! horizontal (geopotential) (require "key_ldfslp" when ln_sco=T) ln_traldf_iso = .true. ! iso-neutral (require "key_ldfslp") ! ! Coefficient aht0 = 1000. ! horizontal eddy diffusivity for tracers [m2/s] ahtb0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] aeiv0 = 1000. ! eddy induced velocity coefficient [m2/s] (require "key_traldf_eiv") / !----------------------------------------------------------------------- &namtdp ! tracer newtonian damping ('key_tradmp') !----------------------------------------------------------------------- ndmp = -1 ! type of damping in temperature and salinity ! ='latitude', damping poleward of 'ndmp' degrees and function ! of the distance-to-coast. Red and Med Seas as ndmp=-1 ! =-1 damping only in Med and Red Seas ndmpf = 1 ! create a damping.coeff NetCDF file (=1) or not (=0) nmldmp = 1 ! type of damping: =0 damping throughout the water column ! =1 no damping in the mixed layer defined by avt >5cm2/s ) ! =2 no damping in the mixed layer defined rho /seconds ; =86400. -> /day) ln_trdmld_restart = .false. ! restart for ML diagnostics ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S / !----------------------------------------------------------------------- &namgap ! level mean model-data gap ('key_diagap') !----------------------------------------------------------------------- ngap = 15 ! time-step frequency of model-data gap computation nprg = 10 ! time-step frequency of gap print in model output / !----------------------------------------------------------------------- &namspr ! surface pressure diagnostic !----------------------------------------------------------------------- nmaxp = 1000 ! maximum of iterations for the solver epsp = 1.e-3 ! absolute precision of the solver niterp = 400 ! number of iteration done by the solver / !----------------------------------------------------------------------- &namflo ! float parameters ("key_float") !----------------------------------------------------------------------- ln_rstflo = .false. ! float restart (T) or not (F) nwritefl = 75 ! frequency of writing in float output file nstockfl = 5475 ! frequency of creation of the float restart file ln_argo = .false. ! Argo type floats (stay at the surface each 10 days) ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T) ! or computed with Blanke' scheme (F) / !----------------------------------------------------------------------- &namptr ! Poleward Transport Diagnostic !----------------------------------------------------------------------- ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F) ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not ! (orca configuration only, need input basins mask file named "subbasins.nc" nf_ptr = 15 ! Frequency of ptr computation [time step] /