- Timestamp:
- 2020-05-14T21:46:00+02:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser
- Property svn:externals
-
old new 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL 8 9 # SETTE 10 ^/utils/CI/sette@HEAD sette
-
- Property svn:externals
-
NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/cfgs/SHARED/namelist_ref
r12314 r12928 5 5 !! namelists 2 - Surface boundary (namsbc, namsbc_flx, namsbc_blk, namsbc_cpl, 6 6 !! namsbc_sas, namtra_qsr, namsbc_rnf, 7 !! nam sbc_isf, namsbc_iscpl, namsbc_apr,7 !! namisf, namsbc_apr, 8 8 !! namsbc_ssr, namsbc_wave, namberg) 9 9 !! 3 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) … … 42 42 nn_leapy = 0 ! Leap year calendar (1) or not (0) 43 43 ln_rstart = .false. ! start from rest (F) or from a restart file (T) 44 nn_euler = 1 ! = 0 : start with forward time step if ln_rstart=T45 nn_rstctl = 0! restart control ==> activated only if ln_rstart=T44 ln_1st_euler = .false. ! =T force a start with forward time step (ln_rstart=T) 45 nn_rstctl = 0 ! restart control ==> activated only if ln_rstart=T 46 46 ! ! = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 47 47 ! ! = 1 nn_date0 read in namelist ; nn_it000 : check consistancy between namelist and restart … … 51 51 cn_ocerst_out = "restart" ! suffix of ocean restart name (output) 52 52 cn_ocerst_outdir = "." ! directory in which to write output ocean restarts 53 ln_iscpl = .false. ! cavity evolution forcing or coupling to ice sheet model54 53 nn_istate = 0 ! output the initial state (1) or not (0) 55 54 ln_rst_list = .false. ! output restarts at list of times using nn_stocklist (T) or at set frequency with nn_stock (F) … … 65 64 ln_clobber = .true. ! clobber (overwrite) an existing file 66 65 nn_chunksz = 0 ! chunksize (bytes) for NetCDF file (works only with iom_nf90 routines) 67 ln_xios_read = . FALSE. ! use XIOS to read restart file (only for a single file restart)66 ln_xios_read = .false. ! use XIOS to read restart file (only for a single file restart) 68 67 nn_wxios = 0 ! use XIOS to write restart file 0 - no, 1 - single file output, 2 - multiple file output 69 68 / … … 72 71 !----------------------------------------------------------------------- 73 72 ln_linssh = .false. ! =T linear free surface ==>> model level are fixed in time 74 rn_isfhmin = 1.00 ! treshold [m] to discriminate grounding ice from floating ice 75 ! 76 rn_rdt = 5400. ! time step for the dynamics and tracer 73 ! 74 rn_Dt = 5400. ! time step for the dynamics and tracer 77 75 rn_atfp = 0.1 ! asselin time filter parameter 78 76 ! 79 77 ln_crs = .false. ! Logical switch for coarsening module (T => fill namcrs) 80 78 ! 81 ln_meshmask = . false. ! =T create a mesh file79 ln_meshmask = .true. ! =T create a mesh file 82 80 / 83 81 !----------------------------------------------------------------------- 84 82 &namcfg ! parameters of the configuration (default: use namusr_def in namelist_cfg) 85 83 !----------------------------------------------------------------------- 86 ln_read_cfg = .false. ! (=T) read the domain configuration file87 ! ! (=F) user defined configuration (F => create/check namusr_def)84 ln_read_cfg = .false. ! (=T) read the domain configuration file 85 ! ! (=F) user defined configuration (F => create/check namusr_def) 88 86 cn_domcfg = "domain_cfg" ! domain configuration filename 89 87 ! 90 ln_closea = .false. ! T => keep closed seas (defined by closea_mask field) in the 91 ! ! domain and apply special treatment of freshwater fluxes. 92 ! ! F => suppress closed seas (defined by closea_mask field) 93 ! ! from the bathymetry at runtime. 94 ! ! If closea_mask field doesn't exist in the domain_cfg file 95 ! ! then this logical does nothing. 96 ln_write_cfg = .false. ! (=T) create the domain configuration file 88 ln_closea = .false. ! (=T => fill namclo) 89 ! ! (=F) no control of net precip/evap over closed sea 90 ! 91 ln_write_cfg = .false. ! (=T) create the domain configuration file 97 92 cn_domcfg_out = "domain_cfg_out" ! newly created domain configuration filename 98 93 ! 99 ln_use_jattr = .false. ! use (T) the file attribute: open_ocean_jstart, if present 100 ! ! in netcdf input files, as the start j-row for reading 94 ln_use_jattr = .false. ! use (T) the file attribute: open_ocean_jstart, if present 95 ! ! in netcdf input files, as the start j-row for reading 96 / 97 !----------------------------------------------------------------------- 98 &namclo ! parameters of the closed sea (cs) behavior (default: OFF) 99 !----------------------------------------------------------------------- 100 ln_maskcs = .false. ! (=T) cs are masked ; So, in this case ln_mask_csundef and ln_clo_rnf have no effect. 101 ! ! (=F => set ln_mask_csundef and ln_clo_rnf) 102 ! ! cs masks are read and net evap/precip over closed sea spread out depending on domain_cfg.nc masks. 103 ! ! See ln_mask_csundef and ln_clo_rnf for specific option related to this case 104 ! 105 ln_mask_csundef = .true. ! (=T) undefined closed seas are masked ; 106 ! ! (=F) undefined closed seas are kept and no specific treatment is done for these closed seas 107 ! 108 ln_clo_rnf = .true. ! (=T) river mouth specified in domain_cfg.nc masks (rnf and emp case) are added to the runoff mask. 109 ! ! allow the treatment of closed sea outflow grid-points to be the same as river mouth grid-points 101 110 / 102 111 !----------------------------------------------------------------------- … … 106 115 ln_tsd_init = .false. ! ocean initialisation 107 116 ln_tsd_dmp = .false. ! T-S restoring (see namtra_dmp) 108 117 109 118 cn_dir = './' ! root directory for the T-S data location 110 119 !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! … … 184 193 !! namsbc_rnf river runoffs (ln_rnf =T) 185 194 !! namsbc_apr Atmospheric Pressure (ln_apr_dyn =T) 186 !! namsbc_isf ice shelf melting/freezing (ln_isfcav =T : read (ln_read_cfg=T) or set or usr_def_zgr )187 !! namsbc_iscpl coupling option between land ice model and ocean (ln_isfcav =T)188 195 !! namsbc_wave external fields from wave model (ln_wave =T) 189 196 !! namberg iceberg floats (ln_icebergs=T) … … 195 202 nn_fsbc = 2 ! frequency of SBC module call 196 203 ! ! (control sea-ice & iceberg model call) 197 ! Type of air-sea fluxes 204 ! Type of air-sea fluxes 198 205 ln_usr = .false. ! user defined formulation (T => check usrdef_sbc) 199 206 ln_flx = .false. ! flux formulation (T => fill namsbc_flx ) 200 207 ln_blk = .false. ! Bulk formulation (T => fill namsbc_blk ) 208 ln_abl = .false. ! ABL formulation (T => fill namsbc_abl ) 201 209 ! ! Type of coupling (Ocean/Ice/Atmosphere) : 202 210 ln_cpl = .false. ! atmosphere coupled formulation ( requires key_oasis3 ) … … 205 213 ! ! =0 no opa-sas OASIS coupling: default single executable config. 206 214 ! ! =1 opa-sas OASIS coupling: multi executable config., OPA component 207 ! ! =2 opa-sas OASIS coupling: multi executable config., SAS component 215 ! ! =2 opa-sas OASIS coupling: multi executable config., SAS component 208 216 ! Sea-ice : 209 nn_ice = 0 ! =0 no ice boundary condition 217 nn_ice = 0 ! =0 no ice boundary condition 210 218 ! ! =1 use observed ice-cover ( => fill namsbc_iif ) 211 ! ! =2 or 3 automatically for SI3 or CICE ("key_si3" or "key_cice") 212 ! ! except in AGRIF zoom where it has to be specified 219 ! ! =2 or 3 for SI3 and CICE, respectively 213 220 ln_ice_embd = .false. ! =T embedded sea-ice (pressure + mass and salt exchanges) 214 221 ! ! =F levitating ice (no pressure, mass and salt exchanges) 215 ! Misc. options of sbc : 222 ! Misc. options of sbc : 216 223 ln_traqsr = .false. ! Light penetration in the ocean (T => fill namtra_qsr) 217 224 ln_dm2dc = .false. ! daily mean to diurnal cycle on short wave … … 222 229 ln_rnf = .false. ! runoffs (T => fill namsbc_rnf) 223 230 ln_apr_dyn = .false. ! Patm gradient added in ocean & ice Eqs. (T => fill namsbc_apr ) 224 ln_isf = .false. ! ice shelf (T => fill namsbc_isf & namsbc_iscpl)225 231 ln_wave = .false. ! Activate coupling with wave (T => fill namsbc_wave) 226 232 ln_cdgw = .false. ! Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave) 227 ln_sdw = .false. ! Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave) 233 ln_sdw = .false. ! Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave) 228 234 nn_sdrift = 0 ! Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift 229 235 ! ! = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)] … … 250 256 / 251 257 !----------------------------------------------------------------------- 252 &namsbc_blk ! namsbc_blk generic Bulk formula 258 &namsbc_blk ! namsbc_blk generic Bulk formula (ln_blk =T) 253 259 !----------------------------------------------------------------------- 254 260 ! ! bulk algorithm : 255 ln_NCAR = .false.! "NCAR" algorithm (Large and Yeager 2008)261 ln_NCAR = .true. ! "NCAR" algorithm (Large and Yeager 2008) 256 262 ln_COARE_3p0 = .false. ! "COARE 3.0" algorithm (Fairall et al. 2003) 257 ln_COARE_3p 5 = .false. ! "COARE 3.5" algorithm (Edson et al. 2013)258 ln_ECMWF = .false. ! "ECMWF" algorithm (IFS cycle 31)263 ln_COARE_3p6 = .false. ! "COARE 3.6" algorithm (Edson et al. 2013) 264 ln_ECMWF = .false. ! "ECMWF" algorithm (IFS cycle 45r1) 259 265 260 266 ln_humi_dpt = .false. ! Supply dewpoint tempearture instead of specific humidity (is true for ERA5) 261 267 262 268 ! 263 rn_zqt = 10. ! Air temperature & humidity reference height (m) 264 rn_zu = 10. ! Wind vector reference height (m) 265 ln_Cd_L12 = .false. ! air-ice drags = F(ice concentration) (Lupkes et al. 2012) 266 ln_Cd_L15 = .false. ! air-ice drags = F(ice concentration) (Lupkes et al. 2015) 267 ln_taudif = .false. ! HF tau contribution: use "mean of stress module - module of the mean stress" data 268 rn_pfac = 1. ! multiplicative factor for precipitation (total & snow) 269 rn_efac = 1. ! multiplicative factor for evaporation (0. or 1.) 270 rn_vfac = 0. ! multiplicative factor for ocean & ice velocity used to 271 ! ! calculate the wind stress (0.=absolute or 1.=relative winds) 272 269 rn_zqt = 10. ! Air temperature & humidity reference height (m) 270 rn_zu = 10. ! Wind vector reference height (m) 271 ln_Cd_L12 = .false. ! air-ice drags = F(ice conc.) (Lupkes et al. 2012) 272 ln_Cd_L15 = .false. ! air-ice drags = F(ice conc.) (Lupkes et al. 2015) 273 ! ! - module of the mean stress" data 274 rn_pfac = 1. ! multipl. factor for precipitation (total & snow) 275 rn_efac = 1. ! multipl. factor for evaporation (0. or 1.) 276 rn_vfac = 0. ! multipl. factor for ocean & ice velocity 277 ! ! used to calculate the wind stress 278 ! ! (0. => absolute or 1. => relative winds) 279 ln_skin_cs = .false. ! use the cool-skin parameterization 280 ln_skin_wl = .false. ! use the warm-layer parameterization 281 ! ! ==> only available in ECMWF and COARE algorithms 282 ln_humi_sph = .true. ! humidity "sn_humi" is specific humidity [kg/kg] 283 ln_humi_dpt = .false. ! humidity "sn_humi" is dew-point temperature [K] 284 ln_humi_rlh = .false. ! humidity "sn_humi" is relative humidity [%] 285 ! 273 286 cn_dir = './' ! root directory for the bulk data location 274 287 !___________!_________________________!___________________!___________!_____________!________!___________!______________________________________!__________!_______________! … … 281 294 sn_tair = 't_10.15JUNE2009_fill' , 6. , 'T_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 282 295 sn_humi = 'q_10.15JUNE2009_fill' , 6. , 'Q_10_MOD', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 296 sn_hpgi = 'NOT USED' , 24. , 'uhpg' , .false. , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'UG' , '' 297 sn_hpgj = 'NOT USED' , 24. , 'vhpg' , .false. , .false., 'monthly' , 'weights_ERAI3D_F128_2_ORCA2_bicubic', 'VG' , '' 283 298 sn_prec = 'ncar_precip.15JUNE2009_fill', -1. , 'PRC_MOD1', .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 284 299 sn_snow = 'ncar_precip.15JUNE2009_fill', -1. , 'SNOW' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 285 300 sn_slp = 'slp.15JUNE2009_fill' , 6. , 'SLP' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 286 sn_tdif = 'taudif_core' , 24 , 'taudif' , .false. , .true. , 'yearly' , 'weights_core_orca2_bilinear_noc.nc' , '' , '' 301 / 302 !----------------------------------------------------------------------- 303 &namsbc_abl ! Atmospheric Boundary Layer formulation (ln_abl = T) 304 !----------------------------------------------------------------------- 305 cn_dir = './' ! root directory for the location of the ABL grid file 306 cn_dom = 'dom_cfg_abl.nc' 307 308 cn_ablrst_in = "restart_abl" ! suffix of abl restart name (input) 309 cn_ablrst_out = "restart_abl" ! suffix of abl restart name (output) 310 cn_ablrst_indir = "." ! directory to read input abl restarts 311 cn_ablrst_outdir = "." ! directory to write output abl restarts 312 313 ln_hpgls_frc = .false. 314 ln_geos_winds = .false. 315 nn_dyn_restore = 2 ! restoring option for dynamical ABL variables: = 0 no restoring 316 ! = 1 equatorial restoring 317 ! = 2 global restoring 318 rn_ldyn_min = 4.5 ! magnitude of the nudging on ABL dynamics at the bottom of the ABL [hour] 319 rn_ldyn_max = 1.5 ! magnitude of the nudging on ABL dynamics at the top of the ABL [hour] 320 rn_ltra_min = 4.5 ! magnitude of the nudging on ABL tracers at the bottom of the ABL [hour] 321 rn_ltra_max = 1.5 ! magnitude of the nudging on ABL tracers at the top of the ABL [hour] 322 nn_amxl = 0 ! mixing length: = 0 Deardorff 80 length-scale 323 ! = 1 length-scale based on the distance to the PBL height 324 ! = 2 Bougeault & Lacarrere 89 length-scale 325 rn_Cm = 0.0667 ! 0.126 in MesoNH 326 rn_Ct = 0.1667 ! 0.143 in MesoNH 327 rn_Ce = 0.4 ! 0.4 in MesoNH 328 rn_Ceps = 0.7 ! 0.85 in MesoNH 329 rn_Rod = 0.15 ! c0 in RMCA17 mixing length formulation (not yet implemented) 330 rn_Ric = 0.139 ! Critical Richardson number (to compute PBL height and diffusivities) 287 331 / 288 332 !----------------------------------------------------------------------- … … 293 337 ! ! -> file cplmask.nc with the float variable called cplmask (jpi,jpj,nn_cplmodel) 294 338 nn_cats_cpl = 5 ! Number of sea ice categories over which coupling is to be carried out (if not 1) 295 296 339 !_____________!__________________________!____________!_____________!______________________!________! 297 340 ! ! description ! multiple ! vector ! vector ! vector ! … … 330 373 sn_rcv_wper = 'none' , 'no' , '' , '' , '' 331 374 sn_rcv_wnum = 'none' , 'no' , '' , '' , '' 332 sn_rcv_w strf= 'none' , 'no' , '' , '' , ''375 sn_rcv_wfreq = 'none' , 'no' , '' , '' , '' 333 376 sn_rcv_wdrag = 'none' , 'no' , '' , '' , '' 334 377 sn_rcv_ts_ice = 'none' , 'no' , '' , '' , '' … … 379 422 nn_chldta = 0 ! RGB : Chl data (=1) or cst value (=0) 380 423 rn_si1 = 23.0 ! 2BD : longest depth of extinction 381 424 382 425 cn_dir = './' ! root directory for the chlorophyl data location 383 426 !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! … … 396 439 ln_sssr_bnd = .true. ! flag to bound erp term (associated with nn_sssr=2) 397 440 rn_sssr_bnd = 4.e0 ! ABS(Max/Min) value of the damping erp term [mm/day] 441 nn_sssr_ice = 1 ! control of sea surface restoring under sea-ice 442 ! 0 = no restoration under ice : * (1-icefrac) 443 ! 1 = restoration everywhere 444 ! >1 = enhanced restoration under ice : 1+(nn_icedmp-1)*icefrac 398 445 399 446 cn_dir = './' ! root directory for the SST/SSS data location … … 443 490 / 444 491 !----------------------------------------------------------------------- 445 &namsbc_isf ! Top boundary layer (ISF) (ln_isfcav =T : read (ln_read_cfg=T) 446 !----------------------------------------------------------------------- or set or usr_def_zgr ) 447 ! ! type of top boundary layer 448 nn_isf = 1 ! ice shelf melting/freezing 449 ! 1 = presence of ISF ; 2 = bg03 parametrisation 450 ! 3 = rnf file for ISF ; 4 = ISF specified freshwater flux 451 ! options 1 and 4 need ln_isfcav = .true. (domzgr) 452 ! ! nn_isf = 1 or 2 cases: 453 rn_gammat0 = 1.e-4 ! gammat coefficient used in blk formula 454 rn_gammas0 = 1.e-4 ! gammas coefficient used in blk formula 455 ! ! nn_isf = 1 or 4 cases: 456 rn_hisf_tbl = 30. ! thickness of the top boundary layer (Losh et al. 2008) 457 ! ! 0 => thickness of the tbl = thickness of the first wet cell 458 ! ! nn_isf = 1 case 459 nn_isfblk = 1 ! 1 ISOMIP like: 2 equations formulation (Hunter et al., 2006) 460 ! ! 2 ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2015) 461 nn_gammablk = 1 ! 0 = cst Gammat (= gammat/s) 462 ! ! 1 = velocity dependend Gamma (u* * gammat/s) (Jenkins et al. 2010) 463 ! ! 2 = velocity and stability dependent Gamma (Holland et al. 1999) 464 465 !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 466 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 467 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 468 !* nn_isf = 4 case 469 sn_fwfisf = 'rnfisf' , -12. ,'sowflisf' , .false. , .true. , 'yearly' , '' , '' , '' 470 !* nn_isf = 3 case 471 sn_rnfisf = 'rnfisf' , -12. ,'sofwfisf' , .false. , .true. , 'yearly' , '' , '' , '' 472 !* nn_isf = 2 and 3 cases 473 sn_depmax_isf ='rnfisf' , -12. ,'sozisfmax', .false. , .true. , 'yearly' , '' , '' , '' 474 sn_depmin_isf ='rnfisf' , -12. ,'sozisfmin', .false. , .true. , 'yearly' , '' , '' , '' 475 !* nn_isf = 2 case 476 sn_Leff_isf = 'rnfisf' , -12. ,'Leff' , .false. , .true. , 'yearly' , '' , '' , '' 477 / 478 !----------------------------------------------------------------------- 479 &namsbc_iscpl ! land ice / ocean coupling option (ln_isfcav =T : read (ln_read_cfg=T) 480 !----------------------------------------------------------------------- or set or usr_def_zgr ) 481 nn_drown = 10 ! number of iteration of the extrapolation loop (fill the new wet cells) 482 ln_hsb = .false. ! activate conservation module (conservation exact after a time of rn_fiscpl) 483 nn_fiscpl = 43800 ! (number of time step) conservation period (maybe should be fix to the coupling frequencey of restart frequency) 492 &namisf ! Top boundary layer (ISF) (default: OFF) 493 !----------------------------------------------------------------------- 494 ! 495 ! ---------------- ice shelf load ------------------------------- 496 ! 497 cn_isfload = 'uniform' ! scheme to compute ice shelf load (ln_isfcav = .true. in domain_cfg.nc) 498 rn_isfload_T = -1.9 499 rn_isfload_S = 34.4 500 ! 501 ! ---------------- ice shelf melt formulation ------------------------------- 502 ! 503 ln_isf = .false. ! activate ice shelf module 504 ln_isfdebug = .false. ! add debug print in ISF code (global min/max/sum of specific variable) 505 cn_isfdir = './' ! directory for all ice shelf input file 506 ! 507 ! ---------------- cavities opened ------------------------------- 508 ! 509 ln_isfcav_mlt = .false. ! ice shelf melting into the cavity (need ln_isfcav = .true. in domain_cfg.nc) 510 cn_isfcav_mlt = '3eq' ! ice shelf melting formulation (spe/2eq/3eq/oasis) 511 ! ! spe = fwfisf is read from a forcing field 512 ! ! 2eq = ISOMIP like: 2 equations formulation (Hunter et al., 2006 for a short description) 513 ! ! 3eq = ISOMIP+ like: 3 equations formulation (Asay-Davis et al., 2016 for a short description) 514 ! ! oasis = fwfisf is given by oasis and pattern by file sn_isfcav_fwf 515 ! ! cn_isfcav_mlt = 2eq or 3eq cases: 516 cn_gammablk = 'vel' ! scheme to compute gammat/s (spe,ad15,hj99) 517 ! ! spe = constant transfert velocity (rn_gammat0, rn_gammas0) 518 ! ! vel = velocity dependent transfert velocity (u* * gammat/s) (Asay-Davis et al. 2016 for a short description) 519 ! ! vel_stab = velocity and stability dependent transfert coeficient (Holland et al. 1999 for a complete description) 520 rn_gammat0 = 1.4e-2 ! gammat coefficient used in spe, vel and vel_stab gamma computation method 521 rn_gammas0 = 4.0e-4 ! gammas coefficient used in spe, vel and vel_stab gamma computation method 522 ! 523 rn_htbl = 30. ! thickness of the top boundary layer (Losh et al. 2008) 524 ! ! 0 => thickness of the tbl = thickness of the first wet cell 525 ! 526 !* 'spe' and 'oasis' case 527 !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 528 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 529 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 530 sn_isfcav_fwf = 'isfmlt_cav', -12. , 'fwflisf' , .false. , .true. , 'yearly' , '' , '' , '' 531 ! 532 ! ---------------- cavities parametrised ------------------------------- 533 ! 534 ln_isfpar_mlt = .false. ! ice shelf melting parametrised 535 cn_isfpar_mlt = 'spe' ! ice shelf melting parametrisation (spe/bg03/oasis) 536 ! ! spe = fwfisf is read from a forcing field 537 ! ! bg03 = melt computed using Beckmann and Goosse parametrisation 538 ! ! oasis = fwfisf is given by oasis and pattern by file sn_isfpar_fwf 539 ! 540 !* all cases 541 !___________!_____________!___________________!___________!_____________!_________!___________!__________!__________!_______________! 542 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights ! rotation ! land/sea mask ! 543 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! filename ! pairing ! filename ! 544 sn_isfpar_zmax = 'isfmlt_par', 0 ,'sozisfmax', .false. , .true. , 'yearly' , '' , '' , '' 545 sn_isfpar_zmin = 'isfmlt_par', 0 ,'sozisfmin', .false. , .true. , 'yearly' , '' , '' , '' 546 !* 'spe' and 'oasis' case 547 sn_isfpar_fwf = 'isfmlt_par' , -12. ,'sofwfisf' , .false. , .true. , 'yearly' , '' , '' , '' 548 !* 'bg03' case 549 sn_isfpar_Leff = 'isfmlt_par', 0. ,'Leff' , .false. , .true. , 'yearly' , '' , '' , '' 550 ! 551 ! ---------------- ice sheet coupling ------------------------------- 552 ! 553 ln_isfcpl = .false. 554 nn_drown = 10 ! number of iteration of the extrapolation loop (fill the new wet cells) 555 ln_isfcpl_cons = .false. 484 556 / 485 557 !----------------------------------------------------------------------- … … 506 578 ln_icebergs = .false. ! activate iceberg floats (force =F with "key_agrif") 507 579 ! 580 ! ! restart 581 cn_icbrst_in = "restart_icb" ! suffix of iceberg restart name (input) 582 cn_icbrst_indir = "./" ! directory from which to read input ocean restarts 583 cn_icbrst_out = "restart_icb" ! suffix of ocean restart name (output) 584 cn_icbrst_outdir = "./" ! directory from which to read output ocean restarts 585 ! 508 586 ! ! diagnostics: 509 587 ln_bergdia = .true. ! Calculate budgets … … 563 641 &namagrif ! AGRIF zoom ("key_agrif") 564 642 !----------------------------------------------------------------------- 643 ln_agrif_2way = .true. ! activate two way nesting 565 644 ln_spc_dyn = .true. ! use 0 as special value for dynamics 566 645 rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [m2/s] 567 646 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [m2/s] 647 rn_trelax_tra = 0.01 ! inverse of relaxation time (in steps) for tracers [] 648 rn_trelax_dyn = 0.01 ! inverse of relaxation time (in steps) for dynamics [] 568 649 ln_chk_bathy = .false. ! =T check the parent bathymetry 569 650 / … … 572 653 !----------------------------------------------------------------------- 573 654 ln_tide = .false. ! Activate tides 574 ln_tide_pot = .true. ! use tidal potential forcing 655 nn_tide_var = 1 ! Variant of tidal parameter set and tide-potential computation 656 ! ! (1: default; 0: compatibility with previous versions) 657 ln_tide_dia = .false. ! Enable tidal diagnostic output 658 ln_tide_pot = .false. ! use tidal potential forcing 659 rn_tide_gamma = 0.7 ! Tidal tilt factor 575 660 ln_scal_load = .false. ! Use scalar approximation for 576 661 rn_scal_load = 0.094 ! load potential 577 662 ln_read_load = .false. ! Or read load potential from file 578 663 cn_tide_load = 'tide_LOAD_grid_T.nc' ! filename for load potential 579 ! 664 ! 580 665 ln_tide_ramp = .false. ! Use linear ramp for tides at startup 581 r dttideramp = 0.! ramp duration in days582 clname(1)= 'DUMMY' ! name of constituent - all tidal components must be set in namelist_cfg666 rn_tide_ramp_dt = 0. ! ramp duration in days 667 sn_tide_cnames(1) = 'DUMMY' ! name of constituent - all tidal components must be set in namelist_cfg 583 668 / 584 669 !----------------------------------------------------------------------- … … 656 741 filtide = 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files 657 742 ln_bdytide_2ddta = .false. ! 658 ln_bdytide_conj = .false. !659 743 / 660 744 … … 683 767 !----------------------------------------------------------------------- 684 768 rn_Cd0 = 1.e-3 ! drag coefficient [-] 685 rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) 769 rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) 686 770 rn_Cdmax = 0.1 ! drag value maximum [-] (logarithmic drag) 687 771 rn_ke0 = 2.5e-3 ! background kinetic energy [m2/s2] (non-linear cases) … … 694 778 !----------------------------------------------------------------------- 695 779 rn_Cd0 = 1.e-3 ! drag coefficient [-] 696 rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) 780 rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) 697 781 rn_Cdmax = 0.1 ! drag value maximum [-] (logarithmic drag) 698 782 rn_ke0 = 2.5e-3 ! background kinetic energy [m2/s2] (non-linear cases) … … 761 845 nn_cen_v = 4 ! =2/4, vertical 2nd order CEN / 4th order COMPACT 762 846 ln_traadv_fct = .false. ! FCT scheme 763 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 764 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 847 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 848 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 765 849 ln_traadv_mus = .false. ! MUSCL scheme 766 850 ln_mus_ups = .false. ! use upstream scheme near river mouths … … 783 867 ln_traldf_triad = .false. ! iso-neutral (triad operator) 784 868 ! 785 ! ! iso-neutral options: 869 ! ! iso-neutral options: 786 870 ln_traldf_msc = .false. ! Method of Stabilizing Correction (both operators) 787 871 rn_slpmax = 0.01 ! slope limit (both operators) … … 793 877 nn_aht_ijk_t = 0 ! space/time variation of eddy coefficient: 794 878 ! ! =-20 (=-30) read in eddy_diffusivity_2D.nc (..._3D.nc) file 795 ! ! = 0 constant 796 ! ! = 10 F(k) =ldf_c1d 797 ! ! = 20 F(i,j) =ldf_c2d 879 ! ! = 0 constant 880 ! ! = 10 F(k) =ldf_c1d 881 ! ! = 20 F(i,j) =ldf_c2d 798 882 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 799 883 ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d 800 884 ! ! = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 801 ! ! time invariant coefficients: aht0 = 1/2 Ud*Ld (lap case) 885 ! ! time invariant coefficients: aht0 = 1/2 Ud*Ld (lap case) 802 886 ! ! or = 1/12 Ud*Ld^3 (blp case) 803 887 rn_Ud = 0.01 ! lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) … … 825 909 nn_aei_ijk_t = 0 ! space/time variation of eddy coefficient: 826 910 ! ! =-20 (=-30) read in eddy_induced_velocity_2D.nc (..._3D.nc) file 827 ! ! = 0 constant 828 ! ! = 10 F(k) =ldf_c1d 829 ! ! = 20 F(i,j) =ldf_c2d 911 ! ! = 0 constant 912 ! ! = 10 F(k) =ldf_c1d 913 ! ! = 20 F(i,j) =ldf_c2d 830 914 ! ! = 21 F(i,j,t) =Treguier et al. JPO 1997 formulation 831 915 ! ! = 30 F(i,j,k) =ldf_c2d * ldf_c1d 832 ! ! time invariant coefficients: aei0 = 1/2 Ue*Le 916 ! ! time invariant coefficients: aei0 = 1/2 Ue*Le 833 917 rn_Ue = 0.02 ! lateral diffusive velocity [m/s] (nn_aht_ijk_t= 0, 10, 20, 30) 834 918 rn_Le = 200.e+3 ! lateral diffusive length [m] (nn_aht_ijk_t= 0, 10) … … 870 954 rn_lf_cutoff = 5.0 ! cutoff frequency for low-pass filter [days] 871 955 rn_zdef_max = 0.9 ! maximum fractional e3t deformation 872 ln_vvl_dbg = . true.! debug prints (T/F)956 ln_vvl_dbg = .false. ! debug prints (T/F) 873 957 / 874 958 !----------------------------------------------------------------------- … … 890 974 ln_dynvor_eeT = .false. ! energy conserving scheme (een using e3t) 891 975 ln_dynvor_een = .false. ! energy & enstrophy scheme 892 nn_een_e3f = 0 ! =0 e3f = mi(mj(e3t))/4 976 nn_een_e3f = 0 ! =0 e3f = mi(mj(e3t))/4 893 977 ! ! =1 e3f = mi(mj(e3t))/mi(mj( tmask)) 894 978 ln_dynvor_msk = .false. ! vorticity multiplied by fmask (=T) ==>>> PLEASE DO NOT ACTIVATE … … 913 997 ln_bt_av = .true. ! Time filtering of barotropic variables 914 998 nn_bt_flt = 1 ! Time filter choice = 0 None 915 ! ! = 1 Boxcar over nn_ barosub-steps916 ! ! = 2 Boxcar over 2*nn_ baro" "999 ! ! = 1 Boxcar over nn_e sub-steps 1000 ! ! = 2 Boxcar over 2*nn_e " " 917 1001 ln_bt_auto = .true. ! Number of sub-step defined from: 918 1002 rn_bt_cmax = 0.8 ! =T : the Maximum Courant Number allowed 919 nn_ baro = 30 ! =F : the number of sub-step in rn_rdt seconds1003 nn_e = 30 ! =F : the number of sub-step in rn_Dt seconds 920 1004 rn_bt_alpha = 0. ! Temporal diffusion parameter (if ln_bt_av=F) 921 1005 / … … 935 1019 ! ! =-30 read in eddy_viscosity_3D.nc file 936 1020 ! ! =-20 read in eddy_viscosity_2D.nc file 937 ! ! = 0 constant 1021 ! ! = 0 constant 938 1022 ! ! = 10 F(k)=c1d 939 1023 ! ! = 20 F(i,j)=F(grid spacing)=c2d … … 941 1025 ! ! = 31 F(i,j,k)=F(grid spacing and local velocity) 942 1026 ! ! = 32 F(i,j,k)=F(local gridscale and deformation rate) 943 ! ! time invariant coefficients : ahm = 1/2 Uv*Lv (lap case) 1027 ! ! time invariant coefficients : ahm = 1/2 Uv*Lv (lap case) 944 1028 ! ! or = 1/12 Uv*Lv^3 (blp case) 945 1029 rn_Uv = 0.1 ! lateral viscous velocity [m/s] (nn_ahm_ijk_t= 0, 10, 20, 30) … … 1065 1149 ! = 0 constant 10 m length scale 1066 1150 ! = 1 0.5m at the equator to 30m poleward of 40 degrees 1067 rn_eice = 4 ! below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4 1151 rn_eice = 4 ! below sea ice: =0 ON ; =4 OFF when ice fraction > 1/4 1068 1152 / 1069 1153 !----------------------------------------------------------------------- … … 1124 1208 ln_mevar = .true. ! variable (T) or constant (F) mixing efficiency 1125 1209 ln_tsdiff = .true. ! account for differential T/S mixing (T) or not (F) 1126 / 1127 1210 1211 cn_dir = './' ! root directory for the iwm data location 1212 !___________!_________________________!___________________!___________!_____________!________!___________!__________________!__________!_______________! 1213 ! ! file name ! frequency (hours) ! variable ! time interp.! clim ! 'yearly'/ ! weights filename ! rotation ! land/sea mask ! 1214 ! ! ! (if <0 months) ! name ! (logical) ! (T/F) ! 'monthly' ! ! pairing ! filename ! 1215 sn_mpb = 'NOT USED' , -12 , 'mixing_power_bot' , .false. , .true. , 'yearly' , '' , '' , '' 1216 sn_mpp = 'NOT USED' , -12 , 'mixing_power_pyc' , .false. , .true. , 'yearly' , '' , '' , '' 1217 sn_mpc = 'NOT USED' , -12 , 'mixing_power_cri' , .false. , .true. , 'yearly' , '' , '' , '' 1218 sn_dsb = 'NOT USED' , -12 , 'decay_scale_bot' , .false. , .true. , 'yearly' , '' , '' , '' 1219 sn_dsc = 'NOT USED' , -12 , 'decay_scale_cri' , .false. , .true. , 'yearly' , '' , '' , '' 1220 / 1128 1221 !!====================================================================== 1129 1222 !! *** Diagnostics namelists *** !! … … 1135 1228 !! namdiu Cool skin and warm layer models (default: OFF) 1136 1229 !! namflo float parameters (default: OFF) 1137 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)1138 1230 !! nam_diadct transports through some sections (default: OFF) 1139 !! nam_diatmb Top Middle Bottom Output (default: OFF)1140 1231 !! nam_dia25h 25h Mean Output (default: OFF) 1141 1232 !! namnc4 netcdf4 chunking and compression settings ("key_netcdf4") … … 1162 1253 !!gm ln_trdmld_instant = .false. ! flag to diagnose trends of instantantaneous or mean ML T/S 1163 1254 !!gm 1164 !-----------------------------------------------------------------------1165 &namptr ! Poleward Transport Diagnostic (default: OFF)1166 !-----------------------------------------------------------------------1167 ln_diaptr = .false. ! Poleward heat and salt transport (T) or not (F)1168 ln_subbas = .false. ! Atlantic/Pacific/Indian basins computation (T) or not1169 1255 / 1170 1256 !----------------------------------------------------------------------- … … 1195 1281 / 1196 1282 !----------------------------------------------------------------------- 1197 &nam_diaharm ! Harmonic analysis of tidal constituents (default: OFF)1198 !-----------------------------------------------------------------------1199 ln_diaharm = .false. ! Choose tidal harmonic output or not1200 nit000_han = 1 ! First time step used for harmonic analysis1201 nitend_han = 75 ! Last time step used for harmonic analysis1202 nstep_han = 15 ! Time step frequency for harmonic analysis1203 tname(1) = 'M2' ! Name of tidal constituents1204 tname(2) = 'K1' ! ---1205 /1206 !-----------------------------------------------------------------------1207 1283 &nam_diadct ! transports through some sections (default: OFF) 1208 1284 !----------------------------------------------------------------------- … … 1213 1289 ! ! -1 : debug all section 1214 1290 ! ! 0 < n : debug section number n 1215 /1216 !-----------------------------------------------------------------------1217 &nam_diatmb ! Top Middle Bottom Output (default: OFF)1218 !-----------------------------------------------------------------------1219 ln_diatmb = .false. ! Choose Top Middle and Bottom output or not1220 1291 / 1221 1292 !----------------------------------------------------------------------- … … 1339 1410 &namctl ! Control prints (default: OFF) 1340 1411 !----------------------------------------------------------------------- 1341 ln_ctl = .FALSE. ! Toggle all report printing on/off (T/F); Ignored if sn_cfctl%l_config is T 1412 sn_cfctl%l_glochk = .FALSE. ! Range sanity checks are local (F) or global (T). Set T for debugging only 1413 sn_cfctl%l_allon = .FALSE. ! IF T activate all options. If F deactivate all unless l_config is T 1342 1414 sn_cfctl%l_config = .TRUE. ! IF .true. then control which reports are written with the following 1343 sn_cfctl%l_runstat = . FALSE.! switches and which areas produce reports with the proc integer settings.1415 sn_cfctl%l_runstat = .TRUE. ! switches and which areas produce reports with the proc integer settings. 1344 1416 sn_cfctl%l_trcstat = .FALSE. ! The default settings for the proc integers should ensure 1345 1417 sn_cfctl%l_oceout = .FALSE. ! that all areas report. 1346 1418 sn_cfctl%l_layout = .FALSE. ! 1347 sn_cfctl%l_mppout = .FALSE. ! 1348 sn_cfctl%l_mpptop = .FALSE. ! 1419 sn_cfctl%l_prtctl = .FALSE. ! 1420 sn_cfctl%l_prttrc = .FALSE. ! 1421 sn_cfctl%l_oasout = .FALSE. ! 1349 1422 sn_cfctl%procmin = 0 ! Minimum area number for reporting [default:0] 1350 1423 sn_cfctl%procmax = 1000000 ! Maximum area number for reporting [default:1000000]
Note: See TracChangeset
for help on using the changeset viewer.