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 5836 for trunk/NEMOGCM/CONFIG – NEMO

Changeset 5836 for trunk/NEMOGCM/CONFIG


Ignore:
Timestamp:
2015-10-26T15:49:40+01:00 (8 years ago)
Author:
cetlod
Message:

merge the simplification branch onto the trunk, see ticket #1612

Location:
trunk/NEMOGCM/CONFIG
Files:
36 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/CONFIG/AMM12/EXP00/namelist_cfg

    r5501 r5836  
    196196/ 
    197197!----------------------------------------------------------------------- 
    198 &namcla        !   cross land advection 
    199 !----------------------------------------------------------------------- 
    200 / 
    201 !----------------------------------------------------------------------- 
    202 &namobc        !   open boundaries parameters                           ("key_obc") 
    203 !----------------------------------------------------------------------- 
    204 / 
    205 !----------------------------------------------------------------------- 
    206198&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    207199!----------------------------------------------------------------------- 
     
    281273&namtra_adv    !   advection scheme for tracer 
    282274!----------------------------------------------------------------------- 
     275   ln_traadv_fct =  .true.   !  FCT scheme 
     276      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     277      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     278      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     279      !                             !        (number of sub-timestep = nn_fct_zts) 
    283280/ 
    284281!----------------------------------------------------------------------- 
     
    289286&namtra_ldf    !   lateral diffusion scheme for tracers 
    290287!---------------------------------------------------------------------------------- 
    291    ln_traldf_hor    =  .true.   !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T) 
    292    ln_traldf_iso    =  .false.  !  iso-neutral                 (needs "key_ldfslp") 
    293    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    294    rn_aht_0         =    50.    !  horizontal eddy diffusivity for tracers [m2/s] 
     288   !                       !  Operator type: 
     289   ln_traldf_lap   =  .true.   !    laplacian operator 
     290   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     291   !                       !  Direction of action: 
     292   ln_traldf_lev   =  .false.  !  iso-level 
     293   ln_traldf_hor   =  .true.   !  horizontal (geopotential) 
     294   ln_traldf_iso   =  .false.  !  iso-neutral 
     295   ln_traldf_triad =  .false.  !  iso-neutral using Griffies triads 
     296   ! 
     297   !                       !  iso-neutral options:         
     298   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     299   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     300   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML     (triad only) 
     301   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     302   ln_botmix_triad =  .false.  !  lateral mixing on bottom         (triad only) 
     303   ! 
     304   !                       !  Coefficients: 
     305   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     306   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     307   !                                !   =  0           constant  
     308   !                                !   = 10 F(k)      =ldf_c1d  
     309   !                                !   = 20 F(i,j)    =ldf_c2d  
     310   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     311   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     312   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     313   rn_aht_0        =   50.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     314   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     315/ 
     316!---------------------------------------------------------------------------------- 
     317&namtra_ldfeiv !   eddy induced velocity param. 
     318!---------------------------------------------------------------------------------- 
     319   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    295320/ 
    296321!----------------------------------------------------------------------- 
     
    306331&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    307332!----------------------------------------------------------------------- 
     333   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     334   ln_dynvor_ens = .false. !  energy conserving scheme 
     335   ln_dynvor_mix = .false. !  mixed scheme 
     336   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     337      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    308338/ 
    309339!----------------------------------------------------------------------- 
     
    324354!----------------------------------------------------------------------- 
    325355   !                       !  Type of the operator : 
    326    ln_dynldf_bilap  =  .true.   !  bilaplacian operator 
    327    ln_dynldf_lap    =  .false.  !  bilaplacian operator 
     356   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     357   ln_dynldf_lap =  .false.    !    laplacian operator 
     358   ln_dynldf_blp =  .true.     !  bilaplacian operator 
    328359   !                       !  Direction of action  : 
    329    ln_dynldf_level  =  .true.   !  iso-level 
    330    ln_dynldf_hor    =  .false.  !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.) 
    331                            !  Coefficient 
    332    rn_ahm_0_lap     = 60.0      !  horizontal laplacian eddy viscosity   [m2/s] 
    333    rn_ahm_0_blp     = -1.0e+10  !  horizontal bilaplacian eddy viscosity [m4/s] 
     360   ln_dynldf_lev =  .true.     !  iso-level 
     361   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     362   ln_dynldf_iso =  .false.    !  iso-neutral 
     363   !                       !  Coefficient 
     364   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     365   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     366   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     367   !                                !  =  0  constant  
     368   !                                !  = 10  F(k)=c1d 
     369   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     370   !                                !  = 30  F(i,j,k)=c2d*c1d 
     371   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     372   rn_ahm_0      =     60.     !  horizontal laplacian eddy viscosity   [m2/s] 
     373   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     374   rn_bhm_0      = 1.0e+10     !  horizontal bilaplacian eddy viscosity [m4/s] 
     375   ! 
     376   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
    334377/ 
    335378!----------------------------------------------------------------------- 
     
    349392!----------------------------------------------------------------------- 
    350393/ 
    351 !------------------------------------------------------------------------ 
    352 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    353 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    354 / 
    355394!----------------------------------------------------------------------- 
    356395&namzdf_gls                !   GLS vertical diffusion                   ("key_zdfgls") 
  • trunk/NEMOGCM/CONFIG/AMM12/cpp_AMM12.fcm

    r4245 r5836  
    1  bld::tool::fppkeys  key_bdy key_tide key_dynspg_ts key_ldfslp key_zdfgls  key_vvl key_diainstant key_mpp_mpi key_iomput 
     1 bld::tool::fppkeys  key_bdy key_tide key_dynspg_ts key_zdfgls  key_vvl key_diainstant key_mpp_mpi key_iomput 
  • trunk/NEMOGCM/CONFIG/C1D_PAPA/EXP00/namelist_cfg

    r5407 r5836  
    172172/ 
    173173!----------------------------------------------------------------------- 
    174 &namcla        !   cross land advection 
    175 !----------------------------------------------------------------------- 
    176 / 
    177 !----------------------------------------------------------------------- 
    178 &namobc        !   open boundaries parameters                           ("key_obc") 
    179 !----------------------------------------------------------------------- 
    180 / 
    181 !----------------------------------------------------------------------- 
    182174&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    183175!----------------------------------------------------------------------- 
     
    226218&namtra_adv    !   advection scheme for tracer 
    227219!----------------------------------------------------------------------- 
     220! C1D : no advection scheme  
    228221/ 
    229222!----------------------------------------------------------------------- 
     
    231224!----------------------------------------------------------------------- 
    232225/ 
    233 !---------------------------------------------------------------------------------- 
     226!----------------------------------------------------------------------- 
    234227&namtra_ldf    !   lateral diffusion scheme for tracers 
    235 !---------------------------------------------------------------------------------- 
    236 !---------------------------------------------------------------------------------- 
    237    ln_traldf_hor    =  .true.   !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T) 
    238    ln_traldf_iso    =  .false.  !  iso-neutral                 (needs "key_ldfslp") 
    239    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    240    rn_aht_0         =     0.    !  horizontal eddy diffusivity for tracers [m2/s] 
     228!----------------------------------------------------------------------- 
     229! C1D : no lateral diffusion   
     230/ 
     231!----------------------------------------------------------------------- 
     232&namtra_ldfeiv !   eddy induced velocity param. 
     233!----------------------------------------------------------------------- 
     234! C1D : no eiv   
    241235/ 
    242236!----------------------------------------------------------------------- 
     
    248242&namdyn_adv    !   formulation of the momentum advection 
    249243!----------------------------------------------------------------------- 
     244! C1D : no advection scheme  
    250245/ 
    251246!----------------------------------------------------------------------- 
     
    285280&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
    286281!----------------------------------------------------------------------- 
    287 / 
    288 !------------------------------------------------------------------------ 
    289 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    290 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    291282/ 
    292283!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/GYRE/EXP00/namelist_cfg

    r5407 r5836  
    156156/ 
    157157!----------------------------------------------------------------------- 
    158 &namcla        !   cross land advection 
    159 !----------------------------------------------------------------------- 
    160 / 
    161 !----------------------------------------------------------------------- 
    162 &namobc        !   open boundaries parameters                           ("key_obc") 
    163 !----------------------------------------------------------------------- 
    164 / 
    165 !----------------------------------------------------------------------- 
    166158&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    167159!----------------------------------------------------------------------- 
     
    223215&namtra_adv    !   advection scheme for tracer 
    224216!----------------------------------------------------------------------- 
    225    ln_traadv_cen2   =  .false.   !  2nd order centered scheme 
    226    ln_traadv_tvd    =  .true.    !  TVD scheme 
    227    ln_traadv_muscl  =  .false.   !  MUSCL scheme 
    228    ln_traadv_muscl2 =  .false.   !  MUSCL2 scheme + cen2 at boundaries 
    229    ln_traadv_ubs    =  .false.   !  UBS scheme 
    230    ln_traadv_qck    =  .false.   !  QUICKEST scheme 
    231    ln_traadv_msc_ups=  .false.   !  use upstream scheme within muscl 
     217   ln_traadv_fct =  .true.   !  FCT scheme 
     218      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     219      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     220      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     221      !                             !        (number of sub-timestep = nn_fct_zts) 
    232222/ 
    233223!----------------------------------------------------------------------- 
     
    238228&namtra_ldf    !   lateral diffusion scheme for tracers 
    239229!---------------------------------------------------------------------------------- 
    240    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    241    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     230   !                       !  Operator type: 
     231   ln_traldf_lap   =  .true.   !    laplacian operator 
     232   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     233   !                       !  Direction of action: 
     234   ln_traldf_lev   =  .false.  !  iso-level 
     235   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     236   ln_traldf_iso   =  .true.   !  iso-neutral 
     237   ln_traldf_triad =  .false.  !  iso-neutral using Griffies triads 
     238   ! 
     239   !                       !  iso-neutral options:         
     240   ln_traldf_msc   =  .false.  !  Method of Stabilizing Correction (both operators) 
     241   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     242   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML     (triad only) 
     243   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     244   ln_botmix_triad =  .false.  !  lateral mixing on bottom         (triad only) 
     245   ! 
     246   !                       !  Coefficients: 
     247   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     248   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     249   !                                !   =  0           constant  
     250   !                                !   = 10 F(k)      =ldf_c1d  
     251   !                                !   = 20 F(i,j)    =ldf_c2d  
     252   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     253   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     254   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     255   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     256   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     257/ 
     258!---------------------------------------------------------------------------------- 
     259&namtra_ldfeiv !   eddy induced velocity param. 
     260!---------------------------------------------------------------------------------- 
     261   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    242262/ 
    243263!----------------------------------------------------------------------- 
     
    253273&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    254274!----------------------------------------------------------------------- 
    255    ln_dynvor_ene = .true.  !  energy    conserving scheme   
    256    ln_dynvor_ens = .false. !  enstrophy conserving scheme     
     275   ln_dynvor_ene = .true.  !  enstrophy conserving scheme 
     276   ln_dynvor_ens = .false. !  energy conserving scheme 
     277   ln_dynvor_mix = .false. !  mixed scheme 
    257278   ln_dynvor_een = .false. !  energy & enstrophy scheme 
     279      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    258280/ 
    259281!----------------------------------------------------------------------- 
     
    270292&namdyn_ldf    !   lateral diffusion on momentum 
    271293!----------------------------------------------------------------------- 
    272    rn_ahm_0_lap     = 100000.   !  horizontal laplacian eddy viscosity   [m2/s] 
     294   !                       !  Type of the operator : 
     295   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     296   ln_dynldf_lap =  .true.     !    laplacian operator 
     297   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     298   !                       !  Direction of action  : 
     299   ln_dynldf_lev =  .true.     !  iso-level 
     300   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     301   ln_dynldf_iso =  .false.    !  iso-neutral 
     302   !                       !  Coefficient 
     303   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     304   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     305   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     306   !                                !  =  0  constant  
     307   !                                !  = 10  F(k)=c1d 
     308   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     309   !                                !  = 30  F(i,j,k)=c2d*c1d 
     310   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     311   rn_ahm_0      = 100000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     312   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     313   rn_bhm_0      =      0.     !  horizontal bilaplacian eddy viscosity [m4/s] 
    273314/ 
    274315!----------------------------------------------------------------------- 
     
    285326!----------------------------------------------------------------------- 
    286327   nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves 
    287 / 
    288 !------------------------------------------------------------------------ 
    289 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    290 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    291328/ 
    292329!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/GYRE/cpp_GYRE.fcm

    r4990 r5836  
    1  bld::tool::fppkeys key_dynspg_flt key_ldfslp key_zdftke key_iomput key_mpp_mpi key_nosignedzero 
     1 bld::tool::fppkeys key_dynspg_flt key_zdftke key_iomput key_mpp_mpi 
  • trunk/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist_cfg

    r5407 r5836  
    161161/ 
    162162!----------------------------------------------------------------------- 
    163 &namcla        !   cross land advection 
    164 !----------------------------------------------------------------------- 
    165 / 
    166 !----------------------------------------------------------------------- 
    167 &namobc        !   open boundaries parameters                           ("key_obc") 
    168 !----------------------------------------------------------------------- 
    169 / 
    170 !----------------------------------------------------------------------- 
    171163&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    172164!----------------------------------------------------------------------- 
     
    228220&namtra_adv    !   advection scheme for tracer 
    229221!----------------------------------------------------------------------- 
    230    ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl  
     222   ln_traadv_fct =  .true.   !  FCT scheme 
     223      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     224      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     225      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     226      !                             !        (number of sub-timestep = nn_fct_zts) 
    231227/ 
    232228!----------------------------------------------------------------------- 
     
    237233&namtra_ldf    !   lateral diffusion scheme for tracers 
    238234!---------------------------------------------------------------------------------- 
    239    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    240    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     235   !                       !  Operator type: 
     236   ln_traldf_lap   =  .true.   !    laplacian operator 
     237   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     238   !                       !  Direction of action: 
     239   ln_traldf_lev   =  .false.  !  iso-level 
     240   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     241   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     242   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     243   ! 
     244   !                       !  iso-neutral options:         
     245   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     246   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     247   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     248   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     249   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     250   ! 
     251   !                       !  Coefficients: 
     252   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     253   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     254   !                                !   =  0           constant  
     255   !                                !   = 10 F(k)      =ldf_c1d  
     256   !                                !   = 20 F(i,j)    =ldf_c2d  
     257   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     258   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     259   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     260   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     261   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     262/ 
     263!---------------------------------------------------------------------------------- 
     264&namtra_ldfeiv !   eddy induced velocity param. 
     265!---------------------------------------------------------------------------------- 
     266   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    241267/ 
    242268!----------------------------------------------------------------------- 
     
    252278&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    253279!----------------------------------------------------------------------- 
    254    ln_dynvor_ene = .true.  !  energy    conserving scheme   
    255    ln_dynvor_ens = .false. !  enstrophy conserving scheme     
    256    ln_dynvor_een = .false. !  energy & enstrophy scheme 
     280   ln_dynvor_ene = .true.  !  enstrophy conserving scheme 
    257281/ 
    258282!----------------------------------------------------------------------- 
     
    268292&namdyn_ldf    !   lateral diffusion on momentum 
    269293!----------------------------------------------------------------------- 
     294   !                       !  Type of the operator : 
     295   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     296   ln_dynldf_lap =  .true.     !    laplacian operator 
     297   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     298   !                       !  Direction of action  : 
     299   ln_dynldf_lev =  .true.     !  iso-level 
     300   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     301   ln_dynldf_iso =  .false.    !  iso-neutral 
     302   !                       !  Coefficient 
     303   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     304   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     305   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     306   !                                !  =  0  constant  
     307   !                                !  = 10  F(k)=c1d 
     308   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     309   !                                !  = 30  F(i,j,k)=c2d*c1d 
     310   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     311   rn_ahm_0      = 100000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     312   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     313   rn_bhm_0      =      0.     !  horizontal bilaplacian eddy viscosity [m4/s] 
     314   ! 
     315   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
     316/ 
    270317   rn_ahm_0_lap     = 100000.   !  horizontal laplacian eddy viscosity   [m2/s] 
    271318/ 
     
    283330!----------------------------------------------------------------------- 
    284331   nn_etau     =   0       !  penetration of tke below the mixed layer (ML) due to internal & intertial waves 
    285 / 
    286 !------------------------------------------------------------------------ 
    287 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    288 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    289332/ 
    290333!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist_top_cfg

    r4152 r5836  
    2323!----------------------------------------------------------------------- 
    2424&namtrc_adv    !   advection scheme for passive tracer  
    25 !-----------------------------------------------------------------------   
    26    ln_trcadv_tvd     =  .true.   !  TVD scheme 
    27    ln_trcadv_muscl   =  .false.  !  MUSCL scheme 
     25!----------------------------------------------------------------------- 
     26   ln_trcadv_fct =  .true.   !  FCT scheme 
     27      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     28      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     29      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     30      !                             !        (number of sub-timestep = nn_fct_zts) 
    2831/ 
    2932!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/GYRE_BFM/cpp_GYRE_BFM.fcm

    r4230 r5836  
    1 bld::tool::fppkeys key_dynspg_flt key_ldfslp key_zdftke key_vectopt_loop key_top key_my_trc key_mpp_mpi key_iomput 
     1bld::tool::fppkeys key_dynspg_flt key_zdftke key_top key_my_trc key_mpp_mpi key_iomput 
    22inc $BFMDIR/src/nemo/bfm.fcm 
  • trunk/NEMOGCM/CONFIG/GYRE_PISCES/EXP00/namelist_cfg

    r5102 r5836  
    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/ 
  • trunk/NEMOGCM/CONFIG/GYRE_PISCES/EXP00/namelist_top_cfg

    r4340 r5836  
    2424&namtrc_adv    !   advection scheme for passive tracer  
    2525!-----------------------------------------------------------------------    
    26    ln_trcadv_tvd     =  .true.   !  TVD scheme 
    27    ln_trcadv_muscl   =  .false.  !  MUSCL scheme 
     26   ln_trcadv_fct =  .true.   !  FCT scheme 
     27      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     28      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     29      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     30      !                             !        (number of sub-timestep = nn_fct_zts) 
    2831/ 
    2932!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/GYRE_PISCES/cpp_GYRE_PISCES.fcm

    r4990 r5836  
    1 bld::tool::fppkeys  key_dynspg_flt key_ldfslp key_zdftke key_top key_pisces_reduced key_iomput key_mpp_mpi 
     1bld::tool::fppkeys  key_dynspg_flt key_zdftke key_top key_pisces_reduced key_mpp_mpi 
  • trunk/NEMOGCM/CONFIG/GYRE_XIOS/EXP00/namelist_cfg

    r5407 r5836  
    150150/ 
    151151!----------------------------------------------------------------------- 
    152 &namcla        !   cross land advection 
    153 !----------------------------------------------------------------------- 
    154 / 
    155 !----------------------------------------------------------------------- 
    156 &namobc        !   open boundaries parameters                           ("key_obc") 
    157 !----------------------------------------------------------------------- 
    158 / 
    159 !----------------------------------------------------------------------- 
    160152&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    161153!----------------------------------------------------------------------- 
     
    200192&namtra_adv    !   advection scheme for tracer 
    201193!----------------------------------------------------------------------- 
    202    ln_traadv_msc_ups=  .false.  !  use upstream scheme within muscl  
     194   ln_traadv_fct =  .true.    !  FCT scheme 
     195      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     196      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     197      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     198      !                             !        (number of sub-timestep = nn_fct_zts) 
    203199/ 
    204200!----------------------------------------------------------------------- 
     
    209205&namtra_ldf    !   lateral diffusion scheme for tracers 
    210206!---------------------------------------------------------------------------------- 
    211    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    212    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     207   !                       !  Operator type: 
     208   ln_traldf_lap   =  .true.   !    laplacian operator 
     209   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     210   !                       !  Direction of action: 
     211   ln_traldf_lev   =  .false.  !  iso-level 
     212   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     213   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     214   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     215   ! 
     216   !                       !  iso-neutral options:         
     217   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     218   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     219   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     220   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     221   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     222   ! 
     223   !                       !  Coefficients: 
     224   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     225   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     226   !                                !   =  0           constant  
     227   !                                !   = 10 F(k)      =ldf_c1d  
     228   !                                !   = 20 F(i,j)    =ldf_c2d  
     229   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     230   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     231   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     232   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     233   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     234/ 
     235!---------------------------------------------------------------------------------- 
     236&namtra_ldfeiv !   eddy induced velocity param. 
     237!---------------------------------------------------------------------------------- 
     238   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    213239/ 
    214240!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/GYRE_XIOS/cpp_GYRE_XIOS.fcm

    r4373 r5836  
    1  bld::tool::fppkeys key_dynspg_flt key_ldfslp key_zdftke key_iomput key_mpp_mpi 
     1 bld::tool::fppkeys key_dynspg_flt key_zdftke key_iomput key_mpp_mpi 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/1_namelist_cfg

    r5656 r5836  
    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!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist_cfg

    r4990 r5836  
    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!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM/cpp_ORCA2_LIM.fcm

    r5385 r5836  
    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 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/1_namelist_cfg

    r5499 r5836  
    116116/ 
    117117!----------------------------------------------------------------------- 
    118 &namcla        !   cross land advection 
    119 !----------------------------------------------------------------------- 
    120 / 
    121 !----------------------------------------------------------------------- 
    122118&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    123119!----------------------------------------------------------------------- 
     
    145141/ 
    146142!----------------------------------------------------------------------- 
    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") 
     143&namtra_adv    !   advection scheme for tracer 
     144!----------------------------------------------------------------------- 
     145   ln_traadv_fct =  .true.    !  FCT scheme 
     146      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     147      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     148      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     149      !                             !        (number of sub-timestep = nn_fct_zts) 
     150/ 
     151!----------------------------------------------------------------------- 
     152&namtra_ldf    !   lateral diffusion scheme for tracers 
     153!---------------------------------------------------------------------------------- 
     154   !                       !  Operator type: 
     155   ln_traldf_lap   =  .true.   !    laplacian operator 
     156   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     157   !                       !  Direction of action: 
     158   ln_traldf_lev   =  .false.  !  iso-level 
     159   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     160   ln_traldf_iso   =  .true.   !  iso-neutral (Standard operator) 
     161   ln_traldf_triad =  .false.  !  iso-neutral (Triads   operator) 
     162   ! 
     163   !                       !  iso-neutral options:         
     164   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     165   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     166   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     167   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     168   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     169   ! 
     170   !                       !  Coefficients: 
     171   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     172   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     173   !                                !   =  0           constant  
     174   !                                !   = 10 F(k)      =ldf_c1d  
     175   !                                !   = 20 F(i,j)    =ldf_c2d  
     176   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     177   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     178   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     179   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     180   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     181/ 
     182!---------------------------------------------------------------------------------- 
     183&namtra_ldfeiv !   eddy induced velocity param. 
     184!---------------------------------------------------------------------------------- 
     185   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    156186/ 
    157187!----------------------------------------------------------------------- 
     
    162192&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    163193!----------------------------------------------------------------------- 
     194   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     195   ln_dynvor_ens = .false. !  energy conserving scheme 
     196   ln_dynvor_mix = .false. !  mixed scheme 
     197   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     198      nn_een_e3f = 0             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
     199   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
    164200/ 
    165201!----------------------------------------------------------------------- 
     
    170206&namdyn_ldf    !   lateral diffusion on momentum 
    171207!----------------------------------------------------------------------- 
    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]  
     208   !                       !  Type of the operator : 
     209   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     210   ln_dynldf_lap =  .false.    !    laplacian operator 
     211   ln_dynldf_blp =  .true.     !  bilaplacian operator 
     212   !                       !  Direction of action  : 
     213   ln_dynldf_lev =  .true.     !  iso-level 
     214   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     215   ln_dynldf_iso =  .false.    !  iso-neutral 
     216   !                       !  Coefficient 
     217   nn_ahm_ijk_t  = 20          !  space/time variation of eddy coef 
     218   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     219   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     220   !                                !  =  0  constant  
     221   !                                !  = 10  F(k)=c1d 
     222   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     223   !                                !  = 30  F(i,j,k)=c2d*c1d 
     224   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     225   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     226   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     227   rn_bhm_0      = 8.5e+11     !  horizontal bilaplacian eddy viscosity [m4/s] 
    176228/ 
    177229!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/namelist_cfg

    r4995 r5836  
    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!----------------------------------------------------------------------- 
     
    123124&namtra_ldf    !   lateral diffusion scheme for tracers 
    124125!---------------------------------------------------------------------------------- 
    125 / 
    126 !----------------------------------------------------------------------- 
     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 (Triads   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 
     167/!----------------------------------------------------------------------- 
    127168&namtra_dmp    !   tracer: T & S newtonian damping 
    128169!----------------------------------------------------------------------- 
     
    135176&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    136177!----------------------------------------------------------------------- 
     178   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     179   ln_dynvor_ens = .false. !  energy conserving scheme 
     180   ln_dynvor_mix = .false. !  mixed scheme 
     181   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     182      nn_een_e3f = 0             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    137183/ 
    138184!----------------------------------------------------------------------- 
     
    140186!----------------------------------------------------------------------- 
    141187   ln_hpg_zco  = .false.   !  z-coordinate - full steps 
    142    ln_hpg_zps  = .false.   !  z-coordinate - partial steps (interpolation) 
     188   ln_hpg_zps  = .false.    !  z-coordinate - partial steps (interpolation) 
    143189   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation) 
    144190   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial) 
    145    ln_hpg_prj  = .true.    !  s-coordinate (Pressure Jacobian scheme) 
     191   ln_hpg_prj  = .true.   !  s-coordinate (Pressure Jacobian scheme) 
    146192   ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T) 
    147193                                 !           centered      time scheme  (F) 
     
    150196&namdyn_ldf    !   lateral diffusion on momentum 
    151197!----------------------------------------------------------------------- 
     198   !                       !  Type of the operator : 
     199   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     200   ln_dynldf_lap =  .true.     !    laplacian operator 
     201   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     202   !                       !  Direction of action  : 
     203   ln_dynldf_lev =  .true.     !  iso-level 
     204   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     205   ln_dynldf_iso =  .false.    !  iso-neutral 
     206   !                       !  Coefficient 
     207   nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
     208   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     209   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     210   !                                !  =  0  constant  
     211   !                                !  = 10  F(k)=c1d 
     212   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     213   !                                !  = 30  F(i,j,k)=c2d*c1d 
     214   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     215   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     216   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     217   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
     218   ! 
     219   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
    152220/ 
    153221!----------------------------------------------------------------------- 
     
    191259!----------------------------------------------------------------------- 
    192260/ 
     261!----------------------------------------------------------------------- 
     262&namobs       !  observation usage                                      ('key_diaobs') 
     263!----------------------------------------------------------------------- 
     264/ 
     265!----------------------------------------------------------------------- 
     266&nam_asminc   !   assimilation increments                               ('key_asminc') 
     267!----------------------------------------------------------------------- 
     268/ 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM3/cpp_ORCA2_LIM3.fcm

    r5385 r5836  
    1  bld::tool::fppkeys key_trabbl key_lim3 key_vvl key_dynspg_ts 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  
     1 bld::tool::fppkeys key_trabbl key_lim3 key_vvl key_dynspg_ts key_zdftke key_zdfddm key_zdftmx  key_mpp_mpi key_diaobs key_asminc key_iomput key_nosignedzero 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/1_namelist_cfg

    r5407 r5836  
    55!!                                    namsbc_cpl, namtra_qsr, namsbc_rnf,  
    66!!                                    namsbc_apr, namsbc_ssr, namsbc_alb) 
    7 !!              4 - lateral boundary (namlbc, namcla, namobc, namagrif, nambdy, nambdy_tide) 
     7!!              4 - lateral boundary (namlbc, namagrif, nambdy, nambdy_tide) 
    88!!              5 - bottom  boundary (nambfr, nambbc, nambbl) 
    99!!              6 - Tracer           (nameos, namtra_adv, namtra_ldf, namtra_dmp) 
     
    302302!!====================================================================== 
    303303!!   namlbc        lateral momentum boundary condition 
    304 !!   namcla        cross land advection 
    305 !!   namobc        open boundaries parameters                           ("key_obc") 
    306304!!   namagrif      agrif nested grid ( read by child model only )       ("key_agrif")  
    307305!!   nambdy        Unstructured open boundaries                         ("key_bdy") 
     
    314312   rn_shlat    =    2.     !  shlat = 0  !  0 < shlat < 2  !  shlat = 2  !  2 < shlat 
    315313                           !  free slip  !   partial slip  !   no slip   ! strong slip 
    316 / 
    317 !----------------------------------------------------------------------- 
    318 &namcla        !   cross land advection 
    319 !----------------------------------------------------------------------- 
    320    nn_cla      =    0      !  advection between 2 ocean pts separates by land 
    321 / 
    322 !----------------------------------------------------------------------- 
    323 &namobc        !   open boundaries parameters                           ("key_obc") 
    324 !----------------------------------------------------------------------- 
    325    ln_obc_clim = .false.   !  climatological obc data files (T) or not (F) 
    326    ln_vol_cst  = .true.    !  impose the total volume conservation (T) or not (F) 
    327    ln_obc_fla  = .false.   !  Flather open boundary condition  
    328    nn_obcdta   =    1      !  = 0 the obc data are equal to the initial state 
    329                            !  = 1 the obc data are read in 'obc.dta' files 
    330    cn_obcdta   = 'annual'  !  set to annual if obc datafile hold 1 year of data 
    331                            !  set to monthly if obc datafile hold 1 month of data 
    332    rn_dpein    =    1.     !  damping time scale for inflow at east  open boundary 
    333    rn_dpwin    =    1.     !     -           -         -       west    -      - 
    334    rn_dpnin    =    1.     !     -           -         -       north   -      - 
    335    rn_dpsin    =    1.     !     -           -         -       south   -      - 
    336    rn_dpeob    = 3000.     !  time relaxation (days) for the east  open boundary 
    337    rn_dpwob    =   15.     !     -           -         -     west    -      - 
    338    rn_dpnob    = 3000.     !     -           -         -     north   -      - 
    339    rn_dpsob    =   15.     !     -           -         -     south   -      - 
    340    rn_volemp   =    1.     !  = 0 the total volume change with the surface flux (E-P-R) 
    341                            !  = 1 the total volume remains constant 
    342314/ 
    343315!----------------------------------------------------------------------- 
     
    394366&nambdy_tide     ! tidal forcing at open boundaries               
    395367!----------------------------------------------------------------------- 
    396    filtide      = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files 
     368    filtide       = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files 
    397369    tide_cpt(1)   ='Q1'  !  names of tidal components used 
    398370    tide_cpt(2)   ='O1'  !  names of tidal components used 
     
    485457/ 
    486458!----------------------------------------------------------------------- 
    487 &namtra_adv    !   advection scheme for tracer  
    488 !----------------------------------------------------------------------- 
    489    ln_traadv_cen2   =  .false.  !  2nd order centered scheme    
    490    ln_traadv_tvd    =  .true.   !  TVD scheme                 
    491    ln_traadv_muscl  =  .false.  !  MUSCL scheme              
    492    ln_traadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries   
    493    ln_traadv_ubs    =  .false.  !  UBS scheme                  
    494    ln_traadv_qck    =  .false.  !  QUCIKEST scheme                  
     459&namtra_adv    !   advection scheme for tracer 
     460!----------------------------------------------------------------------- 
     461   ln_traadv_fct =  .true.    !  FCT scheme 
     462      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     463      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     464      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     465      !                             !        (number of sub-timestep = nn_fct_zts) 
    495466/ 
    496467!----------------------------------------------------------------------- 
    497468&namtra_ldf    !   lateral diffusion scheme for tracer  
    498469!----------------------------------------------------------------------- 
    499    !                       !  Type of the operator :  
    500    ln_traldf_lap    =  .true.   !  laplacian operator        
    501    ln_traldf_bilap  =  .false.  !  bilaplacian operator      
    502    !                       !  Direction of action  : 
    503    ln_traldf_level  =  .false.  !  iso-level                 
    504    ln_traldf_hor    =  .false.  !  horizontal (geopotential)            (require "key_ldfslp" when ln_sco=T) 
    505    ln_traldf_iso    =  .true.   !  iso-neutral                          (require "key_ldfslp") 
    506    ln_traldf_grif   =  .false.  !  griffies skew flux formulation       (require "key_ldfslp")  ! UNDER TEST, DO NOT USE 
    507    ln_traldf_gdia   =  .false.  !  griffies operator strfn diagnostics  (require "key_ldfslp")  ! UNDER TEST, DO NOT USE 
    508    !                       !  Coefficient 
    509    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    510    rn_ahtb_0        =     0.    !  background eddy diffusivity for ldf_iso [m2/s] 
    511    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s]    (require "key_traldf_eiv") 
     470   !                       !  Operator type: 
     471   !                           !  no diffusion: set ln_traldf_lap=..._blp=F  
     472   ln_traldf_lap   =  .true.   !    laplacian operator 
     473   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     474   !                       !  Direction of action: 
     475   ln_traldf_lev   =  .false.  !  iso-level 
     476   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     477   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     478   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     479   ! 
     480   !                       !  iso-neutral options:         
     481   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     482   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     483   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     484   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     485   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     486   ! 
     487   !                       !  Coefficients: 
     488   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
     489   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     490   !                                !   =  0           constant  
     491   !                                !   = 10 F(k)      =ldf_c1d  
     492   !                                !   = 20 F(i,j)    =ldf_c2d  
     493   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     494   !                                !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
     495   !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
     496   rn_aht_0        = 1000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     497   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     498/ 
     499!---------------------------------------------------------------------------------- 
     500&namtra_ldfeiv !   eddy induced velocity param. 
     501!---------------------------------------------------------------------------------- 
     502   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    512503/ 
    513504!----------------------------------------------------------------------- 
     
    546537&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    547538!----------------------------------------------------------------------- 
    548    ln_dynvor_ene = .false. !  enstrophy conserving scheme   
    549    ln_dynvor_ens = .false. !  energy conserving scheme     
    550    ln_dynvor_mix = .false. !  mixed scheme                
    551    ln_dynvor_een = .true.  !  energy & enstrophy scheme   
     539   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     540   ln_dynvor_ens = .false. !  energy conserving scheme 
     541   ln_dynvor_mix = .false. !  mixed scheme 
     542   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     543      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
     544   ln_dynvor_msk = .true.  !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
    552545/ 
    553546!----------------------------------------------------------------------- 
     
    572565&namdyn_ldf    !   lateral diffusion on momentum 
    573566!----------------------------------------------------------------------- 
    574    !                       !  Type of the operator :  
    575    ln_dynldf_lap    =  .false.  !  laplacian operator          
    576    ln_dynldf_bilap  =  .true.   !  bilaplacian operator     
    577    !                       !  Direction of action  :  
    578    ln_dynldf_level  =  .false.  !  iso-level                
    579    ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.) 
    580    ln_dynldf_iso    =  .false.  !  iso-neutral                          (require "key_ldfslp") 
     567   !                       !  Type of the operator : 
     568   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     569   ln_dynldf_lap =  .false.    !    laplacian operator 
     570   ln_dynldf_blp =  .true.     !  bilaplacian operator 
     571   !                       !  Direction of action  : 
     572   ln_dynldf_lev =  .true.     !  iso-level 
     573   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     574   ln_dynldf_iso =  .false.    !  iso-neutral 
    581575   !                       !  Coefficient 
    582    rn_ahm_0_lap     = 40000.    !  horizontal laplacian eddy viscosity   [m2/s] 
    583    rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s] 
    584    rn_ahm_0_blp     = -8.5e+11  !  horizontal bilaplacian eddy viscosity [m4/s]  
     576   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     577   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     578   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     579   !                                !  =  0  constant  
     580   !                                !  = 10  F(k)=c1d 
     581   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     582   !                                !  = 30  F(i,j,k)=c2d*c1d 
     583   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     584   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     585   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     586   rn_bhm_0      = 8.5e+11     !  horizontal bilaplacian eddy viscosity [m4/s] 
    585587/ 
    586588 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/1_namelist_top_cfg

    r4147 r5836  
    3939&namtrc_adv    !   advection scheme for passive tracer  
    4040!----------------------------------------------------------------------- 
    41    ln_trcadv_cen2   =  .false.  !  2nd order centered scheme    
    42    ln_trcadv_tvd    =  .true.  !  TVD scheme 
    43    ln_trcadv_muscl  =  .false.   !  MUSCL scheme 
    44    ln_trcadv_muscl2 =  .false.  !  MUSCL2 scheme + cen2 at boundaries 
    45    ln_trcadv_ubs    =  .false.  !  UBS scheme 
    46    ln_trcadv_qck    =  .false.  !  QUICKEST scheme 
     41   ln_trcadv_fct =  .true.   !  FCT scheme 
     42      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     43      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     44      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     45      !                             !        (number of sub-timestep = nn_fct_zts) 
    4746/ 
    4847!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/namelist_cfg

    r5407 r5836  
    132132/ 
    133133!----------------------------------------------------------------------- 
    134 &namcla        !   cross land advection 
    135 !----------------------------------------------------------------------- 
    136 / 
    137 !----------------------------------------------------------------------- 
    138 &namobc        !   open boundaries parameters                           ("key_obc") 
    139 !----------------------------------------------------------------------- 
    140 / 
    141 !----------------------------------------------------------------------- 
    142134&namagrif      !  AGRIF zoom                                            ("key_agrif") 
    143135!----------------------------------------------------------------------- 
     
    176168/ 
    177169!----------------------------------------------------------------------- 
    178 &namtra_adv    !   advection scheme for tracer  
    179 !-----------------------------------------------------------------------      
     170&namtra_adv    !   advection scheme for tracer 
     171!----------------------------------------------------------------------- 
     172   ln_traadv_fct =  .true.    !  FCT scheme 
     173      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     174      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     175      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     176      !                             !        (number of sub-timestep = nn_fct_zts) 
    180177/ 
    181178!---------------------------------------------------------------------------------- 
    182179&namtra_ldf    !   lateral diffusion scheme for tracers 
    183180!---------------------------------------------------------------------------------- 
     181   !                       !  Operator type: 
     182   !                           !  no diffusion: set ln_traldf_lap=..._blp=F  
     183   ln_traldf_lap   =  .true.   !    laplacian operator 
     184   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     185   !                       !  Direction of action: 
     186   ln_traldf_lev   =  .false.  !  iso-level 
     187   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     188   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     189   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     190   ! 
     191   !                       !  iso-neutral options:         
     192   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     193   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     194   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     195   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     196   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     197   ! 
     198   !                       !  Coefficients: 
     199   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     200   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     201   !                                !   =  0           constant  
     202   !                                !   = 10 F(k)      =ldf_c1d  
     203   !                                !   = 20 F(i,j)    =ldf_c2d  
     204   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     205   !                                !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
     206   !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
     207   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     208   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     209/ 
     210!---------------------------------------------------------------------------------- 
     211&namtra_ldfeiv !   eddy induced velocity param. 
     212!---------------------------------------------------------------------------------- 
     213   ln_ldfeiv     =.true.    ! use eddy induced velocity parameterization 
     214   ln_ldfeiv_dia =.false.   ! diagnose eiv stream function and velocities 
     215   rn_aeiv_0     = 2000.    ! eddy induced velocity coefficient   [m2/s] 
     216   nn_aei_ijk_t  = 21       ! space/time variation of the eiv coeficient 
     217   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     218   !                                !   =  0           constant  
     219   !                                !   = 10 F(k)      =ldf_c1d  
     220   !                                !   = 20 F(i,j)    =ldf_c2d  
     221   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     222   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    184223/ 
    185224!----------------------------------------------------------------------- 
     
    188227/ 
    189228!----------------------------------------------------------------------- 
     229&namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
     230!----------------------------------------------------------------------- 
     231/ 
     232!----------------------------------------------------------------------- 
    190233&namdyn_adv    !   formulation of the momentum advection 
    191234!----------------------------------------------------------------------- 
    192235 
    193236!----------------------------------------------------------------------- 
    194 &namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
    195 !----------------------------------------------------------------------- 
    196 / 
    197 !----------------------------------------------------------------------- 
    198237&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    199238!----------------------------------------------------------------------- 
     239   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     240   ln_dynvor_ens = .false. !  energy conserving scheme 
     241   ln_dynvor_mix = .false. !  mixed scheme 
     242   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     243      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
    200244/ 
    201245!----------------------------------------------------------------------- 
     
    209253&namdyn_ldf    !   lateral diffusion on momentum 
    210254!-----------------------------------------------------------------------       
     255   !                       !  Type of the operator : 
     256   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     257   ln_dynldf_lap =  .true.     !    laplacian operator 
     258   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     259   !                       !  Direction of action  : 
     260   ln_dynldf_lev =  .true.     !  iso-level 
     261   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     262   ln_dynldf_iso =  .false.    !  iso-neutral 
     263   !                       !  Coefficient 
     264   nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
     265   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     266   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     267   !                                !  =  0  constant  
     268   !                                !  = 10  F(k)=c1d 
     269   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     270   !                                !  = 30  F(i,j,k)=c2d*c1d 
     271   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     272   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     273   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     274   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    211275/ 
    212276!----------------------------------------------------------------------- 
     
    221285&namzdf_tke    !   turbulent eddy kinetic dependent vertical diffusion  ("key_zdftke") 
    222286!----------------------------------------------------------------------- 
    223 / 
    224 !------------------------------------------------------------------------ 
    225 &namzdf_kpp    !   K-Profile Parameterization dependent vertical mixing  ("key_zdfkpp", and optionally: 
    226 !------------------------------------------------------------------------ "key_kppcustom" or "key_kpplktb") 
    227287/ 
    228288!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/namelist_top_cfg

    r4340 r5836  
    2424&namtrc_adv    !   advection scheme for passive tracer  
    2525!-----------------------------------------------------------------------    
    26    ln_trcadv_tvd     =  .true.   !  TVD scheme 
    27    ln_trcadv_muscl   =  .false.  !  MUSCL scheme 
     26   ln_trcadv_fct =  .true.   !  FCT scheme 
     27      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     28      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     29      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     30      !                             !        (number of sub-timestep = nn_fct_zts) 
    2831/ 
    2932!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/cpp_ORCA2_LIM_CFC_C14b.fcm

    r4523 r5836  
    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_top key_cfc key_c14b key_iomput key_mpp_mpi 
     1bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_zdftke key_zdfddm key_zdftmx key_top key_cfc key_c14b key_iomput key_mpp_mpi 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_cfg

    r4370 r5836  
    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 = 0             !  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!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_top_cfg

    r5385 r5836  
    6262&namtrc_adv    !   advection scheme for passive tracer  
    6363!----------------------------------------------------------------------- 
    64    ln_trcadv_tvd     =  .false.  !  TVD scheme 
    65    ln_trcadv_muscl   =  .true.   !  MUSCL scheme 
     64   ln_trcadv_mus  =  .true.  !  MUSCL scheme 
     65      ln_mus_ups  =  .false.         !  use upstream scheme near river mouths 
    6666/ 
    6767!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/cpp_ORCA2_LIM_PISCES.fcm

    r4523 r5836  
    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_top key_pisces key_iomput key_mpp_mpi 
     1bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_zdftke key_zdfddm key_zdftmx key_top key_pisces key_mpp_mpi key_iomput 
  • trunk/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_cfg

    r5385 r5836  
    7878!----------------------------------------------------------------------- 
    7979/       
    80  
    81 !----------------------------------------------------------------------- 
    82 &namcla        !   cross land advection 
    83 !----------------------------------------------------------------------- 
    84 / 
    8580!----------------------------------------------------------------------- 
    8681&nambbl        !   bottom boundary layer scheme 
     
    9186!----------------------------------------------------------------------- 
    9287/ 
    93 !----------------------------------------------------------------------- 
    94 &namtra_ldf    !   lateral diffusion scheme for tracer  
    95 !----------------------------------------------------------------------- 
    96    ln_triad_iso     =  .true.  !  griffies operator calculates triads twice => pure lateral mixing in ML (require "key_ldfslp") 
     88!---------------------------------------------------------------------------------- 
     89&namtra_ldf    !   lateral diffusion scheme for tracers 
     90!---------------------------------------------------------------------------------- 
     91   !                       !  Operator type: 
     92   ln_traldf_lap   =  .true.   !    laplacian operator 
     93   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     94   !                       !  Direction of action: 
     95   ln_traldf_lev   =  .false.  !  iso-level 
     96   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     97   ln_traldf_iso   =  .true.   !  iso-neutral 
     98   ln_traldf_triad =  .false.  !  iso-neutral using Griffies triads 
     99   ! 
     100   !                       !  iso-neutral options:         
     101   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     102   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     103   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML     (triad only) 
     104   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     105   ln_botmix_triad =  .false.  !  lateral mixing on bottom         (triad only) 
     106   ! 
     107   !                       !  Coefficients: 
     108   nn_aht_ijk_t    = 0        !  space/time variation of eddy coef 
     109   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     110   !                                !   =  0           constant  
     111   !                                !   = 10 F(k)      =ldf_c1d  
     112   !                                !   = 20 F(i,j)    =ldf_c2d  
     113   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     114   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     115   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     116   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     117   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     118/ 
     119!---------------------------------------------------------------------------------- 
     120&namtra_ldfeiv !   eddy induced velocity param. 
     121!---------------------------------------------------------------------------------- 
     122   ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
     123   ln_ldfeiv_dia =.false.   ! diagnose eiv stream function and velocities 
     124   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
     125   nn_aei_ijk_t  = 0      ! space/time variation of the eiv coeficient 
     126   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     127   !                                !   =  0           constant 
     128   !                                !   = 10 F(k)      =ldf_c1d 
     129   !                                !   = 20 F(i,j)    =ldf_c2d 
     130   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     131   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    97132/ 
    98133!----------------------------------------------------------------------- 
     
    115150   sn_ubl  = 'dyna_grid_U' ,    120            , 'sobblcox' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    116151   sn_vbl  = 'dyna_grid_V' ,    120            , 'sobblcoy' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    117    sn_ahu  = 'dyna_grid_U' ,    120            , 'vozoahtu' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    118    sn_ahv  = 'dyna_grid_V' ,    120            , 'vomeahtv' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    119    sn_ahw  = 'dyna_grid_W' ,    120            , 'voveahtz' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    120    sn_eiu  = 'dyna_grid_U' ,    120            , 'vozoaeiu' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    121    sn_eiv  = 'dyna_grid_V' ,    120            , 'vomeaeiv' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    122    sn_eiw  = 'dyna_grid_W' ,    120            , 'soleaeiw' ,  .true.    , .true. ,   'yearly'  , ''       , ''    , '' 
    123152! 
    124153   cn_dir      = './'      !  root directory for the location of the dynamical files 
    125    ln_degrad   =  .false.  !  flag for degradation -                requires ("key_degrad") 
    126154   ln_dynwzv   =  .true.   !  computation of vertical velocity instead of using the one read in file 
    127155   ln_dynbbl   =  .true.   !  bbl coef are in files, so read them - requires ("key_trabbl") 
  • trunk/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_top_cfg

    r5385 r5836  
    6565&namtrc_adv    !   advection scheme for passive tracer  
    6666!----------------------------------------------------------------------- 
    67    ln_trcadv_tvd     =  .false.  !  TVD scheme 
    68    ln_trcadv_muscl   =  .true.   !  MUSCL scheme 
     67   ln_trcadv_mus =  .true.   !  MUSCL scheme 
     68      ln_mus_ups =  .false.         !  use upstream scheme near river mouths 
    6969/ 
    7070!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/cpp_ORCA2_OFF_PISCES.fcm

    r4523 r5836  
    1 bld::tool::fppkeys key_trabbl key_ldfslp key_traldf_c2d key_traldf_eiv key_top key_offline key_pisces key_iomput key_mpp_mpi 
     1bld::tool::fppkeys key_trabbl key_top key_offline key_pisces key_iomput key_mpp_mpi 
  • trunk/NEMOGCM/CONFIG/ORCA2_SAS_LIM/EXP00/namelist_cfg

    r4370 r5836  
    8787/ 
    8888!----------------------------------------------------------------------- 
    89 &namcla        !   cross land advection 
    90 !----------------------------------------------------------------------- 
    91 / 
    92 !----------------------------------------------------------------------- 
    9389&nameos        !   ocean physical parameters 
    9490!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/ORCA2_SAS_LIM/cpp_ORCA2_SAS_LIM.fcm

    r4523 r5836  
    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 
     1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi 
  • trunk/NEMOGCM/CONFIG/SHARED/domain_def.xml

    r5426 r5836  
    163163 
    164164      <domain id="grid_V" long_name="grid V"/> 
     165       
    165166      <domain_group id="grid_W"> 
    166167         <domain id="grid_W" long_name="grid W"/> 
     
    168169         <domain id="EqW" zoom_ibegin="1" zoom_jbegin="0000" zoom_ni="0000" zoom_nj="1" /> 
    169170      </domain_group> 
     171 
     172      <domain id="grid_F" long_name="grid F"/> 
    170173 
    171174      <domain_group id="scalarpoint">  
  • trunk/NEMOGCM/CONFIG/SHARED/field_def.xml

    r5517 r5836  
    2525 
    2626         <field id="toce"         long_name="temperature"         standard_name="sea_water_potential_temperature"   unit="degC"     grid_ref="grid_T_3D"/> 
    27          <field id="toce_e3t"     long_name="temperature * e3t"                                                     unit="degC*m"   grid_ref="grid_T_3D" > toce * e3t </field > 
     27         <field id="toce_e3t"     long_name="temperature (thickness weighted)"                                      unit="degC"     grid_ref="grid_T_3D" > toce * e3t </field > 
    2828         <field id="soce"         long_name="salinity"            standard_name="sea_water_practical_salinity"      unit="1e-3"     grid_ref="grid_T_3D"/> 
    29          <field id="soce_e3t"     long_name="salinity * e3t"                                                        unit="1e-3*m"   grid_ref="grid_T_3D" > soce * e3t </field > 
     29         <field id="soce_e3t"     long_name="salinity    (thickness weighted)"                                      unit="1e-3"     grid_ref="grid_T_3D" > soce * e3t </field > 
     30 
     31         <!-- t-eddy viscosity coefficients (ldfdyn) --> 
     32         <field id="ahmt_2d"      long_name=" surface t-eddy viscosity coefficient"   unit="m2/s or m4/s" /> 
     33         <field id="ahmt_3d"      long_name=" 3D      t-eddy viscosity coefficient"   unit="m2/s or m4/s"                           grid_ref="grid_T_3D"/> 
    3034 
    3135         <field id="sst"          long_name="sea surface temperature"             standard_name="sea_surface_temperature"             unit="degC"     /> 
     
    367371         <field id="utau"         long_name="Wind Stress along i-axis"                               standard_name="surface_downward_x_stress"   unit="N/m2"                            /> 
    368372         <field id="uoce"         long_name="ocean current along i-axis"                             standard_name="sea_water_x_velocity"        unit="m/s"        grid_ref="grid_U_3D" /> 
    369          <field id="uoce_e3u"     long_name="ocean current along i-axis * e3u"                                                                   unit="m2/s"       grid_ref="grid_U_3D"  > uoce * e3u </field> 
     373         <field id="uoce_e3u"     long_name="ocean current along i-axis  (thickness weighted)"                                                   unit="m/s"        grid_ref="grid_U_3D"  > uoce * e3u </field> 
    370374         <field id="ssu"          long_name="ocean surface current along i-axis"                                                                 unit="m/s"                             /> 
    371375         <field id="sbu"          long_name="ocean bottom current along i-axis"                                                                  unit="m/s"                             /> 
     
    375379         <field id="uoces"        long_name="ocean transport along i-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_U_3D" /> 
    376380 
     381         <!-- u-eddy coefficients (ldftra) --> 
     382         <field id="ahtu_2d"      long_name=" surface u-eddy diffusivity coefficient"   unit="m2/s or m4/s" /> 
     383         <field id="ahtu_3d"      long_name=" 3D u-EIV coefficient"                     unit="m2/s or m4/s"      grid_ref="grid_U_3D"/> 
     384         <field id="aeiu_2d"      long_name=" surface u-EIV coefficient"                unit="m2/s" /> 
     385         <field id="aeiu_3d"      long_name=" 3D u-EIV coefficient"                     unit="m2/s"              grid_ref="grid_U_3D"/> 
     386 
    377387         <!-- variables available with MLE --> 
    378388         <field id="psiu_mle"     long_name="MLE streamfunction along i-axis"   unit="m3/s"   grid_ref="grid_U_3D" /> 
    379389 
    380          <!-- uoce_eiv: available with key_traldf_eiv and key_diaeiv --> 
     390         <!-- uoce_eiv: available EIV --> 
    381391         <field id="uoce_eiv"     long_name="EIV ocean current along i-axis"   standard_name="bolus_sea_water_x_velocity"   unit="m/s"   grid_ref="grid_U_3D" /> 
    382392 
     
    402412         <field id="vtau"         long_name="Wind Stress along j-axis"                               standard_name="surface_downward_y_stress"   unit="N/m2"                            /> 
    403413         <field id="voce"         long_name="ocean current along j-axis"                             standard_name="sea_water_y_velocity"        unit="m/s"        grid_ref="grid_V_3D" /> 
    404          <field id="voce_e3v"     long_name="ocean current along j-axis * e3v"                                                                   unit="m2/s"       grid_ref="grid_V_3D"  > voce * e3v </field> 
     414         <field id="voce_e3v"     long_name="ocean current along j-axis  (thickness weighted)"                                                   unit="m/s"        grid_ref="grid_V_3D"  > voce * e3v </field> 
    405415         <field id="ssv"          long_name="ocean surface current along j-axis"                                                                 unit="m/s"                             /> 
    406416         <field id="sbv"          long_name="ocean bottom current along j-axis"                                                                  unit="m/s"                             /> 
     
    410420         <field id="voces"        long_name="ocean transport along j-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_V_3D" /> 
    411421 
     422         <!-- v-eddy coefficients (ldftra, ldfdyn) --> 
     423         <field id="ahtv_2d"      long_name=" surface v-eddy diffusivity coefficient"     unit="m2/s or (m4/s)^1/2" /> 
     424         <field id="ahtv_3d"      long_name=" 3D v-eddy diffusivity coefficient"          unit="m2/s or (m4/s)^1/2"           grid_ref="grid_V_3D"/> 
     425         <field id="aeiv_2d"      long_name=" surface v-EIV coefficient"                  unit="m2/s" /> 
     426         <field id="aeiv_3d"      long_name=" 3D v-EIV coefficient"                       unit="m2/s"                         grid_ref="grid_V_3D" /> 
     427 
    412428         <!-- variables available with MLE --> 
    413429         <field id="psiv_mle"     long_name="MLE streamfunction along j-axis"   unit="m3/s"   grid_ref="grid_V_3D" /> 
    414430 
    415          <!-- voce_eiv: available with key_traldf_eiv and key_diaeiv --> 
     431         <!-- voce_eiv: available with EIV --> 
    416432         <field id="voce_eiv"     long_name="EIV ocean current along j-axis"   standard_name="bolus_sea_water_y_velocity"   unit="m/s"   grid_ref="grid_V_3D" /> 
    417433 
     
    438454        <field id="wocetr_eff"   long_name="effective ocean vertical transport"                                               unit="m3/s" /> 
    439455 
    440         <!-- woce_eiv: available with key_traldf_eiv and key_diaeiv --> 
     456        <!-- woce_eiv: available with EIV --> 
    441457        <field id="woce_eiv"     long_name="EIV ocean vertical velocity"   standard_name="bolus_upward_sea_water_velocity"   unit="m/s" /> 
    442458 
     
    449465 
    450466        <!-- avt_evd and avm_evd: available with ln_zdfevd --> 
    451         <field id="avt_evd"      long_name="convective enhancement to vertical diffusivity"   standard_name="ocean_vertical_tracer_diffusivity_due_to_convection"     unit="m2/s" /> 
    452         <field id="avm_evd"      long_name="convective enhancement to vertical viscosity"     standard_name="ocean_vertical_momentum_diffusivity_due_to_convection"   unit="m2/s" /> 
     467        <field id="avt_evd"      long_name="convective enhancement of vertical diffusivity"   standard_name="ocean_vertical_tracer_diffusivity_due_to_convection"     unit="m2/s" /> 
     468        <field id="avm_evd"      long_name="convective enhancement of vertical viscosity"     standard_name="ocean_vertical_momentum_diffusivity_due_to_convection"   unit="m2/s" /> 
    453469 
    454470        <!-- avt_tide: available with key_zdftmx --> 
     
    459475        <field id="w_masstr2"    long_name="square of vertical mass transport"   standard_name="square_of_upward_ocean_mass_transport"   unit="kg2/s2" /> 
    460476 
    461         <!-- aht2d and  aht2d_eiv: available with key_traldf_eiv and key_traldf_c2d --> 
     477        <!-- aht2d and  aht2d_eiv --> 
    462478        <field id="aht2d"        long_name="lateral eddy diffusivity"       standard_name="ocean_tracer_xy_laplacian_diffusivity"      unit="m2/s"   grid_ref="grid_W_2D" /> 
    463479        <field id="aht2d_eiv"    long_name="EIV lateral eddy diffusivity"   standard_name="ocean_tracer_bolus_laplacian_diffusivity"   unit="m2/s"   grid_ref="grid_W_2D" /> 
    464480      </field_group> 
    465            
     481         
     482      <!-- F grid --> 
     483         <!-- f-eddy viscosity coefficients (ldfdyn) --> 
     484         <field id="ahmf_2d"      long_name=" surface f-eddy viscosity coefficient"   unit="m2/s or m4/s" /> 
     485         <field id="ahmf_3d"      long_name=" 3D      f-eddy viscosity coefficient"   unit="m2/s or m4/s"                           grid_ref="grid_T_3D"/> 
     486 
    466487      <!-- scalar variables available with key_diaar5 --> 
    467488 
  • trunk/NEMOGCM/CONFIG/SHARED/namelist_ref

    r5656 r5836  
     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!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/SHARED/namelist_top_ref

    r5416 r5836  
    4141&namtrc_adv    !   advection scheme for passive tracer  
    4242!----------------------------------------------------------------------- 
    43    ln_trcadv_cen2    =  .false.  !  2nd order centered scheme    
    44    ln_trcadv_tvd     =  .true.  !  TVD scheme 
    45    ln_trcadv_muscl   =  .false.   !  MUSCL scheme 
    46    ln_trcadv_muscl2  =  .false.  !  MUSCL2 scheme + cen2 at boundaries 
    47    ln_trcadv_ubs     =  .false.  !  UBS scheme 
    48    ln_trcadv_qck     =  .false.  !  QUICKEST scheme 
    49    ln_trcadv_msc_ups =  .false.  !  use upstream scheme within muscl 
     43   ln_trcadv_cen =  .false.  !  2nd order centered scheme 
     44      nn_cen_h   =  4               !  =2/4, horizontal 2nd order CEN / 4th order CEN 
     45      nn_cen_v   =  4               !  =2/4, vertical   2nd order CEN / 4th order COMPACT 
     46   ln_trcadv_fct =  .false.  !  FCT scheme 
     47      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     48      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     49      nn_fct_zts =  0               !  >=1,  2nd order FCT scheme with vertical sub-timestepping 
     50      !                             !        (number of sub-timestep = nn_fct_zts) 
     51   ln_trcadv_mus =  .false.  !  MUSCL scheme 
     52      ln_mus_ups =  .false.         !  use upstream scheme near river mouths 
     53   ln_trcadv_ubs =  .false.  !  UBS scheme 
     54      nn_ubs_v   =  2               !  =2  , vertical 2nd order FCT 
     55   ln_trcadv_qck =  .false.  !  QUICKEST scheme 
    5056/ 
    5157!----------------------------------------------------------------------- 
    5258&namtrc_ldf    !   lateral diffusion scheme for passive tracer  
    5359!----------------------------------------------------------------------- 
    54 !                               !  Type of the operator :  
    55    ln_trcldf_lap    =  .true.   !     laplacian operator        
    56    ln_trcldf_bilap  =  .false.  !     bilaplacian operator      
    57                                 !  Direction of action  : 
    58    ln_trcldf_level  =  .false.  !     iso-level                 
    59    ln_trcldf_hor    =  .false.  !     horizontal (geopotential)         (require "key_ldfslp" when ln_sco=T) 
    60    ln_trcldf_iso    =  .true.   !     iso-neutral                       (require "key_ldfslp") 
    61 !                               !  Coefficient 
    62    rn_ahtrc_0       =  2000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    63    rn_ahtrb_0       =     0.    !     background eddy diffusivity for ldf_iso [m2/s] 
     60!                          !  Type of the operator:   
     61   ln_trcldf_lap   =  .true.   !    laplacian operator 
     62   ln_trcldf_blp   =  .false.  !  bilaplacian operator 
     63   !                       !  Direction of action: 
     64   ln_trcldf_lev   =  .false.  !  iso-level 
     65   ln_trcldf_hor   =  .false.  !  horizontal (geopotential) 
     66   ln_trcldf_iso   =  .true.   !  iso-neutral (standard operator) 
     67   ln_trcldf_triad =  .false.  !  iso-neutral (triad    operator) 
     68   !                       !  Coefficient  
     69   rn_ahtrc_0      = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     70   rn_bhtrc_0      = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
    6471/ 
    6572!----------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.