!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> !! NEMO/OPA Configuration namelist : used to overwrite defaults values defined in SHARED/namelist_ref !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ! !----------------------------------------------------------------------- &namusr_def !----------------------------------------------------------------------- rn_dx = 1000.0 rn_dz = 1.0 nn_wad_test = 1 / !----------------------------------------------------------------------- &namrun ! parameters of the run !----------------------------------------------------------------------- cn_exp = "WAD" ! experience name nn_it000 = 1 ! first time step nn_itend = 3840 ! last time step !nn_itend = 6 ! last time step nn_leapy = 30 ! Leap year calendar (1) or not (0) nn_stock = 48000 ! frequency of creation of a restart file (modulo referenced to 1) ln_clobber = .true. ! clobber (overwrite) an existing file nn_istate = 0 ! output the initial state (1) or not (0) / !----------------------------------------------------------------------- &namcfg ! parameters of the configuration !----------------------------------------------------------------------- ln_read_cfg = .false. ! (=T) read the domain configuration file ! ! (=F) user defined configuration ==>>> see usrdef(_...) modules ln_write_cfg= .true. ! (=T) create the domain configuration file / !----------------------------------------------------------------------- &namdom ! space and time domain (bathymetry, mesh, timestep) !----------------------------------------------------------------------- ln_linssh = .false. ! =T linear free surface ==>> model level are fixed in time ! nn_msh = 1 ! create (>0) a mesh file or not (=0) rn_rdt = 18. ! time step for the dynamics / !----------------------------------------------------------------------- &namcrs ! coarsened grid (for outputs and/or TOP) (ln_crs =T) !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &namtsd ! data : Temperature & Salinity !----------------------------------------------------------------------- cn_dir = './' ! root directory for the location of the runoff files ln_tsd_init = .false. ! Initialisation of ocean T & S with T &S input data (T) or not (F) ln_tsd_tradmp = .false. ! damping of ocean T & S toward T &S input data (T) or not (F) / !----------------------------------------------------------------------- &namsbc ! Surface Boundary Condition (surface module) !----------------------------------------------------------------------- nn_fsbc = 1 ! frequency of surface boundary condition computation ! ! (also = the frequency of sea-ice model call) ln_usr = .true. ! analytical formulation (T => check usrdef_sbc) ln_blk = .false. ! CORE bulk formulation (T => fill namsbc_blk ) nn_ice = 0 ! =0 no ice boundary condition , 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 / !----------------------------------------------------------------------- &namsbc_cpl ! coupled ocean/atmosphere model ("key_oasis3") !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &namtra_qsr ! penetrative solar radiation !----------------------------------------------------------------------- ln_qsr_rgb = .false. ! RGB (Red-Green-Blue) light penetration ln_qsr_2bd = .true. ! 2 bands light penetration nn_chldta = 0 ! RGB : Chl data (=1) or cst value (=0) / !----------------------------------------------------------------------- &namsbc_rnf ! runoffs namelist surface boundary condition !----------------------------------------------------------------------- ln_rnf_mouth = .false. ! specific treatment at rivers mouths / !----------------------------------------------------------------------- &namsbc_apr ! Atmospheric pressure used as ocean forcing or in bulk !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &namsbc_ssr ! surface boundary condition : sea surface restoring !----------------------------------------------------------------------- nn_sssr = 0 ! add a damping term in the surface freshwater flux (=2) rn_deds = -27.7 ! magnitude of the damping on salinity [mm/day] ln_sssr_bnd = .false. ! flag to bound erp term (associated with nn_sssr=2) / !----------------------------------------------------------------------- &namsbc_alb ! albedo parameters !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &namberg ! iceberg parameters !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &namlbc ! lateral momentum boundary condition !----------------------------------------------------------------------- rn_shlat = 0. ! shlat = 0 ! 0 < shlat < 2 ! shlat = 2 ! 2 < shlat / !----------------------------------------------------------------------- &namagrif ! AGRIF zoom ("key_agrif") !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &nam_tide ! tide parameters (#ifdef key_tide) !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &nambdy ! unstructured open boundaries !----------------------------------------------------------------------- ln_bdy = .false. nb_bdy = 0 ! number of open boundary sets ln_coords_file = .false. ! =T : read bdy coordinates from file cn_coords_file = 'coordinates.bdy.nc' ! bdy coordinates files ln_mask_file = .false. ! =T : read mask from file cn_mask_file = '' ! name of mask file (if ln_mask_file=.TRUE.) cn_dyn2d = 'flather' ! nn_dyn2d_dta = 1 ! = 0, bdy data are equal to the initial state ! = 1, bdy data are read in 'bdydata .nc' files ! = 2, use tidal harmonic forcing data from files ! = 3, use external data AND tidal harmonic forcing cn_dyn3d = 'none' ! nn_dyn3d_dta = 0 ! = 0, bdy data are equal to the initial state ! = 1, bdy data are read in 'bdydata .nc' files cn_tra = 'frs' ! nn_tra_dta = 0 ! = 0, bdy data are equal to the initial state ! = 1, bdy data are read in 'bdydata .nc' files cn_ice_lim = 'none' ! nn_ice_lim_dta = 0 ! = 0, bdy data are equal to the initial state ! = 1, bdy data are read in 'bdydata .nc' files rn_ice_tem = 270. ! lim3 only: arbitrary temperature of incoming sea ice rn_ice_sal = 10. ! lim3 only: -- salinity -- rn_ice_age = 30. ! lim3 only: -- age -- ln_tra_dmp =.false. ! open boudaries conditions for tracers ln_dyn3d_dmp =.false. ! open boundary condition for baroclinic velocities rn_time_dmp = 1. ! Damping time scale in days rn_time_dmp_out = 1. ! Outflow damping time scale nn_rimwidth = 10 ! width of the relaxation zone ln_vol = .false. ! total volume correction (see nn_volctl parameter) nn_volctl = 1 ! = 0, the total water flux across open boundaries is zero / !----------------------------------------------------------------------- &nambdy_index !----------------------------------------------------------------------- ctypebdy = 'E' nbdyind = 50 nbdybeg = 1 nbdyend = 34 !ctypebdy = 'W' !nbdyind = 2 !nbdybeg = 1 !nbdyend = 34 / !----------------------------------------------------------------------- &nambdy_dta ! open boundaries - external data !----------------------------------------------------------------------- ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename ! bn_ssh = 'bdyssh_tc7' , 1 , 'sshbdy', .true. , .true. , 'daily' , '' , '' , '' bn_u2d = 'bdyuv_tc7' , 1 , 'ubdy' , .true. , .true. , 'daily' , '' , '' , '' bn_v2d = 'bdyuv_tc7' , 1 , 'vbdy' , .true. , .true. , 'daily' , '' , '' , '' cn_dir = './' ! root directory for the location of the bulk files ln_full_vel = .false. ! / !----------------------------------------------------------------------- &nambdy_tide ! tidal forcing at open boundaries !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &namdrg ! top/bottom drag coefficient (default: NO selection) !----------------------------------------------------------------------- ln_loglayer= .false. ! logarithmic drag: Cd = vkarmn/log(z/z0) |U| / !----------------------------------------------------------------------- &namdrg_bot ! BOTTOM friction !----------------------------------------------------------------------- rn_Cd0 = 1.e-4 ! drag coefficient [-] rn_Uc0 = 0.1 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) rn_Cdmax = 1.e-4 ! drag value maximum [-] (logarithmic drag) rn_ke0 = 2.5e-3 ! background kinetic energy [m2/s2] (non-linear cases) rn_z0 = 3.e-3 ! roughness [m] (ln_loglayer=T) ln_boost = .false. ! =T regional boost of Cd0 ; =F constant rn_boost= 50. ! local boost factor [-] / !----------------------------------------------------------------------- &nambbc ! bottom temperature boundary condition !----------------------------------------------------------------------- ln_trabbc = .false. ! Apply a geothermal heating at the ocean bottom nn_geoflx = 0 ! geothermal heat flux: = 0 no flux / !----------------------------------------------------------------------- &nambbl ! bottom boundary layer scheme !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &nameos ! ocean physical parameters !----------------------------------------------------------------------- ln_teos10 = .false. ! = Use TEOS-10 equation of state ln_eos80 = .false. ! = Use EOS80 equation of state ln_seos = .true. ! = Use simplified equation of state (S-EOS) ! ! ! S-EOS coefficients (ln_seos=T): ! ! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS rn_a0 = 1.6550e-1 ! thermal expension coefficient (nn_eos= 1) rn_b0 = 7.6554e-1 ! saline expension coefficient (nn_eos= 1) rn_lambda1 = 5.9520e-2 ! cabbeling coeff in T^2 (=0 for linear eos) rn_lambda2 = 7.4914e-4 ! cabbeling coeff in S^2 (=0 for linear eos) rn_mu1 = 1.4970e-4 ! thermobaric coeff. in T (=0 for linear eos) rn_mu2 = 1.1090e-5 ! thermobaric coeff. in S (=0 for linear eos) rn_nu = 2.4341e-3 ! cabbeling coeff in T*S (=0 for linear eos) !!org GYRE rn_alpha = 2.0e-4 ! thermal expension coefficient (nn_eos= 1 or 2) !!org GYRE rn_beta = 7.7e-4 ! saline expension coefficient (nn_eos= 2) !!org caution now a0 = alpha / rau0 with rau0 = 1026 / !----------------------------------------------------------------------- &namtra_adv ! advection scheme for tracer (default: No selection) !----------------------------------------------------------------------- ln_traadv_NONE= .false. ! No tracer advection ln_traadv_cen = .false. ! 2nd order centered scheme ln_traadv_mus = .false. ! MUSCL scheme ln_traadv_fct = .true. ! FCT scheme nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order / !----------------------------------------------------------------------- &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) !----------------------------------------------------------------------- / !---------------------------------------------------------------------------------- &namtra_ldf ! lateral diffusion scheme for tracers !---------------------------------------------------------------------------------- ! ! Operator type: ln_traldf_lap = .true. ! laplacian operator ln_traldf_blp = .false. ! bilaplacian operator ! ! Direction of action: ln_traldf_lev = .false. ! iso-level ln_traldf_hor = .false. ! horizontal (geopotential) ln_traldf_iso = .true. ! iso-neutral ln_traldf_triad = .false. ! iso-neutral using Griffies triads ! ! ! iso-neutral options: ln_traldf_msc = .false. ! Method of Stabilizing Correction (both operators) rn_slpmax = 0.01 ! slope limit (both operators) ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only) rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only) ln_botmix_triad = .false. ! lateral mixing on bottom (triad only) ! ! ! Coefficients: nn_aht_ijk_t = 0 ! space/time variation of eddy coef ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file ! ! = 0 constant ! ! = 10 F(k) =ldf_c1d ! ! = 20 F(i,j) =ldf_c2d ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d ! ! = 31 F(i,j,k,t)=F(local velocity) rn_aht_0 = 1000. ! lateral eddy diffusivity (lap. operator) [m2/s] rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s] / !---------------------------------------------------------------------------------- &namtra_ldfeiv ! eddy induced velocity param. !---------------------------------------------------------------------------------- ln_ldfeiv =.false. ! use eddy induced velocity parameterization / !----------------------------------------------------------------------- &namtra_dmp ! tracer: T & S newtonian damping !----------------------------------------------------------------------- ln_tradmp = .false. ! add a damping termn (T) or not (F) / !----------------------------------------------------------------------- &namdyn_adv ! formulation of the momentum advection (default: NO selection) !----------------------------------------------------------------------- ln_dynadv_NONE= .false. ! linear dynamics (no momentum advection) ln_dynadv_vec = .true. ! vector form - 2nd centered scheme nn_dynkeg = 0 ! grad(KE) scheme: =0 C2 ; =1 Hollingsworth correction ln_dynadv_cen2= .false. ! flux form - 2nd order centered scheme ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &namdyn_vor ! option of physics/algorithm (not control by CPP keys) !----------------------------------------------------------------------- ln_dynvor_ene = .false. ! enstrophy conserving scheme ln_dynvor_ens = .false. ! energy conserving scheme ln_dynvor_mix = .false. ! mixed scheme ln_dynvor_een = .true. ! energy & enstrophy scheme nn_een_e3f = 0 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) / !----------------------------------------------------------------------- &namdyn_hpg ! Hydrostatic pressure gradient option !----------------------------------------------------------------------- ln_hpg_zco = .false. ! z-coordinate - full steps ln_hpg_zps = .false. ! z-coordinate - partial steps (interpolation) ln_hpg_sco = .true. ! s-coordinate / !----------------------------------------------------------------------- &namdyn_spg ! surface pressure gradient !----------------------------------------------------------------------- ln_dynspg_ts = .true. ! split-explicit free surface ln_bt_auto = .false. ! Number of sub-step defined from: nn_baro = 12 ! =F : the number of sub-step in rn_rdt seconds / !----------------------------------------------------------------------- &namdyn_ldf ! lateral diffusion on momentum !----------------------------------------------------------------------- ! ! Type of the operator : ! ! no diffusion: set ln_dynldf_lap=..._blp=F ln_dynldf_lap = .true. ! laplacian operator ln_dynldf_blp = .false. ! bilaplacian operator ! ! Direction of action : ln_dynldf_lev = .true. ! iso-level ln_dynldf_hor = .false. ! horizontal (geopotential) ln_dynldf_iso = .false. ! iso-neutral ! ! Coefficient nn_ahm_ijk_t = 0 ! space/time variation of eddy coef ! ! =-30 read in eddy_viscosity_3D.nc file ! ! =-20 read in eddy_viscosity_2D.nc file ! ! = 0 constant ! ! = 10 F(k)=c1d ! ! = 20 F(i,j)=F(grid spacing)=c2d ! ! = 30 F(i,j,k)=c2d*c1d ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) rn_ahm_0 = 1000. ! horizontal laplacian eddy viscosity [m2/s] rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s] rn_bhm_0 = 0. ! horizontal bilaplacian eddy viscosity [m4/s] / !!====================================================================== !! vertical physics namelists !! !!====================================================================== !! namzdf vertical physics !! namzdf_ric richardson number vertical mixing (ln_zdfric=T) !! namzdf_tke TKE vertical mixing (ln_zdftke=T) !! namzdf_gls GLS vertical mixing (ln_zdfgls=T) !! namzdf_iwm tidal mixing parameterization (ln_zdfiwm=T) !!====================================================================== !----------------------------------------------------------------------- &namzdf ! vertical physics (default: NO selection) !----------------------------------------------------------------------- ! ! type of vertical closure ln_zdfcst = .false. ! constant mixing ln_zdfric = .false. ! local Richardson dependent formulation (T => fill namzdf_ric) ln_zdftke = .true. ! Turbulent Kinetic Energy closure (T => fill namzdf_tke) ln_zdfgls = .false. ! Generic Length Scale closure (T => fill namzdf_gls) ln_zdfosm = .false. ! OSMOSIS BL closure (T => fill namzdf_osm) ! ! ! convection ln_zdfevd = .true. ! enhanced vertical diffusion nn_evdm = 1 ! apply on tracer (=0) or on tracer and momentum (=1) rn_evd = 100. ! mixing coefficient [m2/s] ln_zdfnpc = .false. ! Non-Penetrative Convective algorithm nn_npc = 1 ! frequency of application of npc nn_npcp = 365 ! npc control print frequency ! ln_zdfddm = .false. ! double diffusive mixing rn_avts = 1.e-4 ! maximum avs (vertical mixing on salinity) rn_hsbfr = 1.6 ! heat/salt buoyancy flux ratio ! ! ! gravity wave-driven vertical mixing ln_zdfiwm = .false. ! internal wave-induced mixing (T => fill namzdf_iwm) ln_zdfswm = .false. ! surface wave-induced mixing (T => ln_wave=ln_sdw=T ) ! ! ! coefficients rn_avm0 = 1.2e-4 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F) rn_avt0 = 1.2e-5 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F) nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) / !----------------------------------------------------------------------- &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion (ln_zdftke =T) !----------------------------------------------------------------------- nn_etau = 0 ! penetration of tke below the mixed layer (ML) due to internal & intertial waves / !!====================================================================== !! *** Miscellaneous namelists *** !!====================================================================== !----------------------------------------------------------------------- &nammpp ! Massively Parallel Processing ("key_mpp_mpi) !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &namctl ! Control prints & Benchmark !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4") !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &namtrd ! diagnostics on dynamics and/or tracer trends ("key_trddyn" and/or "key_trdtra") ! ! or mixed-layer trends or barotropic vorticity ("key_trdmld" or "key_trdvor") !----------------------------------------------------------------------- ln_glo_trd = .false. ! (T) global domain averaged diag for T, T^2, KE, and PE ln_dyn_trd = .false. ! (T) 3D momentum trend output ln_dyn_mxl = .FALSE. ! (T) 2D momentum trends averaged over the mixed layer (not coded yet) ln_vor_trd = .FALSE. ! (T) 2D barotropic vorticity trends (not coded yet) ln_KE_trd = .false. ! (T) 3D Kinetic Energy trends ln_PE_trd = .false. ! (T) 3D Potential Energy trends ln_tra_trd = .false. ! (T) 3D tracer trend output ln_tra_mxl = .false. ! (T) 2D tracer trends averaged over the mixed layer (not coded yet) nn_trd = 365 ! print frequency (ln_glo_trd=T) (unit=time step) / !!gm nn_ctls = 0 ! control surface type in mixed-layer trends (0,1 or n /seconds ; =86400. -> /day) !!gm cn_trdrst_in = "restart_mld" ! suffix of ocean restart name (input) !!gm cn_trdrst_out = "restart_mld" ! suffix of ocean restart name (output) !!gm ln_trdmld_restart = .false. ! restart for ML diagnostics !!gm ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S !!gm !----------------------------------------------------------------------- &namhsb ! Heat and salt budgets !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &namdct ! transports through sections !----------------------------------------------------------------------- nn_dct = 60 ! time step frequency for transports computing nn_dctwri = 60 ! time step frequency for transports writing nn_secdebug = 0 ! 0 : no section to debug / !----------------------------------------------------------------------- &namobs ! observation usage switch !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &nam_asminc ! assimilation increments ('key_asminc') !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &namsbc_wave ! External fields from wave model !----------------------------------------------------------------------- / !----------------------------------------------------------------------- &namwad ! Wetting and drying !----------------------------------------------------------------------- ln_wd_il = .false ! T/F activation of iterative limiter for wetting and drying scheme ln_wd_dl = .true. ! T/F activation of directional llimiter for wetting drying scheme ln_wd_dl_bc = .true. ! T/F Directional limiteer Baroclinic option ln_wd_dl_rmp = .true. ! T/F Turn on directional limiter ramp rn_wdmin0 = 0.30 ! dpoth at which wetting/drying starts rn_wdmin1 = 0.2 ! Minimum wet depth on dried cells rn_wdmin2 = 0.0001 ! Tolerance of min wet depth on dried cells rn_wdld = 2.5 ! Land elevation below which wetting/drying is allowed nn_wdit = 20 ! Max iterations for W/D limiter /