Changeset 10727 for utils/tools_AGRIF_CMEMS_2020/DOMAINcfg/namelist_ref
- Timestamp:
- 2019-02-27T17:02:02+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
utils/tools_AGRIF_CMEMS_2020/DOMAINcfg/namelist_ref
r9051 r10727 171 171 ppacr2 = 13.000000000000 ! 172 172 / 173 !-----------------------------------------------------------------------174 &namwad ! Wetting and drying (default F)175 !-----------------------------------------------------------------------176 ln_wd = .false. ! T/F activation of wetting and drying177 rn_wdmin1 = 0.1 ! Minimum wet depth on dried cells178 rn_wdmin2 = 0.01 ! Tolerance of min wet depth on dried cells179 rn_wdld = 20.0 ! Land elevation below which wetting/drying is allowed180 nn_wdit = 10 ! Max iterations for W/D limiter181 /182 !-----------------------------------------------------------------------183 &namtsd ! data : Temperature & Salinity184 !-----------------------------------------------------------------------185 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !186 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !187 sn_tem = 'data_1m_potential_temperature_nomask', -1 ,'votemper', .true. , .true. , 'yearly' , '' , '' , ''188 sn_sal = 'data_1m_salinity_nomask' , -1 ,'vosaline', .true. , .true. , 'yearly' , '' , '' , ''189 !190 cn_dir = './' ! root directory for the location of the runoff files191 ln_tsd_init = .true. ! Initialisation of ocean T & S with T & S input data (T) or not (F)192 ln_tsd_tradmp = .true. ! damping of ocean T & S toward T & S input data (T) or not (F)193 /194 !-----------------------------------------------------------------------195 &namcrs ! coarsened grid (for outputs and/or TOP) ("key_crs")196 !-----------------------------------------------------------------------197 nn_factx = 3 ! Reduction factor of x-direction198 nn_facty = 3 ! Reduction factor of y-direction199 nn_binref = 0 ! Bin centering preference: NORTH or EQUAT200 ! 0, coarse grid is binned with preferential treatment of the north fold201 ! 1, coarse grid is binned with centering at the equator202 ! Symmetry with nn_facty being odd-numbered. Asymmetry with even-numbered nn_facty.203 nn_msh_crs = 1 ! create (=1) a mesh file or not (=0)204 nn_crs_kz = 0 ! 0, MEAN of volume boxes205 ! 1, MAX of boxes206 ! 2, MIN of boxes207 ln_crs_wn = .true. ! wn coarsened (T) or computed using horizontal divergence ( F )208 /209 !-----------------------------------------------------------------------210 &namc1d ! 1D configuration options ("key_c1d")211 !-----------------------------------------------------------------------212 rn_lat1d = 50 ! Column latitude (default at PAPA station)213 rn_lon1d = -145 ! Column longitude (default at PAPA station)214 ln_c1d_locpt= .true. ! Localization of 1D config in a grid (T) or independant point (F)215 /216 !-----------------------------------------------------------------------217 &namc1d_dyndmp ! U & V newtonian damping ("key_c1d")218 !-----------------------------------------------------------------------219 ln_dyndmp = .false. ! add a damping term (T) or not (F)220 /221 !-----------------------------------------------------------------------222 &namc1d_uvd ! data: U & V currents ("key_c1d")223 !-----------------------------------------------------------------------224 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !225 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !226 sn_ucur = 'ucurrent' , -1 ,'u_current', .false. , .true. , 'monthly' , '' , 'Ume' , ''227 sn_vcur = 'vcurrent' , -1 ,'v_current', .false. , .true. , 'monthly' , '' , 'Vme' , ''228 !229 cn_dir = './' ! root directory for the location of the files230 ln_uvd_init = .false. ! Initialisation of ocean U & V with U & V input data (T) or not (F)231 ln_uvd_dyndmp = .false. ! damping of ocean U & V toward U & V input data (T) or not (F)232 /233 234 !!======================================================================235 !! *** Surface Boundary Condition namelists ***236 !!======================================================================237 !! namsbc surface boundary condition238 !! namsbc_ana analytical formulation (ln_ana =T)239 !! namsbc_flx flux formulation (ln_flx =T)240 !! namsbc_clio CLIO bulk formulae formulation (ln_blk_clio=T)241 !! namsbc_core CORE bulk formulae formulation (ln_blk_core=T)242 !! namsbc_mfs MFS bulk formulae formulation (ln_blk_mfs =T)243 !! namsbc_cpl CouPLed formulation ("key_oasis3" )244 !! namsbc_sas StAndalone Surface module245 !! namtra_qsr penetrative solar radiation (ln_traqsr =T)246 !! namsbc_rnf river runoffs (ln_rnf =T)247 !! namsbc_isf ice shelf melting/freezing (nn_isf >0)248 !! namsbc_iscpl coupling option between land ice model and ocean249 !! namsbc_apr Atmospheric Pressure (ln_apr_dyn =T)250 !! namsbc_ssr sea surface restoring term (for T and/or S) (ln_ssr =T)251 !! namsbc_alb albedo parameters252 !! namsbc_wave external fields from wave model (ln_wave =T)253 !! namberg iceberg floats (ln_icebergs=T)254 !!======================================================================255 !256 !-----------------------------------------------------------------------257 &namsbc ! Surface Boundary Condition (surface module)258 !-----------------------------------------------------------------------259 nn_fsbc = 5 ! frequency of surface boundary condition computation260 ! (also = the frequency of sea-ice & iceberg model call)261 ! Type of air-sea fluxes262 ln_ana = .false. ! analytical formulation (T => fill namsbc_ana )263 ln_flx = .false. ! flux formulation (T => fill namsbc_flx )264 ln_blk_clio = .false. ! CLIO bulk formulation (T => fill namsbc_clio)265 ln_blk_core = .true. ! CORE bulk formulation (T => fill namsbc_core)266 ln_blk_mfs = .false. ! MFS bulk formulation (T => fill namsbc_mfs )267 ! Type of coupling (Ocean/Ice/Atmosphere) :268 ln_cpl = .false. ! atmosphere coupled formulation ( requires key_oasis3 )269 ln_mixcpl = .false. ! forced-coupled mixed formulation ( requires key_oasis3 )270 nn_components = 0 ! configuration of the opa-sas OASIS coupling271 ! =0 no opa-sas OASIS coupling: default single executable configuration272 ! =1 opa-sas OASIS coupling: multi executable configuration, OPA component273 ! =2 opa-sas OASIS coupling: multi executable configuration, SAS component274 nn_limflx = -1 ! LIM3 Multi-category heat flux formulation (use -1 if LIM3 is not used)275 ! =-1 Use per-category fluxes, bypass redistributor, forced mode only, not yet implemented coupled276 ! = 0 Average per-category fluxes (forced and coupled mode)277 ! = 1 Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled278 ! = 2 Redistribute a single flux over categories (coupled mode only)279 ! Sea-ice :280 nn_ice = 2 ! =0 no ice boundary condition ,281 ! =1 use observed ice-cover ,282 ! =2 ice-model used ("key_lim3", "key_lim2", "key_cice")283 nn_ice_embd = 1 ! =0 levitating ice (no mass exchange, concentration/dilution effect)284 ! =1 levitating ice with mass and salt exchange but no presure effect285 ! =2 embedded sea-ice (full salt and mass exchanges and pressure)286 ! Misc. options of sbc :287 ln_traqsr = .true. ! Light penetration in the ocean (T => fill namtra_qsr )288 ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave289 ln_rnf = .true. ! runoffs (T => fill namsbc_rnf)290 ln_ssr = .true. ! Sea Surface Restoring on T and/or S (T => fill namsbc_ssr)291 nn_fwb = 2 ! FreshWater Budget: =0 unchecked292 ! =1 global mean of e-p-r set to zero at each time step293 ! =2 annual global mean of e-p-r set to zero294 ln_apr_dyn = .false. ! Patm gradient added in ocean & ice Eqs. (T => fill namsbc_apr )295 ln_isf = .false. ! ice shelf (T => fill namsbc_isf)296 ln_wave = .false. ! coupling with surface wave (T => fill namsbc_wave)297 nn_lsm = 0 ! =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) ,298 ! =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field)299 /300 !-----------------------------------------------------------------------301 &namsbc_ana ! analytical surface boundary condition302 !-----------------------------------------------------------------------303 nn_tau000 = 0 ! gently increase the stress over the first ntau_rst time-steps304 rn_utau0 = 0.5 ! uniform value for the i-stress305 rn_vtau0 = 0.e0 ! uniform value for the j-stress306 rn_qns0 = 0.e0 ! uniform value for the total heat flux307 rn_qsr0 = 0.e0 ! uniform value for the solar radiation308 rn_emp0 = 0.e0 ! uniform value for the freswater budget (E-P)309 /310 !-----------------------------------------------------------------------311 &namsbc_flx ! surface boundary condition : flux formulation312 !-----------------------------------------------------------------------313 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !314 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !315 sn_utau = 'utau' , 24 , 'utau' , .false. , .false., 'yearly' , '' , '' , ''316 sn_vtau = 'vtau' , 24 , 'vtau' , .false. , .false., 'yearly' , '' , '' , ''317 sn_qtot = 'qtot' , 24 , 'qtot' , .false. , .false., 'yearly' , '' , '' , ''318 sn_qsr = 'qsr' , 24 , 'qsr' , .false. , .false., 'yearly' , '' , '' , ''319 sn_emp = 'emp' , 24 , 'emp' , .false. , .false., 'yearly' , '' , '' , ''320 321 cn_dir = './' ! root directory for the location of the flux files322 /323 !-----------------------------------------------------------------------324 &namsbc_clio ! namsbc_clio CLIO bulk formulae325 !-----------------------------------------------------------------------326 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !327 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !328 sn_utau = 'taux_1m' , -1 , 'sozotaux', .true. , .true. , 'yearly' , '' , '' , ''329 sn_vtau = 'tauy_1m' , -1 , 'sometauy', .true. , .true. , 'yearly' , '' , '' , ''330 sn_wndm = 'flx' , -1 , 'socliowi', .true. , .true. , 'yearly' , '' , '' , ''331 sn_tair = 'flx' , -1 , 'socliot2', .true. , .true. , 'yearly' , '' , '' , ''332 sn_humi = 'flx' , -1 , 'socliohu', .true. , .true. , 'yearly' , '' , '' , ''333 sn_ccov = 'flx' , -1 , 'socliocl', .false. , .true. , 'yearly' , '' , '' , ''334 sn_prec = 'flx' , -1 , 'socliopl', .false. , .true. , 'yearly' , '' , '' , ''335 336 cn_dir = './' ! root directory for the location of the bulk files are337 /338 !-----------------------------------------------------------------------339 &namsbc_core ! namsbc_core CORE bulk formulae340 !-----------------------------------------------------------------------341 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !342 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !343 sn_wndi = 'u_10.15JUNE2009_fill' , 6 , 'U_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bicubic_noc.nc' , 'Uwnd' , ''344 sn_wndj = 'v_10.15JUNE2009_fill' , 6 , 'V_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bicubic_noc.nc' , 'Vwnd' , ''345 sn_qsr = 'ncar_rad.15JUNE2009_fill' , 24 , 'SWDN_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''346 sn_qlw = 'ncar_rad.15JUNE2009_fill' , 24 , 'LWDN_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''347 sn_tair = 't_10.15JUNE2009_fill' , 6 , 'T_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''348 sn_humi = 'q_10.15JUNE2009_fill' , 6 , 'Q_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''349 sn_prec = 'ncar_precip.15JUNE2009_fill', -1 , 'PRC_MOD1', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''350 sn_snow = 'ncar_precip.15JUNE2009_fill', -1 , 'SNOW' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''351 sn_tdif = 'taudif_core' , 24 , 'taudif' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , ''352 353 cn_dir = './' ! root directory for the location of the bulk files354 ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data355 rn_zqt = 10. ! Air temperature and humidity reference height (m)356 rn_zu = 10. ! Wind vector reference height (m)357 rn_pfac = 1. ! multiplicative factor for precipitation (total & snow)358 rn_efac = 1. ! multiplicative factor for evaporation (0. or 1.)359 rn_vfac = 0. ! multiplicative factor for ocean/ice velocity360 ! in the calculation of the wind stress (0.=absolute winds or 1.=relative winds)361 /362 !-----------------------------------------------------------------------363 &namsbc_mfs ! namsbc_mfs MFS bulk formulae364 !-----------------------------------------------------------------------365 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !366 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !367 sn_wndi = 'ecmwf' , 6 , 'u10' , .true. , .false., 'daily' ,'bicubic.nc' , '' , ''368 sn_wndj = 'ecmwf' , 6 , 'v10' , .true. , .false., 'daily' ,'bicubic.nc' , '' , ''369 sn_clc = 'ecmwf' , 6 , 'clc' , .true. , .false., 'daily' ,'bilinear.nc', '' , ''370 sn_msl = 'ecmwf' , 6 , 'msl' , .true. , .false., 'daily' ,'bicubic.nc' , '' , ''371 sn_tair = 'ecmwf' , 6 , 't2' , .true. , .false., 'daily' ,'bicubic.nc' , '' , ''372 sn_rhm = 'ecmwf' , 6 , 'rh' , .true. , .false., 'daily' ,'bilinear.nc', '' , ''373 sn_prec = 'ecmwf' , 6 , 'precip' , .true. , .true. , 'daily' ,'bicubic.nc' , '' , ''374 375 cn_dir = './ECMWF/' ! root directory for the location of the bulk files376 /377 !-----------------------------------------------------------------------378 &namsbc_cpl ! coupled ocean/atmosphere model ("key_oasis3")379 !-----------------------------------------------------------------------380 ! ! description ! multiple ! vector ! vector ! vector !381 ! ! ! categories ! reference ! orientation ! grids !382 ! send383 sn_snd_temp = 'weighted oce and ice' , 'no' , '' , '' , ''384 sn_snd_alb = 'weighted ice' , 'no' , '' , '' , ''385 sn_snd_thick = 'none' , 'no' , '' , '' , ''386 sn_snd_crt = 'none' , 'no' , 'spherical' , 'eastward-northward' , 'T'387 sn_snd_co2 = 'coupled' , 'no' , '' , '' , ''388 ! receive389 sn_rcv_w10m = 'none' , 'no' , '' , '' , ''390 sn_rcv_taumod = 'coupled' , 'no' , '' , '' , ''391 sn_rcv_tau = 'oce only' , 'no' , 'cartesian' , 'eastward-northward', 'U,V'392 sn_rcv_dqnsdt = 'coupled' , 'no' , '' , '' , ''393 sn_rcv_qsr = 'oce and ice' , 'no' , '' , '' , ''394 sn_rcv_qns = 'oce and ice' , 'no' , '' , '' , ''395 sn_rcv_emp = 'conservative' , 'no' , '' , '' , ''396 sn_rcv_rnf = 'coupled' , 'no' , '' , '' , ''397 sn_rcv_cal = 'coupled' , 'no' , '' , '' , ''398 sn_rcv_co2 = 'coupled' , 'no' , '' , '' , ''399 !400 nn_cplmodel = 1 ! Maximum number of models to/from which NEMO is potentialy sending/receiving data401 ln_usecplmask = .false. ! use a coupling mask file to merge data received from several models402 ! ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel)403 /404 !-----------------------------------------------------------------------405 &namsbc_sas ! analytical surface boundary condition406 !-----------------------------------------------------------------------407 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !408 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !409 sn_usp = 'sas_grid_U', 120 , 'vozocrtx', .true. , .true. , 'yearly' , '' , '' , ''410 sn_vsp = 'sas_grid_V', 120 , 'vomecrty', .true. , .true. , 'yearly' , '' , '' , ''411 sn_tem = 'sas_grid_T', 120 , 'sosstsst', .true. , .true. , 'yearly' , '' , '' , ''412 sn_sal = 'sas_grid_T', 120 , 'sosaline', .true. , .true. , 'yearly' , '' , '' , ''413 sn_ssh = 'sas_grid_T', 120 , 'sossheig', .true. , .true. , 'yearly' , '' , '' , ''414 sn_e3t = 'sas_grid_T', 120 , 'e3t_m' , .true. , .true. , 'yearly' , '' , '' , ''415 sn_frq = 'sas_grid_T', 120 , 'frq_m' , .true. , .true. , 'yearly' , '' , '' , ''416 417 ln_3d_uve = .true. ! specify whether we are supplying a 3D u,v and e3 field418 ln_read_frq = .false. ! specify whether we must read frq or not419 cn_dir = './' ! root directory for the location of the bulk files are420 /421 !-----------------------------------------------------------------------422 &namtra_qsr ! penetrative solar radiation (ln_traqsr=T)423 !-----------------------------------------------------------------------424 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !425 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !426 sn_chl ='chlorophyll', -1 , 'CHLA' , .true. , .true. , 'yearly' , '' , '' , ''427 428 cn_dir = './' ! root directory for the location of the runoff files429 ln_qsr_rgb = .true. ! RGB (Red-Green-Blue) light penetration430 ln_qsr_2bd = .false. ! 2 bands light penetration431 ln_qsr_bio = .false. ! bio-model light penetration432 nn_chldta = 1 ! RGB : Chl data (=1) or cst value (=0)433 rn_abs = 0.58 ! RGB & 2 bands: fraction of light (rn_si1)434 rn_si0 = 0.35 ! RGB & 2 bands: shortess depth of extinction435 rn_si1 = 23.0 ! 2 bands: longest depth of extinction436 ln_qsr_ice = .true. ! light penetration for ice-model LIM3437 /438 !-----------------------------------------------------------------------439 &namsbc_rnf ! runoffs namelist surface boundary condition (ln_rnf=T)440 !-----------------------------------------------------------------------441 ! ! file name ! frequency (hours) ! variable ! time interp. ! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !442 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !443 sn_rnf = 'runoff_core_monthly', -1 , 'sorunoff', .true. , .true. , 'yearly' , '' , '' , ''444 sn_cnf = 'runoff_core_monthly', 0 , 'socoefr0', .false. , .true. , 'yearly' , '' , '' , ''445 sn_s_rnf = 'runoffs' , 24 , 'rosaline', .true. , .true. , 'yearly' , '' , '' , ''446 sn_t_rnf = 'runoffs' , 24 , 'rotemper', .true. , .true. , 'yearly' , '' , '' , ''447 sn_dep_rnf = 'runoffs' , 0 , 'rodepth' , .false. , .true. , 'yearly' , '' , '' , ''448 449 cn_dir = './' ! root directory for the location of the runoff files450 ln_rnf_mouth= .true. ! specific treatment at rivers mouths451 rn_hrnf = 15.e0 ! depth over which enhanced vertical mixing is used (ln_rnf_mouth=T)452 rn_avt_rnf = 1.e-3 ! value of the additional vertical mixing coef. [m2/s] (ln_rnf_mouth=T)453 rn_rfact = 1.e0 ! multiplicative factor for runoff454 ln_rnf_depth= .false. ! read in depth information for runoff455 ln_rnf_tem = .false. ! read in temperature information for runoff456 ln_rnf_sal = .false. ! read in salinity information for runoff457 ln_rnf_depth_ini = .false. ! compute depth at initialisation from runoff file458 rn_rnf_max = 5.735e-4 ! max value of the runoff climatologie over global domain ( ln_rnf_depth_ini = .true )459 rn_dep_max = 150. ! depth over which runoffs is spread ( ln_rnf_depth_ini = .true )460 nn_rnf_depth_file = 0 ! create (=1) a runoff depth file or not (=0)461 /462 !-----------------------------------------------------------------------463 &namsbc_isf ! Top boundary layer (ISF) (nn_isf >0)464 !-----------------------------------------------------------------------465 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !466 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !467 ! nn_isf == 4468 sn_fwfisf = 'rnfisf' , -12 ,'sowflisf', .false. , .true. , 'yearly' , '' , '' , ''469 ! nn_isf == 3470 sn_rnfisf = 'rnfisf' , -12 ,'sofwfisf', .false. , .true. , 'yearly' , '' , '' , ''471 ! nn_isf == 2 and 3472 sn_depmax_isf='rnfisf' , -12 ,'sozisfmax', .false. , .true. , 'yearly' , '' , '' , ''473 sn_depmin_isf='rnfisf' , -12 ,'sozisfmin', .false. , .true. , 'yearly' , '' , '' , ''474 ! nn_isf == 2475 sn_Leff_isf = 'rnfisf' , -12 ,'Leff' , .false. , .true. , 'yearly' , '' , '' , ''476 !477 ! for all case478 nn_isf = 1 ! ice shelf melting/freezing479 ! 1 = presence of ISF 2 = bg03 parametrisation480 ! 3 = rnf file for isf 4 = ISF fwf specified481 ! option 1 and 4 need ln_isfcav = .true. (domzgr)482 ! only for nn_isf = 1 or 2483 rn_gammat0 = 1.e-4 ! gammat coefficient used in blk formula484 rn_gammas0 = 1.e-4 ! gammas coefficient used in blk formula485 ! only for nn_isf = 1 or 4486 rn_hisf_tbl = 30. ! thickness of the top boundary layer (Losh et al. 2008)487 ! ! 0 => thickness of the tbl = thickness of the first wet cell488 ! only for nn_isf = 1489 nn_isfblk = 1 ! 1 ISOMIP like: 2 equations formulation (Hunter et al., 2006)490 ! ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015)491 nn_gammablk = 1 ! 0 = cst Gammat (= gammat/s)492 ! ! 1 = velocity dependend Gamma (u* * gammat/s) (Jenkins et al. 2010)493 ! ! 2 = velocity and stability dependent Gamma (Holland et al. 1999)494 /495 !-----------------------------------------------------------------------496 &namsbc_iscpl ! land ice / ocean coupling option497 !-----------------------------------------------------------------------498 nn_drown = 10 ! number of iteration of the extrapolation loop (fill the new wet cells)499 ln_hsb = .false. ! activate conservation module (conservation exact after a time of rn_fiscpl)500 nn_fiscpl = 43800 ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency)501 /502 !-----------------------------------------------------------------------503 &namsbc_apr ! Atmospheric pressure used as ocean forcing or in bulk504 !-----------------------------------------------------------------------505 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !506 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !507 sn_apr = 'patm' , -1 ,'somslpre', .true. , .true. , 'yearly' , '' , '' , ''508 509 cn_dir = './' ! root directory for the location of the bulk files510 rn_pref = 101000. ! reference atmospheric pressure [N/m2]/511 ln_ref_apr = .false. ! ref. pressure: global mean Patm (T) or a constant (F)512 ln_apr_obc = .false. ! inverse barometer added to OBC ssh data513 /514 !-----------------------------------------------------------------------515 &namsbc_ssr ! surface boundary condition : sea surface restoring (ln_ssr=T)516 !-----------------------------------------------------------------------517 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !518 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !519 sn_sst = 'sst_data', 24 , 'sst' , .false. , .false., 'yearly' , '' , '' , ''520 sn_sss = 'sss_data', -1 , 'sss' , .true. , .true. , 'yearly' , '' , '' , ''521 522 cn_dir = './' ! root directory for the location of the runoff files523 nn_sstr = 0 ! add a retroaction term in the surface heat flux (=1) or not (=0)524 nn_sssr = 2 ! add a damping term in the surface freshwater flux (=2)525 ! or to SSS only (=1) or no damping term (=0)526 rn_dqdt = -40. ! magnitude of the retroaction on temperature [W/m2/K]527 rn_deds = -166.67 ! magnitude of the damping on salinity [mm/day]528 ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2)529 rn_sssr_bnd = 4.e0 ! ABS(Max/Min) value of the damping erp term [mm/day]530 /531 !-----------------------------------------------------------------------532 &namsbc_alb ! albedo parameters533 !-----------------------------------------------------------------------534 nn_ice_alb = 0 ! parameterization of ice/snow albedo535 ! 0: Shine & Henderson-Sellers (JGR 1985)536 ! 1: "home made" based on Brandt et al. (J. Climate 2005)537 ! and Grenfell & Perovich (JGR 2004)538 rn_albice = 0.53 ! albedo of bare puddled ice (values from 0.49 to 0.58)539 ! 0.53 (default) => if nn_ice_alb=0540 ! 0.50 (default) => if nn_ice_alb=1541 /542 !-----------------------------------------------------------------------543 &namsbc_wave ! External fields from wave model (ln_wave=T)544 !-----------------------------------------------------------------------545 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !546 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename !547 sn_cdg = 'cdg_wave', 1 , 'drag_coeff', .true. , .false., 'daily' , '' , '' , ''548 sn_usd = 'sdw_wave', 1 , 'u_sd2d' , .true. , .false., 'daily' , '' , '' , ''549 sn_vsd = 'sdw_wave', 1 , 'v_sd2d' , .true. , .false., 'daily' , '' , '' , ''550 sn_wn = 'sdw_wave', 1 , 'wave_num' , .true. , .false., 'daily' , '' , '' , ''551 !552 cn_dir_cdg = './' ! root directory for the location of drag coefficient files553 ln_cdgw = .false. ! Neutral drag coefficient read from wave model554 ln_sdw = .false. ! Computation of 3D stokes drift555 /556 !-----------------------------------------------------------------------557 &namberg ! iceberg parameters (default: No iceberg)558 !-----------------------------------------------------------------------559 ln_icebergs = .false. ! iceberg floats or not560 ln_bergdia = .true. ! Calculate budgets561 nn_verbose_level = 1 ! Turn on more verbose output if level > 0562 nn_verbose_write = 15 ! Timesteps between verbose messages563 nn_sample_rate = 1 ! Timesteps between sampling for trajectory storage564 ! Initial mass required for an iceberg of each class565 rn_initial_mass = 8.8e7, 4.1e8, 3.3e9, 1.8e10, 3.8e10, 7.5e10, 1.2e11, 2.2e11, 3.9e11, 7.4e11566 ! Proportion of calving mass to apportion to each class567 rn_distribution = 0.24, 0.12, 0.15, 0.18, 0.12, 0.07, 0.03, 0.03, 0.03, 0.02568 ! Ratio between effective and real iceberg mass (non-dim)569 ! i.e. number of icebergs represented at a point570 rn_mass_scaling = 2000, 200, 50, 20, 10, 5, 2, 1, 1, 1571 ! thickness of newly calved bergs (m)572 rn_initial_thickness = 40., 67., 133., 175., 250., 250., 250., 250., 250., 250.573 rn_rho_bergs = 850. ! Density of icebergs574 rn_LoW_ratio = 1.5 ! Initial ratio L/W for newly calved icebergs575 ln_operator_splitting = .true. ! Use first order operator splitting for thermodynamics576 rn_bits_erosion_fraction = 0. ! Fraction of erosion melt flux to divert to bergy bits577 rn_sicn_shift = 0. ! Shift of sea-ice concn in erosion flux (0<sicn_shift<1)578 ln_passive_mode = .false. ! iceberg - ocean decoupling579 nn_test_icebergs = 10 ! Create test icebergs of this class (-1 = no)580 ! Put a test iceberg at each gridpoint in box (lon1,lon2,lat1,lat2)581 rn_test_box = 108.0, 116.0, -66.0, -58.0582 rn_speed_limit = 0. ! CFL speed limit for a berg583 584 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !585 ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename !586 sn_icb = 'calving', -1 , 'calvingmask', .true. , .true. , 'yearly' , '' , '' , ''587 588 cn_dir = './'589 /590 591 173 !!====================================================================== 592 174 !! *** Lateral boundary condition *** … … 615 197 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [m2/s] 616 198 ln_chk_bathy = .FALSE. ! 617 /618 !-----------------------------------------------------------------------619 &nam_tide ! tide parameters ("key_tide")620 !-----------------------------------------------------------------------621 ln_tide_pot = .true. ! use tidal potential forcing622 ln_tide_ramp= .false. !623 rdttideramp = 0. !624 clname(1) = 'DUMMY' ! name of constituent - all tidal components must be set in namelist_cfg625 199 / 626 200 !----------------------------------------------------------------------- … … 659 233 / 660 234 !----------------------------------------------------------------------- 661 &nambdy_dta ! open boundaries - external data ("key_bdy")662 !-----------------------------------------------------------------------663 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !664 ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename !665 bn_ssh = 'amm12_bdyT_u2d', 24 , 'sossheig', .true. , .false. , 'daily' , '' , '' , ''666 bn_u2d = 'amm12_bdyU_u2d', 24 , 'vobtcrtx', .true. , .false. , 'daily' , '' , '' , ''667 bn_v2d = 'amm12_bdyV_u2d', 24 , 'vobtcrty', .true. , .false. , 'daily' , '' , '' , ''668 bn_u3d = 'amm12_bdyU_u3d', 24 , 'vozocrtx', .true. , .false. , 'daily' , '' , '' , ''669 bn_v3d = 'amm12_bdyV_u3d', 24 , 'vomecrty', .true. , .false. , 'daily' , '' , '' , ''670 bn_tem = 'amm12_bdyT_tra', 24 , 'votemper', .true. , .false. , 'daily' , '' , '' , ''671 bn_sal = 'amm12_bdyT_tra', 24 , 'vosaline', .true. , .false. , 'daily' , '' , '' , ''672 ! for lim2673 ! bn_frld = 'amm12_bdyT_ice', 24 , 'ileadfra', .true. , .false. , 'daily' , '' , '' , ''674 ! bn_hicif = 'amm12_bdyT_ice', 24 , 'iicethic', .true. , .false. , 'daily' , '' , '' , ''675 ! bn_hsnif = 'amm12_bdyT_ice', 24 , 'isnowthi', .true. , .false. , 'daily' , '' , '' , ''676 ! for lim3677 ! bn_a_i = 'amm12_bdyT_ice', 24 , 'ileadfra', .true. , .false. , 'daily' , '' , '' , ''678 ! bn_ht_i = 'amm12_bdyT_ice', 24 , 'iicethic', .true. , .false. , 'daily' , '' , '' , ''679 ! bn_ht_s = 'amm12_bdyT_ice', 24 , 'isnowthi', .true. , .false. , 'daily' , '' , '' , ''680 681 cn_dir = 'bdydta/' ! root directory for the location of the bulk files682 ln_full_vel = .false. !683 /684 !-----------------------------------------------------------------------685 &nambdy_tide ! tidal forcing at open boundaries686 !-----------------------------------------------------------------------687 filtide = 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files688 ln_bdytide_2ddta = .false. !689 ln_bdytide_conj = .false. !690 /691 692 !!======================================================================693 !! *** Bottom boundary condition ***694 !!======================================================================695 !! nambfr bottom friction696 !! nambbc bottom temperature boundary condition697 !! nambbl bottom boundary layer scheme ("key_trabbl")698 !!======================================================================699 !700 !-----------------------------------------------------------------------701 &nambfr ! bottom friction (default: linear)702 !-----------------------------------------------------------------------703 nn_bfr = 1 ! type of bottom friction : = 0 : free slip, = 1 : linear friction704 ! = 2 : nonlinear friction705 rn_bfri1 = 4.e-4 ! bottom drag coefficient (linear case)706 rn_bfri2 = 1.e-3 ! bottom drag coefficient (non linear case). Minimum coeft if ln_loglayer=T707 rn_bfri2_max= 1.e-1 ! max. bottom drag coefficient (non linear case and ln_loglayer=T)708 rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m2/s2)709 rn_bfrz0 = 3.e-3 ! bottom roughness [m] if ln_loglayer=T710 ln_bfr2d = .false. ! horizontal variation of the bottom friction coef (read a 2D mask file )711 rn_bfrien = 50. ! local multiplying factor of bfr (ln_bfr2d=T)712 rn_tfri1 = 4.e-4 ! top drag coefficient (linear case)713 rn_tfri2 = 2.5e-3 ! top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T714 rn_tfri2_max= 1.e-1 ! max. top drag coefficient (non linear case and ln_loglayer=T)715 rn_tfeb2 = 0.0 ! top turbulent kinetic energy background (m2/s2)716 rn_tfrz0 = 3.e-3 ! top roughness [m] if ln_loglayer=T717 ln_tfr2d = .false. ! horizontal variation of the top friction coef (read a 2D mask file )718 rn_tfrien = 50. ! local multiplying factor of tfr (ln_tfr2d=T)719 720 ln_bfrimp = .true. ! implicit bottom friction (requires ln_zdfexp = .false. if true)721 ln_loglayer = .false. ! logarithmic formulation (non linear case)722 /723 !-----------------------------------------------------------------------724 &nambbc ! bottom temperature boundary condition (default: NO)725 !-----------------------------------------------------------------------726 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask !727 ! ! ! (if <0 months) ! name ! (logical) ! (T/F ) ! 'monthly' ! filename ! pairing ! filename !728 sn_qgh ='geothermal_heating.nc', -12. , 'heatflow', .false. , .true. , 'yearly' , '' , '' , ''729 !730 ln_trabbc = .false. ! Apply a geothermal heating at the ocean bottom731 nn_geoflx = 2 ! geothermal heat flux: = 0 no flux732 ! = 1 constant flux733 ! = 2 variable flux (read in geothermal_heating.nc in mW/m2)734 rn_geoflx_cst = 86.4e-3 ! Constant value of geothermal heat flux [W/m2]735 cn_dir = './' ! root directory for the location of the runoff files736 /737 !-----------------------------------------------------------------------738 &nambbl ! bottom boundary layer scheme ("key_trabbl")739 !-----------------------------------------------------------------------740 nn_bbl_ldf = 1 ! diffusive bbl (=1) or not (=0)741 nn_bbl_adv = 0 ! advective bbl (=1/2) or not (=0)742 rn_ahtbbl = 1000. ! lateral mixing coefficient in the bbl [m2/s]743 rn_gambbl = 10. ! advective bbl coefficient [s]744 /745 746 !!======================================================================747 !! Tracer (T & S ) namelists748 !!======================================================================749 !! nameos equation of state750 !! namtra_adv advection scheme751 !! namtra_adv_mle mixed layer eddy param. (Fox-Kemper param.)752 !! namtra_ldf lateral diffusion scheme753 !! namtra_ldfeiv eddy induced velocity param.754 !! namtra_dmp T & S newtonian damping755 !!======================================================================756 !757 !-----------------------------------------------------------------------758 &nameos ! ocean physical parameters759 !-----------------------------------------------------------------------760 ln_teos10 = .false. ! = Use TEOS-10 equation of state761 ln_eos80 = .false. ! = Use EOS80 equation of state762 ln_seos = .false. ! = Use simplified equation of state (S-EOS)763 !764 ! ! S-EOS coefficients (ln_seos=T):765 ! ! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS766 rn_a0 = 1.6550e-1 ! thermal expension coefficient (nn_eos= 1)767 rn_b0 = 7.6554e-1 ! saline expension coefficient (nn_eos= 1)768 rn_lambda1 = 5.9520e-2 ! cabbeling coeff in T^2 (=0 for linear eos)769 rn_lambda2 = 7.4914e-4 ! cabbeling coeff in S^2 (=0 for linear eos)770 rn_mu1 = 1.4970e-4 ! thermobaric coeff. in T (=0 for linear eos)771 rn_mu2 = 1.1090e-5 ! thermobaric coeff. in S (=0 for linear eos)772 rn_nu = 2.4341e-3 ! cabbeling coeff in T*S (=0 for linear eos)773 /774 !-----------------------------------------------------------------------775 &namtra_adv ! advection scheme for tracer (default: NO advection)776 !-----------------------------------------------------------------------777 ln_traadv_cen = .false. ! 2nd order centered scheme778 nn_cen_h = 4 ! =2/4, horizontal 2nd order CEN / 4th order CEN779 nn_cen_v = 4 ! =2/4, vertical 2nd order CEN / 4th order COMPACT780 ln_traadv_fct = .false. ! FCT scheme781 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order782 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order783 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping784 ! ! (number of sub-timestep = nn_fct_zts)785 ln_traadv_mus = .false. ! MUSCL scheme786 ln_mus_ups = .false. ! use upstream scheme near river mouths787 ln_traadv_ubs = .false. ! UBS scheme788 nn_ubs_v = 2 ! =2 , vertical 2nd order FCT / COMPACT 4th order789 ln_traadv_qck = .false. ! QUICKEST scheme790 /791 !-----------------------------------------------------------------------792 &namtra_adv_mle ! mixed layer eddy parametrisation (Fox-Kemper param) (default: NO)793 !-----------------------------------------------------------------------794 ln_mle = .false. ! (T) use the Mixed Layer Eddy (MLE) parameterisation795 rn_ce = 0.06 ! magnitude of the MLE (typical value: 0.06 to 0.08)796 nn_mle = 1 ! MLE type: =0 standard Fox-Kemper ; =1 new formulation797 rn_lf = 5.e+3 ! typical scale of mixed layer front (meters) (case rn_mle=0)798 rn_time = 172800. ! time scale for mixing momentum across the mixed layer (seconds) (case rn_mle=0)799 rn_lat = 20. ! reference latitude (degrees) of MLE coef. (case rn_mle=1)800 nn_mld_uv = 0 ! space interpolation of MLD at u- & v-pts (0=min,1=averaged,2=max)801 nn_conv = 0 ! =1 no MLE in case of convection ; =0 always MLE802 rn_rho_c_mle= 0.01 ! delta rho criterion used to calculate MLD for FK803 /804 !-----------------------------------------------------------------------805 &namtra_ldf ! lateral diffusion scheme for tracers (default: NO diffusion)806 !-----------------------------------------------------------------------807 ! ! Operator type:808 ! ! no diffusion: set ln_traldf_lap=..._blp=F809 ln_traldf_lap = .false. ! laplacian operator810 ln_traldf_blp = .false. ! bilaplacian operator811 !812 ! ! Direction of action:813 ln_traldf_lev = .false. ! iso-level814 ln_traldf_hor = .false. ! horizontal (geopotential)815 ln_traldf_iso = .false. ! iso-neutral (standard operator)816 ln_traldf_triad = .false. ! iso-neutral (triad operator)817 !818 ! ! iso-neutral options:819 ln_traldf_msc = .false. ! Method of Stabilizing Correction (both operators)820 rn_slpmax = 0.01 ! slope limit (both operators)821 ln_triad_iso = .false. ! pure horizontal mixing in ML (triad only)822 rn_sw_triad = 1 ! =1 switching triad ; =0 all 4 triads used (triad only)823 ln_botmix_triad = .false. ! lateral mixing on bottom (triad only)824 !825 ! ! Coefficients:826 nn_aht_ijk_t = 0 ! space/time variation of eddy coef827 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file828 ! ! = 0 constant829 ! ! = 10 F(k) =ldf_c1d830 ! ! = 20 F(i,j) =ldf_c2d831 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation832 ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d833 ! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing)834 rn_aht_0 = 2000. ! lateral eddy diffusivity (lap. operator) [m2/s]835 rn_bht_0 = 1.e+12 ! lateral eddy diffusivity (bilap. operator) [m4/s]836 /837 !-----------------------------------------------------------------------838 &namtra_ldfeiv ! eddy induced velocity param. (default: NO)839 !-----------------------------------------------------------------------840 ln_ldfeiv =.false. ! use eddy induced velocity parameterization841 ln_ldfeiv_dia =.false. ! diagnose eiv stream function and velocities842 rn_aeiv_0 = 2000. ! eddy induced velocity coefficient [m2/s]843 nn_aei_ijk_t = 21 ! space/time variation of the eiv coeficient844 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file845 ! ! = 0 constant846 ! ! = 10 F(k) =ldf_c1d847 ! ! = 20 F(i,j) =ldf_c2d848 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation849 ! ! = 30 F(i,j,k) =ldf_c2d + ldf_c1d850 /851 !-----------------------------------------------------------------------852 &namtra_dmp ! tracer: T & S newtonian damping (default: NO)853 !-----------------------------------------------------------------------854 ln_tradmp = .true. ! add a damping termn (T) or not (F)855 nn_zdmp = 0 ! vertical shape =0 damping throughout the water column856 ! =1 no damping in the mixing layer (kz criteria)857 ! =2 no damping in the mixed layer (rho crieria)858 cn_resto ='resto.nc' ! Name of file containing restoration coeff. field (use dmp_tools to create this)859 /860 861 !!======================================================================862 !! *** Dynamics namelists ***863 !!======================================================================864 !! namdyn_adv formulation of the momentum advection865 !! namdyn_vor advection scheme866 !! namdyn_hpg hydrostatic pressure gradient867 !! namdyn_spg surface pressure gradient868 !! namdyn_ldf lateral diffusion scheme869 !!======================================================================870 !871 !-----------------------------------------------------------------------872 &namdyn_adv ! formulation of the momentum advection (default: vector form)873 !-----------------------------------------------------------------------874 ln_dynadv_vec = .true. ! vector form (T) or flux form (F)875 nn_dynkeg = 0 ! scheme for grad(KE): =0 C2 ; =1 Hollingsworth correction876 ln_dynadv_cen2= .false. ! flux form - 2nd order centered scheme877 ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme878 ln_dynzad_zts = .false. ! Use (T) sub timestepping for vertical momentum advection879 /880 !-----------------------------------------------------------------------881 235 &nam_vvl ! vertical coordinate options (default: zstar) 882 236 !----------------------------------------------------------------------- … … 892 246 ln_vvl_dbg = .true. ! debug prints (T/F) 893 247 / 894 !-----------------------------------------------------------------------895 &namdyn_vor ! Vorticity / Coriolis scheme (default: NO)896 !-----------------------------------------------------------------------897 ln_dynvor_ene = .false. ! enstrophy conserving scheme898 ln_dynvor_ens = .false. ! energy conserving scheme899 ln_dynvor_mix = .false. ! mixed scheme900 ln_dynvor_een = .false. ! energy & enstrophy scheme901 nn_een_e3f = 1 ! e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1)902 ln_dynvor_msk = .false. ! vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) ! PLEASE DO NOT ACTIVATE903 /904 !-----------------------------------------------------------------------905 &namdyn_hpg ! Hydrostatic pressure gradient option (default: zps)906 !-----------------------------------------------------------------------907 ln_hpg_zco = .false. ! z-coordinate - full steps908 ln_hpg_zps = .true. ! z-coordinate - partial steps (interpolation)909 ln_hpg_sco = .false. ! s-coordinate (standard jacobian formulation)910 ln_hpg_isf = .false. ! s-coordinate (sco ) adapted to isf911 ln_hpg_djc = .false. ! s-coordinate (Density Jacobian with Cubic polynomial)912 ln_hpg_prj = .false. ! s-coordinate (Pressure Jacobian scheme)913 /914 !-----------------------------------------------------------------------915 &namdyn_spg ! surface pressure gradient (default: NO)916 !-----------------------------------------------------------------------917 ln_dynspg_exp = .false. ! explicit free surface918 ln_dynspg_ts = .false. ! split-explicit free surface919 ln_bt_fw = .true. ! Forward integration of barotropic Eqs.920 ln_bt_av = .true. ! Time filtering of barotropic variables921 nn_bt_flt = 1 ! Time filter choice = 0 None922 ! ! = 1 Boxcar over nn_baro sub-steps923 ! ! = 2 Boxcar over 2*nn_baro " "924 ln_bt_auto = .true. ! Number of sub-step defined from:925 rn_bt_cmax = 0.8 ! =T : the Maximum Courant Number allowed926 nn_baro = 30 ! =F : the number of sub-step in rn_rdt seconds927 /928 !-----------------------------------------------------------------------929 &namdyn_ldf ! lateral diffusion on momentum (default: NO)930 !-----------------------------------------------------------------------931 ! ! Type of the operator :932 ! ! no diffusion: set ln_dynldf_lap=..._blp=F933 ln_dynldf_lap = .false. ! laplacian operator934 ln_dynldf_blp = .false. ! bilaplacian operator935 ! ! Direction of action :936 ln_dynldf_lev = .false. ! iso-level937 ln_dynldf_hor = .false. ! horizontal (geopotential)938 ln_dynldf_iso = .false. ! iso-neutral939 ! ! Coefficient940 nn_ahm_ijk_t = 0 ! space/time variation of eddy coef941 ! ! =-30 read in eddy_viscosity_3D.nc file942 ! ! =-20 read in eddy_viscosity_2D.nc file943 ! ! = 0 constant944 ! ! = 10 F(k)=c1d945 ! ! = 20 F(i,j)=F(grid spacing)=c2d946 ! ! = 30 F(i,j,k)=c2d*c1d947 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity)948 rn_ahm_0 = 40000. ! horizontal laplacian eddy viscosity [m2/s]949 rn_ahm_b = 0. ! background eddy viscosity for ldf_iso [m2/s]950 rn_bhm_0 = 1.e+12 ! horizontal bilaplacian eddy viscosity [m4/s]951 !952 ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km)953 /954 955 !!======================================================================956 !! Tracers & Dynamics vertical physics namelists957 !!======================================================================958 !! namzdf vertical physics959 !! namzdf_ric richardson number dependent vertical mixing ("key_zdfric")960 !! namzdf_tke TKE dependent vertical mixing ("key_zdftke")961 !! namzdf_gls GLS vertical mixing ("key_zdfgls")962 !! namzdf_ddm double diffusive mixing parameterization ("key_zdfddm")963 !! namzdf_tmx tidal mixing parameterization ("key_zdftmx")964 !!======================================================================965 !966 !-----------------------------------------------------------------------967 &namzdf ! vertical physics968 !-----------------------------------------------------------------------969 rn_avm0 = 1.2e-4 ! vertical eddy viscosity [m2/s] (background Kz if not "key_zdfcst")970 rn_avt0 = 1.2e-5 ! vertical eddy diffusivity [m2/s] (background Kz if not "key_zdfcst")971 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0)972 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0)973 ln_zdfevd = .true. ! enhanced vertical diffusion (evd) (T) or not (F)974 nn_evdm = 0 ! evd apply on tracer (=0) or on tracer and momentum (=1)975 rn_avevd = 100. ! evd mixing coefficient [m2/s]976 ln_zdfnpc = .false. ! Non-Penetrative Convective algorithm (T) or not (F)977 nn_npc = 1 ! frequency of application of npc978 nn_npcp = 365 ! npc control print frequency979 ln_zdfexp = .false. ! time-stepping: split-explicit (T) or implicit (F) time stepping980 nn_zdfexp = 3 ! number of sub-timestep for ln_zdfexp=T981 /982 !-----------------------------------------------------------------------983 &namzdf_ric ! richardson number dependent vertical diffusion ("key_zdfric" )984 !-----------------------------------------------------------------------985 rn_avmri = 100.e-4 ! maximum value of the vertical viscosity986 rn_alp = 5. ! coefficient of the parameterization987 nn_ric = 2 ! coefficient of the parameterization988 rn_ekmfc = 0.7 ! Factor in the Ekman depth Equation989 rn_mldmin = 1.0 ! minimum allowable mixed-layer depth estimate (m)990 rn_mldmax = 1000.0 ! maximum allowable mixed-layer depth estimate (m)991 rn_wtmix = 10.0 ! vertical eddy viscosity coeff [m2/s] in the mixed-layer992 rn_wvmix = 10.0 ! vertical eddy diffusion coeff [m2/s] in the mixed-layer993 ln_mldw = .true. ! Flag to use or not the mixed layer depth param.994 /995 !-----------------------------------------------------------------------996 &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke")997 !-----------------------------------------------------------------------998 rn_ediff = 0.1 ! coef. for vertical eddy coef. (avt=rn_ediff*mxl*sqrt(e) )999 rn_ediss = 0.7 ! coef. of the Kolmogoroff dissipation1000 rn_ebb = 67.83 ! coef. of the surface input of tke (=67.83 suggested when ln_mxl0=T)1001 rn_emin = 1.e-6 ! minimum value of tke [m2/s2]1002 rn_emin0 = 1.e-4 ! surface minimum value of tke [m2/s2]1003 rn_bshear = 1.e-20 ! background shear (>0) currently a numerical threshold (do not change it)1004 nn_mxl = 2 ! mixing length: = 0 bounded by the distance to surface and bottom1005 ! = 1 bounded by the local vertical scale factor1006 ! = 2 first vertical derivative of mixing length bounded by 11007 ! = 3 as =2 with distinct disspipative an mixing length scale1008 nn_pdl = 1 ! Prandtl number function of richarson number (=1, avt=pdl(Ri)*avm) or not (=0, avt=avm)1009 ln_mxl0 = .true. ! surface mixing length scale = F(wind stress) (T) or not (F)1010 rn_mxl0 = 0.04 ! surface buoyancy lenght scale minimum value1011 ln_lc = .true. ! Langmuir cell parameterisation (Axell 2002)1012 rn_lc = 0.15 ! coef. associated to Langmuir cells1013 nn_etau = 1 ! penetration of tke below the mixed layer (ML) due to near intertial waves1014 ! = 0 no penetration1015 ! = 1 add a tke source below the ML1016 ! = 2 add a tke source just at the base of the ML1017 ! = 3 as = 1 applied on HF part of the stress (ln_cpl=T)1018 rn_efr = 0.05 ! fraction of surface tke value which penetrates below the ML (nn_etau=1 or 2)1019 nn_htau = 1 ! type of exponential decrease of tke penetration below the ML1020 ! = 0 constant 10 m length scale1021 ! = 1 0.5m at the equator to 30m poleward of 40 degrees1022 /1023 !-----------------------------------------------------------------------1024 &namzdf_gls ! GLS vertical diffusion ("key_zdfgls")1025 !-----------------------------------------------------------------------1026 rn_emin = 1.e-7 ! minimum value of e [m2/s2]1027 rn_epsmin = 1.e-12 ! minimum value of eps [m2/s3]1028 ln_length_lim = .true. ! limit on the dissipation rate under stable stratification (Galperin et al., 1988)1029 rn_clim_galp = 0.267 ! galperin limit1030 ln_sigpsi = .true. ! Activate or not Burchard 2001 mods on psi schmidt number in the wb case1031 rn_crban = 100. ! Craig and Banner 1994 constant for wb tke flux1032 rn_charn = 70000. ! Charnock constant for wb induced roughness length1033 rn_hsro = 0.02 ! Minimum surface roughness1034 rn_frac_hs = 1.3 ! Fraction of wave height as roughness (if nn_z0_met=2)1035 nn_z0_met = 2 ! Method for surface roughness computation (0/1/2)1036 nn_bc_surf = 1 ! surface condition (0/1=Dir/Neum)1037 nn_bc_bot = 1 ! bottom condition (0/1=Dir/Neum)1038 nn_stab_func = 2 ! stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB)1039 nn_clos = 1 ! predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)1040 /1041 !-----------------------------------------------------------------------1042 &namzdf_ddm ! double diffusive mixing parameterization ("key_zdfddm")1043 !-----------------------------------------------------------------------1044 rn_avts = 1.e-4 ! maximum avs (vertical mixing on salinity)1045 rn_hsbfr = 1.6 ! heat/salt buoyancy flux ratio1046 /1047 !-----------------------------------------------------------------------1048 &namzdf_tmx ! tidal mixing parameterization ("key_zdftmx")1049 !-----------------------------------------------------------------------1050 rn_htmx = 500. ! vertical decay scale for turbulence (meters)1051 rn_n2min = 1.e-8 ! threshold of the Brunt-Vaisala frequency (s-1)1052 rn_tfe = 0.333 ! tidal dissipation efficiency1053 rn_me = 0.2 ! mixing efficiency1054 ln_tmx_itf = .true. ! ITF specific parameterisation1055 rn_tfe_itf = 1. ! ITF tidal dissipation efficiency1056 /1057 !-----------------------------------------------------------------------1058 &namzdf_tmx_new ! internal wave-driven mixing parameterization ("key_zdftmx_new" & "key_zdfddm")1059 !-----------------------------------------------------------------------1060 nn_zpyc = 1 ! pycnocline-intensified dissipation scales as N (=1) or N^2 (=2)1061 ln_mevar = .true. ! variable (T) or constant (F) mixing efficiency1062 ln_tsdiff = .true. ! account for differential T/S mixing (T) or not (F)1063 /1064 1065 1066 248 !!====================================================================== 1067 249 !! *** Miscellaneous namelists *** … … 1086 268 &namctl ! Control prints & Benchmark 1087 269 !----------------------------------------------------------------------- 1088 ln_ctl = .false. ! trends control print (expensive!) 270 ln_ctl = .FALSE. ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T 271 sn_cfctl%l_config = .TRUE. ! IF .true. then control which reports are written with the following 272 sn_cfctl%l_runstat = .FALSE. ! switches and which areas produce reports with the proc integer settings. 273 sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 274 sn_cfctl%l_oceout = .FALSE. ! that all areas report. 275 sn_cfctl%l_layout = .FALSE. ! 276 sn_cfctl%l_mppout = .FALSE. ! 277 sn_cfctl%l_mpptop = .FALSE. ! 278 sn_cfctl%procmin = 0 ! Minimum area number for reporting [default:0] 279 sn_cfctl%procmax = 1000000 ! Maximum area number for reporting [default:1000000] 280 sn_cfctl%procincr = 1 ! Increment for optional subsetting of areas [default:1] 281 sn_cfctl%ptimincr = 1 ! Timestep increment for writing time step progress info 1089 282 nn_print = 0 ! level of print (0 no extra print) 1090 283 nn_ictls = 0 ! start i indice of control sum (use to compare mono versus … … 1094 287 nn_isplt = 1 ! number of processors in i-direction 1095 288 nn_jsplt = 1 ! number of processors in j-direction 1096 nn_bench = 0 ! Bench mode (1/0): CAUTION use zero except for bench 1097 ! (no physical validity of the results) 1098 nn_timing = 0 ! timing by routine activated (=1) creates timing.output file, or not (=0) 1099 nn_diacfl = 0 ! Write out CFL diagnostics (=1) in cfl_diagnostics.ascii, or not (=0) 1100 / 1101 !----------------------------------------------------------------------- 1102 &namsto ! Stochastic parametrization of EOS (default: NO) 1103 !----------------------------------------------------------------------- 1104 ln_sto_eos = .false. ! stochastic equation of state 1105 nn_sto_eos = 1 ! number of independent random walks 1106 rn_eos_stdxy= 1.4 ! random walk horz. standard deviation (in grid points) 1107 rn_eos_stdz = 0.7 ! random walk vert. standard deviation (in grid points) 1108 rn_eos_tcor = 1440. ! random walk time correlation (in timesteps) 1109 nn_eos_ord = 1 ! order of autoregressive processes 1110 nn_eos_flt = 0 ! passes of Laplacian filter 1111 rn_eos_lim = 2.0 ! limitation factor (default = 3.0) 1112 ln_rststo = .false. ! start from mean parameter (F) or from restart file (T) 1113 ln_rstseed = .true. ! read seed of RNG from restart file 1114 cn_storst_in = "restart_sto" ! suffix of stochastic parameter restart file (input) 1115 cn_storst_out = "restart_sto" ! suffix of stochastic parameter restart file (output) 1116 / 1117 1118 !!====================================================================== 1119 !! *** Diagnostics namelists *** 1120 !!====================================================================== 1121 !! namtrd dynamics and/or tracer trends (default F) 1122 !! namptr Poleward Transport Diagnostics (default F) 1123 !! namhsb Heat and salt budgets (default F) 1124 !! namdiu Cool skin and warm layer models (default F) 1125 !! namflo float parameters ("key_float") 1126 !! nam_diaharm Harmonic analysis of tidal constituents ("key_diaharm") 1127 !! namdct transports through some sections ("key_diadct") 1128 !! nam_diatmb Top Middle Bottom Output (default F) 1129 !! nam_dia25h 25h Mean Output (default F) 1130 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") 1131 !!====================================================================== 1132 ! 1133 !----------------------------------------------------------------------- 1134 &namtrd ! trend diagnostics (default F) 1135 !----------------------------------------------------------------------- 1136 ln_glo_trd = .false. ! (T) global domain averaged diag for T, T^2, KE, and PE 1137 ln_dyn_trd = .false. ! (T) 3D momentum trend output 1138 ln_dyn_mxl = .false. ! (T) 2D momentum trends averaged over the mixed layer (not coded yet) 1139 ln_vor_trd = .false. ! (T) 2D barotropic vorticity trends (not coded yet) 1140 ln_KE_trd = .false. ! (T) 3D Kinetic Energy trends 1141 ln_PE_trd = .false. ! (T) 3D Potential Energy trends 1142 ln_tra_trd = .false. ! (T) 3D tracer trend output 1143 ln_tra_mxl = .false. ! (T) 2D tracer trends averaged over the mixed layer (not coded yet) 1144 nn_trd = 365 ! print frequency (ln_glo_trd=T) (unit=time step) 1145 / 1146 !!gm nn_ctls = 0 ! control surface type in mixed-layer trends (0,1 or n<jpk) 1147 !!gm rn_ucf = 1. ! unit conversion factor (=1 -> /seconds ; =86400. -> /day) 1148 !!gm cn_trdrst_in = "restart_mld" ! suffix of ocean restart name (input) 1149 !!gm cn_trdrst_out = "restart_mld" ! suffix of ocean restart name (output) 1150 !!gm ln_trdmld_restart = .false. ! restart for ML diagnostics 1151 !!gm ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S 1152 !!gm 1153 !----------------------------------------------------------------------- 1154 &namptr ! Poleward Transport Diagnostic (default F) 1155 !----------------------------------------------------------------------- 1156 ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F) 1157 ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not 1158 / 1159 !----------------------------------------------------------------------- 1160 &namhsb ! Heat and salt budgets (default F) 1161 !----------------------------------------------------------------------- 1162 ln_diahsb = .false. ! check the heat and salt budgets (T) or not (F) 1163 / 1164 !----------------------------------------------------------------------- 1165 &namdiu ! Cool skin and warm layer models (default F) 1166 !----------------------------------------------------------------------- 1167 ln_diurnal = .false. ! 1168 ln_diurnal_only = .false. ! 1169 / 1170 !----------------------------------------------------------------------- 1171 &namflo ! float parameters ("key_float") 1172 !----------------------------------------------------------------------- 1173 jpnfl = 1 ! total number of floats during the run 1174 jpnnewflo = 0 ! number of floats for the restart 1175 ln_rstflo = .false. ! float restart (T) or not (F) 1176 nn_writefl = 75 ! frequency of writing in float output file 1177 nn_stockfl = 5475 ! frequency of creation of the float restart file 1178 ln_argo = .false. ! Argo type floats (stay at the surface each 10 days) 1179 ln_flork4 = .false. ! trajectories computed with a 4th order Runge-Kutta (T) 1180 ! ! or computed with Blanke' scheme (F) 1181 ln_ariane = .true. ! Input with Ariane tool convention(T) 1182 ln_flo_ascii= .true. ! Output with Ariane tool netcdf convention(F) or ascii file (T) 1183 / 1184 !----------------------------------------------------------------------- 1185 &nam_diaharm ! Harmonic analysis of tidal constituents ("key_diaharm") 1186 !----------------------------------------------------------------------- 1187 nit000_han = 1 ! First time step used for harmonic analysis 1188 nitend_han = 75 ! Last time step used for harmonic analysis 1189 nstep_han = 15 ! Time step frequency for harmonic analysis 1190 tname(1) = 'M2' ! Name of tidal constituents 1191 tname(2) = 'K1' 1192 / 1193 !----------------------------------------------------------------------- 1194 &namdct ! transports through some sections ("key_diadct") 1195 !----------------------------------------------------------------------- 1196 nn_dct = 15 ! time step frequency for transports computing 1197 nn_dctwri = 15 ! time step frequency for transports writing 1198 nn_secdebug= 112 ! 0 : no section to debug 1199 ! ! -1 : debug all section 1200 ! ! 0 < n : debug section number n 1201 / 1202 !----------------------------------------------------------------------- 1203 &nam_diatmb ! Top Middle Bottom Output (default F) 1204 !----------------------------------------------------------------------- 1205 ln_diatmb = .false. ! Choose Top Middle and Bottom output or not 1206 / 1207 !----------------------------------------------------------------------- 1208 &nam_dia25h ! 25h Mean Output (default F) 1209 !----------------------------------------------------------------------- 1210 ln_dia25h = .false. ! Choose 25h mean output or not 1211 / 1212 !----------------------------------------------------------------------- 1213 &namnc4 ! netcdf4 chunking and compression settings ("key_netcdf4") 1214 !----------------------------------------------------------------------- 1215 nn_nchunks_i= 4 ! number of chunks in i-dimension 1216 nn_nchunks_j= 4 ! number of chunks in j-dimension 1217 nn_nchunks_k= 31 ! number of chunks in k-dimension 1218 ! ! setting nn_nchunks_k = jpk will give a chunk size of 1 in the vertical which 1219 ! ! is optimal for postprocessing which works exclusively with horizontal slabs 1220 ln_nc4zip = .true. ! (T) use netcdf4 chunking and compression 1221 ! ! (F) ignore chunking information and produce netcdf3-compatible files 1222 / 1223 1224 !!====================================================================== 1225 !! *** Observation & Assimilation *** 1226 !!====================================================================== 1227 !! namobs observation and model comparison 1228 !! nam_asminc assimilation increments ('key_asminc') 1229 !!====================================================================== 1230 ! 1231 !----------------------------------------------------------------------- 1232 &namobs ! observation usage switch 1233 !----------------------------------------------------------------------- 1234 ln_diaobs = .false. ! Logical switch for the observation operator 1235 ln_t3d = .false. ! Logical switch for T profile observations 1236 ln_s3d = .false. ! Logical switch for S profile observations 1237 ln_sla = .false. ! Logical switch for SLA observations 1238 ln_sst = .false. ! Logical switch for SST observations 1239 ln_sic = .false. ! Logical switch for Sea Ice observations 1240 ln_vel3d = .false. ! Logical switch for velocity observations 1241 ln_altbias = .false. ! Logical switch for altimeter bias correction 1242 ln_nea = .false. ! Logical switch for rejection of observations near land 1243 ln_grid_global = .true. ! Logical switch for global distribution of observations 1244 ln_grid_search_lookup = .false. ! Logical switch for obs grid search w/lookup table 1245 ln_ignmis = .true. ! Logical switch for ignoring missing files 1246 ln_s_at_t = .false. ! Logical switch for computing model S at T obs if not there 1247 ln_sstnight = .false. ! Logical switch for calculating night-time average for SST obs 1248 ! All of the *files* variables below are arrays. Use namelist_cfg to add more files 1249 cn_profbfiles = 'profiles_01.nc' ! Profile feedback input observation file names 1250 cn_slafbfiles = 'sla_01.nc' ! SLA feedback input observation file names 1251 cn_sstfbfiles = 'sst_01.nc' ! SST feedback input observation file names 1252 cn_sicfbfiles = 'sic_01.nc' ! SIC feedback input observation file names 1253 cn_velfbfiles = 'vel_01.nc' ! Velocity feedback input observation file names 1254 cn_altbiasfile = 'altbias.nc' ! Altimeter bias input file name 1255 cn_gridsearchfile='gridsearch.nc' ! Grid search file name 1256 rn_gridsearchres = 0.5 ! Grid search resolution 1257 rn_dobsini = 00010101.000000 ! Initial date in window YYYYMMDD.HHMMSS 1258 rn_dobsend = 00010102.000000 ! Final date in window YYYYMMDD.HHMMSS 1259 nn_1dint = 0 ! Type of vertical interpolation method 1260 nn_2dint = 0 ! Type of horizontal interpolation method 1261 nn_msshc = 0 ! MSSH correction scheme 1262 rn_mdtcorr = 1.61 ! MDT correction 1263 rn_mdtcutoff = 65.0 ! MDT cutoff for computed correction 1264 nn_profdavtypes = -1 ! Profile daily average types - array 1265 ln_sstbias = .false. ! 1266 cn_sstbias_files = 'sstbias.nc' ! 1267 / 1268 !----------------------------------------------------------------------- 1269 &nam_asminc ! assimilation increments ('key_asminc') 1270 !----------------------------------------------------------------------- 1271 ln_bkgwri = .false. ! Logical switch for writing out background state 1272 ln_trainc = .false. ! Logical switch for applying tracer increments 1273 ln_dyninc = .false. ! Logical switch for applying velocity increments 1274 ln_sshinc = .false. ! Logical switch for applying SSH increments 1275 ln_asmdin = .false. ! Logical switch for Direct Initialization (DI) 1276 ln_asmiau = .false. ! Logical switch for Incremental Analysis Updating (IAU) 1277 nitbkg = 0 ! Timestep of background in [0,nitend-nit000-1] 1278 nitdin = 0 ! Timestep of background for DI in [0,nitend-nit000-1] 1279 nitiaustr = 1 ! Timestep of start of IAU interval in [0,nitend-nit000-1] 1280 nitiaufin = 15 ! Timestep of end of IAU interval in [0,nitend-nit000-1] 1281 niaufn = 0 ! Type of IAU weighting function 1282 ln_salfix = .false. ! Logical switch for ensuring that the sa > salfixmin 1283 salfixmin = -9999 ! Minimum salinity after applying the increments 1284 nn_divdmp = 0 ! Number of iterations of divergence damping operator 1285 / 289 ln_timing = .false. ! timing by routine write out in timing.output file 290 ln_diacfl = .false. ! CFL diagnostics write out in cfl_diagnostics.ascii 291 /
Note: See TracChangeset
for help on using the changeset viewer.