!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> !! 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 = "ORCA2" ! 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 = 5475 ! last time step ndate0 = 010101 ! initial calendar date yymmdd (used if nrstdt=1) nleapy = 0 ! Leap year calendar (1) or not (0) ninist = 0 ! output the initial state (1) or not (0) nstock = 5475 ! frequency of creation of a restart file nwrite = 5475 ! 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 = .false. ! z-coordinate - full steps (T/F) ("key_zco" may also be defined) ln_zps = .true. ! 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 = .true. ! 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 = 2 ! =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_ssr) 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 = 0 ! gently increase the stress over the first ntau_rst time-steps rn_utau0 = 0.5 ! 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 ! starting ! ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! sn_utau = 'utau.nc' , 24. , 'utau' , .false. , 0 , 0 sn_vtau = 'vtau.nc' , 24. , 'vtau' , .false. , 0 , 0 sn_qtot = 'qtot.nc' , 24. , 'qtot' , .false. , 0 , 0 sn_qsr = 'qsr.nc' , 24. , 'qsr' , .false. , 0 , 0 sn_emp = 'emp.nc' , 24. , 'emp' , .false. , 0 , 0 ! 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 ! starting ! ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! sn_utau = 'taux_1m' , -12. , 'sozotaux' , .false. , 1 , 0 sn_vtau = 'tauy_1m' , -12. , 'sometauy' , .false. , 1 , 0 sn_wndm = 'flx' , -12. , 'socliowi' , .false. , 1 , 0 sn_tair = 'flx' , -12. , 'socliot2' , .false. , 1 , 0 sn_humi = 'flx' , -12. , 'socliohu' , .false. , 1 , 0 sn_ccov = 'flx' , -12. , 'socliocl' , .false. , 1 , 0 sn_prec = 'flx' , -12. , 'socliopl' , .false. , 1 , 0 ! 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 ! starting ! ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! sn_wndi = 'u10' , 24. , 'U_10_MOD' , .false. , 1 , 0 sn_wndj = 'v10' , 24. , 'V_10_MOD' , .false. , 1 , 0 sn_qsr = 'rad' , 24. , 'SWDN_MOD' , .false. , 1 , 0 sn_qlw = 'rad' , 24. , 'LWDN_MOD' , .false. , 1 , 0 sn_tair = 't10' , 24. , 'T_10_MOD' , .false. , 1 , 0 sn_humi = 'q10' , 24. , 'Q_10_MOD' , .false. , 1 , 0 sn_prec = 'precip' , -12. , 'RAIN' , .false. , 1 , 0 sn_snow = 'precip' , -12. , 'SNOW' , .false. , 1 , 0 ! 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 ! starting ! ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! sn_rnf = 'runoff_1m_nomask.nc', -12. , 'sorunoff' , .true. , 1 , 0 sn_cnf = 'runoff_1m_nomask.nc', 0. , 'socoefr' , .false. , 1 , 0 ! cn_dir = './' ! root directory for the location of the runoff files nn_runoff = 2 ! no runoff (0), runoff (1), runoff+river mouth ups adv with cen2 (2) 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 ! starting ! ! ! ! (if <0 months) ! name ! (logical) ! (0/1) ! record ! sn_sst = 'sst_data.nc' , 24. , 'sst' , .false. , 0 , 0 sn_sss = 'sss_data.nc' , -12. , 'sss' , .true. , 0 , 0 ! nn_sstr =0/1 add a retroaction term in the surface heat flux ! nn_sssr =O/1/2 add a damping term in the surface freshwater flux ! (=1, salt flux, concentration/dillution emps only) ! (=2, volume flux, both emp and emps are updated) ! dqdt magnitude of the retroaction on temperature [W/m2/K] ! deds magnitude of the damping on salinity [mm/day] ! 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 = 2. ! shlat = 0 : free slip ! 0 < shlat < 2 : partial slip ! shlat = 2 : no slip ! 2 < shlat : strong slip / !----------------------------------------------------------------------- &namcla ! cross land advection !----------------------------------------------------------------------- n_cla = 1 ! advection between 2 ocean pts separates by land / !----------------------------------------------------------------------- &namobc ! open boundaries parameters ("key_obc") !----------------------------------------------------------------------- nobc_dta = 0 ! = 0 the obc data are equal to the initial state ! = 1 the obc data are read in 'obc.dta' files rdpein = 1. ! ??? rdpwin = 1. ! ??? rdpnin = 30. ! ??? rdpsin = 1. ! ??? rdpeob = 1500. ! time relaxation (days) for the east open boundary rdpwob = 15. ! " " " west " rdpnob = 150. ! " " " 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= .true. ! climatological obc data files (T) or not (F) ln_vol_cst = .false. ! 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 = 'bdymask_grid_T.nc' ! name of mask file (at T-points) filbdy_data_T = 'bdydata_grid_T.nc' ! name of data file (at T-points) filbdy_data_U = 'bdydata_grid_U.nc' ! name of data file (at U-points) filbdy_data_V = 'bdydata_grid_V.nc' ! name of data file (at V-points) ln_bdy_clim = .false. ! must contain 1 or 12 time dumps and be cyclic ln_bdy_vol = .false. ! total volume correction (see volbdy parameter) ln_bdy_fla = .true. ! flather boundary conditions nbdy_dta = 1 ! = 0, bdy data are equal to the initial state ! = 1, bdy data are read in 'bdydata .nc' files nb_rimwidth = 1 ! width of the relaxation zone volbdy = 1 ! = 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 ("key_bdy_tides") !----------------------------------------------------------------------- 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 = .true. ! 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 = 1 ! 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 = 2 ! 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 = 0 ! 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 = 0.001 ! saline expension coefficient (neos= 2) / !----------------------------------------------------------------------- &nam_traadv ! advection scheme for tracer !----------------------------------------------------------------------- ln_traadv_cen2 = .true. ! 2nd order centered scheme ln_traadv_tvd = .false. ! 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 = 2000. ! horizontal eddy diffusivity for tracers [m2/s] ahtb0 = 0. ! background eddy diffusivity for ldf_iso [m2/s] aeiv0 = 2000. ! 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] /