- Timestamp:
- 2017-12-13T15:58:53+01:00 (7 years ago)
- Location:
- branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/ISOMIP
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/ISOMIP/EXP00/namelist_cfg
r8599 r9019 2 2 !! NEMO/OPA : ISOMIP Configuration namelist to overwrite reference dynamical namelist 3 3 !!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 4 !-----------------------------------------------------------------------5 &namusr_def ! ISOMIP user defined namelist6 !-----------------------------------------------------------------------7 ln_zps = .true. ! z-partial-step coordinate8 rn_lam0 = 0.0 ! longitude of first raw and column T-point (jphgr_msh = 1)9 rn_phi0 = -80.0 ! latitude of first raw and column T-point (jphgr_msh = 1)10 rn_e1deg = 0.3 ! zonal grid-spacing (degrees)11 rn_e2deg = 0.1 ! meridional grid-spacing (degrees)12 rn_e3 = 30. ! vertical resolution13 /14 4 !----------------------------------------------------------------------- 15 5 &namrun ! parameters of the run … … 29 19 ln_read_cfg = .false. ! (=T) read the domain configuration file 30 20 ! ! (=F) user defined configuration ==>>> see usrdef(_...) modules 31 !32 ln_write_cfg= .true. ! (=T) create the domain configuration file33 cn_domcfg_out = "ISOMIP_cfg_out" ! newly created domain configuration filename34 21 / 35 22 !----------------------------------------------------------------------- … … 38 25 ln_linssh = .false. ! =T linear free surface ==>> model level are fixed in time 39 26 ! 40 nn_msh = 1! create (>0) a mesh file or not (=0)27 nn_msh = 0 ! create (>0) a mesh file or not (=0) 41 28 ! 42 29 rn_rdt = 1800. ! time step for the dynamics (and tracer if nn_acc=0) … … 64 51 ! =1 use observed ice-cover , 65 52 ! =2 ice-model used 66 nn_ice_embd = 0 ! =0 levitating ice (no mass exchange, concentration/dilution effect)67 ! = 1 levitating ice with mass and salt exchange but no presure effect68 ! =2 embedded sea-ice (full salt and mass exchanges and pressure)53 ln_ice_embd = .false. ! =F levitating ice with mass and salt exchange but no presure effect 54 ! =T embedded sea-ice (full salt and mass exchanges and pressure) 55 ! Misc. options of sbc : 69 56 ln_traqsr = .false. ! Light penetration (T) or not (F) 70 57 ln_rnf = .false. ! runoffs (T => fill namsbc_rnf) … … 179 166 / 180 167 !----------------------------------------------------------------------- 181 &nambfr ! bottom friction 182 !----------------------------------------------------------------------- 183 nn_bfr = 2 ! type of bottom friction : = 0 : free slip, = 1 : linear friction 184 ! = 2 : nonlinear friction 185 rn_bfri1 = 4.e-4 ! bottom drag coefficient (linear case) 186 rn_bfri2 = 1.e-3 ! bottom drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 187 rn_bfri2_max = 1.e-1 ! max. bottom drag coefficient (non linear case and ln_loglayer=T) 188 rn_bfeb2 = 2.5e-3 ! bottom turbulent kinetic energy background (m2/s2) 189 rn_bfrz0 = 3.e-3 ! bottom roughness [m] if ln_loglayer=T 190 ln_bfr2d = .false. ! horizontal variation of the bottom friction coef (read a 2D mask file ) 191 rn_bfrien = 50. ! local multiplying factor of bfr (ln_bfr2d=T) 192 rn_tfri1 = 4.e-4 ! top drag coefficient (linear case) 193 rn_tfri2 = 2.5e-3 ! top drag coefficient (non linear case). Minimum coeft if ln_loglayer=T 194 rn_tfri2_max = 1.e-1 ! max. top drag coefficient (non linear case and ln_loglayer=T) 195 rn_tfeb2 = 0.0 ! top turbulent kinetic energy background (m2/s2) 196 rn_tfrz0 = 3.e-3 ! top roughness [m] if ln_loglayer=T 197 ln_tfr2d = .false. ! horizontal variation of the top friction coef (read a 2D mask file ) 198 rn_tfrien = 50. ! local multiplying factor of tfr (ln_tfr2d=T) 199 200 ln_bfrimp = .true. ! implicit bottom friction (requires ln_zdfexp = .false. if true) 201 ln_loglayer = .false. ! logarithmic formulation (non linear case) 168 &namdrg ! top/bottom drag coefficient (default: NO selection) 169 !----------------------------------------------------------------------- 170 ln_non_lin = .true. ! non-linear drag: Cd = Cd0 |U| 171 / 172 !----------------------------------------------------------------------- 173 &namdrg_top ! TOP friction (ln_isfcav=T) 174 !----------------------------------------------------------------------- 175 rn_Cd0 = 2.5e-3 ! drag coefficient [-] 176 rn_Uc0 = 0.16 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) 177 rn_Cdmax = 0.1 ! drag value maximum [-] (logarithmic drag) 178 rn_ke0 = 0.0e-0 ! background kinetic energy [m2/s2] (non-linear cases) 179 rn_z0 = 3.0e-3 ! roughness [m] (ln_loglayer=T) 180 ln_boost = .false. ! =T regional boost of Cd0 ; =F constant 181 rn_boost= 50. ! local boost factor [-] 182 / 183 !----------------------------------------------------------------------- 184 &namdrg_bot ! BOTTOM friction 185 !----------------------------------------------------------------------- 186 rn_Cd0 = 1.e-3 ! drag coefficient [-] 187 rn_Uc0 = 0.4 ! ref. velocity [m/s] (linear drag=Cd0*Uc0) 188 rn_Cdmax = 0.1 ! drag value maximum [-] (logarithmic drag) 189 rn_ke0 = 2.5e-3 ! background kinetic energy [m2/s2] (non-linear cases) 190 rn_z0 = 3.e-3 ! roughness [m] (ln_loglayer=T) 191 ln_boost = .false. ! =T regional boost of Cd0 ; =F constant 192 rn_boost= 50. ! local boost factor [-] 202 193 / 203 194 !----------------------------------------------------------------------- … … 217 208 ln_teos10 = .false. ! = Use TEOS-10 equation of state 218 209 ln_eos80 = .true. ! = Use EOS80 equation of state 219 ! ! rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS220 210 / 221 211 !----------------------------------------------------------------------- … … 225 215 nn_fct_h = 2 ! =2/4, horizontal 2nd / 4th order 226 216 nn_fct_v = 2 ! =2/4, vertical 2nd / COMPACT 4th order 227 nn_fct_zts = 0 ! >=1, 2nd order FCT scheme with vertical sub-timestepping228 ! ! (number of sub-timestep = nn_fct_zts)229 217 / 230 218 !----------------------------------------------------------------------- … … 275 263 !----------------------------------------------------------------------- 276 264 &namdyn_adv ! formulation of the momentum advection 265 !----------------------------------------------------------------------- 266 ln_dynadv_vec = .true. ! vector form (T) or flux form (F) 267 nn_dynkeg = 0 ! scheme for grad(KE): =0 C2 ; =1 Hollingsworth correction 268 ln_dynadv_cen2= .false. ! flux form - 2nd order centered scheme 269 ln_dynadv_ubs = .false. ! flux form - 3rd order UBS scheme 270 / 277 271 !----------------------------------------------------------------------- 278 272 / … … 324 318 / 325 319 !----------------------------------------------------------------------- 326 &namzdf ! vertical physics 327 !----------------------------------------------------------------------- 328 rn_avm0 = 1.0e-3 ! vertical eddy viscosity [m2/s] (background Kz if not "key_zdfcst") 329 rn_avt0 = 5.0e-5 ! vertical eddy diffusivity [m2/s] (background Kz if not "key_zdfcst") 330 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 331 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) 332 ln_zdfevd = .true. ! enhanced vertical diffusion (evd) (T) or not (F) 333 nn_evdm = 1 ! evd apply on tracer (=0) or on tracer and momentum (=1) 334 rn_avevd = 0.1 ! evd mixing coefficient [m2/s] 335 ln_zdfnpc = .false. ! Non-Penetrative Convective algorithm (T) or not (F) 336 nn_npc = 1 ! frequency of application of npc 337 nn_npcp = 365 ! npc control print frequency 338 ln_zdfexp = .false. ! time-stepping: split-explicit (T) or implicit (F) time stepping 339 nn_zdfexp = 3 ! number of sub-timestep for ln_zdfexp=T 340 / 341 !----------------------------------------------------------------------- 342 &namzdf_ric ! richardson number dependent vertical diffusion ("key_zdfric" ) 343 !----------------------------------------------------------------------- 344 / 345 !----------------------------------------------------------------------- 346 &namzdf_tke ! turbulent eddy kinetic dependent vertical diffusion ("key_zdftke") 347 !----------------------------------------------------------------------- 348 / 349 !----------------------------------------------------------------------- 350 &namzdf_gls ! GLS vertical diffusion ("key_zdfgls") 351 !----------------------------------------------------------------------- 352 / 353 !----------------------------------------------------------------------- 354 &namzdf_ddm ! double diffusive mixing parameterization ("key_zdfddm") 355 !----------------------------------------------------------------------- 356 / 357 !----------------------------------------------------------------------- 358 &namzdf_tmx ! tidal mixing parameterization ("key_zdftmx") 359 !----------------------------------------------------------------------- 360 ln_tmx_itf = .false. ! ITF specific parameterisation 320 &namzdf ! vertical physics (default: NO selection) 321 !----------------------------------------------------------------------- 322 ! ! type of vertical closure 323 ln_zdfcst = .true. ! constant mixing 324 ln_zdfric = .false. ! local Richardson dependent formulation (T => fill namzdf_ric) 325 ln_zdftke = .false. ! Turbulent Kinetic Energy closure (T => fill namzdf_tke) 326 ln_zdfgls = .false. ! Generic Length Scale closure (T => fill namzdf_gls) 327 ln_zdfosm = .false. ! OSMOSIS BL closure (T => fill namzdf_osm) 328 ! 329 ! ! convection 330 ln_zdfevd = .true. ! enhanced vertical diffusion 331 nn_evdm = 1 ! apply on tracer (=0) or on tracer and momentum (=1) 332 rn_evd = 0.1 ! mixing coefficient [m2/s] 333 ln_zdfnpc = .false. ! Non-Penetrative Convective algorithm 334 nn_npc = 1 ! frequency of application of npc 335 nn_npcp = 365 ! npc control print frequency 336 ! 337 ln_zdfddm = .false. ! double diffusive mixing 338 ! 339 ! ! gravity wave-driven vertical mixing 340 ln_zdfiwm = .false. ! internal wave-induced mixing (T => fill namzdf_iwm) 341 ln_zdfswm = .false. ! surface wave-induced mixing (T => ln_wave=ln_sdw=T ) 342 ! 343 ! ! coefficients 344 rn_avm0 = 1.e-3 ! vertical eddy viscosity [m2/s] (background Kz if ln_zdfcst=F) 345 rn_avt0 = 5.e-5 ! vertical eddy diffusivity [m2/s] (background Kz if ln_zdfcst=F) 346 nn_avb = 0 ! profile for background avt & avm (=1) or not (=0) 347 nn_havtb = 0 ! horizontal shape for avtb (=1) or not (=0) 361 348 / 362 349 !----------------------------------------------------------------------- … … 378 365 / 379 366 !----------------------------------------------------------------------- 380 &namflo ! float parameters ("key_float")381 !-----------------------------------------------------------------------382 /383 !-----------------------------------------------------------------------384 &namptr ! Poleward Transport Diagnostic385 !-----------------------------------------------------------------------386 /387 !-----------------------------------------------------------------------388 367 &namhsb ! Heat and salt budgets 389 368 !----------------------------------------------------------------------- 390 369 / 391 370 !----------------------------------------------------------------------- 392 &namdct ! transports through sections 393 !----------------------------------------------------------------------- 394 / 395 !----------------------------------------------------------------------- 396 &namobs ! observation usage switch ('key_diaobs') 371 &namobs ! observation usage switch (ln_diaobs =T) 397 372 !----------------------------------------------------------------------- 398 373 / -
branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/ISOMIP/MY_SRC/usrdef_hgr.F90
r8018 r9019 68 68 IF( nn_timing == 1 ) CALL timing_start('usr_def_hgr') 69 69 ! 70 IF(lwp) WRITE(numout,*) 71 IF(lwp) WRITE(numout,*) 'usr_def_hgr : ISOMIP configuration' 72 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 70 IF(lwp) THEN 71 WRITE(numout,*) 72 WRITE(numout,*) 'usr_def_hgr : ISOMIP configuration' 73 WRITE(numout,*) '~~~~~~~~~~~' 74 WRITE(numout,*) 75 WRITE(numout,*) ' ===>> geographical mesh on the sphere with regular grid-spacing' 76 WRITE(numout,*) ' given by rn_e1deg and rn_e2deg' 77 ENDIF 73 78 ! 74 ! !== grid point position ==! (in kilometers) 75 76 IF(lwp) WRITE(numout,*) 77 IF(lwp) WRITE(numout,*) ' geographical mesh on the sphere with regular grid-spacing' 78 IF(lwp) WRITE(numout,*) ' given by rn_e1deg and rn_e2deg' 79 79 ! !== grid point position ==! (in degrees) 80 80 DO jj = 1, jpj 81 81 DO ji = 1, jpi … … 94 94 pphiv(ji,jj) = rn_phi0 + rn_e2deg * zvj 95 95 pphif(ji,jj) = rn_phi0 + rn_e2deg * zfj 96 96 97 ! !== Horizontal scale factors ==! (in meters) 97 98 ! e1 -
branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/ISOMIP/MY_SRC/usrdef_sbc.F90
r7715 r9019 17 17 USE dom_oce ! ocean space and time domain 18 18 USE sbc_oce ! Surface boundary condition: ocean fields 19 USE sbc_ice ! Surface boundary condition: ice fields 19 20 USE phycst ! physical constants 20 21 ! -
branches/2017/dev_merge_2017/NEMOGCM/CONFIG/TEST_CASES/ISOMIP/cpp_ISOMIP.fcm
r7715 r9019 1 bld::tool::fppkeys key_zdfcstkey_iomput key_mpp_mpi key_nosignedzero1 bld::tool::fppkeys key_iomput key_mpp_mpi key_nosignedzero
Note: See TracChangeset
for help on using the changeset viewer.