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/GYRE_PISCES/EXP00/namelist_cfg – 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/GYRE_PISCES/EXP00/namelist_cfg

    r5102 r5870  
    104104/ 
    105105!----------------------------------------------------------------------- 
    106 &namcla        !   cross land advection 
    107 !----------------------------------------------------------------------- 
    108 / 
    109 !----------------------------------------------------------------------- 
    110106&nambfr        !   bottom friction 
    111107!----------------------------------------------------------------------- 
    … …  
    143139&namtra_adv    !   advection scheme for tracer 
    144140!----------------------------------------------------------------------- 
    145    ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl  
     141   ln_traadv_fct =  .true.   !  FCT scheme 
     142      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     143      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     144      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     145      !                             !        (number of sub-timestep = nn_fct_zts) 
    146146/ 
    147147!---------------------------------------------------------------------------------- 
    148148&namtra_ldf    !   lateral diffusion scheme for tracers 
    149149!---------------------------------------------------------------------------------- 
    150    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    151    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     150   !                       !  Operator type: 
     151   ln_traldf_lap   =  .true.   !    laplacian operator 
     152   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     153   !                       !  Direction of action: 
     154   ln_traldf_lev   =  .false.  !  iso-level 
     155   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     156   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     157   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     158   ! 
     159   !                       !  iso-neutral options:         
     160   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     161   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     162   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     163   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     164   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     165   ! 
     166   !                       !  Coefficients: 
     167   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     168   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     169   !                                !   =  0           constant  
     170   !                                !   = 10 F(k)      =ldf_c1d  
     171   !                                !   = 20 F(i,j)    =ldf_c2d  
     172   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     173   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     174   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     175   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     176   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     177/ 
     178!---------------------------------------------------------------------------------- 
     179&namtra_ldfeiv !   eddy induced velocity param. 
     180!---------------------------------------------------------------------------------- 
     181   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    152182/ 
    153183!----------------------------------------------------------------------- 
    … …  
    163193&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    164194!----------------------------------------------------------------------- 
    165    ln_dynvor_ene = .true.  !  energy    conserving scheme   
    166    ln_dynvor_ens = .false. !  enstrophy conserving scheme     
     195   ln_dynvor_ene = .true.  !  enstrophy conserving scheme 
     196   ln_dynvor_ens = .false. !  energy conserving scheme 
     197   ln_dynvor_mix = .false. !  mixed scheme 
    167198   ln_dynvor_een = .false. !  energy & enstrophy scheme 
     199      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    168200/ 
    169201!----------------------------------------------------------------------- 
    … …  
    176208&namdyn_ldf    !   lateral diffusion on momentum 
    177209!----------------------------------------------------------------------- 
     210   !                       !  Type of the operator : 
     211   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     212   ln_dynldf_lap =  .true.    !    laplacian operator 
     213   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     214   !                       !  Direction of action  : 
     215   ln_dynldf_lev =  .true.    !  iso-level 
     216   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     217   ln_dynldf_iso =  .false.    !  iso-neutral 
     218   !                       !  Coefficient 
     219   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     220   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     221   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     222   !                                !  =  0  constant  
     223   !                                !  = 10  F(k)=c1d 
     224   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     225   !                                !  = 30  F(i,j,k)=c2d*c1d 
     226   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     227   rn_ahm_0      = 100000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     228   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     229   rn_bhm_0      =      0.      !  horizontal bilaplacian eddy viscosity [m4/s] 
     230   ! 
     231   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
     232/ 
    178233   rn_ahm_0_lap     = 100000.   !  horizontal laplacian eddy viscosity   [m2/s] 
    179234/ 
Note: See TracChangeset for help on using the changeset viewer.