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/ORCA2_LIM – NEMO

Ignore:
Timestamp:
2015-11-09T18:33:54+01:00 (9 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.

Location:
branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_LIM
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/1_namelist_cfg

    r5656 r5870  
    44!----------------------------------------------------------------------- 
    55&namrun        !   parameters of the run 
     6 nn_it000=1 
    67!----------------------------------------------------------------------- 
    78   cn_exp      = "Agulhas" !  experience name  
    8    nn_itend    =   10950   !  last  time step 
     9   nn_itend    =   10950 
    910   nn_stock    =   10950   !  frequency of creation of a restart file (modulo referenced to 1) 
    1011   nn_write    =   10950   !  frequency of write in the output file   (modulo referenced to nn_it000) 
    11    ln_clobber  = .true.    !  clobber (overwrite) an existing file 
     12   ln_clobber  =  .true. 
    1213/ 
    1314!----------------------------------------------------------------------- 
     
    116117/ 
    117118!----------------------------------------------------------------------- 
    118 &namcla        !   cross land advection 
    119 !----------------------------------------------------------------------- 
    120 / 
    121 !----------------------------------------------------------------------- 
    122119&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    123120!----------------------------------------------------------------------- 
     
    145142/ 
    146143!----------------------------------------------------------------------- 
    147 &namtra_adv    !   advection scheme for tracer  
    148 !----------------------------------------------------------------------- 
    149 / 
    150 !----------------------------------------------------------------------- 
    151 &namtra_ldf    !   lateral diffusion scheme for tracers  
    152 !----------------------------------------------------------------------- 
    153  
    154    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    155    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s]    (require "key_traldf_eiv") 
     144&namtra_adv    !   advection scheme for tracer 
     145!----------------------------------------------------------------------- 
     146   ln_traadv_fct =  .true.    !  FCT scheme 
     147      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order 
     148      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order 
     149      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     150      !                             !        (number of sub-timestep = nn_fct_zts) 
     151/ 
     152!----------------------------------------------------------------------- 
     153&namtra_ldf    !   lateral diffusion scheme for tracers 
     154!----------------------------------------------------------------------- 
     155   !                       !  Operator type: 
     156   ln_traldf_lap   =  .true.   !    laplacian operator 
     157   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     158   !                       !  Direction of action: 
     159   ln_traldf_lev   =  .false.  !  iso-level 
     160   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     161   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     162   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     163   ! 
     164   !                       !  iso-neutral options:         
     165   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     166   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     167   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     168   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     169   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     170   ! 
     171   !                       !  Coefficients: 
     172   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     173   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     174   !                                !   =  0           constant 
     175   !                                !   = 10 F(k)      =ldf_c1d 
     176   !                                !   = 20 F(i,j)    =ldf_c2d 
     177   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     178   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     179   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     180   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     181   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     182/ 
     183!---------------------------------------------------------------------------------- 
     184&namtra_ldfeiv !   eddy induced velocity param. 
     185!---------------------------------------------------------------------------------- 
     186   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    156187/ 
    157188!----------------------------------------------------------------------- 
     
    162193&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    163194!----------------------------------------------------------------------- 
     195   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     196   ln_dynvor_ens = .false. !  energy conserving scheme 
     197   ln_dynvor_mix = .false. !  mixed scheme 
     198   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     199      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
     200   ln_dynvor_msk = .true.  !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
    164201/ 
    165202!----------------------------------------------------------------------- 
     
    170207&namdyn_ldf    !   lateral diffusion on momentum 
    171208!----------------------------------------------------------------------- 
    172    !                       !  Type of the operator :  
    173    ln_dynldf_lap    =  .false.  !  laplacian operator          
    174    ln_dynldf_bilap  =  .true.   !  bilaplacian operator     
    175    rn_ahm_0_blp     = -8.5e+11  !  horizontal bilaplacian eddy viscosity [m4/s]  
     209   !                       !  Type of the operator : 
     210   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F 
     211   ln_dynldf_lap =  .false.    !    laplacian operator 
     212   ln_dynldf_blp =  .true.     !  bilaplacian operator 
     213   !                       !  Direction of action  : 
     214   ln_dynldf_lev =  .true.     !  iso-level 
     215   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     216   ln_dynldf_iso =  .false.    !  iso-neutral 
     217   !                       !  Coefficient 
     218   nn_ahm_ijk_t  =   0         !  space/time variation of eddy coef 
     219   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     220   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     221   !                                !  =  0  constant 
     222   !                                !  = 10  F(k)=c1d 
     223   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     224   !                                !  = 30  F(i,j,k)=c2d*c1d 
     225   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     226   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     227   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     228   rn_bhm_0      = 8.5e+11      !  horizontal bilaplacian eddy viscosity [m4/s] 
    176229/ 
    177230!----------------------------------------------------------------------- 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist_cfg

    r4990 r5870  
    9393/ 
    9494!----------------------------------------------------------------------- 
    95 &namcla        !   cross land advection 
    96 !----------------------------------------------------------------------- 
    97 / 
    98 !----------------------------------------------------------------------- 
    9995&nambfr        !   bottom friction 
    10096!----------------------------------------------------------------------- 
     
    115111&namtra_adv    !   advection scheme for tracer 
    116112!----------------------------------------------------------------------- 
     113   ln_traadv_fct =  .true.    !  FCT scheme 
     114      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     115      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     116      nn_fct_zts =  0               !  > 1 , 2nd order FCT scheme with vertical sub-timestepping 
     117      !                             !        (number of sub-timestep = nn_fct_zts) 
    117118/ 
    118119!----------------------------------------------------------------------- 
     
    120121!----------------------------------------------------------------------- 
    121122/ 
    122 !----------------------------------------------------------------------- 
    123 &namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
    124 !----------------------------------------------------------------------- 
    125 / 
    126123!---------------------------------------------------------------------------------- 
    127124&namtra_ldf    !   lateral diffusion scheme for tracers 
    128125!---------------------------------------------------------------------------------- 
     126   !                       !  Operator type: 
     127   ln_traldf_lap   =  .true.   !    laplacian operator 
     128   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     129   !                       !  Direction of action: 
     130   ln_traldf_lev   =  .false.  !  iso-level 
     131   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     132   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     133   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     134   ! 
     135   !                       !  iso-neutral options:         
     136   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     137   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     138   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     139   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     140   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     141   ! 
     142   !                       !  Coefficients: 
     143   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     144   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     145   !                                !   =  0           constant 
     146   !                                !   = 10 F(k)      =ldf_c1d 
     147   !                                !   = 20 F(i,j)    =ldf_c2d 
     148   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     149   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     150   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     151   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     152   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     153/ 
     154!---------------------------------------------------------------------------------- 
     155&namtra_ldfeiv !   eddy induced velocity param. 
     156!---------------------------------------------------------------------------------- 
     157   ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
     158   ln_ldfeiv_dia =.true.   ! diagnose eiv stream function and velocities 
     159   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
     160   nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
     161   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     162   !                                !   =  0           constant 
     163   !                                !   = 10 F(k)      =ldf_c1d 
     164   !                                !   = 20 F(i,j)    =ldf_c2d 
     165   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     166   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    129167/ 
    130168!----------------------------------------------------------------------- 
     
    139177&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    140178!----------------------------------------------------------------------- 
     179   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     180   ln_dynvor_ens = .false. !  energy conserving scheme 
     181   ln_dynvor_mix = .false. !  mixed scheme 
     182   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     183      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    141184/ 
    142185!----------------------------------------------------------------------- 
     
    147190&namdyn_ldf    !   lateral diffusion on momentum 
    148191!----------------------------------------------------------------------- 
     192   !                       !  Type of the operator : 
     193   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     194   ln_dynldf_lap =  .true.     !    laplacian operator 
     195   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     196   !                       !  Direction of action  : 
     197   ln_dynldf_lev =  .true.     !  iso-level 
     198   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     199   ln_dynldf_iso =  .false.    !  iso-neutral 
     200   !                       !  Coefficient 
     201   nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
     202   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     203   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     204   !                                !  =  0  constant  
     205   !                                !  = 10  F(k)=c1d 
     206   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     207   !                                !  = 30  F(i,j,k)=c2d*c1d 
     208   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     209   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     210   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     211   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    149212/ 
    150213!----------------------------------------------------------------------- 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_LIM/cpp_ORCA2_LIM.fcm

    r5385 r5870  
    1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi key_diaobs key_asminc  
     1bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_zdftke key_zdfddm key_zdftmx  key_mpp_mpi key_iomput 
Note: See TracChangeset for help on using the changeset viewer.