!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> !! NEMO/OCE : Reference namelist_ref !! !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> !! NEMO/OPA : 1 - run manager (namrun) !! namelists 2 - Domain (namcfg, namdom, namtsd, namcrs, namc1d, namc1d_uvd) !! 3 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_cpl, !! namsbc_sas, namtra_qsr, namsbc_rnf, !! namsbc_isf, namsbc_iscpl, namsbc_apr, !! namsbc_ssr, namsbc_wave, namberg) !! 4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) !! 5 - top/bot boundary (namdrg, namdrg_top, namdrg_bot, nambbc, nambbl) !! 6 - Tracer (nameos, namtra_adv, namtra_ldf, namtra_eiv, namtra_dmp) !! 7 - dynamics (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) !! 8 - Vertical physics (namzdf, namzdf_ric, namzdf_tke, namzdf_gls, namzdf_iwm) !! 9 - miscellaneous (nammpp, namctl) !! 10 - diagnostics (namnc4, namtrd, namspr, namflo, namhsb, namsto) !! 11 - Obs & Assim (namobs, nam_asminc) !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> !!====================================================================== !! *** Domain & Run management namelists *** !! !! !! !! namrun parameters of the run !! namdom space and time domain !! namcfg parameters of the configuration (default: user defined GYRE) !! namwad Wetting and drying (default NO) !! namtsd data: temperature & salinity (default NO) !! namcrs coarsened grid (for outputs and/or TOP) (ln_crs =T) !! namc1d 1D configuration options ("key_c1d") !! namc1d_dyndmp 1D newtonian damping applied on currents ("key_c1d") !! namc1d_uvd 1D data (currents) ("key_c1d") !!====================================================================== ! !----------------------------------------------------------------------- &namrun ! parameters of the run !----------------------------------------------------------------------- nn_no = 0 ! Assimilation cycle index cn_exp = "ORCA2" ! experience name nn_it000 = 1 ! first time step nn_itend = 5475 ! last time step (std 5475) nn_date0 = 010101 ! date at nit_0000 (format yyyymmdd) used if ln_rstart=F or (ln_rstart=T and nn_rstctl=0 or 1) nn_time0 = 0 ! initial time of day in hhmm nn_leapy = 0 ! Leap year calendar (1) or not (0) ln_rstart = .false. ! start from rest (F) or from a restart file (T) nn_euler = 1 ! = 0 : start with forward time step if ln_rstart=T nn_rstctl = 0 ! restart control ==> activated only if ln_rstart=T ! ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist ! ! = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart ! ! = 2 nn_date0 read in restart ; nn_it000 : check consistancy between namelist and restart cn_ocerst_in = "restart" ! suffix of ocean restart name (input) cn_ocerst_indir = "." ! directory from which to read input ocean restarts cn_ocerst_out = "restart" ! suffix of ocean restart name (output) cn_ocerst_outdir= "." ! directory in which to write output ocean restarts ln_iscpl = .false. ! cavity evolution forcing or coupling to ice sheet model nn_istate = 0 ! output the initial state (1) or not (0) ln_rst_list = .false. ! output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) nn_stock = 5475 ! frequency of creation of a restart file (modulo referenced to 1) nn_stocklist = 0,0,0,0,0,0,0,0,0,0 ! List of timesteps when a restart file is to be written nn_write = 5475 ! frequency of write in the output file (modulo referenced to nn_it000) ln_mskland = .false. ! mask land points in NetCDF outputs (costly: + ~15%) ln_cfmeta = .false. ! output additional data to netCDF files required for compliance with the CF metadata standard ln_clobber = .true. ! clobber (overwrite) an existing file nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) ln_xios_read = .FALSE. ! use XIOS to read restart file (only for a single file restart) nn_wxios = 0 ! use XIOS to write restart file 0 - no, 1 - single file output, 2 - multiple file output / !----------------------------------------------------------------------- &namdom ! time and space domain !----------------------------------------------------------------------- ln_linssh = .false. ! =T linear free surface ==>> model level are fixed in time rn_isfhmin = 1.00 ! treshold [m] to discriminate grounding ice from floating ice ! rn_rdt = 5760. ! time step for the dynamics and tracer rn_atfp = 0.1 ! asselin time filter parameter ! ln_crs = .false. ! Logical switch for coarsening module (T => fill namcrs) ! ln_meshmask = .false. ! =T create a mesh file / !----------------------------------------------------------------------- &namcfg ! parameters of the configuration (default: use namusr_def in namelist_cfg) !----------------------------------------------------------------------- ln_read_cfg = .false. ! (=T) read the domain configuration file ! ! (=F) user defined configuration (F => create/check namusr_def) cn_domcfg = "domain_cfg" ! domain configuration filename ! ln_closea = .false. ! T => keep closed seas (defined by closea_mask field) in the ! ! domain and apply special treatment of freshwater fluxes. ! ! F => suppress closed seas (defined by closea_mask field) ! ! from the bathymetry at runtime. ! ! If closea_mask field doesn't exist in the domain_cfg file ! ! then this logical does nothing. ln_write_cfg= .false. ! (=T) create the domain configuration file cn_domcfg_out = "domain_cfg_out" ! newly created domain configuration filename ! ln_use_jattr = .false. ! use (T) the file attribute: open_ocean_jstart, if present ! ! in netcdf input files, as the start j-row for reading / !----------------------------------------------------------------------- &namtsd ! Temperature & Salinity Data (init/dmp) (default: NO) !----------------------------------------------------------------------- ! ! =T read T-S fields for: ln_tsd_init = .false. ! ocean initialisation ln_tsd_dmp = .false. ! T-S restoring (see namtra_dmp) cn_dir = './' ! root directory for the T-S data location !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! sn_tem = 'data_1m_potential_temperature_nomask', -1 , 'votemper', .true. , .true. , 'yearly' , '' , '' , '' sn_sal = 'data_1m_salinity_nomask' , -1 , 'vosaline', .true. , .true. , 'yearly' , '' , '' , '' / !----------------------------------------------------------------------- &namwad ! Wetting and Drying (WaD) (default: NO) !----------------------------------------------------------------------- ln_wd_il = .false ! T/F activation of iterative limiter ln_wd_dl = .false. ! T/F activation of directional limiter ln_wd_dl_bc = .false. ! T/F Directional limiteer Baroclinic option ln_wd_dl_rmp= .false. ! T/F Turn on directional limiter ramp rn_wdmin0 = 0.30 ! depth at which WaD 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 WaD is allowed nn_wdit = 20 ! Max iterations for WaD limiter / !----------------------------------------------------------------------- &namcrs ! coarsened grid (for outputs and/or TOP) (ln_crs =T) !----------------------------------------------------------------------- nn_factx = 3 ! Reduction factor of x-direction nn_facty = 3 ! Reduction factor of y-direction nn_binref = 0 ! Bin centering preference: NORTH or EQUAT ! ! 0, coarse grid is binned with preferential treatment of the north fold ! ! 1, coarse grid is binned with centering at the equator ! ! Symmetry with nn_facty being odd-numbered. Asymmetry with even-numbered nn_facty. ln_msh_crs = .false. ! =T create a mesh & mask file nn_crs_kz = 0 ! 0, MEAN of volume boxes ! ! 1, MAX of boxes ! ! 2, MIN of boxes ln_crs_wn = .true. ! wn coarsened (T) or computed using horizontal divergence ( F ) / !----------------------------------------------------------------------- &namc1d ! 1D configuration options ("key_c1d" default: PAPA station) !----------------------------------------------------------------------- rn_lat1d = 50 ! Column latitude (default at PAPA station) rn_lon1d = -145 ! Column longitude (default at PAPA station) ln_c1d_locpt= .true. ! Localization of 1D config in a grid (T) or independant point (F) / !----------------------------------------------------------------------- &namc1d_dyndmp ! U & V newtonian damping ("key_c1d" default: NO) !----------------------------------------------------------------------- ln_dyndmp = .false. ! add a damping term (T) or not (F) / !----------------------------------------------------------------------- &namc1d_uvd ! data: U & V currents ("key_c1d" default: NO) !----------------------------------------------------------------------- ! ! =T read U-V fields for: ln_uvd_init = .false. ! ocean initialisation ln_uvd_dyndmp = .false. ! U-V restoring cn_dir = './' ! root directory for the U-V data location !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! sn_ucur = 'ucurrent' , -1 ,'u_current', .false. , .true. , 'monthly' , '' , 'Ume' , '' sn_vcur = 'vcurrent' , -1 ,'v_current', .false. , .true. , 'monthly' , '' , 'Vme' , '' / !!====================================================================== !! *** Surface Boundary Condition namelists *** !! !! !! !! namsbc surface boundary condition manager (default: NO selection) !! namsbc_flx flux formulation (ln_flx =T) !! namsbc_blk Bulk formulae formulation (ln_blk =T) !! namsbc_cpl CouPLed formulation ("key_oasis3" ) !! namsbc_sas Stand-Alone Surface module (SAS_SRC only) !! namsbc_iif Ice-IF: use observed ice cover (nn_ice = 1 ) !! namtra_qsr penetrative solar radiation (ln_traqsr =T) !! namsbc_ssr sea surface restoring term (for T and/or S) (ln_ssr =T) !! namsbc_rnf river runoffs (ln_rnf =T) !! namsbc_apr Atmospheric Pressure (ln_apr_dyn =T) !! namsbc_isf ice shelf melting/freezing (ln_isfcav =T : read (ln_read_cfg=T) or set or usr_def_zgr ) !! namsbc_iscpl coupling option between land ice model and ocean (ln_isfcav =T) !! namsbc_wave external fields from wave model (ln_wave =T) !! namberg iceberg floats (ln_icebergs=T) !!====================================================================== ! !----------------------------------------------------------------------- &namsbc ! Surface Boundary Condition manager (default: NO selection) !----------------------------------------------------------------------- nn_fsbc = 5 ! frequency of SBC module call ! ! (control sea-ice & iceberg model call) ! Type of air-sea fluxes ln_usr = .false. ! user defined formulation (T => check usrdef_sbc) ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) ln_blk = .false. ! Bulk formulation (T => fill namsbc_blk ) ! ! Type of coupling (Ocean/Ice/Atmosphere) : ln_cpl = .false. ! atmosphere coupled formulation ( requires key_oasis3 ) ln_mixcpl = .false. ! forced-coupled mixed formulation ( requires key_oasis3 ) nn_components = 0 ! configuration of the opa-sas OASIS coupling ! ! =0 no opa-sas OASIS coupling: default single executable config. ! ! =1 opa-sas OASIS coupling: multi executable config., OPA component ! ! =2 opa-sas OASIS coupling: multi executable config., SAS component ! Sea-ice : nn_ice = 0 ! =0 no ice boundary condition ! ! =1 use observed ice-cover ( => fill namsbc_iif ) ! ! =2 or 3 automatically for LIM3 or CICE ("key_lim3" or "key_cice") ! ! except in AGRIF zoom where it has to be specified ln_ice_embd = .false. ! =T embedded sea-ice (pressure + mass and salt exchanges) ! ! =F levitating ice (no pressure, mass and salt exchanges) ! Misc. options of sbc : ln_traqsr = .false. ! Light penetration in the ocean (T => fill namtra_qsr) ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave ln_ssr = .false. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr) nn_fwb = 0 ! FreshWater Budget: =0 unchecked ! ! =1 global mean of e-p-r set to zero at each time step ! ! =2 annual global mean of e-p-r set to zero ln_rnf = .false. ! runoffs (T => fill namsbc_rnf) ln_apr_dyn = .false. ! Patm gradient added in ocean & ice Eqs. (T => fill namsbc_apr ) ln_isf = .false. ! ice shelf (T => fill namsbc_isf & namsbc_iscpl) ln_wave = .false. ! Activate coupling with wave (T => fill namsbc_wave) ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave) ln_sdw = .false. ! Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave) nn_sdrift = 0 ! Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift ! ! = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)] ! ! = 1 Phillips: v_z=v_o*[exp(2*k*z)-beta*sqrt(-2*k*pi*z)*erfc(sqrt(-2*k*z))] ! ! = 2 Phillips as (1) but using the wave frequency from a wave model ln_tauwoc = .false. ! Activate ocean stress modified by external wave induced stress (T => ln_wave=.true. & fill namsbc_wave) ln_tauw = .false. ! Activate ocean stress components from wave model ln_stcor = .false. ! Activate Stokes Coriolis term (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave) nn_lsm = 0 ! =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) , ! =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field) / !----------------------------------------------------------------------- &namsbc_flx ! surface boundary condition : flux formulation (ln_flx =T) !----------------------------------------------------------------------- cn_dir = './' ! root directory for the fluxes data location !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 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' , '' , '' , '' / !----------------------------------------------------------------------- &namsbc_blk ! namsbc_blk generic Bulk formula (ln_blk =T) !----------------------------------------------------------------------- ! ! bulk algorithm : ln_NCAR = .false. ! "NCAR" algorithm (Large and Yeager 2008) ln_COARE_3p0= .false. ! "COARE 3.0" algorithm (Fairall et al. 2003) ln_COARE_3p5= .false. ! "COARE 3.5" algorithm (Edson et al. 2013) ln_ECMWF = .false. ! "ECMWF" algorithm (IFS cycle 31) ! rn_zqt = 10. ! Air temperature & humidity reference height (m) rn_zu = 10. ! Wind vector reference height (m) ln_Cd_L12 = .false. ! air-ice drags = F(ice concentration) (Lupkes et al. 2012) ln_Cd_L15 = .false. ! air-ice drags = F(ice concentration) (Lupkes et al. 2015) ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data rn_pfac = 1. ! multiplicative factor for precipitation (total & snow) rn_efac = 1. ! multiplicative factor for evaporation (0. or 1.) rn_vfac = 0. ! multiplicative factor for ocean & ice velocity used to ! ! calculate the wind stress (0.=absolute or 1.=relative winds) cn_dir = './' ! root directory for the bulk data location !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! sn_wndi = 'u_10.15JUNE2009_fill' , 6 , 'U_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bicubic_noc.nc' , 'Uwnd' , '' sn_wndj = 'v_10.15JUNE2009_fill' , 6 , 'V_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bicubic_noc.nc' , 'Vwnd' , '' sn_qsr = 'ncar_rad.15JUNE2009_fill' , 24 , 'SWDN_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' sn_qlw = 'ncar_rad.15JUNE2009_fill' , 24 , 'LWDN_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' sn_tair = 't_10.15JUNE2009_fill' , 6 , 'T_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' sn_humi = 'q_10.15JUNE2009_fill' , 6 , 'Q_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' sn_prec = 'ncar_precip.15JUNE2009_fill', -1 , 'PRC_MOD1', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' sn_snow = 'ncar_precip.15JUNE2009_fill', -1 , 'SNOW' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' sn_slp = 'slp.15JUNE2009_fill' , 6 , 'SLP' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' sn_tdif = 'taudif_core' , 24 , 'taudif' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' / !----------------------------------------------------------------------- &namsbc_cpl ! coupled ocean/atmosphere model ("key_oasis3") !----------------------------------------------------------------------- nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentially sending/receiving data ln_usecplmask = .false. ! use a coupling mask file to merge data received from several models ! ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) nn_cats_cpl = 5 ! Number of sea ice categories over which coupling is to be carried out (if not 1) !_____________!__________________________!____________!_____________!______________________!________! ! ! description ! multiple ! vector ! vector ! vector ! ! ! ! categories ! reference ! orientation ! grids ! !*** send *** sn_snd_temp = 'weighted oce and ice' , 'no' , '' , '' , '' sn_snd_alb = 'weighted ice' , 'no' , '' , '' , '' sn_snd_thick = 'none' , 'no' , '' , '' , '' sn_snd_crt = 'none' , 'no' , 'spherical' , 'eastward-northward' , 'T' sn_snd_co2 = 'coupled' , 'no' , '' , '' , '' sn_snd_crtw = 'none' , 'no' , '' , '' , 'U,V' sn_snd_ifrac = 'none' , 'no' , '' , '' , '' sn_snd_wlev = 'coupled' , 'no' , '' , '' , '' sn_snd_cond = 'weighted ice' , 'no' , '' , '' , '' sn_snd_thick1 = 'ice and snow' , 'no' , '' , '' , '' sn_snd_mpnd = 'weighted ice' , 'no' , '' , '' , '' sn_snd_sstfrz = 'coupled' , 'no' , '' , '' , '' sn_snd_ttilyr = 'weighted ice' , 'no' , '' , '' , '' !*** receive *** sn_rcv_w10m = 'none' , 'no' , '' , '' , '' sn_rcv_taumod = 'coupled' , 'no' , '' , '' , '' sn_rcv_tau = 'oce only' , 'no' , 'cartesian' , 'eastward-northward' , 'U,V' sn_rcv_dqnsdt = 'coupled' , 'no' , '' , '' , '' sn_rcv_qsr = 'oce and ice' , 'no' , '' , '' , '' sn_rcv_qns = 'oce and ice' , 'no' , '' , '' , '' sn_rcv_emp = 'conservative' , 'no' , '' , '' , '' sn_rcv_rnf = 'coupled' , 'no' , '' , '' , '' sn_rcv_cal = 'coupled' , 'no' , '' , '' , '' sn_rcv_co2 = 'coupled' , 'no' , '' , '' , '' sn_rcv_hsig = 'none' , 'no' , '' , '' , '' sn_rcv_iceflx = 'none' , 'no' , '' , '' , '' sn_rcv_mslp = 'none' , 'no' , '' , '' , '' sn_rcv_phioc = 'none' , 'no' , '' , '' , '' sn_rcv_sdrfx = 'none' , 'no' , '' , '' , '' sn_rcv_sdrfy = 'none' , 'no' , '' , '' , '' sn_rcv_wper = 'none' , 'no' , '' , '' , '' sn_rcv_wnum = 'none' , 'no' , '' , '' , '' sn_rcv_wstrf = 'none' , 'no' , '' , '' , '' sn_rcv_wdrag = 'none' , 'no' , '' , '' , '' sn_rcv_ts_ice = 'none' , 'no' , '' , '' , '' sn_rcv_isf = 'none' , 'no' , '' , '' , '' sn_rcv_icb = 'none' , 'no' , '' , '' , '' sn_rcv_tauwoc = 'none' , 'no' , '' , '' , '' sn_rcv_tauw = 'none' , 'no' , '' , '' , '' sn_rcv_wdrag = 'none' , 'no' , '' , '' , '' / !----------------------------------------------------------------------- &namsbc_sas ! Stand-Alone Surface module: ocean data (SAS_SRC only) !----------------------------------------------------------------------- l_sasread = .true. ! =T Read in file ; =F set all to 0. (see sbcssm) ln_3d_uve = .false. ! specify whether we are supplying a 3D u,v and e3 field ln_read_frq = .false. ! specify whether we must read frq or not cn_dir = './' ! root directory for the ocean data location !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! sn_usp = 'sas_grid_U' , 120 , 'uos' , .true. , .true. , 'yearly' , '' , '' , '' sn_vsp = 'sas_grid_V' , 120 , 'vos' , .true. , .true. , 'yearly' , '' , '' , '' sn_tem = 'sas_grid_T' , 120 , 'sosstsst', .true. , .true. , 'yearly' , '' , '' , '' sn_sal = 'sas_grid_T' , 120 , 'sosaline', .true. , .true. , 'yearly' , '' , '' , '' sn_ssh = 'sas_grid_T' , 120 , 'sossheig', .true. , .true. , 'yearly' , '' , '' , '' sn_e3t = 'sas_grid_T' , 120 , 'e3t_m' , .true. , .true. , 'yearly' , '' , '' , '' sn_frq = 'sas_grid_T' , 120 , 'frq_m' , .true. , .true. , 'yearly' , '' , '' , '' / !----------------------------------------------------------------------- &namsbc_iif ! Ice-IF : use observed ice cover (nn_ice = 1) !----------------------------------------------------------------------- cn_dir = './' ! root directory for the ice cover data location !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! sn_ice ='ice_cover_clim.nc' , -12. ,'ice_cover', .true. , .true. , 'yearly' , '' , '' , '' / !----------------------------------------------------------------------- &namtra_qsr ! penetrative solar radiation (ln_traqsr =T) !----------------------------------------------------------------------- ! ! type of penetration (default: NO selection) ln_qsr_rgb = .false. ! RGB light penetration (Red-Green-Blue) ln_qsr_2bd = .false. ! 2BD light penetration (two bands) ln_qsr_bio = .false. ! bio-model light penetration ! ! RGB & 2BD choices: rn_abs = 0.58 ! RGB & 2BD: fraction absorbed in the very near surface rn_si0 = 0.35 ! RGB & 2BD: shortess depth of extinction nn_chldta = 0 ! RGB : Chl data (=1) or cst value (=0) rn_si1 = 23.0 ! 2BD : longest depth of extinction cn_dir = './' ! root directory for the chlorophyl data location !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! sn_chl ='chlorophyll' , -1 , 'CHLA' , .true. , .true. , 'yearly' , '' , '' , '' / !----------------------------------------------------------------------- &namsbc_ssr ! surface boundary condition : sea surface restoring (ln_ssr =T) !----------------------------------------------------------------------- nn_sstr = 0 ! add a retroaction term to the surface heat flux (=1) or not (=0) rn_dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K] nn_sssr = 0 ! add a damping term to the surface freshwater flux (=2) ! ! or to SSS only (=1) or no damping term (=0) rn_deds = -166.67 ! magnitude of the damping on salinity [mm/day] ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2) rn_sssr_bnd = 4.e0 ! ABS(Max/Min) value of the damping erp term [mm/day] cn_dir = './' ! root directory for the SST/SSS data location !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! sn_sst = 'sst_data' , 24 , 'sst' , .false. , .false., 'yearly' , '' , '' , '' sn_sss = 'sss_data' , -1 , 'sss' , .true. , .true. , 'yearly' , '' , '' , '' / !----------------------------------------------------------------------- &namsbc_rnf ! runoffs (ln_rnf =T) !----------------------------------------------------------------------- ln_rnf_mouth= .false. ! specific treatment at rivers mouths rn_hrnf = 15.e0 ! depth over which enhanced vertical mixing is used (ln_rnf_mouth=T) rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T) rn_rfact = 1.e0 ! multiplicative factor for runoff ln_rnf_depth= .false. ! read in depth information for runoff ln_rnf_tem = .false. ! read in temperature information for runoff ln_rnf_sal = .false. ! read in salinity information for runoff ln_rnf_depth_ini = .false. ! compute depth at initialisation from runoff file rn_rnf_max = 5.735e-4 ! max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true ) rn_dep_max = 150. ! depth over which runoffs is spread ( ln_rnf_depth_ini = .true ) nn_rnf_depth_file = 0 ! create (=1) a runoff depth file or not (=0) cn_dir = './' ! root directory for the runoff data location !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! sn_rnf = 'runoff_core_monthly' , -1 , 'sorunoff', .true. , .true. , 'yearly' , '' , '' , '' sn_cnf = 'runoff_core_monthly' , 0 , 'socoefr0', .false. , .true. , 'yearly' , '' , '' , '' sn_s_rnf = 'runoffs' , 24 , 'rosaline', .true. , .true. , 'yearly' , '' , '' , '' sn_t_rnf = 'runoffs' , 24 , 'rotemper', .true. , .true. , 'yearly' , '' , '' , '' sn_dep_rnf = 'runoffs' , 0 , 'rodepth' , .false. , .true. , 'yearly' , '' , '' , '' / !----------------------------------------------------------------------- &namsbc_apr ! Atmospheric pressure used as ocean forcing (ln_apr_dyn =T) !----------------------------------------------------------------------- rn_pref = 101000. ! reference atmospheric pressure [N/m2]/ ln_ref_apr = .false. ! ref. pressure: global mean Patm (T) or a constant (F) ln_apr_obc = .false. ! inverse barometer added to OBC ssh data cn_dir = './' ! root directory for the Patm data location !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! sn_apr = 'patm' , -1 ,'somslpre' , .true. , .true. , 'yearly' , '' , '' , '' / !----------------------------------------------------------------------- &namsbc_isf ! Top boundary layer (ISF) (ln_isfcav =T : read (ln_read_cfg=T) !----------------------------------------------------------------------- or set or usr_def_zgr ) ! ! type of top boundary layer nn_isf = 1 ! ice shelf melting/freezing ! 1 = presence of ISF ; 2 = bg03 parametrisation ! 3 = rnf file for ISF ; 4 = ISF specified freshwater flux ! option 1 and 4 need ln_isfcav = .true. (domzgr) ! ! nn_isf = 1 or 2 cases: rn_gammat0 = 1.e-4 ! gammat coefficient used in blk formula rn_gammas0 = 1.e-4 ! gammas coefficient used in blk formula ! ! nn_isf = 1 or 4 cases: rn_hisf_tbl = 30. ! thickness of the top boundary layer (Losh et al. 2008) ! ! 0 => thickness of the tbl = thickness of the first wet cell ! ! nn_isf = 1 case nn_isfblk = 1 ! 1 ISOMIP like: 2 equations formulation (Hunter et al., 2006) ! ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) nn_gammablk = 1 ! 0 = cst Gammat (= gammat/s) ! ! 1 = velocity dependend Gamma (u* * gammat/s) (Jenkins et al. 2010) ! ! 2 = velocity and stability dependent Gamma (Holland et al. 1999) !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! ! ! 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 ! !* nn_isf = 4 case sn_fwfisf = 'rnfisf' , -12 ,'sowflisf' , .false. , .true. , 'yearly' , '' , '' , '' !* nn_isf = 3 case sn_rnfisf = 'rnfisf' , -12 ,'sofwfisf' , .false. , .true. , 'yearly' , '' , '' , '' !* nn_isf = 2 and 3 cases sn_depmax_isf='rnfisf' , -12 ,'sozisfmax', .false. , .true. , 'yearly' , '' , '' , '' sn_depmin_isf='rnfisf' , -12 ,'sozisfmin', .false. , .true. , 'yearly' , '' , '' , '' !* nn_isf = 2 case sn_Leff_isf = 'rnfisf' , -12 ,'Leff' , .false. , .true. , 'yearly' , '' , '' , '' / !----------------------------------------------------------------------- &namsbc_iscpl ! land ice / ocean coupling option (ln_isfcav =T : read (ln_read_cfg=T) !----------------------------------------------------------------------- or set or usr_def_zgr ) nn_drown = 10 ! number of iteration of the extrapolation loop (fill the new wet cells) ln_hsb = .false. ! activate conservation module (conservation exact after a time of rn_fiscpl) nn_fiscpl = 43800 ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency) / !----------------------------------------------------------------------- &namsbc_wave ! External fields from wave model (ln_wave=T) !----------------------------------------------------------------------- cn_dir = './' ! root directory for the waves data location !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! sn_cdg = 'sdw_wave' , 1 , 'drag_coeff' , .true. , .false., 'daily' , '' , '' , '' sn_usd = 'sdw_wave' , 1 , 'u_sd2d' , .true. , .false., 'daily' , '' , '' , '' sn_vsd = 'sdw_wave' , 1 , 'v_sd2d' , .true. , .false., 'daily' , '' , '' , '' sn_hsw = 'sdw_wave' , 1 , 'hs' , .true. , .false., 'daily' , '' , '' , '' sn_wmp = 'sdw_wave' , 1 , 'wmp' , .true. , .false., 'daily' , '' , '' , '' sn_wfr = 'sdw_wave' , 1 , 'wfr' , .true. , .false., 'daily' , '' , '' , '' sn_wnum = 'sdw_wave' , 1 , 'wave_num' , .true. , .false., 'daily' , '' , '' , '' sn_tauwoc = 'sdw_wave' , 1 , 'wave_stress', .true. , .false., 'daily' , '' , '' , '' sn_tauwx = 'sdw_wave' , 1 , 'wave_stress', .true. , .false., 'daily' , '' , '' , '' sn_tauwy = 'sdw_wave' , 1 , 'wave_stress', .true. , .false., 'daily' , '' , '' , '' / !----------------------------------------------------------------------- &namberg ! iceberg parameters (default: NO) !----------------------------------------------------------------------- ln_icebergs = .false. ! activate iceberg floats (force =F with "key_agrif") ! ! ! diagnostics: ln_bergdia = .true. ! Calculate budgets nn_verbose_level = 1 ! Turn on more verbose output if level > 0 nn_verbose_write = 15 ! Timesteps between verbose messages nn_sample_rate = 1 ! Timesteps between sampling for trajectory storage ! ! ! iceberg setting: ! ! Initial mass required for an iceberg of each class rn_initial_mass = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11 ! ! Proportion of calving mass to apportion to each class rn_distribution = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02 ! ! Ratio between effective and real iceberg mass (non-dim) ! ! i.e. number of icebergs represented at a point rn_mass_scaling = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1 ! thickness of newly calved bergs (m) rn_initial_thickness = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250. ! rn_rho_bergs = 850. ! Density of icebergs rn_LoW_ratio = 1.5 ! Initial ratio L/W for newly calved icebergs ln_operator_splitting = .true. ! Use first order operator splitting for thermodynamics rn_bits_erosion_fraction = 0. ! Fraction of erosion melt flux to divert to bergy bits rn_sicn_shift = 0. ! Shift of sea-ice concn in erosion flux (0 fill namdrg_bot ln_lin = .false. ! linear drag: Cd = Cd0 Uc0 & namdrg_top) ln_non_lin = .false. ! non-linear drag: Cd = Cd0 |U| ln_loglayer= .false. ! logarithmic drag: Cd = vkarmn/log(z/z0) |U| ! ln_drgimp = .true. ! implicit top/bottom friction flag / !----------------------------------------------------------------------- &namdrg_top ! TOP friction (ln_NONE=F & ln_isfcav=T) !----------------------------------------------------------------------- rn_Cd0 = 1.e-3 ! drag coefficient [-] rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) rn_Cdmax = 0.1 ! drag value maximum [-] (logarithmic drag) rn_ke0 = 2.5e-3 ! background kinetic energy [m2/s2] (non-linear cases) rn_z0 = 3.0e-3 ! roughness [m] (ln_loglayer=T) ln_boost = .false. ! =T regional boost of Cd0 ; =F constant rn_boost= 50. ! local boost factor [-] / !----------------------------------------------------------------------- &namdrg_bot ! BOTTOM friction (ln_NONE=F) !----------------------------------------------------------------------- rn_Cd0 = 1.e-3 ! drag coefficient [-] rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) rn_Cdmax = 0.1 ! 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 (default: NO) !----------------------------------------------------------------------- ln_trabbc = .false. ! Apply a geothermal heating at the ocean bottom nn_geoflx = 2 ! geothermal heat flux: = 1 constant flux ! ! = 2 read variable flux [mW/m2] rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [mW/m2] cn_dir = './' ! root directory for the geothermal data location !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! sn_qgh ='geothermal_heating.nc' , -12. , 'heatflow', .false. , .true. , 'yearly' , '' , '' , '' / !----------------------------------------------------------------------- &nambbl ! bottom boundary layer scheme (default: NO) !----------------------------------------------------------------------- ln_trabbl = .false. ! Bottom Boundary Layer parameterisation flag nn_bbl_ldf = 1 ! diffusive bbl (=1) or not (=0) nn_bbl_adv = 0 ! advective bbl (=1/2) or not (=0) rn_ahtbbl = 1000. ! lateral mixing coefficient in the bbl [m2/s] rn_gambbl = 10. ! advective bbl coefficient [s] / !!====================================================================== !! Tracer (T & S) namelists !! !! !! !! nameos equation of state (default: NO selection) !! namtra_adv advection scheme (default: NO selection) !! namtra_ldf lateral diffusion scheme (default: NO selection) !! namtra_adv_mle mixed layer eddy param. (Fox-Kemper param.) (default: NO) !! namtra_eiv eddy induced velocity param. (default: NO) !! namtra_dmp T & S newtonian damping (default: NO) !!====================================================================== ! !----------------------------------------------------------------------- &nameos ! ocean Equation Of Seawater (default: NO selection) !----------------------------------------------------------------------- ln_teos10 = .false. ! = Use TEOS-10 ln_eos80 = .false. ! = Use EOS80 ln_seos = .false. ! = Use S-EOS (simplified Eq.) ! ! ! 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 rn_b0 = 7.6554e-1 ! saline expension coefficient 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) / !----------------------------------------------------------------------- &namtra_adv ! advection scheme for tracer (default: NO selection) !----------------------------------------------------------------------- ln_traadv_NONE= .false. ! No tracer advection ln_traadv_cen = .false. ! 2nd order centered scheme nn_cen_h = 4 ! =2/4, horizontal 2nd order CEN / 4th order CEN nn_cen_v = 4 ! =2/4, vertical 2nd order CEN / 4th order COMPACT ln_traadv_fct = .false. ! FCT scheme nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order ln_traadv_mus = .false. ! MUSCL scheme ln_mus_ups = .false. ! use upstream scheme near river mouths ln_traadv_ubs = .false. ! UBS scheme nn_ubs_v = 2 ! =2 , vertical 2nd order FCT / COMPACT 4th order ln_traadv_qck = .false. ! QUICKEST scheme / !----------------------------------------------------------------------- &namtra_ldf ! lateral diffusion scheme for tracers (default: NO selection) !----------------------------------------------------------------------- ! ! Operator type: ln_traldf_NONE = .false. ! No explicit diffusion ln_traldf_lap = .false. ! 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 = .false. ! iso-neutral (standard operator) ln_traldf_triad = .false. ! iso-neutral (triad operator) ! ! ! 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 coefficient: ! ! =-20 (=-30) read in eddy_diffusivity_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 and grid-spacing) ! ! time invariant coefficients: aht0 = 1/2 Ud*Ld (lap case) ! ! or = 1/12 Ud*Ld^3 (blp case) rn_Ud = 0.01 ! lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) rn_Ld = 200.e+3 ! lateral diffusive length [m] (nn_aht_ijk_t= 0, 10) / !----------------------------------------------------------------------- &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper) (default: NO) !----------------------------------------------------------------------- ln_mle = .false. ! (T) use the Mixed Layer Eddy (MLE) parameterisation rn_ce = 0.06 ! magnitude of the MLE (typical value: 0.06 to 0.08) nn_mle = 1 ! MLE type: =0 standard Fox-Kemper ; =1 new formulation rn_lf = 5.e+3 ! typical scale of mixed layer front (meters) (case rn_mle=0) rn_time = 172800. ! time scale for mixing momentum across the mixed layer (seconds) (case rn_mle=0) rn_lat = 20. ! reference latitude (degrees) of MLE coef. (case rn_mle=1) nn_mld_uv = 0 ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max) nn_conv = 0 ! =1 no MLE in case of convection ; =0 always MLE rn_rho_c_mle= 0.01 ! delta rho criterion used to calculate MLD for FK / !----------------------------------------------------------------------- &namtra_eiv ! eddy induced velocity param. (default: NO) !----------------------------------------------------------------------- ln_ldfeiv = .false. ! use eddy induced velocity parameterization ! ! ! Coefficients: nn_aei_ijk_t = 0 ! space/time variation of eddy coefficient: ! ! =-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 ! ! time invariant coefficients: aei0 = 1/2 Ue*Le rn_Ue = 0.02 ! lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) rn_Le = 200.e+3 ! lateral diffusive length [m] (nn_aht_ijk_t= 0, 10) ! ln_ldfeiv_dia =.false. ! diagnose eiv stream function and velocities / !----------------------------------------------------------------------- &namtra_dmp ! tracer: T & S newtonian damping (default: NO) !----------------------------------------------------------------------- ln_tradmp = .false. ! add a damping term (using resto.nc coef.) nn_zdmp = 0 ! vertical shape =0 damping throughout the water column ! ! =1 no damping in the mixing layer (kz criteria) ! ! =2 no damping in the mixed layer (rho crieria) cn_resto = 'resto.nc' ! Name of file containing restoration coeff. field (use dmp_tools to create this) / !!====================================================================== !! *** Dynamics namelists *** !! !! !! !! nam_vvl vertical coordinate options (default: z-star) !! namdyn_adv formulation of the momentum advection (default: NO selection) !! namdyn_vor advection scheme (default: NO selection) !! namdyn_hpg hydrostatic pressure gradient (default: NO selection) !! namdyn_spg surface pressure gradient (default: NO selection) !! namdyn_ldf lateral diffusion scheme (default: NO selection) !! namdta_dyn offline TOP: dynamics read in files (OFF_SRC only) !!====================================================================== ! !----------------------------------------------------------------------- &nam_vvl ! vertical coordinate options (default: z-star) !----------------------------------------------------------------------- ln_vvl_zstar = .true. ! z-star vertical coordinate ln_vvl_ztilde = .false. ! z-tilde vertical coordinate: only high frequency variations ln_vvl_layer = .false. ! full layer vertical coordinate ln_vvl_ztilde_as_zstar = .false. ! ztilde vertical coordinate emulating zstar ln_vvl_zstar_at_eqtor = .false. ! ztilde near the equator rn_ahe3 = 0.0 ! thickness diffusion coefficient rn_rst_e3t = 30.0 ! ztilde to zstar restoration timescale [days] rn_lf_cutoff = 5.0 ! cutoff frequency for low-pass filter [days] rn_zdef_max = 0.9 ! maximum fractional e3t deformation ln_vvl_dbg = .true. ! debug prints (T/F) / !----------------------------------------------------------------------- &namdyn_adv ! formulation of the momentum advection (default: NO selection) !----------------------------------------------------------------------- ln_dynadv_NONE= .false. ! linear dynamics (no momentum advection) ln_dynadv_vec = .false. ! 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 ! Vorticity / Coriolis scheme (default: NO selection) !----------------------------------------------------------------------- ln_dynvor_ene = .false. ! energy conserving scheme ln_dynvor_ens = .false. ! enstrophy conserving scheme ln_dynvor_mix = .false. ! mixed scheme ln_dynvor_een = .false. ! energy & enstrophy scheme nn_een_e3f = 1 ! =0 e3f = mean masked e3t divided by 4 ! ! =1 e3f = mean masked e3t divided by the sum of mask ln_dynvor_msk = .false. ! vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) ! PLEASE DO NOT ACTIVATE / !----------------------------------------------------------------------- &namdyn_hpg ! Hydrostatic pressure gradient option (default: NO selection) !----------------------------------------------------------------------- ln_hpg_zco = .false. ! z-coordinate - full steps ln_hpg_zps = .false. ! z-coordinate - partial steps (interpolation) ln_hpg_sco = .false. ! s-coordinate (standard jacobian formulation) ln_hpg_isf = .false. ! s-coordinate (sco ) adapted to isf ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial) ln_hpg_prj = .false. ! s-coordinate (Pressure Jacobian scheme) / !----------------------------------------------------------------------- &namdyn_spg ! surface pressure gradient (default: NO selection) !----------------------------------------------------------------------- ln_dynspg_exp = .false. ! explicit free surface ln_dynspg_ts = .false. ! split-explicit free surface ln_bt_fw = .true. ! Forward integration of barotropic Eqs. ln_bt_av = .true. ! Time filtering of barotropic variables nn_bt_flt = 1 ! Time filter choice = 0 None ! ! = 1 Boxcar over nn_baro sub-steps ! ! = 2 Boxcar over 2*nn_baro " " ln_bt_auto = .true. ! Number of sub-step defined from: rn_bt_cmax = 0.8 ! =T : the Maximum Courant Number allowed nn_baro = 30 ! =F : the number of sub-step in rn_rdt seconds rn_bt_alpha = 0. ! Temporal diffusion parameter (if ln_bt_av=F) / !----------------------------------------------------------------------- &namdyn_ldf ! lateral diffusion on momentum (default: NO selection) !----------------------------------------------------------------------- ! ! Type of the operator : ln_dynldf_NONE= .false. ! No operator (i.e. no explicit diffusion) ln_dynldf_lap = .false. ! laplacian operator ln_dynldf_blp = .false. ! bilaplacian operator ! ! Direction of action : ln_dynldf_lev = .false. ! iso-level ln_dynldf_hor = .false. ! horizontal (geopotential) ln_dynldf_iso = .false. ! iso-neutral (lap only) ! ! Coefficient nn_ahm_ijk_t = 0 ! space/time variation of eddy coefficient : ! ! =-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) ! ! = 32 F(i,j,k)=F(local gridscale and deformation rate) ! ! time invariant coefficients : ahm = 1/2 Uv*Lv (lap case) ! ! or = 1/12 Uv*Lv^3 (blp case) rn_Uv = 0.1 ! lateral viscous velocity [m/s] (nn_ahm_ijk_t= 0, 10, 20, 30) rn_Lv = 10.e+3 ! lateral viscous length [m] (nn_ahm_ijk_t= 0, 10) ! ! Smagorinsky settings (nn_ahm_ijk_t= 32) : rn_csmc = 3.5 ! Smagorinsky constant of proportionality rn_minfac = 1.0 ! multiplier of theorectical lower limit rn_maxfac = 1.0 ! multiplier of theorectical upper limit ! ! iso-neutral laplacian operator (ln_dynldf_iso=T) : rn_ahm_b = 0.0 ! background eddy viscosity [m2/s] / !----------------------------------------------------------------------- &namdta_dyn ! offline ocean input files (OFF_SRC only) !----------------------------------------------------------------------- ln_dynrnf = .false. ! runoffs option enabled (T) or not (F) ln_dynrnf_depth = .false. ! runoffs is spread in vertical (T) or not (F) ! fwbcorr = 3.786e-06 ! annual global mean of empmr for ssh correction cn_dir = './' ! root directory for the ocean data location !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! sn_tem = 'dyna_grid_T' , 120 , 'votemper' , .true. , .true. , 'yearly' , '' , '' , '' sn_sal = 'dyna_grid_T' , 120 , 'vosaline' , .true. , .true. , 'yearly' , '' , '' , '' sn_mld = 'dyna_grid_T' , 120 , 'somixhgt' , .true. , .true. , 'yearly' , '' , '' , '' sn_emp = 'dyna_grid_T' , 120 , 'sowaflup' , .true. , .true. , 'yearly' , '' , '' , '' sn_fmf = 'dyna_grid_T' , 120 , 'iowaflup' , .true. , .true. , 'yearly' , '' , '' , '' sn_ice = 'dyna_grid_T' , 120 , 'soicecov' , .true. , .true. , 'yearly' , '' , '' , '' sn_qsr = 'dyna_grid_T' , 120 , 'soshfldo' , .true. , .true. , 'yearly' , '' , '' , '' sn_wnd = 'dyna_grid_T' , 120 , 'sowindsp' , .true. , .true. , 'yearly' , '' , '' , '' sn_uwd = 'dyna_grid_U' , 120 , 'uocetr_eff', .true. , .true. , 'yearly' , '' , '' , '' sn_vwd = 'dyna_grid_V' , 120 , 'vocetr_eff', .true. , .true. , 'yearly' , '' , '' , '' sn_wwd = 'dyna_grid_W' , 120 , 'wocetr_eff', .true. , .true. , 'yearly' , '' , '' , '' sn_avt = 'dyna_grid_W' , 120 , 'voddmavs' , .true. , .true. , 'yearly' , '' , '' , '' sn_ubl = 'dyna_grid_U' , 120 , 'sobblcox' , .true. , .true. , 'yearly' , '' , '' , '' sn_vbl = 'dyna_grid_V' , 120 , 'sobblcoy' , .true. , .true. , 'yearly' , '' , '' , '' / !!====================================================================== !! vertical physics namelists !! !! !! !! namzdf vertical physics manager (default: NO selection) !! 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_osm OSM vertical diffusion (ln_zdfosm=T) !! namzdf_iwm tidal mixing parameterization (ln_zdfiwm=T) !!====================================================================== ! !----------------------------------------------------------------------- &namzdf ! vertical physics manager (default: NO selection) !----------------------------------------------------------------------- ! ! type of vertical closure (required) ln_zdfcst = .false. ! constant mixing ln_zdfric = .false. ! local Richardson dependent formulation (T => fill namzdf_ric) ln_zdftke = .false. ! 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 = .false. ! enhanced vertical diffusion nn_evdm = 0 ! 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_ric ! richardson number dependent vertical diffusion (ln_zdfric =T) !----------------------------------------------------------------------- rn_avmri = 100.e-4 ! maximum value of the vertical viscosity rn_alp = 5. ! coefficient of the parameterization nn_ric = 2 ! coefficient of the parameterization ln_mldw = .false. ! enhanced mixing in the Ekman layer rn_ekmfc = 0.7 ! Factor in the Ekman depth Equation rn_mldmin = 1.0 ! minimum allowable mixed-layer depth estimate (m) rn_mldmax = 1000.0 ! maximum allowable mixed-layer depth estimate (m) rn_wtmix = 10.0 ! vertical eddy viscosity coeff [m2/s] in the mixed-layer rn_wvmix = 10.0 ! vertical eddy diffusion coeff [m2/s] in the mixed-layer / !----------------------------------------------------------------------- &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion (ln_zdftke =T) !----------------------------------------------------------------------- rn_ediff = 0.1 ! coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) ) rn_ediss = 0.7 ! coef. of the Kolmogoroff dissipation rn_ebb = 67.83 ! coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T) rn_emin = 1.e-6 ! minimum value of tke [m2/s2] rn_emin0 = 1.e-4 ! surface minimum value of tke [m2/s2] rn_bshear = 1.e-20 ! background shear (>0) currently a numerical threshold (do not change it) nn_pdl = 1 ! Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm) nn_mxl = 2 ! mixing length: = 0 bounded by the distance to surface and bottom ! ! = 1 bounded by the local vertical scale factor ! ! = 2 first vertical derivative of mixing length bounded by 1 ! ! = 3 as =2 with distinct disspipative an mixing length scale ln_mxl0 = .true. ! surface mixing length scale = F(wind stress) (T) or not (F) rn_mxl0 = 0.04 ! surface buoyancy lenght scale minimum value ln_drg = .false. ! top/bottom friction added as boundary condition of TKE ln_lc = .true. ! Langmuir cell parameterisation (Axell 2002) rn_lc = 0.15 ! coef. associated to Langmuir cells nn_etau = 1 ! penetration of tke below the mixed layer (ML) due to NIWs ! = 0 none ; = 1 add a tke source below the ML ! = 2 add a tke source just at the base of the ML ! = 3 as = 1 applied on HF part of the stress (ln_cpl=T) rn_efr = 0.05 ! fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2) nn_htau = 1 ! type of exponential decrease of tke penetration below the ML ! = 0 constant 10 m length scale ! = 1 0.5m at the equator to 30m poleward of 40 degrees / !----------------------------------------------------------------------- &namzdf_gls ! GLS vertical diffusion (ln_zdfgls =T) !----------------------------------------------------------------------- rn_emin = 1.e-7 ! minimum value of e [m2/s2] rn_epsmin = 1.e-12 ! minimum value of eps [m2/s3] ln_length_lim = .true. ! limit on the dissipation rate under stable stratification (Galperin et al., 1988) rn_clim_galp = 0.267 ! galperin limit ln_sigpsi = .true. ! Activate or not Burchard 2001 mods on psi schmidt number in the wb case rn_crban = 100. ! Craig and Banner 1994 constant for wb tke flux rn_charn = 70000. ! Charnock constant for wb induced roughness length rn_hsro = 0.02 ! Minimum surface roughness rn_frac_hs = 1.3 ! Fraction of wave height as roughness (if nn_z0_met>1) nn_z0_met = 2 ! Method for surface roughness computation (0/1/2/3) ! ! =3 requires ln_wave=T nn_bc_surf = 1 ! surface condition (0/1=Dir/Neum) nn_bc_bot = 1 ! bottom condition (0/1=Dir/Neum) nn_stab_func = 2 ! stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB) nn_clos = 1 ! predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen) / !----------------------------------------------------------------------- &namzdf_osm ! OSM vertical diffusion (ln_zdfosm =T) !----------------------------------------------------------------------- ln_use_osm_la = .false. ! Use namelist rn_osm_la rn_osm_la = 0.3 ! Turbulent Langmuir number rn_osm_dstokes = 5. ! Depth scale of Stokes drift (m) nn_ave = 0 ! choice of horizontal averaging on avt, avmu, avmv ln_dia_osm = .true. ! output OSMOSIS-OBL variables rn_osm_hbl0 = 10. ! initial hbl value ln_kpprimix = .true. ! Use KPP-style Ri# mixing below BL rn_riinfty = 0.7 ! Highest local Ri_g permitting shear instability rn_difri = 0.005 ! max Ri# diffusivity at Ri_g = 0 (m^2/s) ln_convmix = .true. ! Use convective instability mixing below BL rn_difconv = 1. ! diffusivity when unstable below BL (m2/s) nn_osm_wave = 0 ! Method used to calculate Stokes drift ! = 2: Use ECMWF wave fields ! = 1: Pierson Moskowitz wave spectrum ! = 0: Constant La# = 0.3 / !----------------------------------------------------------------------- &namzdf_iwm ! internal wave-driven mixing parameterization (ln_zdfiwm =T) !----------------------------------------------------------------------- nn_zpyc = 1 ! pycnocline-intensified dissipation scales as N (=1) or N^2 (=2) ln_mevar = .true. ! variable (T) or constant (F) mixing efficiency ln_tsdiff = .true. ! account for differential T/S mixing (T) or not (F) / !!====================================================================== !! *** Diagnostics namelists *** !! !! !! !! namtrd dynamics and/or tracer trends (default NO) !! namptr Poleward Transport Diagnostics (default NO) !! namhsb Heat and salt budgets (default NO) !! namdiu Cool skin and warm layer models (default NO) !! namdiu Cool skin and warm layer models (default NO) !! namflo float parameters ("key_float") !! nam_diaharm Harmonic analysis of tidal constituents ("key_diaharm") !! namdct transports through some sections ("key_diadct") !! nam_diatmb Top Middle Bottom Output (default NO) !! nam_dia25h 25h Mean Output (default NO) !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") !!====================================================================== ! !----------------------------------------------------------------------- &namtrd ! trend diagnostics (default NO) !----------------------------------------------------------------------- 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 !----------------------------------------------------------------------- &namptr ! Poleward Transport Diagnostic (default NO) !----------------------------------------------------------------------- ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F) ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not / !----------------------------------------------------------------------- &namhsb ! Heat and salt budgets (default NO) !----------------------------------------------------------------------- ln_diahsb = .false. ! check the heat and salt budgets (T) or not (F) / !----------------------------------------------------------------------- &namdiu ! Cool skin and warm layer models (default NO) !----------------------------------------------------------------------- ln_diurnal = .false. ! ln_diurnal_only = .false. ! / !----------------------------------------------------------------------- &namflo ! float parameters ("key_float") !----------------------------------------------------------------------- jpnfl = 1 ! total number of floats during the run jpnnewflo = 0 ! number of floats for the restart ln_rstflo = .false. ! float restart (T) or not (F) nn_writefl = 75 ! frequency of writing in float output file nn_stockfl = 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) ln_ariane = .true. ! Input with Ariane tool convention(T) ln_flo_ascii= .true. ! Output with Ariane tool netcdf convention(F) or ascii file (T) / !----------------------------------------------------------------------- &nam_diaharm ! Harmonic analysis of tidal constituents ("key_diaharm") !----------------------------------------------------------------------- nit000_han = 1 ! First time step used for harmonic analysis nitend_han = 75 ! Last time step used for harmonic analysis nstep_han = 15 ! Time step frequency for harmonic analysis tname(1) = 'M2' ! Name of tidal constituents tname(2) = 'K1' / !----------------------------------------------------------------------- &namdct ! transports through some sections ("key_diadct") !----------------------------------------------------------------------- nn_dct = 15 ! time step frequency for transports computing nn_dctwri = 15 ! time step frequency for transports writing nn_secdebug= 112 ! 0 : no section to debug ! ! -1 : debug all section ! ! 0 < n : debug section number n / !----------------------------------------------------------------------- &nam_diatmb ! Top Middle Bottom Output (default NO) !----------------------------------------------------------------------- ln_diatmb = .false. ! Choose Top Middle and Bottom output or not / !----------------------------------------------------------------------- &nam_dia25h ! 25h Mean Output (default NO) !----------------------------------------------------------------------- ln_dia25h = .false. ! Choose 25h mean output or not / !----------------------------------------------------------------------- &namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4") !----------------------------------------------------------------------- nn_nchunks_i= 4 ! number of chunks in i-dimension nn_nchunks_j= 4 ! number of chunks in j-dimension nn_nchunks_k= 31 ! number of chunks in k-dimension ! ! setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which ! ! is optimal for postprocessing which works exclusively with horizontal slabs ln_nc4zip = .true. ! (T) use netcdf4 chunking and compression ! ! (F) ignore chunking information and produce netcdf3-compatible files / !!====================================================================== !! *** Miscellaneous namelists *** !! !! !! !! nammpp Massively Parallel Processing ("key_mpp_mpi") !! namctl Control prints (default NO) !! namsto Stochastic parametrization of EOS (default NO) !!====================================================================== ! !----------------------------------------------------------------------- &nammpp ! Massively Parallel Processing ("key_mpp_mpi") !----------------------------------------------------------------------- cn_mpi_send = 'I' ! mpi send/recieve type ='S', 'B', or 'I' for standard send, ! ! buffer blocking send or immediate non-blocking sends, resp. nn_buffer = 0 ! size in bytes of exported buffer ('B' case), 0 no exportation ln_nnogather= .false. ! activate code to avoid mpi_allgather use at the northfold jpni = 0 ! jpni number of processors following i (set automatically if < 1) jpnj = 0 ! jpnj number of processors following j (set automatically if < 1) jpnij = 0 ! jpnij number of local domains (set automatically if < 1) / !----------------------------------------------------------------------- &namctl ! Control prints (default: NO) !----------------------------------------------------------------------- ln_ctl = .false. ! trends control print (expensive!) nn_print = 0 ! level of print (0 no extra print) nn_ictls = 0 ! start i indice of control sum (use to compare mono versus nn_ictle = 0 ! end i indice of control sum multi processor runs nn_jctls = 0 ! start j indice of control over a subdomain) nn_jctle = 0 ! end j indice of control nn_isplt = 1 ! number of processors in i-direction nn_jsplt = 1 ! number of processors in j-direction ln_timing = .false. ! timing by routine write out in timing.output file ln_diacfl = .false. ! CFL diagnostics write out in cfl_diagnostics.ascii / !----------------------------------------------------------------------- &namsto ! Stochastic parametrization of EOS (default: NO) !----------------------------------------------------------------------- ln_sto_eos = .false. ! stochastic equation of state nn_sto_eos = 1 ! number of independent random walks rn_eos_stdxy= 1.4 ! random walk horz. standard deviation (in grid points) rn_eos_stdz = 0.7 ! random walk vert. standard deviation (in grid points) rn_eos_tcor = 1440. ! random walk time correlation (in timesteps) nn_eos_ord = 1 ! order of autoregressive processes nn_eos_flt = 0 ! passes of Laplacian filter rn_eos_lim = 2.0 ! limitation factor (default = 3.0) ln_rststo = .false. ! start from mean parameter (F) or from restart file (T) ln_rstseed = .true. ! read seed of RNG from restart file cn_storst_in = "restart_sto" ! suffix of stochastic parameter restart file (input) cn_storst_out = "restart_sto" ! suffix of stochastic parameter restart file (output) / !!====================================================================== !! *** Observation & Assimilation *** !! !! !! !! namobs observation and model comparison (default: NO) !! nam_asminc assimilation increments ('key_asminc') !!====================================================================== ! !----------------------------------------------------------------------- &namobs ! observation usage switch (default: NO) !----------------------------------------------------------------------- ln_diaobs = .false. ! Logical switch for the observation operator ! ln_t3d = .false. ! Logical switch for T profile observations ln_s3d = .false. ! Logical switch for S profile observations ln_sla = .false. ! Logical switch for SLA observations ln_sst = .false. ! Logical switch for SST observations ln_sss = .false. ! Logical swithc for SSS observations ln_sic = .false. ! Logical switch for Sea Ice observations ln_vel3d = .false. ! Logical switch for velocity observations ln_altbias = .false. ! Logical switch for altimeter bias correction ln_sstbias = .false. ! Logical switch for SST bias correction ln_nea = .false. ! Logical switch for rejection of observations near land ln_grid_global = .true. ! Logical switch for global distribution of observations ln_grid_search_lookup = .false. ! Logical switch for obs grid search w/lookup table ln_ignmis = .true. ! Logical switch for ignoring missing files ln_s_at_t = .false. ! Logical switch for computing model S at T obs if not there ln_sstnight = .false. ! Logical switch for calculating night-time average for SST obs ln_sla_fp_indegs = .true. ! Logical for SLA: T=> averaging footprint is in degrees, F=> in metres ln_sst_fp_indegs = .true. ! Logical for SST: T=> averaging footprint is in degrees, F=> in metres ln_sss_fp_indegs = .true. ! Logical for SSS: T=> averaging footprint is in degrees, F=> in metres ln_sic_fp_indegs = .true. ! Logical for SIC: T=> averaging footprint is in degrees, F=> in metres ! All of the *files* variables below are arrays. Use namelist_cfg to add more files cn_profbfiles = 'profiles_01.nc' ! Profile feedback input observation file names cn_slafbfiles = 'sla_01.nc' ! SLA feedback input observation file names cn_sstfbfiles = 'sst_01.nc' ! SST feedback input observation file names cn_sssfbfiles = 'sss_01.nc' ! SSS feedback input observation file names cn_sicfbfiles = 'sic_01.nc' ! SIC feedback input observation file names cn_velfbfiles = 'vel_01.nc' ! Velocity feedback input observation file names cn_altbiasfile = 'altbias.nc' ! Altimeter bias input file name cn_sstbiasfiles = 'sstbias.nc' ! SST bias input file name cn_gridsearchfile='gridsearch.nc' ! Grid search file name rn_gridsearchres = 0.5 ! Grid search resolution rn_mdtcorr = 1.61 ! MDT correction rn_mdtcutoff = 65.0 ! MDT cutoff for computed correction rn_dobsini = 00010101.000000 ! Initial date in window YYYYMMDD.HHMMSS rn_dobsend = 00010102.000000 ! Final date in window YYYYMMDD.HHMMSS rn_sla_avglamscl = 0. ! E/W diameter of SLA observation footprint (metres/degrees) rn_sla_avgphiscl = 0. ! N/S diameter of SLA observation footprint (metres/degrees) rn_sst_avglamscl = 0. ! E/W diameter of SST observation footprint (metres/degrees) rn_sst_avgphiscl = 0. ! N/S diameter of SST observation footprint (metres/degrees) rn_sss_avglamscl = 0. ! E/W diameter of SSS observation footprint (metres/degrees) rn_sss_avgphiscl = 0. ! N/S diameter of SSS observation footprint (metres/degrees) rn_sic_avglamscl = 0. ! E/W diameter of SIC observation footprint (metres/degrees) rn_sic_avgphiscl = 0. ! N/S diameter of SIC observation footprint (metres/degrees) nn_1dint = 0 ! Type of vertical interpolation method nn_2dint = 0 ! Default horizontal interpolation method nn_2dint_sla = 0 ! Horizontal interpolation method for SLA nn_2dint_sst = 0 ! Horizontal interpolation method for SST nn_2dint_sss = 0 ! Horizontal interpolation method for SSS nn_2dint_sic = 0 ! Horizontal interpolation method for SIC nn_msshc = 0 ! MSSH correction scheme nn_profdavtypes = -1 ! Profile daily average types - array / !----------------------------------------------------------------------- &nam_asminc ! assimilation increments ('key_asminc') !----------------------------------------------------------------------- ln_bkgwri = .false. ! Logical switch for writing out background state ln_trainc = .false. ! Logical switch for applying tracer increments ln_dyninc = .false. ! Logical switch for applying velocity increments ln_sshinc = .false. ! Logical switch for applying SSH increments ln_asmdin = .false. ! Logical switch for Direct Initialization (DI) ln_asmiau = .false. ! Logical switch for Incremental Analysis Updating (IAU) nitbkg = 0 ! Timestep of background in [0,nitend-nit000-1] nitdin = 0 ! Timestep of background for DI in [0,nitend-nit000-1] nitiaustr = 1 ! Timestep of start of IAU interval in [0,nitend-nit000-1] nitiaufin = 15 ! Timestep of end of IAU interval in [0,nitend-nit000-1] niaufn = 0 ! Type of IAU weighting function ln_salfix = .false. ! Logical switch for ensuring that the sa > salfixmin salfixmin = -9999 ! Minimum salinity after applying the increments nn_divdmp = 0 ! Number of iterations of divergence damping operator /