New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 5870 for branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/SHARED/namelist_ref – NEMO

Ignore:
Timestamp:
2015-11-09T18:33:54+01:00 (8 years ago)
Author:
acc
Message:

Branch 2015/dev_r5803_NOC_WAD. Merge in trunk changes from 5803 to 5869 in preparation for merge. Also tidied and reorganised some wetting and drying code. Renamed wadlmt.F90 to wetdry.F90. Wetting drying code changes restricted to domzgr.F90, domvvl.F90 nemogcm.F90 sshwzv.F90, dynspg_ts.F90, wetdry.F90 and dynhpg.F90. Code passes full SETTE tests with ln_wd=.false.. Still awaiting test case for checking with ln_wd=.false.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/SHARED/namelist_ref

    r5656 r5870  
     1!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
     2!!                            namelist_ref 
    13!!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
    24!! NEMO/OPA  :  1 - run manager      (namrun) 
     
    57!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf, 
    68!!                                    namsbc_apr, namsbc_ssr, namsbc_alb) 
    7 !!              4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 
     9!!              4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 
    810!!              5 - bottom  boundary (nambfr, nambbc, nambbl) 
    9 !!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_dmp) 
     11!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_ldfeiv, namtra_dmp) 
    1012!!              7 - dynamics         (namdyn_adv, namdyn_vor, namdyn_hpg, namdyn_spg, namdyn_ldf) 
    11 !!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_kpp, namzdf_ddm, namzdf_tmx) 
     13!!              8 - Verical physics  (namzdf, namzdf_ric, namzdf_tke, namzdf_ddm, namzdf_tmx) 
    1214!!              9 - diagnostics      (namnc4, namtrd, namspr, namflo, namhsb, namsto) 
    1315!!             10 - miscellaneous    (namsol, nammpp, namctl) 
     
    200202!----------------------------------------------------------------------- 
    201203&namtsd    !   data : Temperature  & Salinity 
    202 !----------------------------------------------------------------------- 
    203204!----------------------------------------------------------------------- 
    204205!          !  file name                            ! frequency (hours) ! variable  ! time interp. !  clim  ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     
    545546!!====================================================================== 
    546547!!   namlbc        lateral momentum boundary condition 
    547 !!   namcla        cross land advection 
    548548!!   namobc        open boundaries parameters                           ("key_obc") 
    549549!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif") 
     
    558558                           !  free slip  !   partial slip  !   no slip   ! strong slip 
    559559   ln_vorlat   = .false.   !  consistency of vorticity boundary condition with analytical eqs. 
    560 / 
    561 !----------------------------------------------------------------------- 
    562 &namcla        !   cross land advection 
    563 !----------------------------------------------------------------------- 
    564    nn_cla      =    0      !  advection between 2 ocean pts separates by land 
    565 / 
    566 !----------------------------------------------------------------------- 
    567 &namobc        !   open boundaries parameters                           ("key_obc") 
    568 !----------------------------------------------------------------------- 
    569    ln_obc_clim = .false.   !  climatological obc data files (T) or not (F) 
    570    ln_vol_cst  = .true.    !  impose the total volume conservation (T) or not (F) 
    571    ln_obc_fla  = .false.   !  Flather open boundary condition 
    572    nn_obcdta   =    1      !  = 0 the obc data are equal to the initial state 
    573                            !  = 1 the obc data are read in 'obc.dta' files 
    574    cn_obcdta   = 'annual'  !  set to annual if obc datafile hold 1 year of data 
    575                            !  set to monthly if obc datafile hold 1 month of data 
    576    rn_dpein    =    1.     !  damping time scale for inflow at east  open boundary 
    577    rn_dpwin    =    1.     !     -           -         -       west    -      - 
    578    rn_dpnin    =    1.     !     -           -         -       north   -      - 
    579    rn_dpsin    =    1.     !     -           -         -       south   -      - 
    580    rn_dpeob    = 3000.     !  time relaxation (days) for the east  open boundary 
    581    rn_dpwob    =   15.     !     -           -         -     west    -      - 
    582    rn_dpnob    = 3000.     !     -           -         -     north   -      - 
    583    rn_dpsob    =   15.     !     -           -         -     south   -      - 
    584    rn_volemp   =    1.     !  = 0 the total volume change with the surface flux (E-P-R) 
    585                            !  = 1 the total volume remains constant 
    586560/ 
    587561!----------------------------------------------------------------------- 
     
    711685                           !     = 2 variable flux (read in geothermal_heating.nc in mW/m2) 
    712686   rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2] 
    713  
    714687/ 
    715688!----------------------------------------------------------------------- 
     
    725698!!                        Tracer (T & S ) namelists 
    726699!!====================================================================== 
    727 !!   nameos        equation of state 
    728 !!   namtra_adv    advection scheme 
     700!!   nameos           equation of state 
     701!!   namtra_adv       advection scheme 
    729702!!   namtra_adv_mle   mixed layer eddy param. (Fox-Kemper param.) 
    730 !!   namtra_ldf    lateral diffusion scheme 
    731 !!   namtra_dmp    T & S newtonian damping 
     703!!   namtra_ldf       lateral diffusion scheme 
     704!!   namtra_ldfeiv    eddy induced velocity param. 
     705!!   namtra_dmp       T & S newtonian damping 
    732706!!====================================================================== 
    733707! 
     
    740714                                 !  = 1, S-EOS   (simplified eos) 
    741715   ln_useCT    = .true.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 
    742    !                             ! 
     716                                 ! 
    743717   !                     ! S-EOS coefficients : 
    744    !                             !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
     718                                 !  rd(T,S,Z)*rau0 = -a0*(1+.5*lambda*dT+mu*Z+nu*dS)*dT+b0*dS 
    745719   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1) 
    746720   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1) 
     
    754728&namtra_adv    !   advection scheme for tracer 
    755729!----------------------------------------------------------------------- 
    756    ln_traadv_cen2   =  .false.   !  2nd order centered scheme 
    757    ln_traadv_tvd    =  .true.    !  TVD scheme 
    758    ln_traadv_muscl  =  .false.   !  MUSCL scheme 
    759    ln_traadv_muscl2 =  .false.   !  MUSCL2 scheme + cen2 at boundaries 
    760    ln_traadv_ubs    =  .false.   !  UBS scheme 
    761    ln_traadv_qck    =  .false.   !  QUICKEST scheme 
    762    ln_traadv_msc_ups=  .false.   !  use upstream scheme within muscl 
    763    ln_traadv_tvd_zts=  .false.  !  TVD scheme with sub-timestepping of vertical tracer advection 
     730   ln_traadv_cen =  .false.  !  2nd order centered scheme 
     731      nn_cen_h   =  4               !  =2/4, horizontal 2nd order CEN / 4th order CEN 
     732      nn_cen_v   =  4               !  =2/4, vertical   2nd order CEN / 4th order COMPACT 
     733   ln_traadv_fct =  .false.  !  FCT scheme 
     734      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     735      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     736      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     737      !                             !        (number of sub-timestep = nn_fct_zts) 
     738   ln_traadv_mus =  .false.  !  MUSCL scheme 
     739      ln_mus_ups =  .false.         !  use upstream scheme near river mouths 
     740   ln_traadv_ubs =  .false.  !  UBS scheme 
     741      nn_ubs_v   =  2               !  =2  , vertical 2nd order FCT 
     742   ln_traadv_qck =  .false.  !  QUICKEST scheme 
    764743/ 
    765744!----------------------------------------------------------------------- 
    766745&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param) 
    767746!----------------------------------------------------------------------- 
    768    ln_mle    = .true.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
     747   ln_mle    = .false.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
    769748   rn_ce     = 0.06        ! magnitude of the MLE (typical value: 0.06 to 0.08) 
    770749   nn_mle    = 1           ! MLE type: =0 standard Fox-Kemper ; =1 new formulation 
     
    780759!---------------------------------------------------------------------------------- 
    781760   !                       !  Operator type: 
    782    ln_traldf_lap    =  .true.   !  laplacian operator 
    783    ln_traldf_bilap  =  .false.  !  bilaplacian operator 
     761   !                           !  no diffusion: set ln_traldf_lap=..._blp=F  
     762   ln_traldf_lap   =  .false.  !    laplacian operator 
     763   ln_traldf_blp   =  .false.  !  bilaplacian operator 
    784764   !                       !  Direction of action: 
    785    ln_traldf_level  =  .false.  !  iso-level 
    786    ln_traldf_hor    =  .false.  !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T) 
    787    ln_traldf_iso    =  .true.   !  iso-neutral                 (needs "key_ldfslp") 
    788    !                 !  Griffies parameters              (all need "key_ldfslp") 
    789    ln_traldf_grif   =  .false.  !  use griffies triads 
    790    ln_traldf_gdia   =  .false.  !  output griffies eddy velocities 
    791    ln_triad_iso     =  .false.  !  pure lateral mixing in ML 
    792    ln_botmix_grif   =  .false.  !  lateral mixing on bottom 
    793    !                       !  Coefficients 
    794    ! Eddy-induced (GM) advection always used with Griffies; otherwise needs "key_traldf_eiv" 
    795    ! Value rn_aeiv_0 is ignored unless = 0 with Held-Larichev spatially varying aeiv 
    796    !                                  (key_traldf_c2d & key_traldf_eiv & key_orca_r2, _r1 or _r05) 
    797    rn_aeiv_0        =  2000.    !  eddy induced velocity coefficient [m2/s] 
    798    rn_aht_0         =  2000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    799    rn_ahtb_0        =     0.    !  background eddy diffusivity for ldf_iso [m2/s] 
    800    !                                           (normally=0; not used with Griffies) 
    801    rn_slpmax        =     0.01  !  slope limit 
    802    rn_chsmag        =     1.    !  multiplicative factor in Smagorinsky diffusivity 
    803    rn_smsh          =     1.    !  Smagorinsky diffusivity: = 0 - use only sheer 
    804    rn_aht_m         =  2000.    !  upper limit or stability criteria for lateral eddy diffusivity (m2/s) 
     765   ln_traldf_lev   =  .false.  !  iso-level 
     766   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     767   ln_traldf_iso   =  .false.  !  iso-neutral (standard operator) 
     768   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     769   ! 
     770   !                       !  iso-neutral options:         
     771   ln_traldf_msc   =  .false.  !  Method of Stabilizing Correction (both operators) 
     772   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     773   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     774   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     775   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     776   ! 
     777   !                       !  Coefficients: 
     778   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     779   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     780   !                                !   =  0           constant  
     781   !                                !   = 10 F(k)      =ldf_c1d  
     782   !                                !   = 20 F(i,j)    =ldf_c2d  
     783   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     784   !                                !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
     785   !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
     786   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     787   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     788/ 
     789!---------------------------------------------------------------------------------- 
     790&namtra_ldfeiv !   eddy induced velocity param. 
     791!---------------------------------------------------------------------------------- 
     792   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
     793   ln_ldfeiv_dia =.false.   ! diagnose eiv stream function and velocities 
     794   rn_aeiv_0     = 2000.    ! eddy induced velocity coefficient   [m2/s] 
     795   nn_aei_ijk_t  = 21       ! space/time variation of the eiv coeficient 
     796   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     797   !                                !   =  0           constant  
     798   !                                !   = 10 F(k)      =ldf_c1d  
     799   !                                !   = 20 F(i,j)    =ldf_c2d  
     800   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     801   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    805802/ 
    806803!----------------------------------------------------------------------- 
     
    853850   ln_dynvor_ens = .false. !  energy conserving scheme 
    854851   ln_dynvor_mix = .false. !  mixed scheme 
    855    ln_dynvor_een = .true.  !  energy & enstrophy scheme 
    856    ln_dynvor_een_old = .false.  !  energy & enstrophy scheme - original formulation 
     852   ln_dynvor_een = .false. !  energy & enstrophy scheme 
     853      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
     854   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes)  ! PLEASE DO NOT USE 
    857855/ 
    858856!----------------------------------------------------------------------- 
     
    879877!----------------------------------------------------------------------- 
    880878   !                       !  Type of the operator : 
    881    ln_dynldf_lap    =  .true.   !  laplacian operator 
    882    ln_dynldf_bilap  =  .false.  !  bilaplacian operator 
     879   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     880   ln_dynldf_lap =  .false.    !    laplacian operator 
     881   ln_dynldf_blp =  .false.    !  bilaplacian operator 
    883882   !                       !  Direction of action  : 
    884    ln_dynldf_level  =  .false.  !  iso-level 
    885    ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.) 
    886    ln_dynldf_iso    =  .false.  !  iso-neutral                          (require "key_ldfslp") 
     883   ln_dynldf_lev =  .false.    !  iso-level 
     884   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     885   ln_dynldf_iso =  .false.    !  iso-neutral 
    887886   !                       !  Coefficient 
    888    rn_ahm_0_lap     = 40000.    !  horizontal laplacian eddy viscosity   [m2/s] 
    889    rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s] 
    890    rn_ahm_0_blp     =     0.    !  horizontal bilaplacian eddy viscosity [m4/s] 
    891    rn_cmsmag_1      =     3.    !  constant in laplacian Smagorinsky viscosity 
    892    rn_cmsmag_2      =     3     !  constant in bilaplacian Smagorinsky viscosity 
    893    rn_cmsh          =     1.    !  1 or 0 , if 0 -use only shear for Smagorinsky viscosity 
    894    rn_ahm_m_blp     =    -1.e12 !  upper limit for bilap  abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 
    895    rn_ahm_m_lap     = 40000.    !  upper limit for lap  ahm < min(dx^2/16rdt, rn_ahm_m_lap) 
     887   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     888   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     889   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     890   !                                !  =  0  constant  
     891   !                                !  = 10  F(k)=c1d 
     892   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     893   !                                !  = 30  F(i,j,k)=c2d*c1d 
     894   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     895   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     896   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     897   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
     898   ! 
     899   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
    896900/ 
    897901 
     
    902906!!    namzdf_ric    richardson number dependent vertical mixing         ("key_zdfric") 
    903907!!    namzdf_tke    TKE dependent vertical mixing                       ("key_zdftke") 
    904 !!    namzdf_kpp    KPP dependent vertical mixing                       ("key_zdfkpp") 
    905908!!    namzdf_ddm    double diffusive mixing parameterization            ("key_zdfddm") 
    906909!!    namzdf_tmx    tidal mixing parameterization                       ("key_zdftmx") 
     
    963966                           !        = 0  constant 10 m length scale 
    964967                           !        = 1  0.5m at the equator to 30m poleward of 40 degrees 
    965 / 
    966 !------------------------------------------------------------------------ 
    967 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    968 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    969    ln_kpprimix = .true.    !  shear instability mixing 
    970    rn_difmiw   =  1.0e-04  !  constant internal wave viscosity [m2/s] 
    971    rn_difsiw   =  0.1e-04  !  constant internal wave diffusivity [m2/s] 
    972    rn_riinfty  =  0.8      !  local Richardson Number limit for shear instability 
    973    rn_difri    =  0.0050   !  maximum shear mixing at Rig = 0    [m2/s] 
    974    rn_bvsqcon  = -0.01e-07 !  Brunt-Vaisala squared for maximum convection [1/s2] 
    975    rn_difcon   =  1.       !  maximum mixing in interior convection [m2/s] 
    976    nn_avb      =  0        !  horizontal averaged (=1) or not (=0) on avt and amv 
    977    nn_ave      =  1        !  constant (=0) or profile (=1) background on avt 
    978968/ 
    979969!----------------------------------------------------------------------- 
     
    12931283   rn_htrmax         =  200.0   ! max. depth of transition range 
    12941284/ 
     1285!----------------------------------------------------------------------- 
     1286&namwad  !   Wetting and drying 
     1287!----------------------------------------------------------------------- 
     1288   ln_wd             = .false.  ! T/F activation of wetting and drying 
     1289   rn_wdmin1         =  0.1     ! Minimum wet depth on dried cells 
     1290   rn_wdmin2         =  0.01    ! Tolerance of min wet depth on dried cells 
     1291   rn_wdld           =  20.0    ! Land elevation below which wetting/dryins is allowed 
     1292   nn_wdit           =  10      ! Max iterations for W/D limiter 
     1293/ 
Note: See TracChangeset for help on using the changeset viewer.