Changeset 5778


Ignore:
Timestamp:
2015-10-06T13:46:54+02:00 (5 years ago)
Author:
gm
Message:

#1593: LDF-ADV, III. Phasing of the improvements/simplifications of ADV & LDF momentum trends: associated namelist & cpp

Location:
branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG
Files:
22 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/AMM12/EXP00/namelist_cfg

    r5771 r5778  
    304304   !                       !  Coefficients: 
    305305   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
    306    !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     306   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    307307   !                                !   =  0           constant  
    308308   !                                !   = 10 F(k)      =ldf_c1d  
     
    331331&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    332332!----------------------------------------------------------------------- 
     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) 
     338   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
    333339/ 
    334340!----------------------------------------------------------------------- 
     
    349355!----------------------------------------------------------------------- 
    350356   !                       !  Type of the operator : 
    351    ln_dynldf_bilap  =  .true.   !  bilaplacian operator 
    352    ln_dynldf_lap    =  .false.  !  bilaplacian operator 
     357   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     358   ln_dynldf_lap =  .false.    !    laplacian operator 
     359   ln_dynldf_blp =  .true.     !  bilaplacian operator 
    353360   !                       !  Direction of action  : 
    354    ln_dynldf_level  =  .true.   !  iso-level 
    355    ln_dynldf_hor    =  .false.  !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.) 
    356                            !  Coefficient 
    357    rn_ahm_0_lap     = 60.0      !  horizontal laplacian eddy viscosity   [m2/s] 
    358    rn_ahm_0_blp     = -1.0e+10  !  horizontal bilaplacian eddy viscosity [m4/s] 
     361   ln_dynldf_lev =  .true.     !  iso-level 
     362   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     363   ln_dynldf_iso =  .false.    !  iso-neutral 
     364   !                       !  Coefficient 
     365   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     366   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     367   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     368   !                                !  =  0  constant  
     369   !                                !  = 10  F(k)=c1d 
     370   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     371   !                                !  = 30  F(i,j,k)=c2d*c1d 
     372   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     373   rn_ahm_0      =     60.     !  horizontal laplacian eddy viscosity   [m2/s] 
     374   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     375   rn_bhm_0      = 1.0e+10     !  horizontal bilaplacian eddy viscosity [m4/s] 
     376   ! 
     377   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
    359378/ 
    360379!----------------------------------------------------------------------- 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/C1D_PAPA/EXP00/namelist_cfg

    r5771 r5778  
    224224!----------------------------------------------------------------------- 
    225225/ 
    226 !---------------------------------------------------------------------------------- 
     226!----------------------------------------------------------------------- 
    227227&namtra_ldf    !   lateral diffusion scheme for tracers 
    228 !---------------------------------------------------------------------------------- 
    229 !---------------------------------------------------------------------------------- 
    230    ln_traldf_hor    =  .true.   !  horizontal (geopotential)   (needs "key_ldfslp" when ln_sco=T) 
    231    ln_traldf_iso    =  .false.  !  iso-neutral                 (needs "key_ldfslp") 
    232    rn_aeiv_0        =     0.    !  eddy induced velocity coefficient [m2/s] 
    233    rn_aht_0         =     0.    !  horizontal eddy diffusivity for tracers [m2/s] 
    234 / 
    235 !---------------------------------------------------------------------------------- 
    236 &namtra_ldf    !   lateral diffusion scheme for tracers 
    237 !---------------------------------------------------------------------------------- 
    238    !                       !  Operator type: 
    239    ln_traldf_lap   =  .false.  !    laplacian operator 
    240    ln_traldf_blp   =  .false.  !  bilaplacian operator 
    241    !                       !  Direction of action: 
    242    ln_traldf_lev   =  .false.  !  iso-level 
    243    ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
    244    ln_traldf_iso   =  .false.  !  iso-neutral 
    245    ln_traldf_triad =  .false.  !  iso-neutral using Griffies triads 
    246 / 
    247 !---------------------------------------------------------------------------------- 
     228!----------------------------------------------------------------------- 
     229! C1D : no lateral diffusion   
     230/ 
     231!----------------------------------------------------------------------- 
    248232&namtra_ldfeiv !   eddy induced velocity param. 
    249 !---------------------------------------------------------------------------------- 
    250    ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
     233!----------------------------------------------------------------------- 
     234! C1D : no eiv   
    251235/ 
    252236!----------------------------------------------------------------------- 
     
    258242&namdyn_adv    !   formulation of the momentum advection 
    259243!----------------------------------------------------------------------- 
     244! C1D : no advection scheme  
    260245/ 
    261246!----------------------------------------------------------------------- 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/GYRE/EXP00/namelist_cfg

    r5771 r5778  
    273273&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    274274!----------------------------------------------------------------------- 
    275    ln_dynvor_ene = .true.  !  energy    conserving scheme   
    276    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 
    277278   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) 
     280   ln_dynvor_msk = .true.  !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
    278281/ 
    279282!----------------------------------------------------------------------- 
     
    290293&namdyn_ldf    !   lateral diffusion on momentum 
    291294!----------------------------------------------------------------------- 
    292    rn_ahm_0_lap     = 100000.   !  horizontal laplacian eddy viscosity   [m2/s] 
     295   !                       !  Type of the operator : 
     296   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     297   ln_dynldf_lap =  .true.     !    laplacian operator 
     298   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     299   !                       !  Direction of action  : 
     300   ln_dynldf_lev =  .true.     !  iso-level 
     301   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     302   ln_dynldf_iso =  .false.    !  iso-neutral 
     303   !                       !  Coefficient 
     304   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     305   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     306   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     307   !                                !  =  0  constant  
     308   !                                !  = 10  F(k)=c1d 
     309   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     310   !                                !  = 30  F(i,j,k)=c2d*c1d 
     311   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     312   rn_ahm_0      = 100000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     313   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     314   rn_bhm_0      =      0.     !  horizontal bilaplacian eddy viscosity [m4/s] 
    293315/ 
    294316!----------------------------------------------------------------------- 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist_cfg

    r5771 r5778  
    251251   !                       !  Coefficients: 
    252252   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
    253    !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     253   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    254254   !                                !   =  0           constant  
    255255   !                                !   = 10 F(k)      =ldf_c1d  
     
    278278&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    279279!----------------------------------------------------------------------- 
    280    ln_dynvor_ene = .true.  !  energy    conserving scheme   
    281    ln_dynvor_ens = .false. !  enstrophy conserving scheme     
    282    ln_dynvor_een = .false. !  energy & enstrophy scheme 
     280   ln_dynvor_ene = .true.  !  enstrophy conserving scheme 
     281   ln_dynvor_msk = .true.  !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
    283282/ 
    284283!----------------------------------------------------------------------- 
     
    294293&namdyn_ldf    !   lateral diffusion on momentum 
    295294!----------------------------------------------------------------------- 
     295   !                       !  Type of the operator : 
     296   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     297   ln_dynldf_lap =  .true.     !    laplacian operator 
     298   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     299   !                       !  Direction of action  : 
     300   ln_dynldf_lev =  .true.     !  iso-level 
     301   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     302   ln_dynldf_iso =  .false.    !  iso-neutral 
     303   !                       !  Coefficient 
     304   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     305   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     306   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     307   !                                !  =  0  constant  
     308   !                                !  = 10  F(k)=c1d 
     309   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     310   !                                !  = 30  F(i,j,k)=c2d*c1d 
     311   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     312   rn_ahm_0      = 100000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     313   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     314   rn_bhm_0      =      0.     !  horizontal bilaplacian eddy viscosity [m4/s] 
     315   ! 
     316   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
     317/ 
    296318   rn_ahm_0_lap     = 100000.   !  horizontal laplacian eddy viscosity   [m2/s] 
    297319/ 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/GYRE_BFM/cpp_GYRE_BFM.fcm

    r5759 r5778  
    1 bld::tool::fppkeys key_dynspg_flt 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 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/GYRE_PISCES/EXP00/namelist_cfg

    r5771 r5778  
    166166   !                       !  Coefficients: 
    167167   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
    168    !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     168   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    169169   !                                !   =  0           constant  
    170170   !                                !   = 10 F(k)      =ldf_c1d  
     
    193193&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    194194!----------------------------------------------------------------------- 
    195    ln_dynvor_ene = .true.  !  energy    conserving scheme   
    196    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 
    197198   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) 
     200   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
    198201/ 
    199202!----------------------------------------------------------------------- 
     
    206209&namdyn_ldf    !   lateral diffusion on momentum 
    207210!----------------------------------------------------------------------- 
     211   !                       !  Type of the operator : 
     212   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     213   ln_dynldf_lap =  .true.    !    laplacian operator 
     214   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     215   !                       !  Direction of action  : 
     216   ln_dynldf_lev =  .true.    !  iso-level 
     217   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     218   ln_dynldf_iso =  .false.    !  iso-neutral 
     219   !                       !  Coefficient 
     220   nn_ahm_ijk_t  = 0           !  space/time variation of eddy coef 
     221   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     222   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     223   !                                !  =  0  constant  
     224   !                                !  = 10  F(k)=c1d 
     225   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     226   !                                !  = 30  F(i,j,k)=c2d*c1d 
     227   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     228   rn_ahm_0      = 100000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     229   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     230   rn_bhm_0      =      0.      !  horizontal bilaplacian eddy viscosity [m4/s] 
     231   ! 
     232   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
     233/ 
    208234   rn_ahm_0_lap     = 100000.   !  horizontal laplacian eddy viscosity   [m2/s] 
    209235/ 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/1_namelist_cfg

    r5774 r5778  
    152152&namtra_ldf    !   lateral diffusion scheme for tracers 
    153153!----------------------------------------------------------------------- 
    154    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
     154   !                       !  Operator type: 
     155   !                           !  no diffusion: set ln_traldf_lap=..._blp=F  
     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    = 0         !  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 and grid-spacing) 
     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] 
    155182/ 
    156183!---------------------------------------------------------------------------------- 
    157184&namtra_ldfeiv !   eddy induced velocity param. 
    158185!---------------------------------------------------------------------------------- 
    159    rn_aeiv_0     = 0.   ! eddy induced velocity coefficient   [m2/s] 
     186   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    160187/ 
    161188!----------------------------------------------------------------------- 
     
    166193&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    167194!----------------------------------------------------------------------- 
     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) 
    168201/ 
    169202!----------------------------------------------------------------------- 
     
    177210   ln_dynldf_lap    =  .false.  !  laplacian operator          
    178211   ln_dynldf_bilap  =  .true.   !  bilaplacian operator     
    179    rn_ahm_0_blp     = -8.5e+11  !  horizontal bilaplacian eddy viscosity [m4/s]  
     212   rn_ahm_0_blp     =  8.5e+11  !  horizontal bilaplacian eddy viscosity [m4/s]  
    180213/ 
    181214!----------------------------------------------------------------------- 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist_cfg

    r5771 r5778  
    180180   !                       !  Coefficients: 
    181181   nn_aht_ijk_t    = 21        !  space/time variation of eddy coef 
    182    !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     182   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    183183   !                                !   =  0           constant  
    184184   !                                !   = 10 F(k)      =ldf_c1d  
     
    214214&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    215215!----------------------------------------------------------------------- 
     216   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     217   ln_dynvor_ens = .false. !  energy conserving scheme 
     218   ln_dynvor_mix = .false. !  mixed scheme 
     219   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     220      nn_een_e3f = 1             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
     221   ln_dynvor_msk = .true.  !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
    216222/ 
    217223!----------------------------------------------------------------------- 
     
    222228&namdyn_ldf    !   lateral diffusion on momentum 
    223229!----------------------------------------------------------------------- 
     230   !                       !  Type of the operator : 
     231   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     232   ln_dynldf_lap =  .true.     !    laplacian operator 
     233   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     234   !                       !  Direction of action  : 
     235   ln_dynldf_lev =  .true.     !  iso-level 
     236   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     237   ln_dynldf_iso =  .false.    !  iso-neutral 
     238   !                       !  Coefficient 
     239   nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
     240   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     241   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     242   !                                !  =  0  constant  
     243   !                                !  = 10  F(k)=c1d 
     244   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     245   !                                !  = 30  F(i,j,k)=c2d*c1d 
     246   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     247   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     248   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     249   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    224250/ 
    225251!----------------------------------------------------------------------- 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_LIM/cpp_ORCA2_LIM.fcm

    r5766 r5778  
    1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt  key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx  key_mpp_mpi key_iomput 
     1bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_zdftke key_zdfddm key_zdftmx  key_mpp_mpi key_iomput 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/1_namelist_cfg

    r5771 r5778  
    150150/ 
    151151!----------------------------------------------------------------------- 
    152 &namtra_ldf    !   lateral diffusion scheme for tracers  
    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") 
     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!----------------------------------------------------------------------- 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/namelist_cfg

    r5771 r5778  
    162162   ! 
    163163   !                       !  Coefficients: 
    164    nn_aht_ijk_t    = 21        !  space/time variation of eddy coef 
    165    !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     164   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     165   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    166166   !                                !   =  0           constant  
    167167   !                                !   = 10 F(k)      =ldf_c1d  
     
    197197&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    198198!----------------------------------------------------------------------- 
     199   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     200   ln_dynvor_ens = .false. !  energy conserving scheme 
     201   ln_dynvor_mix = .false. !  mixed scheme 
     202   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     203      nn_een_e3f = 0             !  e3f = masked averaging of e3t divided by 4 (=0) or by the sum of mask (=1) 
     204   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
    199205/ 
    200206!----------------------------------------------------------------------- 
     
    202208!----------------------------------------------------------------------- 
    203209   ln_hpg_zco  = .false.   !  z-coordinate - full steps 
    204    ln_hpg_zps  = .false.   !  z-coordinate - partial steps (interpolation) 
     210   ln_hpg_zps  = .true.    !  z-coordinate - partial steps (interpolation) 
    205211   ln_hpg_sco  = .false.   !  s-coordinate (standard jacobian formulation) 
    206212   ln_hpg_djc  = .false.   !  s-coordinate (Density Jacobian with Cubic polynomial) 
    207    ln_hpg_prj  = .true.    !  s-coordinate (Pressure Jacobian scheme) 
     213   ln_hpg_prj  = .false.   !  s-coordinate (Pressure Jacobian scheme) 
    208214   ln_dynhpg_imp = .false. !  time stepping: semi-implicit time scheme  (T) 
    209215                                 !           centered      time scheme  (F) 
     
    212218&namdyn_ldf    !   lateral diffusion on momentum 
    213219!----------------------------------------------------------------------- 
     220   !                       !  Type of the operator : 
     221   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     222   ln_dynldf_lap =  .true.     !    laplacian operator 
     223   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     224   !                       !  Direction of action  : 
     225   ln_dynldf_lev =  .true.     !  iso-level 
     226   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     227   ln_dynldf_iso =  .false.    !  iso-neutral 
     228   !                       !  Coefficient 
     229   nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
     230   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     231   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     232   !                                !  =  0  constant  
     233   !                                !  = 10  F(k)=c1d 
     234   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     235   !                                !  = 30  F(i,j,k)=c2d*c1d 
     236   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     237   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     238   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     239   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
     240   ! 
     241   ! Caution in 20 and 30 cases the coefficient have to be given for a 1 degree grid (~111km) 
    214242/ 
    215243!----------------------------------------------------------------------- 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_LIM3/cpp_ORCA2_LIM3.fcm

    r5759 r5778  
    1  bld::tool::fppkeys key_trabbl key_lim3 key_vvl key_dynspg_ts key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx  key_mpp_mpi key_diaobs key_asminc key_iomput key_nosignedzero 
     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 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/1_namelist_cfg

    r5771 r5778  
    366366&nambdy_tide     ! tidal forcing at open boundaries               
    367367!----------------------------------------------------------------------- 
    368    filtide      = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files 
     368    filtide       = 'bdydta/amm12_bdytide_'         !  file name root of tidal forcing files 
    369369    tide_cpt(1)   ='Q1'  !  names of tidal components used 
    370370    tide_cpt(2)   ='O1'  !  names of tidal components used 
     
    468468&namtra_ldf    !   lateral diffusion scheme for tracer  
    469469!----------------------------------------------------------------------- 
    470    !                       !  Type of the operator :  
    471    ln_traldf_lap    =  .true.   !  laplacian operator        
    472    ln_traldf_bilap  =  .false.  !  bilaplacian operator      
    473    !                       !  Direction of action  : 
    474    ln_traldf_level  =  .false.  !  iso-level                 
    475    ln_traldf_hor    =  .false.  !  horizontal (geopotential)            (require "key_ldfslp" when ln_sco=T) 
    476    ln_traldf_iso    =  .true.   !  iso-neutral                          (require "key_ldfslp") 
    477    ln_traldf_grif   =  .false.  !  griffies skew flux formulation       (require "key_ldfslp")  ! UNDER TEST, DO NOT USE 
    478    ln_traldf_gdia   =  .false.  !  griffies operator strfn diagnostics  (require "key_ldfslp")  ! UNDER TEST, DO NOT USE 
    479    !                       !  Coefficient 
    480    rn_aht_0         =  1000.    !  horizontal eddy diffusivity for tracers [m2/s] 
    481    rn_ahtb_0        =     0.    !  background eddy diffusivity for ldf_iso [m2/s] 
    482    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 
    483503/ 
    484504!----------------------------------------------------------------------- 
     
    517537&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    518538!----------------------------------------------------------------------- 
    519    ln_dynvor_ene = .false. !  enstrophy conserving scheme   
    520    ln_dynvor_ens = .false. !  energy conserving scheme     
    521    ln_dynvor_mix = .false. !  mixed scheme                
    522    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) 
    523545/ 
    524546!----------------------------------------------------------------------- 
     
    543565&namdyn_ldf    !   lateral diffusion on momentum 
    544566!----------------------------------------------------------------------- 
    545    !                       !  Type of the operator :  
    546    ln_dynldf_lap    =  .false.  !  laplacian operator          
    547    ln_dynldf_bilap  =  .true.   !  bilaplacian operator     
    548    !                       !  Direction of action  :  
    549    ln_dynldf_level  =  .false.  !  iso-level                
    550    ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.) 
    551    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 
    552575   !                       !  Coefficient 
    553    rn_ahm_0_lap     = 40000.    !  horizontal laplacian eddy viscosity   [m2/s] 
    554    rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s] 
    555    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] 
    556587/ 
    557588 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/namelist_cfg

    r5771 r5778  
    179179&namtra_ldf    !   lateral diffusion scheme for tracers 
    180180!---------------------------------------------------------------------------------- 
     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 
    181223/ 
    182224!----------------------------------------------------------------------- 
     
    185227/ 
    186228!----------------------------------------------------------------------- 
     229&namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
     230!----------------------------------------------------------------------- 
     231/ 
     232!----------------------------------------------------------------------- 
    187233&namdyn_adv    !   formulation of the momentum advection 
    188234!----------------------------------------------------------------------- 
    189235 
    190236!----------------------------------------------------------------------- 
    191 &namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
    192 !----------------------------------------------------------------------- 
    193 / 
    194 !----------------------------------------------------------------------- 
    195237&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    196238!----------------------------------------------------------------------- 
     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) 
     244   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
    197245/ 
    198246!----------------------------------------------------------------------- 
     
    206254&namdyn_ldf    !   lateral diffusion on momentum 
    207255!-----------------------------------------------------------------------       
     256   !                       !  Type of the operator : 
     257   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     258   ln_dynldf_lap =  .true.     !    laplacian operator 
     259   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     260   !                       !  Direction of action  : 
     261   ln_dynldf_lev =  .true.     !  iso-level 
     262   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     263   ln_dynldf_iso =  .false.    !  iso-neutral 
     264   !                       !  Coefficient 
     265   nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
     266   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     267   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     268   !                                !  =  0  constant  
     269   !                                !  = 10  F(k)=c1d 
     270   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     271   !                                !  = 30  F(i,j,k)=c2d*c1d 
     272   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     273   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     274   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     275   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    208276/ 
    209277!----------------------------------------------------------------------- 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/cpp_ORCA2_LIM_CFC_C14b.fcm

    r5759 r5778  
    1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt 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 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_cfg

    r5771 r5778  
    141141   ! 
    142142   !                       !  Coefficients: 
    143    nn_aht_ijk_t    = 21        !  space/time variation of eddy coef 
    144    !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     143   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     144   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    145145   !                                !   =  0           constant 
    146146   !                                !   = 10 F(k)      =ldf_c1d 
     
    177177&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    178178!----------------------------------------------------------------------- 
     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) 
     184   ln_dynvor_msk = .false. !  vorticity multiplied by fmask (=T) or not (=F) (all vorticity schemes) 
    179185/ 
    180186!----------------------------------------------------------------------- 
     
    185191&namdyn_ldf    !   lateral diffusion on momentum 
    186192!----------------------------------------------------------------------- 
     193   !                       !  Type of the operator : 
     194   !                           !  no diffusion: set ln_dynldf_lap=..._blp=F  
     195   ln_dynldf_lap =  .true.     !    laplacian operator 
     196   ln_dynldf_blp =  .false.    !  bilaplacian operator 
     197   !                       !  Direction of action  : 
     198   ln_dynldf_lev =  .true.     !  iso-level 
     199   ln_dynldf_hor =  .false.    !  horizontal (geopotential) 
     200   ln_dynldf_iso =  .false.    !  iso-neutral 
     201   !                       !  Coefficient 
     202   nn_ahm_ijk_t  = -30         !  space/time variation of eddy coef 
     203   !                                !  =-30  read in eddy_viscosity_3D.nc file 
     204   !                                !  =-20  read in eddy_viscosity_2D.nc file 
     205   !                                !  =  0  constant  
     206   !                                !  = 10  F(k)=c1d 
     207   !                                !  = 20  F(i,j)=F(grid spacing)=c2d 
     208   !                                !  = 30  F(i,j,k)=c2d*c1d 
     209   !                                !  = 31  F(i,j,k)=F(grid spacing and local velocity) 
     210   rn_ahm_0      =  40000.     !  horizontal laplacian eddy viscosity   [m2/s] 
     211   rn_ahm_b      =      0.     !  background eddy viscosity for ldf_iso [m2/s] 
     212   rn_bhm_0      = 1.e+12      !  horizontal bilaplacian eddy viscosity [m4/s] 
    187213/ 
    188214!----------------------------------------------------------------------- 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/cpp_ORCA2_LIM_PISCES.fcm

    r5766 r5778  
    1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_top key_pisces key_mpp_mpi key_iomput 
     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 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_cfg

    r5771 r5778  
    106106   ! 
    107107   !                       !  Coefficients: 
    108    nn_aht_ijk_t    = 21        !  space/time variation of eddy coef 
    109    !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     108   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     109   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    110110   !                                !   =  0           constant  
    111111   !                                !   = 10 F(k)      =ldf_c1d  
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/ORCA2_SAS_LIM/cpp_ORCA2_SAS_LIM.fcm

    r5759 r5778  
    1  bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt 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 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/SHARED/domain_def.xml

    r5426 r5778  
    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">  
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/SHARED/field_def.xml

    r5765 r5778  
    2828         <field id="soce"         long_name="salinity"            standard_name="sea_water_practical_salinity"      unit="1e-3"     grid_ref="grid_T_3D"/> 
    2929         <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"     /> 
     
    475479        <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" /> 
    476480      </field_group> 
    477            
     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 
    478487      <!-- scalar variables available with key_diaar5 --> 
    479488 
  • branches/2015/dev_r5721_CNRS9_NOC3_LDF/NEMOGCM/CONFIG/SHARED/namelist_ref

    r5771 r5778  
    714714                                 !  = 1, S-EOS   (simplified eos) 
    715715   ln_useCT    = .true.  ! use of Conservative Temp. ==> surface CT converted in Pot. Temp. in sbcssm 
    716    !                             ! 
     716                                 ! 
    717717   !                     ! S-EOS coefficients : 
    718    !                             !  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 
    719719   rn_a0       =  1.6550e-1      !  thermal expension coefficient (nn_eos= 1) 
    720720   rn_b0       =  7.6554e-1      !  saline  expension coefficient (nn_eos= 1) 
     
    745745&namtra_adv_mle !   mixed layer eddy parametrisation (Fox-Kemper param) 
    746746!----------------------------------------------------------------------- 
    747    ln_mle    = .true.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
     747   ln_mle    = .false.      ! (T) use the Mixed Layer Eddy (MLE) parameterisation 
    748748   rn_ce     = 0.06        ! magnitude of the MLE (typical value: 0.06 to 0.08) 
    749749   nn_mle    = 1           ! MLE type: =0 standard Fox-Kemper ; =1 new formulation 
     
    759759!---------------------------------------------------------------------------------- 
    760760   !                       !  Operator type: 
    761    ln_traldf_lap   =  .true.   !    laplacian operator 
     761   !                           !  no diffusion: set ln_traldf_lap=..._blp=F  
     762   ln_traldf_lap   =  .false.  !    laplacian operator 
    762763   ln_traldf_blp   =  .false.  !  bilaplacian operator 
    763764   !                       !  Direction of action: 
    764765   ln_traldf_lev   =  .false.  !  iso-level 
    765766   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
    766    ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     767   ln_traldf_iso   =  .false.  !  iso-neutral (standard operator) 
    767768   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
    768769   ! 
    769770   !                       !  iso-neutral options:         
    770    ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     771   ln_traldf_msc   =  .false.  !  Method of Stabilizing Correction (both operators) 
    771772   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
    772773   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     
    776777   !                       !  Coefficients: 
    777778   nn_aht_ijk_t    = 0         !  space/time variation of eddy coef 
    778    !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     779   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
    779780   !                                !   =  0           constant  
    780781   !                                !   = 10 F(k)      =ldf_c1d  
    781782   !                                !   = 20 F(i,j)    =ldf_c2d  
    782783   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
    783    !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    784    !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     784   !                                !   = 30 F(i,j,k)  =ldf_c2d * ldf_c1d 
     785   !                                !   = 31 F(i,j,k,t)=F(local velocity and grid-spacing) 
    785786   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
    786787   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     
    791792   ln_ldfeiv     =.false.   ! use eddy induced velocity parameterization 
    792793   ln_ldfeiv_dia =.false.   ! diagnose eiv stream function and velocities 
    793    rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
    794    nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
     794   rn_aeiv_0     = 2000.    ! eddy induced velocity coefficient   [m2/s] 
     795   nn_aei_ijk_t  = 21       ! space/time variation of the eiv coeficient 
    795796   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
    796797   !                                !   =  0           constant  
     
    849850   ln_dynvor_ens = .false. !  energy conserving scheme 
    850851   ln_dynvor_mix = .false. !  mixed scheme 
    851    ln_dynvor_een = .true.  !  energy & enstrophy scheme 
    852    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) 
    853855/ 
    854856!----------------------------------------------------------------------- 
     
    875877!----------------------------------------------------------------------- 
    876878   !                       !  Type of the operator : 
    877    ln_dynldf_lap    =  .true.   !  laplacian operator 
    878    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 
    879882   !                       !  Direction of action  : 
    880    ln_dynldf_level  =  .false.  !  iso-level 
    881    ln_dynldf_hor    =  .true.   !  horizontal (geopotential)            (require "key_ldfslp" in s-coord.) 
    882    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 
    883886   !                       !  Coefficient 
    884    rn_ahm_0_lap     = 40000.    !  horizontal laplacian eddy viscosity   [m2/s] 
    885    rn_ahmb_0        =     0.    !  background eddy viscosity for ldf_iso [m2/s] 
    886    rn_ahm_0_blp     =     0.    !  horizontal bilaplacian eddy viscosity [m4/s] 
    887    rn_cmsmag_1      =     3.    !  constant in laplacian Smagorinsky viscosity 
    888    rn_cmsmag_2      =     3     !  constant in bilaplacian Smagorinsky viscosity 
    889    rn_cmsh          =     1.    !  1 or 0 , if 0 -use only shear for Smagorinsky viscosity 
    890    rn_ahm_m_blp     =    -1.e12 !  upper limit for bilap  abs(ahm) < min( dx^4/128rdt, rn_ahm_m_blp) 
    891    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) 
    892900/ 
    893901 
Note: See TracChangeset for help on using the changeset viewer.