Changeset 5836 for trunk


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

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

Location:
trunk/NEMOGCM
Files:
31 deleted
217 edited
8 copied

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!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/NEMO/LIM_SRC_2/limdmp_2.F90

    r4624 r5836  
    7171         CALL fld_read( kt, nn_fsbc, sf_icedmp ) 
    7272         ! 
    73 !CDIR COLLAPSE 
    7473         hicif(:,:) = MAX( 0._wp,                     &        ! h >= 0         avoid spurious out of physical range 
    7574            &         hicif(:,:) - rdt_ice * resto_ice(:,:,1) * ( hicif(:,:) - sf_icedmp(jp_hicif)%fnow(:,:,1) )  )  
    76 !CDIR COLLAPSE 
    7775         frld (:,:) = MAX( 0._wp, MIN( 1._wp,         &        ! 0<= frld<=1    values which blow the run up 
    7876            &         frld (:,:) - rdt_ice * resto_ice(:,:,1) * ( frld (:,:) - sf_icedmp(jp_frld )%fnow(:,:,1) )  )  ) 
  • trunk/NEMOGCM/NEMO/LIM_SRC_2/limrhg_2.F90

    r5123 r5836  
    160160      !------------------------------------------------------------------- 
    161161 
    162 !CDIR NOVERRCHK 
    163162      DO jj = k_j1 , k_jpj-1 
    164 !CDIR NOVERRCHK 
    165163         DO ji = 1 , jpi 
    166164            ! only the sinus changes its sign with the hemisphere 
     
    245243         ! Computation of free drift field for free slip boundary conditions. 
    246244 
    247 !CDIR NOVERRCHK 
    248245         DO jj = k_j1, k_jpj-1 
    249 !CDIR NOVERRCHK 
    250246            DO ji = 1, fs_jpim1 
    251247               !- Rate of strain tensor. 
     
    401397iflag:   DO jter = 1 , nbitdr                                   !    Relaxation    ! 
    402398            !                                                   ! ================ ! 
    403 !CDIR NOVERRCHK 
    404399            DO jj = k_j1+1, k_jpj-1 
    405 !CDIR NOVERRCHK 
    406400               DO ji = 2, fs_jpim1   ! NO vector opt. 
    407401                  ! 
  • trunk/NEMOGCM/NEMO/LIM_SRC_2/limsbc_2.F90

    r5407 r5836  
    319319         ! 
    320320         IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN     !==  Ice time-step only  ==! (i.e. surface module time-step) 
    321 !CDIR NOVERRCHK 
     321            ! 
    322322            DO jj = 1, jpj                               !* modulus of ice-ocean relative velocity at I-point 
    323 !CDIR NOVERRCHK 
    324323               DO ji = 1, jpi 
    325324                  zu_i  = u_ice(ji,jj) - u_oce(ji,jj)                   ! ice-ocean relative velocity at I-point 
     
    328327               END DO 
    329328            END DO 
    330 !CDIR NOVERRCHK 
    331329            DO jj = 1, jpjm1                             !* update the modulus of stress at ocean surface (T-point) 
    332 !CDIR NOVERRCHK 
    333330               DO ji = 1, jpim1   ! NO vector opt. 
    334331                  !                                               ! modulus of U_ice-U_oce at T-point 
     
    383380         ! 
    384381         IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN     !==  Ice time-step only  ==! (i.e. surface module time-step) 
    385 !CDIR NOVERRCHK 
     382            ! 
    386383            DO jj = 2, jpjm1                          !* modulus of the ice-ocean velocity at T-point 
    387 !CDIR NOVERRCHK 
    388384               DO ji = fs_2, fs_jpim1 
    389385                  zu_t  = u_ice(ji,jj) + u_ice(ji-1,jj) - u_oce(ji,jj) - u_oce(ji-1,jj)   ! 2*(U_ice-U_oce) at T-point 
  • trunk/NEMOGCM/NEMO/LIM_SRC_2/limthd_2.F90

    r5407 r5836  
    196196      !-------------------------------------------------------------------------- 
    197197 
    198       !CDIR NOVERRCHK 
    199198      DO jj = 1, jpj 
    200          !CDIR NOVERRCHK 
    201199         DO ji = 1, jpi 
    202200            zthsnice       = hsnif(ji,jj) + hicif(ji,jj) 
  • trunk/NEMOGCM/NEMO/LIM_SRC_2/limthd_lac_2.F90

    r3625 r5836  
    134134      !--------------------------------------------------------------------- 
    135135       
    136 !CDIR NOVERRCHK 
    137136      DO ji = kideb , kiut 
    138137         iicefr       = 1 - MAX( 0, INT( SIGN( 1.5 * zone , zfrl_old(ji) - 1.0 + epsi13 ) ) ) 
  • trunk/NEMOGCM/NEMO/LIM_SRC_3/limadv.F90

    r5429 r5836  
    9797 
    9898      !  Initialize volumes of boxes  (=area if adv_x first called, =psm otherwise)                                      
    99       psm (:,:)  = MAX( pcrh * e12t(:,:) + ( 1.0 - pcrh ) * psm(:,:) , epsi20 ) 
     99      psm (:,:)  = MAX( pcrh * e1e2t(:,:) + ( 1.0 - pcrh ) * psm(:,:) , epsi20 ) 
    100100 
    101101      !  Calculate fluxes and moments between boxes i<-->i+1               
     
    282282 
    283283      !  Initialize volumes of boxes (=area if adv_x first called, =psm otherwise) 
    284       psm(:,:)  = MAX(  pcrh * e12t(:,:) + ( 1.0 - pcrh ) * psm(:,:) , epsi20  ) 
     284      psm(:,:)  = MAX(  pcrh * e1e2t(:,:) + ( 1.0 - pcrh ) * psm(:,:) , epsi20  ) 
    285285 
    286286      !  Calculate fluxes and moments between boxes j<-->j+1               
  • trunk/NEMOGCM/NEMO/LIM_SRC_3/limcons.F90

    r5183 r5836  
    185185         zfs_b  = glob_sum(  ( sfx_bri(:,:) + sfx_bog(:,:) + sfx_bom(:,:) + sfx_sum(:,:) + sfx_sni(:,:) +  & 
    186186            &                  sfx_opw(:,:) + sfx_res(:,:) + sfx_dyn(:,:)                                  & 
    187             &                ) *  e12t(:,:) * tmask(:,:,1) * zconv ) 
     187            &                ) *  e1e2t(:,:) * tmask(:,:,1) * zconv ) 
    188188 
    189189         ! water flux 
    190190         zfw_b  = glob_sum( -( wfx_bog(:,:) + wfx_bom(:,:) + wfx_sum(:,:) + wfx_sni(:,:) + wfx_opw(:,:) +  & 
    191191            &                  wfx_res(:,:) + wfx_dyn(:,:) + wfx_snw(:,:) + wfx_sub(:,:) + wfx_spr(:,:)    & 
    192             &                ) *  e12t(:,:) * tmask(:,:,1) * zconv ) 
     192            &                ) *  e1e2t(:,:) * tmask(:,:,1) * zconv ) 
    193193 
    194194         ! heat flux 
    195195         zft_b  = glob_sum(  ( hfx_sum(:,:) + hfx_bom(:,:) + hfx_bog(:,:) + hfx_dif(:,:) + hfx_opw(:,:) + hfx_snw(:,:)  &  
    196196            &                - hfx_thd(:,:) - hfx_dyn(:,:) - hfx_res(:,:) - hfx_sub(:,:) - hfx_spr(:,:)   & 
    197             &                ) *  e12t(:,:) * tmask(:,:,1) * zconv ) 
    198  
    199          zvi_b  = glob_sum( SUM( v_i * rhoic + v_s * rhosn, dim=3 ) * e12t * tmask(:,:,1) * zconv ) 
    200  
    201          zsmv_b = glob_sum( SUM( smv_i * rhoic            , dim=3 ) * e12t * tmask(:,:,1) * zconv ) 
     197            &                ) *  e1e2t(:,:) * tmask(:,:,1) * zconv ) 
     198 
     199         zvi_b  = glob_sum( SUM( v_i * rhoic + v_s * rhosn, dim=3 ) * e1e2t * tmask(:,:,1) * zconv ) 
     200 
     201         zsmv_b = glob_sum( SUM( smv_i * rhoic            , dim=3 ) * e1e2t * tmask(:,:,1) * zconv ) 
    202202 
    203203         zei_b  = glob_sum( ( SUM( SUM( e_i(:,:,1:nlay_i,:), dim=4 ), dim=3 ) +  & 
    204204            &                 SUM( SUM( e_s(:,:,1:nlay_s,:), dim=4 ), dim=3 )    & 
    205                             ) * e12t * tmask(:,:,1) * zconv ) 
     205                            ) * e1e2t * tmask(:,:,1) * zconv ) 
    206206 
    207207      ELSEIF( icount == 1 ) THEN 
     
    210210         zfs  = glob_sum(  ( sfx_bri(:,:) + sfx_bog(:,:) + sfx_bom(:,:) + sfx_sum(:,:) + sfx_sni(:,:) +  & 
    211211            &                sfx_opw(:,:) + sfx_res(:,:) + sfx_dyn(:,:)                                  &  
    212             &              ) * e12t(:,:) * tmask(:,:,1) * zconv ) - zfs_b 
     212            &              ) * e1e2t(:,:) * tmask(:,:,1) * zconv ) - zfs_b 
    213213 
    214214         ! water flux 
    215215         zfw  = glob_sum( -( wfx_bog(:,:) + wfx_bom(:,:) + wfx_sum(:,:) + wfx_sni(:,:) + wfx_opw(:,:) +  & 
    216216            &                wfx_res(:,:) + wfx_dyn(:,:) + wfx_snw(:,:) + wfx_sub(:,:) + wfx_spr(:,:)    & 
    217             &              ) * e12t(:,:) * tmask(:,:,1) * zconv ) - zfw_b 
     217            &              ) * e1e2t(:,:) * tmask(:,:,1) * zconv ) - zfw_b 
    218218 
    219219         ! heat flux 
    220220         zft  = glob_sum(  ( hfx_sum(:,:) + hfx_bom(:,:) + hfx_bog(:,:) + hfx_dif(:,:) + hfx_opw(:,:) + hfx_snw(:,:)  &  
    221221            &              - hfx_thd(:,:) - hfx_dyn(:,:) - hfx_res(:,:) - hfx_sub(:,:) - hfx_spr(:,:)   & 
    222             &              ) * e12t(:,:) * tmask(:,:,1) * zconv ) - zft_b 
     222            &              ) * e1e2t(:,:) * tmask(:,:,1) * zconv ) - zft_b 
    223223  
    224224         ! outputs 
    225225         zvi  = ( ( glob_sum( SUM( v_i * rhoic + v_s * rhosn, dim=3 )  & 
    226             &                    * e12t * tmask(:,:,1) * zconv ) - zvi_b ) * r1_rdtice - zfw ) * rday 
     226            &                    * e1e2t * tmask(:,:,1) * zconv ) - zvi_b ) * r1_rdtice - zfw ) * rday 
    227227 
    228228         zsmv = ( ( glob_sum( SUM( smv_i * rhoic            , dim=3 )  & 
    229             &                    * e12t * tmask(:,:,1) * zconv ) - zsmv_b ) * r1_rdtice + zfs ) * rday 
     229            &                    * e1e2t * tmask(:,:,1) * zconv ) - zsmv_b ) * r1_rdtice + zfs ) * rday 
    230230 
    231231         zei  =   glob_sum( ( SUM( SUM( e_i(:,:,1:nlay_i,:), dim=4 ), dim=3 ) +  & 
    232232            &                 SUM( SUM( e_s(:,:,1:nlay_s,:), dim=4 ), dim=3 )    & 
    233             &                ) * e12t * tmask(:,:,1) * zconv ) * r1_rdtice - zei_b * r1_rdtice + zft 
     233            &                ) * e1e2t * tmask(:,:,1) * zconv ) * r1_rdtice - zei_b * r1_rdtice + zft 
    234234 
    235235         ! zvtrp and zetrp must be close to 0 if the advection scheme is conservative 
    236          zvtrp = glob_sum( ( diag_trp_vi * rhoic + diag_trp_vs * rhosn ) * e12t * tmask(:,:,1) * zconv ) * rday  
    237          zetrp = glob_sum( ( diag_trp_ei         + diag_trp_es         ) * e12t * tmask(:,:,1) * zconv ) 
     236         zvtrp = glob_sum( ( diag_trp_vi * rhoic + diag_trp_vs * rhosn ) * e1e2t * tmask(:,:,1) * zconv ) * rday  
     237         zetrp = glob_sum( ( diag_trp_ei         + diag_trp_es         ) * e1e2t * tmask(:,:,1) * zconv ) 
    238238 
    239239         zvmin = glob_min( v_i ) 
     
    242242 
    243243         ! set threshold values and calculate the ice area (+epsi10 to set a threshold > 0 when there is no ice)  
    244          zarea   = glob_sum( SUM( a_i + epsi10, dim=3 ) * e12t * zconv ) ! in 1.e9 m2 
     244         zarea   = glob_sum( SUM( a_i + epsi10, dim=3 ) * e1e2t * zconv ) ! in 1.e9 m2 
    245245         zv_sill = zarea * 2.5e-5 
    246246         zs_sill = zarea * 25.e-5 
     
    286286#if ! defined key_bdy 
    287287      ! heat flux 
    288       zhfx  = glob_sum( ( hfx_in - hfx_out - diag_heat - diag_trp_ei - diag_trp_es - hfx_sub ) * e12t * tmask(:,:,1) * zconv )  
     288      zhfx  = glob_sum( ( hfx_in - hfx_out - diag_heat - diag_trp_ei - diag_trp_es - hfx_sub ) * e1e2t * tmask(:,:,1) * zconv )  
    289289      ! salt flux 
    290       zsfx  = glob_sum( ( sfx + diag_smvi ) * e12t * tmask(:,:,1) * zconv ) * rday 
     290      zsfx  = glob_sum( ( sfx + diag_smvi ) * e1e2t * tmask(:,:,1) * zconv ) * rday 
    291291      ! water flux 
    292       zvfx  = glob_sum( ( wfx_ice + wfx_snw + wfx_spr + wfx_sub + diag_vice + diag_vsnw ) * e12t * tmask(:,:,1) * zconv ) * rday 
     292      zvfx  = glob_sum( ( wfx_ice + wfx_snw + wfx_spr + wfx_sub + diag_vice + diag_vsnw ) * e1e2t * tmask(:,:,1) * zconv ) * rday 
    293293 
    294294      ! set threshold values and calculate the ice area (+epsi10 to set a threshold > 0 when there is no ice)  
    295       zarea   = glob_sum( SUM( a_i + epsi10, dim=3 ) * e12t * zconv ) ! in 1.e9 m2 
     295      zarea   = glob_sum( SUM( a_i + epsi10, dim=3 ) * e1e2t * zconv ) ! in 1.e9 m2 
    296296      zv_sill = zarea * 2.5e-5 
    297297      zs_sill = zarea * 25.e-5 
  • trunk/NEMOGCM/NEMO/LIM_SRC_3/limctl.F90

    r5167 r5836  
    306306               WRITE(numout,*) ' - Cell values ' 
    307307               WRITE(numout,*) '   ~~~~~~~~~~~ ' 
    308                WRITE(numout,*) ' cell area     : ', e12t(ji,jj) 
     308               WRITE(numout,*) ' cell area     : ', e1e2t(ji,jj) 
    309309               WRITE(numout,*) ' at_i          : ', at_i(ji,jj)        
    310310               WRITE(numout,*) ' vt_i          : ', vt_i(ji,jj)        
     
    350350               WRITE(numout,*) ' - Cell values ' 
    351351               WRITE(numout,*) '   ~~~~~~~~~~~ ' 
    352                WRITE(numout,*) ' cell area     : ', e12t(ji,jj) 
     352               WRITE(numout,*) ' cell area     : ', e1e2t(ji,jj) 
    353353               WRITE(numout,*) ' at_i          : ', at_i(ji,jj)        
    354354               WRITE(numout,*) ' vt_i          : ', vt_i(ji,jj)        
  • trunk/NEMOGCM/NEMO/LIM_SRC_3/limdiahsb.F90

    r5215 r5836  
    7171 
    7272      ! 1/area 
    73       z1_area = 1._wp / MAX( glob_sum( e12t(:,:) * tmask(:,:,1) ), epsi06 ) 
    74  
    75       rswitch = MAX( 0._wp , SIGN( 1._wp , glob_sum( e12t(:,:) * tmask(:,:,1) ) - epsi06 ) ) 
     73      z1_area = 1._wp / MAX( glob_sum( e1e2t(:,:) * tmask(:,:,1) ), epsi06 ) 
     74 
     75      rswitch = MAX( 0._wp , SIGN( 1._wp , glob_sum( e1e2t(:,:) * tmask(:,:,1) ) - epsi06 ) ) 
    7676      ! ----------------------- ! 
    7777      ! 1 -  Content variations ! 
    7878      ! ----------------------- ! 
    79       zbg_ivo = glob_sum( vt_i(:,:) * e12t(:,:) * tmask(:,:,1) ) ! volume ice  
    80       zbg_svo = glob_sum( vt_s(:,:) * e12t(:,:) * tmask(:,:,1) ) ! volume snow 
    81       zbg_are = glob_sum( at_i(:,:) * e12t(:,:) * tmask(:,:,1) ) ! area 
    82       zbg_sal = glob_sum( SUM( smv_i(:,:,:), dim=3 ) * e12t(:,:) * tmask(:,:,1) )       ! mean salt content 
    83       zbg_tem = glob_sum( ( tm_i(:,:) - rt0 ) * vt_i(:,:) * e12t(:,:) * tmask(:,:,1) )  ! mean temp content 
    84  
    85       !zbg_ihc = glob_sum( et_i(:,:) * e12t(:,:) * tmask(:,:,1) ) / MAX( zbg_ivo,epsi06 ) ! ice heat content 
    86       !zbg_shc = glob_sum( et_s(:,:) * e12t(:,:) * tmask(:,:,1) ) / MAX( zbg_svo,epsi06 ) ! snow heat content 
     79      zbg_ivo = glob_sum( vt_i(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! volume ice  
     80      zbg_svo = glob_sum( vt_s(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! volume snow 
     81      zbg_are = glob_sum( at_i(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! area 
     82      zbg_sal = glob_sum( SUM( smv_i(:,:,:), dim=3 ) * e1e2t(:,:) * tmask(:,:,1) )       ! mean salt content 
     83      zbg_tem = glob_sum( ( tm_i(:,:) - rt0 ) * vt_i(:,:) * e1e2t(:,:) * tmask(:,:,1) )  ! mean temp content 
     84 
     85      !zbg_ihc = glob_sum( et_i(:,:) * e1e2t(:,:) * tmask(:,:,1) ) / MAX( zbg_ivo,epsi06 ) ! ice heat content 
     86      !zbg_shc = glob_sum( et_s(:,:) * e1e2t(:,:) * tmask(:,:,1) ) / MAX( zbg_svo,epsi06 ) ! snow heat content 
    8787 
    8888      ! Volume 
    8989      ztmp = rswitch * z1_area * r1_rau0 * rday 
    90       zbg_vfx     = ztmp * glob_sum(     emp(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    91       zbg_vfx_bog = ztmp * glob_sum( wfx_bog(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    92       zbg_vfx_opw = ztmp * glob_sum( wfx_opw(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    93       zbg_vfx_sni = ztmp * glob_sum( wfx_sni(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    94       zbg_vfx_dyn = ztmp * glob_sum( wfx_dyn(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    95       zbg_vfx_bom = ztmp * glob_sum( wfx_bom(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    96       zbg_vfx_sum = ztmp * glob_sum( wfx_sum(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    97       zbg_vfx_res = ztmp * glob_sum( wfx_res(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    98       zbg_vfx_spr = ztmp * glob_sum( wfx_spr(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    99       zbg_vfx_snw = ztmp * glob_sum( wfx_snw(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    100       zbg_vfx_sub = ztmp * glob_sum( wfx_sub(:,:) * e12t(:,:) * tmask(:,:,1) ) 
     90      zbg_vfx     = ztmp * glob_sum(     emp(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     91      zbg_vfx_bog = ztmp * glob_sum( wfx_bog(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     92      zbg_vfx_opw = ztmp * glob_sum( wfx_opw(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     93      zbg_vfx_sni = ztmp * glob_sum( wfx_sni(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     94      zbg_vfx_dyn = ztmp * glob_sum( wfx_dyn(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     95      zbg_vfx_bom = ztmp * glob_sum( wfx_bom(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     96      zbg_vfx_sum = ztmp * glob_sum( wfx_sum(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     97      zbg_vfx_res = ztmp * glob_sum( wfx_res(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     98      zbg_vfx_spr = ztmp * glob_sum( wfx_spr(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     99      zbg_vfx_snw = ztmp * glob_sum( wfx_snw(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     100      zbg_vfx_sub = ztmp * glob_sum( wfx_sub(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
    101101 
    102102      ! Salt 
    103       zbg_sfx     = ztmp * glob_sum(     sfx(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    104       zbg_sfx_bri = ztmp * glob_sum( sfx_bri(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    105       zbg_sfx_res = ztmp * glob_sum( sfx_res(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    106       zbg_sfx_dyn = ztmp * glob_sum( sfx_dyn(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    107  
    108       zbg_sfx_bog = ztmp * glob_sum( sfx_bog(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    109       zbg_sfx_opw = ztmp * glob_sum( sfx_opw(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    110       zbg_sfx_sni = ztmp * glob_sum( sfx_sni(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    111       zbg_sfx_bom = ztmp * glob_sum( sfx_bom(:,:) * e12t(:,:) * tmask(:,:,1) ) 
    112       zbg_sfx_sum = ztmp * glob_sum( sfx_sum(:,:) * e12t(:,:) * tmask(:,:,1) ) 
     103      zbg_sfx     = ztmp * glob_sum(     sfx(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     104      zbg_sfx_bri = ztmp * glob_sum( sfx_bri(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     105      zbg_sfx_res = ztmp * glob_sum( sfx_res(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     106      zbg_sfx_dyn = ztmp * glob_sum( sfx_dyn(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     107 
     108      zbg_sfx_bog = ztmp * glob_sum( sfx_bog(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     109      zbg_sfx_opw = ztmp * glob_sum( sfx_opw(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     110      zbg_sfx_sni = ztmp * glob_sum( sfx_sni(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     111      zbg_sfx_bom = ztmp * glob_sum( sfx_bom(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
     112      zbg_sfx_sum = ztmp * glob_sum( sfx_sum(:,:) * e1e2t(:,:) * tmask(:,:,1) ) 
    113113 
    114114      ! Heat budget 
    115       zbg_ihc      = glob_sum( et_i(:,:) * e12t(:,:) * 1.e-20 ) ! ice heat content  [1.e20 J] 
    116       zbg_shc      = glob_sum( et_s(:,:) * e12t(:,:) * 1.e-20 ) ! snow heat content [1.e20 J] 
    117       zbg_hfx_dhc  = glob_sum( diag_heat(:,:) * e12t(:,:) * tmask(:,:,1) ) ! [in W] 
    118       zbg_hfx_spr  = glob_sum( hfx_spr(:,:) * e12t(:,:) * tmask(:,:,1) ) ! [in W] 
    119  
    120       zbg_hfx_thd  = glob_sum( hfx_thd(:,:) * e12t(:,:) * tmask(:,:,1) ) ! [in W] 
    121       zbg_hfx_dyn  = glob_sum( hfx_dyn(:,:) * e12t(:,:) * tmask(:,:,1) ) ! [in W] 
    122       zbg_hfx_res  = glob_sum( hfx_res(:,:) * e12t(:,:) * tmask(:,:,1) ) ! [in W] 
    123       zbg_hfx_sub  = glob_sum( hfx_sub(:,:) * e12t(:,:) * tmask(:,:,1) ) ! [in W] 
    124       zbg_hfx_snw  = glob_sum( hfx_snw(:,:) * e12t(:,:) * tmask(:,:,1) ) ! [in W] 
    125       zbg_hfx_sum  = glob_sum( hfx_sum(:,:) * e12t(:,:) * tmask(:,:,1) ) ! [in W] 
    126       zbg_hfx_bom  = glob_sum( hfx_bom(:,:) * e12t(:,:) * tmask(:,:,1) ) ! [in W] 
    127       zbg_hfx_bog  = glob_sum( hfx_bog(:,:) * e12t(:,:) * tmask(:,:,1) ) ! [in W] 
    128       zbg_hfx_dif  = glob_sum( hfx_dif(:,:) * e12t(:,:) * tmask(:,:,1) ) ! [in W] 
    129       zbg_hfx_opw  = glob_sum( hfx_opw(:,:) * e12t(:,:) * tmask(:,:,1) ) ! [in W] 
    130       zbg_hfx_out  = glob_sum( hfx_out(:,:) * e12t(:,:) * tmask(:,:,1) ) ! [in W] 
    131       zbg_hfx_in   = glob_sum(  hfx_in(:,:) * e12t(:,:) * tmask(:,:,1) ) ! [in W] 
     115      zbg_ihc      = glob_sum( et_i(:,:) * e1e2t(:,:) ) * 1.e-20 ! ice heat content  [1.e20 J] 
     116      zbg_shc      = glob_sum( et_s(:,:) * e1e2t(:,:) ) * 1.e-20 ! snow heat content [1.e20 J] 
     117      zbg_hfx_dhc  = glob_sum( diag_heat(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 
     118      zbg_hfx_spr  = glob_sum( hfx_spr(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 
     119 
     120      zbg_hfx_thd  = glob_sum( hfx_thd(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 
     121      zbg_hfx_dyn  = glob_sum( hfx_dyn(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 
     122      zbg_hfx_res  = glob_sum( hfx_res(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 
     123      zbg_hfx_sub  = glob_sum( hfx_sub(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 
     124      zbg_hfx_snw  = glob_sum( hfx_snw(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 
     125      zbg_hfx_sum  = glob_sum( hfx_sum(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 
     126      zbg_hfx_bom  = glob_sum( hfx_bom(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 
     127      zbg_hfx_bog  = glob_sum( hfx_bog(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 
     128      zbg_hfx_dif  = glob_sum( hfx_dif(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 
     129      zbg_hfx_opw  = glob_sum( hfx_opw(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 
     130      zbg_hfx_out  = glob_sum( hfx_out(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 
     131      zbg_hfx_in   = glob_sum(  hfx_in(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! [in W] 
    132132     
    133133      ! --------------------------------------------- ! 
    134134      ! 2 - Trends due to forcing and ice growth/melt ! 
    135135      ! --------------------------------------------- ! 
    136       z_frc_vol = r1_rau0 * glob_sum( - emp(:,:) * e12t(:,:) * tmask(:,:,1) ) ! volume fluxes 
    137       z_frc_sal = r1_rau0 * glob_sum(   sfx(:,:) * e12t(:,:) * tmask(:,:,1) ) ! salt fluxes 
     136      z_frc_vol = r1_rau0 * glob_sum( - emp(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! volume fluxes 
     137      z_frc_sal = r1_rau0 * glob_sum(   sfx(:,:) * e1e2t(:,:) * tmask(:,:,1) ) ! salt fluxes 
    138138      z_bg_grme = glob_sum( - ( wfx_bog(:,:) + wfx_opw(:,:) + wfx_sni(:,:) + wfx_dyn(:,:) + & 
    139139                          &     wfx_bom(:,:) + wfx_sum(:,:) + wfx_res(:,:) + wfx_snw(:,:) + & 
    140                           &     wfx_sub(:,:) ) * e12t(:,:) * tmask(:,:,1) ) ! volume fluxes 
     140                          &     wfx_sub(:,:) ) * e1e2t(:,:) * tmask(:,:,1) ) ! volume fluxes 
    141141      ! 
    142142      frc_vol  = frc_vol  + z_frc_vol  * rdt_ice 
  • trunk/NEMOGCM/NEMO/LIM_SRC_3/limdyn.F90

    r5123 r5836  
    191191         CALL prt_ctl(tab2d_1=delta_i   , clinfo1=' lim_dyn  : delta_i   :') 
    192192         CALL prt_ctl(tab2d_1=strength  , clinfo1=' lim_dyn  : strength  :') 
    193          CALL prt_ctl(tab2d_1=e12t      , clinfo1=' lim_dyn  : cell area :') 
     193         CALL prt_ctl(tab2d_1=e1e2t     , clinfo1=' lim_dyn  : cell area :') 
    194194         CALL prt_ctl(tab2d_1=at_i      , clinfo1=' lim_dyn  : at_i      :') 
    195195         CALL prt_ctl(tab2d_1=vt_i      , clinfo1=' lim_dyn  : vt_i      :') 
  • trunk/NEMOGCM/NEMO/LIM_SRC_3/limhdf.F90

    r5429 r5836  
    7676         DO jj = 2, jpjm1   
    7777            DO ji = fs_2 , fs_jpim1   ! vector opt. 
    78                efact(ji,jj) = ( e2u(ji,jj) + e2u(ji-1,jj) + e1v(ji,jj) + e1v(ji,jj-1) ) * r1_e12t(ji,jj) 
     78               efact(ji,jj) = ( e2u(ji,jj) + e2u(ji-1,jj) + e1v(ji,jj) + e1v(ji,jj-1) ) * r1_e1e2t(ji,jj) 
    7979            END DO 
    8080         END DO 
     
    107107         DO jj= 2, jpjm1                                 ! diffusive trend : divergence of the fluxes 
    108108            DO ji = fs_2 , fs_jpim1   ! vector opt.  
    109                zdiv(ji,jj) = ( zflu(ji,jj) - zflu(ji-1,jj) + zflv(ji,jj) - zflv(ji,jj-1) ) * r1_e12t(ji,jj) 
     109               zdiv(ji,jj) = ( zflu(ji,jj) - zflu(ji-1,jj) + zflv(ji,jj) - zflv(ji,jj-1) ) * r1_e1e2t(ji,jj) 
    110110            END DO 
    111111         END DO 
     
    149149      DO jj= 2, jpjm1                                 ! diffusive trend : divergence of the fluxes 
    150150         DO ji = fs_2 , fs_jpim1   ! vector opt.  
    151             zdiv(ji,jj) = ( zflu(ji,jj) - zflu(ji-1,jj) + zflv(ji,jj) - zflv(ji,jj-1) ) * r1_e12t(ji,jj) 
     151            zdiv(ji,jj) = ( zflu(ji,jj) - zflu(ji-1,jj) + zflv(ji,jj) - zflv(ji,jj-1) ) * r1_e1e2t(ji,jj) 
    152152            ptab(ji,jj) = ztab0(ji,jj) + 0.5 * ( zdiv(ji,jj) + zdiv0(ji,jj) ) 
    153153         END DO 
  • trunk/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90

    r5202 r5836  
    377377         CALL prt_ctl_info(' - Cell values : ') 
    378378         CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
    379          CALL prt_ctl(tab2d_1=e12t , clinfo1=' lim_itd_me  : cell area :') 
     379         CALL prt_ctl(tab2d_1=e1e2t, clinfo1=' lim_itd_me  : cell area :') 
    380380         CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_itd_me  : at_i      :') 
    381381         CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_itd_me  : vt_i      :') 
  • trunk/NEMOGCM/NEMO/LIM_SRC_3/limrhg.F90

    r5429 r5836  
    355355               divu_i(ji,jj) = (  e2u(ji,jj) * u_ice(ji,jj) - e2u(ji-1,jj) * u_ice(ji-1,jj)   & 
    356356                  &             + e1v(ji,jj) * v_ice(ji,jj) - e1v(ji,jj-1) * v_ice(ji,jj-1)   & 
    357                   &            ) * r1_e12t(ji,jj) 
     357                  &            ) * r1_e1e2t(ji,jj) 
    358358 
    359359               zdt(ji,jj) = ( ( u_ice(ji,jj) * r1_e2u(ji,jj) - u_ice(ji-1,jj) * r1_e2u(ji-1,jj) ) * e2t(ji,jj) * e2t(ji,jj)   & 
    360360                  &         - ( v_ice(ji,jj) * r1_e1v(ji,jj) - v_ice(ji,jj-1) * r1_e1v(ji,jj-1) ) * e1t(ji,jj) * e1t(ji,jj)   & 
    361                   &         ) * r1_e12t(ji,jj) 
     361                  &         ) * r1_e1e2t(ji,jj) 
    362362 
    363363               ! 
    364364               zds(ji,jj) = ( ( u_ice(ji,jj+1) * r1_e1u(ji,jj+1) - u_ice(ji,jj) * r1_e1u(ji,jj) ) * e1f(ji,jj) * e1f(ji,jj)   & 
    365365                  &         + ( v_ice(ji+1,jj) * r1_e2v(ji+1,jj) - v_ice(ji,jj) * r1_e2v(ji,jj) ) * e2f(ji,jj) * e2f(ji,jj)   & 
    366                   &         ) * r1_e12f(ji,jj) * ( 2._wp - fmask(ji,jj,1) )   & 
     366                  &         ) * r1_e1e2f(ji,jj) * ( 2._wp - fmask(ji,jj,1) )   & 
    367367                  &         * zmask(ji,jj) * zmask(ji,jj+1) * zmask(ji+1,jj) * zmask(ji+1,jj+1) 
    368368 
     
    386386               zdst          = ( e2u(ji,jj) * v_ice1(ji,jj) - e2u(ji-1,jj  ) * v_ice1(ji-1,jj  )   & 
    387387                  &            + e1v(ji,jj) * u_ice2(ji,jj) - e1v(ji  ,jj-1) * u_ice2(ji  ,jj-1)   & 
    388                   &            ) * r1_e12t(ji,jj) 
     388                  &            ) * r1_e1e2t(ji,jj) 
    389389 
    390390               delta          = SQRT( divu_i(ji,jj)**2 + ( zdt(ji,jj)**2 + zdst**2 ) * usecc2 )   
     
    394394               zddc  = (  ( v_ice1(ji,jj+1) * r1_e1u(ji,jj+1) - v_ice1(ji,jj) * r1_e1u(ji,jj) ) * e1f(ji,jj) * e1f(ji,jj)  & 
    395395                  &     + ( u_ice2(ji+1,jj) * r1_e2v(ji+1,jj) - u_ice2(ji,jj) * r1_e2v(ji,jj) ) * e2f(ji,jj) * e2f(ji,jj)  & 
    396                   &    ) * r1_e12f(ji,jj) 
     396                  &    ) * r1_e1e2f(ji,jj) 
    397397 
    398398               zdtc  = (- ( v_ice1(ji,jj+1) * r1_e1u(ji,jj+1) - v_ice1(ji,jj) * r1_e1u(ji,jj) ) * e1f(ji,jj) * e1f(ji,jj)  & 
    399399                  &     + ( u_ice2(ji+1,jj) * r1_e2v(ji+1,jj) - u_ice2(ji,jj) * r1_e2v(ji,jj) ) * e2f(ji,jj) * e2f(ji,jj)  & 
    400                   &    ) * r1_e12f(ji,jj) 
     400                  &    ) * r1_e1e2f(ji,jj) 
    401401 
    402402               zddc = SQRT( zddc**2 + ( zdtc**2 + zds(ji,jj)**2 ) * usecc2 ) + rn_creepl 
     
    423423                  &             + ( zs2(ji+1,jj) * e2t(ji+1,jj)**2 - zs2(ji,jj) * e2t(ji,jj)**2 ) * r1_e2u(ji,jj)          & 
    424424                  &             + 2.0 * ( zs12(ji,jj) * e1f(ji,jj)**2 - zs12(ji,jj-1) * e1f(ji,jj-1)**2 ) * r1_e1u(ji,jj)  & 
    425                   &                ) * r1_e12u(ji,jj) 
     425                  &                ) * r1_e1e2u(ji,jj) 
    426426               ! contribution of zs1, zs2 and zs12 to zf2 
    427427               zf2(ji,jj) = 0.5 * ( ( zs1(ji,jj+1) - zs1(ji,jj) ) * e1v(ji,jj)  & 
    428428                  &             - ( zs2(ji,jj+1) * e1t(ji,jj+1)**2 - zs2(ji,jj) * e1t(ji,jj)**2 ) * r1_e1v(ji,jj)          & 
    429429                  &             + 2.0 * ( zs12(ji,jj) * e2f(ji,jj)**2 - zs12(ji-1,jj) * e2f(ji-1,jj)**2 ) * r1_e2v(ji,jj)  & 
    430                   &               )  * r1_e12v(ji,jj) 
     430                  &               )  * r1_e1e2v(ji,jj) 
    431431            END DO 
    432432         END DO 
     
    607607               divu_i(ji,jj) = (  e2u(ji,jj) * u_ice(ji,jj) - e2u(ji-1,jj  ) * u_ice(ji-1,jj  )   & 
    608608                  &             + e1v(ji,jj) * v_ice(ji,jj) - e1v(ji  ,jj-1) * v_ice(ji  ,jj-1)   & 
    609                   &            ) * r1_e12t(ji,jj) 
     609                  &            ) * r1_e1e2t(ji,jj) 
    610610 
    611611               zdt(ji,jj) = ( ( u_ice(ji,jj) * r1_e2u(ji,jj) - u_ice(ji-1,jj) * r1_e2u(ji-1,jj) ) * e2t(ji,jj) * e2t(ji,jj)  & 
    612612                  &          -( v_ice(ji,jj) * r1_e1v(ji,jj) - v_ice(ji,jj-1) * r1_e1v(ji,jj-1) ) * e1t(ji,jj) * e1t(ji,jj)  & 
    613                   &         ) * r1_e12t(ji,jj) 
     613                  &         ) * r1_e1e2t(ji,jj) 
    614614               ! 
    615615               ! SB modif because ocean has no slip boundary condition  
    616616               zds(ji,jj) = ( ( u_ice(ji,jj+1) * r1_e1u(ji,jj+1) - u_ice(ji,jj) * r1_e1u(ji,jj) ) * e1f(ji,jj) * e1f(ji,jj)  & 
    617617                  &          +( v_ice(ji+1,jj) * r1_e2v(ji+1,jj) - v_ice(ji,jj) * r1_e2v(ji,jj) ) * e2f(ji,jj) * e2f(ji,jj)  & 
    618                   &         ) * r1_e12f(ji,jj) * ( 2.0 - fmask(ji,jj,1) )                                     & 
     618                  &         ) * r1_e1e2f(ji,jj) * ( 2.0 - fmask(ji,jj,1) )                                     & 
    619619                  &         * zmask(ji,jj) * zmask(ji,jj+1) * zmask(ji+1,jj) * zmask(ji+1,jj+1) 
    620620 
    621621               zdst = ( e2u(ji,jj) * v_ice1(ji,jj) - e2u(ji-1,jj  ) * v_ice1(ji-1,jj  )    & 
    622                   &   + e1v(ji,jj) * u_ice2(ji,jj) - e1v(ji  ,jj-1) * u_ice2(ji  ,jj-1) ) * r1_e12t(ji,jj) 
     622                  &   + e1v(ji,jj) * u_ice2(ji,jj) - e1v(ji  ,jj-1) * u_ice2(ji  ,jj-1) ) * r1_e1e2t(ji,jj) 
    623623 
    624624               delta = SQRT( divu_i(ji,jj)**2 + ( zdt(ji,jj)**2 + zdst**2 ) * usecc2 )   
     
    637637         DO ji = fs_2, fs_jpim1 
    638638            zdst           = (  e2u(ji,jj) * v_ice1(ji,jj) - e2u( ji-1, jj   ) * v_ice1(ji-1,jj)  &    
    639                &              + e1v(ji,jj) * u_ice2(ji,jj) - e1v( ji  , jj-1 ) * u_ice2(ji,jj-1) ) * r1_e12t(ji,jj)  
     639               &              + e1v(ji,jj) * u_ice2(ji,jj) - e1v( ji  , jj-1 ) * u_ice2(ji,jj-1) ) * r1_e1e2t(ji,jj)  
    640640            shear_i(ji,jj) = SQRT( zdt(ji,jj) * zdt(ji,jj) + zdst * zdst ) 
    641641         END DO 
  • trunk/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90

    r5407 r5836  
    325325         CALL prt_ctl_info(' - Cell values : ') 
    326326         CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
    327          CALL prt_ctl(tab2d_1=e12t , clinfo1=' lim_thd  : cell area :') 
     327         CALL prt_ctl(tab2d_1=e1e2t, clinfo1=' lim_thd  : cell area :') 
    328328         CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_thd  : at_i      :') 
    329329         CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_thd  : vt_i      :') 
     
    382382         CALL prt_ctl_info(' - Cell values : ') 
    383383         CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
    384          CALL prt_ctl(tab2d_1=e12t , clinfo1=' lim_itd_th  : cell area :') 
     384         CALL prt_ctl(tab2d_1=e1e2t, clinfo1=' lim_itd_th  : cell area :') 
    385385         CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_itd_th  : at_i      :') 
    386386         CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_itd_th  : vt_i      :') 
  • trunk/NEMOGCM/NEMO/LIM_SRC_3/limtrp.F90

    r5202 r5836  
    9595      ENDIF 
    9696 
    97       zsm(:,:) = e12t(:,:) 
     97      zsm(:,:) = e1e2t(:,:) 
    9898       
    9999      !                             !-------------------------------------! 
     
    162162         ! transported fields                                         
    163163         !------------------------- 
    164          z0opw(:,:,1) = ato_i(:,:) * e12t(:,:)             ! Open water area  
    165          DO jl = 1, jpl 
    166             z0snw (:,:,jl)  = v_s  (:,:,jl) * e12t(:,:)    ! Snow volume 
    167             z0ice(:,:,jl)   = v_i  (:,:,jl) * e12t(:,:)    ! Ice  volume 
    168             z0ai  (:,:,jl)  = a_i  (:,:,jl) * e12t(:,:)    ! Ice area 
    169             z0smi (:,:,jl)  = smv_i(:,:,jl) * e12t(:,:)    ! Salt content 
    170             z0oi (:,:,jl)   = oa_i (:,:,jl) * e12t(:,:)    ! Age content 
    171             z0es (:,:,jl)   = e_s  (:,:,1,jl) * e12t(:,:)  ! Snow heat content 
     164         z0opw(:,:,1) = ato_i(:,:) * e1e2t(:,:)             ! Open water area  
     165         DO jl = 1, jpl 
     166            z0snw (:,:,jl)  = v_s  (:,:,  jl) * e1e2t(:,:)  ! Snow volume 
     167            z0ice(:,:,jl)   = v_i  (:,:,  jl) * e1e2t(:,:)  ! Ice  volume 
     168            z0ai  (:,:,jl)  = a_i  (:,:,  jl) * e1e2t(:,:)  ! Ice area 
     169            z0smi (:,:,jl)  = smv_i(:,:,  jl) * e1e2t(:,:)  ! Salt content 
     170            z0oi (:,:,jl)   = oa_i (:,:,  jl) * e1e2t(:,:)  ! Age content 
     171            z0es (:,:,jl)   = e_s  (:,:,1,jl) * e1e2t(:,:)  ! Snow heat content 
    172172            DO jk = 1, nlay_i 
    173                z0ei  (:,:,jk,jl) = e_i  (:,:,jk,jl) * e12t(:,:) ! Ice  heat content 
     173               z0ei  (:,:,jk,jl) = e_i  (:,:,jk,jl) * e1e2t(:,:) ! Ice  heat content 
    174174            END DO 
    175175         END DO 
     
    263263         ! Recover the properties from their contents 
    264264         !------------------------------------------- 
    265          ato_i(:,:) = z0opw(:,:,1) * r1_e12t(:,:) 
    266          DO jl = 1, jpl 
    267             v_i  (:,:,jl)   = z0ice(:,:,jl) * r1_e12t(:,:) 
    268             v_s  (:,:,jl)   = z0snw(:,:,jl) * r1_e12t(:,:) 
    269             smv_i(:,:,jl)   = z0smi(:,:,jl) * r1_e12t(:,:) 
    270             oa_i (:,:,jl)   = z0oi (:,:,jl) * r1_e12t(:,:) 
    271             a_i  (:,:,jl)   = z0ai (:,:,jl) * r1_e12t(:,:) 
    272             e_s  (:,:,1,jl) = z0es (:,:,jl) * r1_e12t(:,:) 
     265         ato_i(:,:) = z0opw(:,:,1) * r1_e1e2t(:,:) 
     266         DO jl = 1, jpl 
     267            v_i  (:,:,  jl) = z0ice(:,:,jl) * r1_e1e2t(:,:) 
     268            v_s  (:,:,  jl) = z0snw(:,:,jl) * r1_e1e2t(:,:) 
     269            smv_i(:,:,  jl) = z0smi(:,:,jl) * r1_e1e2t(:,:) 
     270            oa_i (:,:,  jl) = z0oi (:,:,jl) * r1_e1e2t(:,:) 
     271            a_i  (:,:,  jl) = z0ai (:,:,jl) * r1_e1e2t(:,:) 
     272            e_s  (:,:,1,jl) = z0es (:,:,jl) * r1_e1e2t(:,:) 
    273273            DO jk = 1, nlay_i 
    274                e_i(:,:,jk,jl) = z0ei(:,:,jk,jl) * r1_e12t(:,:) 
     274               e_i(:,:,jk,jl) = z0ei(:,:,jk,jl) * r1_e1e2t(:,:) 
    275275            END DO 
    276276         END DO 
  • trunk/NEMOGCM/NEMO/LIM_SRC_3/limupdate1.F90

    r5215 r5836  
    146146         CALL prt_ctl_info(' - Cell values : ') 
    147147         CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
    148          CALL prt_ctl(tab2d_1=e12t       , clinfo1=' lim_update1  : cell area   :') 
     148         CALL prt_ctl(tab2d_1=e1e2t      , clinfo1=' lim_update1  : cell area   :') 
    149149         CALL prt_ctl(tab2d_1=at_i       , clinfo1=' lim_update1  : at_i        :') 
    150150         CALL prt_ctl(tab2d_1=vt_i       , clinfo1=' lim_update1  : vt_i        :') 
  • trunk/NEMOGCM/NEMO/LIM_SRC_3/limupdate2.F90

    r5410 r5836  
    191191         CALL prt_ctl_info(' - Cell values : ') 
    192192         CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
    193          CALL prt_ctl(tab2d_1=e12t       , clinfo1=' lim_update2  : cell area   :') 
     193         CALL prt_ctl(tab2d_1=e1e2t      , clinfo1=' lim_update2  : cell area   :') 
    194194         CALL prt_ctl(tab2d_1=at_i       , clinfo1=' lim_update2  : at_i        :') 
    195195         CALL prt_ctl(tab2d_1=vt_i       , clinfo1=' lim_update2  : vt_i        :') 
  • trunk/NEMOGCM/NEMO/NST_SRC/agrif_opa_sponge.F90

    r5656 r5836  
    210210               DO jj = j1,j2-1 
    211211                  DO ji = i1,i2-1 
    212                      zabe1 = fsaht_spu(ji,jj) * umask(ji,jj,jk) * re2u_e1u(ji,jj) * fse3u_n(ji,jj,jk) 
    213                      zabe2 = fsaht_spv(ji,jj) * vmask(ji,jj,jk) * re1v_e2v(ji,jj) * fse3v_n(ji,jj,jk) 
     212                     zabe1 = fsaht_spu(ji,jj) * umask(ji,jj,jk) * e2_e1u(ji,jj) * fse3u_n(ji,jj,jk) 
     213                     zabe2 = fsaht_spv(ji,jj) * vmask(ji,jj,jk) * e1_e2v(ji,jj) * fse3v_n(ji,jj,jk) 
    214214                     ztu(ji,jj,jk) = zabe1 * ( tsbdiff(ji+1,jj  ,jk,jn) - tsbdiff(ji,jj,jk,jn) )  
    215215                     ztv(ji,jj,jk) = zabe2 * ( tsbdiff(ji  ,jj+1,jk,jn) - tsbdiff(ji,jj,jk,jn) ) 
     
    239239 
    240240                     IF (.NOT. tabspongedone_tsn(ji,jj)) THEN  
    241                         zbtr = r1_e12t(ji,jj) / fse3t_n(ji,jj,jk) 
     241                        zbtr = r1_e1e2t(ji,jj) / fse3t_n(ji,jj,jk) 
    242242                        ! horizontal diffusive trends 
    243243                        ztsa = zbtr * (  ztu(ji,jj,jk) - ztu(ji-1,jj,jk) + ztv(ji,jj,jk) - ztv(ji  ,jj-1,jk)  ) 
     
    290290            DO jj = j1,j2 
    291291               DO ji = i1+1,i2   ! vector opt. 
    292                   zbtr = r1_e12t(ji,jj) / fse3t_n(ji,jj,jk) * fsahm_spt(ji,jj) 
     292                  zbtr = r1_e1e2t(ji,jj) / fse3t_n(ji,jj,jk) * fsahm_spt(ji,jj) 
    293293                  hdivdiff(ji,jj,jk) = (  e2u(ji  ,jj)*fse3u_n(ji  ,jj,jk) * ubdiff(ji  ,jj,jk) & 
    294294                                     &   -e2u(ji-1,jj)*fse3u_n(ji-1,jj,jk) * ubdiff(ji-1,jj,jk) ) * zbtr 
     
    298298            DO jj = j1,j2-1 
    299299               DO ji = i1,i2   ! vector opt. 
    300                   zbtr = r1_e12f(ji,jj) * fse3f_n(ji,jj,jk) * fsahm_spf(ji,jj) 
     300                  zbtr = r1_e1e2f(ji,jj) * fse3f_n(ji,jj,jk) * fsahm_spf(ji,jj) 
    301301                  rotdiff(ji,jj,jk) = (-e1u(ji,jj+1) * ubdiff(ji,jj+1,jk) & 
    302302                                       +e1u(ji,jj  ) * ubdiff(ji,jj  ,jk) &  
     
    396396            DO jj = j1+1,j2 
    397397               DO ji = i1,i2   ! vector opt. 
    398                   zbtr = r1_e12t(ji,jj) / fse3t_n(ji,jj,jk) * fsahm_spt(ji,jj) 
     398                  zbtr = r1_e1e2t(ji,jj) / fse3t_n(ji,jj,jk) * fsahm_spt(ji,jj) 
    399399                  hdivdiff(ji,jj,jk) = ( e1v(ji,jj  ) * fse3v(ji,jj  ,jk) * vbdiff(ji,jj  ,jk)  & 
    400400                                     &  -e1v(ji,jj-1) * fse3v(ji,jj-1,jk) * vbdiff(ji,jj-1,jk)  ) * zbtr 
     
    403403            DO jj = j1,j2 
    404404               DO ji = i1,i2-1   ! vector opt. 
    405                   zbtr = r1_e12f(ji,jj) * fse3f_n(ji,jj,jk) * fsahm_spf(ji,jj) 
     405                  zbtr = r1_e1e2f(ji,jj) * fse3f_n(ji,jj,jk) * fsahm_spf(ji,jj) 
    406406                  rotdiff(ji,jj,jk) = ( e2v(ji+1,jj) * vbdiff(ji+1,jj,jk) &  
    407407                                    &  -e2v(ji  ,jj) * vbdiff(ji  ,jj,jk) & 
  • trunk/NEMOGCM/NEMO/NST_SRC/agrif_top_sponge.F90

    r5656 r5836  
    7474               DO jj = j1,j2-1 
    7575                  DO ji = i1,i2-1 
    76                      zabe1 = fsaht_spu(ji,jj) * umask(ji,jj,jk) * re2u_e1u(ji,jj) * fse3u_n(ji,jj,jk) 
    77                      zabe2 = fsaht_spv(ji,jj) * vmask(ji,jj,jk) * re1v_e2v(ji,jj) * fse3v_n(ji,jj,jk) 
     76                     zabe1 = fsaht_spu(ji,jj) * umask(ji,jj,jk) * e2_e1u(ji,jj) * fse3u_n(ji,jj,jk) 
     77                     zabe2 = fsaht_spv(ji,jj) * vmask(ji,jj,jk) * e1_e2v(ji,jj) * fse3v_n(ji,jj,jk) 
    7878                     ztu(ji,jj) = zabe1 * ( trbdiff(ji+1,jj  ,jk,jn) - trbdiff(ji,jj,jk,jn) ) 
    7979                     ztv(ji,jj) = zabe2 * ( trbdiff(ji  ,jj+1,jk,jn) - trbdiff(ji,jj,jk,jn) ) 
     
    8585 
    8686                     IF (.NOT. tabspongedone_trn(ji,jj)) THEN  
    87                         zbtr = r1_e12t(ji,jj) / fse3t(ji,jj,jk) 
     87                        zbtr = r1_e1e2t(ji,jj) / fse3t(ji,jj,jk) 
    8888                        ! horizontal diffusive trends 
    8989                        ztra = zbtr * (  ztu(ji,jj) - ztu(ji-1,jj  ) + ztv(ji,jj) - ztv(ji  ,jj-1)  ) 
  • trunk/NEMOGCM/NEMO/OFF_SRC/domrea.F90

    r5504 r5836  
    44   !! Ocean initialization : domain initialization 
    55   !!============================================================================== 
     6   !! History :  OPA  ! 1990-10  (C. Levy - G. Madec)  Original code 
     7   !!                 ! 1992-01  (M. Imbard) insert time step initialization 
     8   !!                 ! 1996-06  (G. Madec) generalized vertical coordinate  
     9   !!                 ! 1997-02  (G. Madec) creation of domwri.F 
     10   !!                 ! 2001-05  (E.Durand - G. Madec) insert closed sea 
     11   !!  NEMO      1.0  ! 2002-08  (G. Madec)  F90: Free form and module 
     12   !!---------------------------------------------------------------------- 
    613 
    714   !!---------------------------------------------------------------------- 
     
    1017   !!   dom_ctl        : control print for the ocean domain 
    1118   !!---------------------------------------------------------------------- 
    12    !! * Modules used 
    1319   USE oce             !  
     20   USE trc_oce         ! shared ocean/biogeochemical variables 
    1421   USE dom_oce         ! ocean space and time domain 
    1522   USE phycst          ! physical constants 
     23   USE domstp          ! domain: set the time-step 
     24   ! 
    1625   USE in_out_manager  ! I/O manager 
    1726   USE lib_mpp         ! distributed memory computing library 
    18  
    19    USE domstp          ! domain: set the time-step 
    20  
    2127   USE lbclnk          ! lateral boundary condition - MPP exchanges 
    22    USE trc_oce         ! shared ocean/biogeochemical variables 
    2328   USE wrk_nemo   
    2429    
     
    2631   PRIVATE 
    2732 
    28    !! * Routine accessibility 
    29    PUBLIC dom_rea       ! called by opa.F90 
     33   PUBLIC   dom_rea    ! called by nemogcm.F90 
    3034 
    3135   !! * Substitutions 
     
    3337#  include "vectopt_loop_substitute.h90" 
    3438   !!---------------------------------------------------------------------- 
    35    !! NEMO/OFF 3.3 , NEMO Consortium (2010) 
     39   !! NEMO/OFF 3.7 , NEMO Consortium (2015) 
    3640   !! $Id$ 
    3741   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3842   !!---------------------------------------------------------------------- 
    39  
    4043CONTAINS 
    4144 
     
    5154      !!      - dom_stp: defined the model time step 
    5255      !!      - dom_rea: read the meshmask file if nmsh=1 
    53       !! 
    54       !! History : 
    55       !!        !  90-10  (C. Levy - G. Madec)  Original code 
    56       !!        !  91-11  (G. Madec) 
    57       !!        !  92-01  (M. Imbard) insert time step initialization 
    58       !!        !  96-06  (G. Madec) generalized vertical coordinate  
    59       !!        !  97-02  (G. Madec) creation of domwri.F 
    60       !!        !  01-05  (E.Durand - G. Madec) insert closed sea 
    61       !!   8.5  !  02-08  (G. Madec)  F90: Free form and module 
    62       !!---------------------------------------------------------------------- 
    63       !! * Local declarations 
    64       INTEGER ::   jk                ! dummy loop argument 
    65       INTEGER ::   iconf = 0         ! temporary integers 
    66       !!---------------------------------------------------------------------- 
    67  
     56      !!---------------------------------------------------------------------- 
     57      INTEGER ::   jk          ! dummy loop index 
     58      INTEGER ::   iconf = 0   ! local integers 
     59      !!---------------------------------------------------------------------- 
     60      ! 
    6861      IF(lwp) THEN 
    6962         WRITE(numout,*) 
     
    7164         WRITE(numout,*) '~~~~~~~~' 
    7265      ENDIF 
    73  
    74       CALL dom_nam      ! read namelist ( namrun, namdom, namcla ) 
     66      ! 
     67      CALL dom_nam      ! read namelist ( namrun, namdom ) 
    7568      CALL dom_zgr      ! Vertical mesh and bathymetry option 
    7669      CALL dom_grd      ! Create a domain file 
    77  
    78      ! 
    79       ! - ML - Used in dom_vvl_sf_nxt and lateral diffusion routines 
    80       !        but could be usefull in many other routines 
    81       e12t    (:,:) = e1t(:,:) * e2t(:,:) 
    82       e1e2t   (:,:) = e1t(:,:) * e2t(:,:) 
    83       e12u    (:,:) = e1u(:,:) * e2u(:,:) 
    84       e12v    (:,:) = e1v(:,:) * e2v(:,:) 
    85       e12f    (:,:) = e1f(:,:) * e2f(:,:) 
    86       r1_e12t (:,:) = 1._wp    / e12t(:,:) 
    87       r1_e12u (:,:) = 1._wp    / e12u(:,:) 
    88       r1_e12v (:,:) = 1._wp    / e12v(:,:) 
    89       r1_e12f (:,:) = 1._wp    / e12f(:,:) 
    90       re2u_e1u(:,:) = e2u(:,:) / e1u(:,:) 
    91       re1v_e2v(:,:) = e1v(:,:) / e2v(:,:) 
    92       ! 
    93       hu(:,:) = 0._wp                          ! Ocean depth at U- and V-points 
     70      ! 
     71      !                                      ! associated horizontal metrics 
     72      ! 
     73      r1_e1t(:,:) = 1._wp / e1t(:,:)   ;   r1_e2t (:,:) = 1._wp / e2t(:,:) 
     74      r1_e1u(:,:) = 1._wp / e1u(:,:)   ;   r1_e2u (:,:) = 1._wp / e2u(:,:) 
     75      r1_e1v(:,:) = 1._wp / e1v(:,:)   ;   r1_e2v (:,:) = 1._wp / e2v(:,:) 
     76      r1_e1f(:,:) = 1._wp / e1f(:,:)   ;   r1_e2f (:,:) = 1._wp / e2f(:,:) 
     77      ! 
     78      e1e2t (:,:) = e1t(:,:) * e2t(:,:)   ;   r1_e1e2t(:,:) = 1._wp / e1e2t(:,:) 
     79      e1e2u (:,:) = e1u(:,:) * e2u(:,:)   ;   r1_e1e2u(:,:) = 1._wp / e1e2u(:,:) 
     80      e1e2v (:,:) = e1v(:,:) * e2v(:,:)   ;   r1_e1e2v(:,:) = 1._wp / e1e2v(:,:) 
     81      e1e2f (:,:) = e1f(:,:) * e2f(:,:)   ;   r1_e1e2f(:,:) = 1._wp / e1e2f(:,:) 
     82      !    
     83      e2_e1u(:,:) = e2u(:,:) / e1u(:,:) 
     84      e1_e2v(:,:) = e1v(:,:) / e2v(:,:) 
     85      ! 
     86      hu(:,:) = 0._wp                        ! Ocean depth at U- and V-points 
    9487      hv(:,:) = 0._wp 
    9588      DO jk = 1, jpk 
     
    10093      hur(:,:) = 1._wp / ( hu(:,:) + 1._wp - umask(:,:,1) ) * umask(:,:,1) 
    10194      hvr(:,:) = 1._wp / ( hv(:,:) + 1._wp - vmask(:,:,1) ) * vmask(:,:,1) 
    102  
     95      ! 
    10396      CALL dom_stp      ! Time step 
    10497      CALL dom_msk      ! Masks 
    10598      CALL dom_ctl      ! Domain control 
    106  
     99      ! 
    107100   END SUBROUTINE dom_rea 
     101 
    108102 
    109103   SUBROUTINE dom_nam 
     
    115109      !! ** input   : - namrun namelist 
    116110      !!              - namdom namelist 
    117       !!              - namcla namelist 
    118111      !!---------------------------------------------------------------------- 
    119112      USE ioipsl 
    120       INTEGER  ::   ios                 ! Local integer output status for namelist read 
     113      INTEGER  ::   ios   ! Local integer output status for namelist read 
     114      ! 
    121115      NAMELIST/namrun/ cn_ocerst_indir, cn_ocerst_outdir, nn_stocklist, ln_rst_list,               & 
    122116         &             nn_no   , cn_exp    , cn_ocerst_in, cn_ocerst_out, ln_rstart , nn_rstctl,   & 
     
    130124         &             ppsur, ppa0, ppa1, ppkth, ppacr, ppdzmin, pphmax, ldbletanh, & 
    131125         &             ppa2, ppkth2, ppacr2 
    132       NAMELIST/namcla/ nn_cla 
    133126#if defined key_netcdf4 
    134127      NAMELIST/namnc4/ nn_nchunks_i, nn_nchunks_j, nn_nchunks_k, ln_nc4zip 
     
    178171      nstocklist = nn_stocklist 
    179172      nwrite = nn_write 
    180  
    181  
     173      ! 
    182174      !                             ! control of output frequency 
    183175      IF ( nstock == 0 .OR. nstock > nitend ) THEN 
     
    275267      rdth      = rn_rdth 
    276268 
    277       REWIND( numnam_ref )              ! Namelist namcla in reference namelist : Cross land advection 
    278       READ  ( numnam_ref, namcla, IOSTAT = ios, ERR = 905) 
    279 905   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in reference namelist', lwp ) 
    280  
    281       REWIND( numnam_cfg )              ! Namelist namcla in configuration namelist : Cross land advection 
    282       READ  ( numnam_cfg, namcla, IOSTAT = ios, ERR = 906 ) 
    283 906   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namcla in configuration namelist', lwp ) 
    284       IF(lwm) WRITE( numond, namcla ) 
    285  
    286       IF(lwp) THEN 
    287          WRITE(numout,*) 
    288          WRITE(numout,*) '   Namelist namcla' 
    289          WRITE(numout,*) '      cross land advection                 nn_cla    = ', nn_cla 
    290       ENDIF 
    291  
    292269#if defined key_netcdf4 
    293270      !                             ! NetCDF 4 case   ("key_netcdf4" defined) 
     
    321298   END SUBROUTINE dom_nam 
    322299 
     300 
    323301   SUBROUTINE dom_zgr 
    324302      !!---------------------------------------------------------------------- 
     
    374352   END SUBROUTINE dom_zgr 
    375353 
     354 
    376355   SUBROUTINE dom_ctl 
    377356      !!---------------------------------------------------------------------- 
     
    382361      !! ** Method  :   compute and print extrema of masked scale factors 
    383362      !! 
    384       !! History : 
    385       !!   8.5  !  02-08  (G. Madec)    Original code 
    386       !!---------------------------------------------------------------------- 
    387       !! * Local declarations 
     363      !!---------------------------------------------------------------------- 
    388364      INTEGER ::   iimi1, ijmi1, iimi2, ijmi2, iima1, ijma1, iima2, ijma2 
    389365      INTEGER, DIMENSION(2) ::   iloc      !  
     
    421397         ijma2 = iloc(2) + njmpp - 1 
    422398      ENDIF 
    423  
     399      ! 
    424400      IF(lwp) THEN 
    425401         WRITE(numout,"(14x,'e1t maxi: ',1f10.2,' at i = ',i5,' j= ',i5)") ze1max, iima1, ijma1 
     
    428404         WRITE(numout,"(14x,'e2t mini: ',1f10.2,' at i = ',i5,' j= ',i5)") ze2min, iimi2, ijmi2 
    429405      ENDIF 
    430  
     406      ! 
    431407   END SUBROUTINE dom_ctl 
     408 
    432409 
    433410   SUBROUTINE dom_grd 
     
    538515         CALL iom_get( inum2, jpdom_data, 'facvolt', facvol ) 
    539516#endif 
    540  
    541517         !                                                         ! horizontal mesh (inum3) 
    542518         CALL iom_get( inum3, jpdom_data, 'glamt', glamt ) 
     
    756732      !!                                     (min value = 1 over land) 
    757733      !!---------------------------------------------------------------------- 
    758       ! 
    759734      INTEGER ::   ji, jj   ! dummy loop indices 
    760735      REAL(wp), POINTER, DIMENSION(:,:) :: zmbk 
     
    785760   END SUBROUTINE zgr_bot_level 
    786761 
     762 
    787763   SUBROUTINE dom_msk 
    788764      !!--------------------------------------------------------------------- 
     
    799775      !!               tpol     : ??? 
    800776      !!---------------------------------------------------------------------- 
    801       ! 
    802       INTEGER  ::   ji, jj, jk                   ! dummy loop indices 
    803       INTEGER  ::   iif, iil, ijf, ijl       ! local integers 
     777      INTEGER  ::   ji, jj, jk           ! dummy loop indices 
     778      INTEGER  ::   iif, iil, ijf, ijl   ! local integers 
    804779      INTEGER, POINTER, DIMENSION(:,:) ::  imsk  
    805       ! 
    806780      !!--------------------------------------------------------------------- 
    807781       
     
    853827      ! 3. Ocean/land mask at wu-, wv- and w points  
    854828      !---------------------------------------------- 
    855       wmask (:,:,1) = tmask(:,:,1) ! ???????? 
    856       wumask(:,:,1) = umask(:,:,1) ! ???????? 
    857       wvmask(:,:,1) = vmask(:,:,1) ! ???????? 
    858       DO jk=2,jpk 
    859          wmask (:,:,jk)=tmask(:,:,jk) * tmask(:,:,jk-1) 
    860          wumask(:,:,jk)=umask(:,:,jk) * umask(:,:,jk-1)    
    861          wvmask(:,:,jk)=vmask(:,:,jk) * vmask(:,:,jk-1) 
     829      wmask (:,:,1) = tmask(:,:,1)        ! surface value 
     830      wumask(:,:,1) = umask(:,:,1)  
     831      wvmask(:,:,1) = vmask(:,:,1) 
     832      DO jk = 2, jpk                      ! deeper value 
     833         wmask (:,:,jk) = tmask(:,:,jk) * tmask(:,:,jk-1) 
     834         wumask(:,:,jk) = umask(:,:,jk) * umask(:,:,jk-1)    
     835         wvmask(:,:,jk) = vmask(:,:,jk) * vmask(:,:,jk-1) 
    862836      END DO 
    863837      ! 
  • trunk/NEMOGCM/NEMO/OFF_SRC/dtadyn.F90

    r5768 r5836  
    2626   USE trc_oce         ! share ocean/biogeo variables 
    2727   USE phycst          ! physical constants 
     28   USE ldftra          ! lateral diffusivity coefficients 
    2829   USE trabbl          ! active tracer: bottom boundary layer 
    2930   USE ldfslp          ! lateral diffusion: iso-neutral slopes 
    30    USE ldfeiv          ! eddy induced velocity coef.  
    31    USE ldftra_oce      ! ocean tracer   lateral physics 
    3231   USE zdfmxl          ! vertical physics: mixed layer depth 
    3332   USE eosbn2          ! equation of state - Brunt Vaisala frequency 
     
    4039   USE fldread         ! read input fields  
    4140   USE timing          ! Timing 
     41   USE wrk_nemo 
    4242 
    4343   IMPLICIT NONE 
     
    5050   LOGICAL            ::   ln_dynwzv    !: vertical velocity read in a file (T) or computed from u/v (F) 
    5151   LOGICAL            ::   ln_dynbbl    !: bbl coef read in a file (T) or computed (F) 
    52    LOGICAL            ::   ln_degrad    !: degradation option enabled or not 
    5352   LOGICAL            ::   ln_dynrnf    !: read runoff data in file (T) or set to zero (F) 
    5453 
    55    INTEGER  , PARAMETER ::   jpfld = 21     ! maximum number of fields to read 
     54   INTEGER  , PARAMETER ::   jpfld = 15     ! maximum number of fields to read 
    5655   INTEGER  , SAVE      ::   jf_tem         ! index of temperature 
    5756   INTEGER  , SAVE      ::   jf_sal         ! index of salinity 
     
    6867   INTEGER  , SAVE      ::   jf_ubl         ! index of u-bbl coef 
    6968   INTEGER  , SAVE      ::   jf_vbl         ! index of v-bbl coef 
    70    INTEGER  , SAVE      ::   jf_ahu         ! index of u-diffusivity coef 
    71    INTEGER  , SAVE      ::   jf_ahv         ! index of v-diffusivity coef  
    72    INTEGER  , SAVE      ::   jf_ahw         ! index of w-diffusivity coef 
    73    INTEGER  , SAVE      ::   jf_eiu         ! index of u-eiv 
    74    INTEGER  , SAVE      ::   jf_eiv         ! index of v-eiv 
    75    INTEGER  , SAVE      ::   jf_eiw         ! index of w-eiv 
    7669   INTEGER  , SAVE      ::   jf_fmf         ! index of downward salt flux 
    7770 
     
    112105      !!             - interpolates data if needed 
    113106      !!---------------------------------------------------------------------- 
    114       ! 
    115       USE oce, ONLY:  zts    => tsa  
     107      USE oce, ONLY:  zts    => tsa 
    116108      USE oce, ONLY:  zuslp  => ua   , zvslp  => va 
    117       USE oce, ONLY:  zwslpi => rotb , zwslpj => rotn 
    118       USE oce, ONLY:  zu     => ub   , zv     => vb,  zw => hdivb 
     109      USE oce, ONLY:  zwslpi => ua_sv , zwslpj => va_sv 
     110      USE oce, ONLY:  zu     => ub   , zv     => vb,  zw => rke 
    119111      ! 
    120112      INTEGER, INTENT(in) ::   kt   ! ocean time-step index 
     113      ! 
     114!      REAL(wp), DIMENSION(jpi,jpj,jpk,jpts)  :: zts 
     115!      REAL(wp), DIMENSION(jpi,jpj,jpk     )  :: zuslp, zvslp, zwslpi, zwslpj 
     116!      REAL(wp), DIMENSION(jpi,jpj,jpk     )  :: zu, zv, zw 
     117      ! 
    121118      ! 
    122119      INTEGER  ::   ji, jj     ! dummy loop indices 
     
    138135         CALL fld_read( kt, 1, sf_dyn )      !==   read data at kt time step   ==! 
    139136         ! 
    140          IF( lk_ldfslp .AND. .NOT.lk_c1d .AND. sf_dyn(jf_tem)%ln_tint ) THEN    ! Computes slopes (here avt is used as workspace)                        
     137         IF( l_ldfslp .AND. .NOT.lk_c1d .AND. sf_dyn(jf_tem)%ln_tint ) THEN    ! Computes slopes (here avt is used as workspace)                        
    141138            zts(:,:,:,jp_tem) = sf_dyn(jf_tem)%fdta(:,:,:,1) * tmask(:,:,:)   ! temperature 
    142139            zts(:,:,:,jp_sal) = sf_dyn(jf_sal)%fdta(:,:,:,1) * tmask(:,:,:)   ! salinity  
     
    162159      ENDIF 
    163160      !  
    164       IF( lk_ldfslp .AND. .NOT.lk_c1d ) THEN    ! Computes slopes (here avt is used as workspace)                        
     161      IF( l_ldfslp .AND. .NOT.lk_c1d ) THEN    ! Computes slopes (here avt is used as workspace)                        
    165162         iswap_tem = 0 
    166163         IF(  kt /= nit000 .AND. ( sf_dyn(jf_tem)%nrec_a(2) - nrecprev_tem ) /= 0 )  iswap_tem = 1 
     
    267264      rnf (:,:)        = sf_dyn(jf_rnf)%fnow(:,:,1) * tmask(:,:,1)    ! river runoffs  
    268265 
     266      !                                               ! update eddy diffusivity coeff. and/or eiv coeff. at kt 
     267      IF( l_ldftra_time .OR. l_ldfeiv_time )   CALL ldf_tra( kt )  
    269268      !                                                      ! bbl diffusive coef 
    270269#if defined key_trabbl && ! defined key_c1d 
     
    276275         CALL bbl( kt, nit000, 'TRC') 
    277276      END IF 
    278 #endif 
    279 #if ( ! defined key_degrad && defined key_traldf_c2d && defined key_traldf_eiv ) && ! defined key_c1d  
    280       aeiw(:,:)        = sf_dyn(jf_eiw)%fnow(:,:,1) * tmask(:,:,1)    ! w-eiv 
    281       !                                                           ! Computes the horizontal values from the vertical value 
    282       DO jj = 2, jpjm1 
    283          DO ji = fs_2, fs_jpim1   ! vector opt. 
    284             aeiu(ji,jj) = .5 * ( aeiw(ji,jj) + aeiw(ji+1,jj  ) )  ! Average the diffusive coefficient at u- v- points 
    285             aeiv(ji,jj) = .5 * ( aeiw(ji,jj) + aeiw(ji  ,jj+1) )  ! at u- v- points 
    286          END DO 
    287       END DO 
    288       CALL lbc_lnk( aeiu, 'U', 1. )   ;   CALL lbc_lnk( aeiv, 'V', 1. )    ! lateral boundary condition 
    289 #endif 
    290        
    291 #if defined key_degrad && ! defined key_c1d  
    292       !                                          ! degrad option : diffusive and eiv coef are 3D 
    293       ahtu(:,:,:) = sf_dyn(jf_ahu)%fnow(:,:,:) * umask(:,:,:) 
    294       ahtv(:,:,:) = sf_dyn(jf_ahv)%fnow(:,:,:) * vmask(:,:,:) 
    295       ahtw(:,:,:) = sf_dyn(jf_ahw)%fnow(:,:,:) * tmask(:,:,:) 
    296 #  if defined key_traldf_eiv  
    297       aeiu(:,:,:) = sf_dyn(jf_eiu)%fnow(:,:,:) * umask(:,:,:) 
    298       aeiv(:,:,:) = sf_dyn(jf_eiv)%fnow(:,:,:) * vmask(:,:,:) 
    299       aeiw(:,:,:) = sf_dyn(jf_eiw)%fnow(:,:,:) * tmask(:,:,:) 
    300 #  endif 
    301277#endif 
    302278      ! 
     
    339315      TYPE(FLD_N), DIMENSION(jpfld) ::  slf_d    ! array of namelist informations on the fields to read 
    340316      TYPE(FLD_N) :: sn_tem, sn_sal, sn_mld, sn_emp, sn_ice, sn_qsr, sn_wnd, sn_rnf  ! informations about the fields to be read 
    341       TYPE(FLD_N) :: sn_uwd, sn_vwd, sn_wwd, sn_avt, sn_ubl, sn_vbl          !   "                                 " 
    342       TYPE(FLD_N) :: sn_ahu, sn_ahv, sn_ahw, sn_eiu, sn_eiv, sn_eiw, sn_fmf  !   "                                 " 
    343       !!---------------------------------------------------------------------- 
    344       ! 
    345       NAMELIST/namdta_dyn/cn_dir, ln_dynwzv, ln_dynbbl, ln_degrad, ln_dynrnf,    & 
     317      TYPE(FLD_N) :: sn_uwd, sn_vwd, sn_wwd, sn_avt, sn_ubl, sn_vbl, sn_fmf          !   "                                 " 
     318      NAMELIST/namdta_dyn/cn_dir, ln_dynwzv, ln_dynbbl, ln_dynrnf,    & 
    346319         &                sn_tem, sn_sal, sn_mld, sn_emp, sn_ice, sn_qsr, sn_wnd, sn_rnf,  & 
    347          &                sn_uwd, sn_vwd, sn_wwd, sn_avt, sn_ubl, sn_vbl,          & 
    348          &                sn_ahu, sn_ahv, sn_ahw, sn_eiu, sn_eiv, sn_eiw, sn_fmf 
     320         &                sn_uwd, sn_vwd, sn_wwd, sn_avt, sn_ubl, sn_vbl, sn_fmf   
     321      !!---------------------------------------------------------------------- 
    349322      ! 
    350323      REWIND( numnam_ref )              ! Namelist namdta_dyn in reference namelist : Offline: init. of dynamical data 
     
    365338         WRITE(numout,*) '      vertical velocity read from file (T) or computed (F) ln_dynwzv  = ', ln_dynwzv 
    366339         WRITE(numout,*) '      bbl coef read from file (T) or computed (F)          ln_dynbbl  = ', ln_dynbbl 
    367          WRITE(numout,*) '      degradation option enabled (T) or not (F)            ln_degrad  = ', ln_degrad 
    368340         WRITE(numout,*) '      river runoff option enabled (T) or not (F)           ln_dynrnf  = ', ln_dynrnf 
    369341         WRITE(numout,*) 
    370342      ENDIF 
    371343      !  
    372       IF( ln_degrad .AND. .NOT.lk_degrad ) THEN 
    373          CALL ctl_warn( 'dta_dyn_init: degradation option requires key_degrad activated ; force ln_degrad to false' ) 
    374          ln_degrad = .FALSE. 
    375       ENDIF 
    376344      IF( ln_dynbbl .AND. ( .NOT.lk_trabbl .OR. lk_c1d ) ) THEN 
    377345         CALL ctl_warn( 'dta_dyn_init: bbl option requires key_trabbl activated ; force ln_dynbbl to false' ) 
     
    395363      ENDIF 
    396364 
    397       ! 
    398       IF( .NOT.ln_degrad ) THEN     ! no degrad option 
    399          IF( lk_traldf_eiv .AND. ln_dynbbl ) THEN        ! eiv & bbl 
    400                  jf_ubl  = jfld + 1 ;        jf_vbl  = jfld + 2 ;        jf_eiw  = jfld + 3   ;   jfld = jf_eiw 
    401            slf_d(jf_ubl) = sn_ubl  ;   slf_d(jf_vbl) = sn_vbl  ;   slf_d(jf_eiw) = sn_eiw 
    402          ENDIF 
    403          IF( .NOT.lk_traldf_eiv .AND. ln_dynbbl ) THEN   ! no eiv & bbl 
     365      IF( ln_dynbbl ) THEN         ! eiv & bbl 
    404366                 jf_ubl  = jfld + 1 ;        jf_vbl  = jfld + 2 ;  jfld = jf_vbl 
    405367           slf_d(jf_ubl) = sn_ubl  ;   slf_d(jf_vbl) = sn_vbl 
    406          ENDIF 
    407          IF( lk_traldf_eiv .AND. .NOT.ln_dynbbl ) THEN   ! eiv & no bbl 
    408            jf_eiw = jfld + 1 ; jfld = jf_eiw ; slf_d(jf_eiw) = sn_eiw 
    409          ENDIF 
    410       ELSE 
    411               jf_ahu  = jfld + 1 ;        jf_ahv  = jfld + 2 ;        jf_ahw  = jfld + 3  ;  jfld = jf_ahw 
    412         slf_d(jf_ahu) = sn_ahu  ;   slf_d(jf_ahv) = sn_ahv  ;   slf_d(jf_ahw) = sn_ahw 
    413         IF( lk_traldf_eiv .AND. ln_dynbbl ) THEN         ! eiv & bbl 
    414                  jf_ubl  = jfld + 1 ;        jf_vbl  = jfld + 2 ; 
    415            slf_d(jf_ubl) = sn_ubl  ;   slf_d(jf_vbl) = sn_vbl 
    416                  jf_eiu  = jfld + 3 ;        jf_eiv  = jfld + 4 ;    jf_eiw  = jfld + 5   ;  jfld = jf_eiw  
    417            slf_d(jf_eiu) = sn_eiu  ;   slf_d(jf_eiv) = sn_eiv  ;    slf_d(jf_eiw) = sn_eiw 
    418         ENDIF 
    419         IF( .NOT.lk_traldf_eiv .AND. ln_dynbbl ) THEN    ! no eiv & bbl 
    420                  jf_ubl  = jfld + 1 ;        jf_vbl  = jfld + 2 ;  jfld = jf_vbl 
    421            slf_d(jf_ubl) = sn_ubl  ;   slf_d(jf_vbl) = sn_vbl 
    422         ENDIF 
    423         IF( lk_traldf_eiv .AND. .NOT.ln_dynbbl ) THEN    ! eiv & no bbl 
    424                  jf_eiu  = jfld + 1 ;         jf_eiv  = jfld + 2 ;    jf_eiw  = jfld + 3   ; jfld = jf_eiw  
    425            slf_d(jf_eiu) = sn_eiu  ;   slf_d(jf_eiv) = sn_eiv  ;   slf_d(jf_eiw) = sn_eiw 
    426         ENDIF 
    427       ENDIF 
    428    
     368      ENDIF 
     369 
     370 
    429371      ALLOCATE( sf_dyn(jfld), STAT=ierr )         ! set sf structure 
    430372      IF( ierr > 0 ) THEN 
     
    452394      END DO 
    453395      ! 
    454       IF( lk_ldfslp .AND. .NOT.lk_c1d ) THEN                  ! slopes  
     396      IF( l_ldfslp .AND. .NOT.lk_c1d ) THEN                  ! slopes  
    455397         IF( sf_dyn(jf_tem)%ln_tint ) THEN      ! time interpolation 
    456398            ALLOCATE( uslpdta (jpi,jpj,jpk,2), vslpdta (jpi,jpj,jpk,2),    & 
     
    511453               zv  = pv(ji  ,jj  ,jk) * vmask(ji  ,jj  ,jk) * e1v(ji  ,jj  ) * fse3v(ji  ,jj  ,jk) 
    512454               zv1 = pv(ji  ,jj-1,jk) * vmask(ji  ,jj-1,jk) * e1v(ji  ,jj-1) * fse3v(ji  ,jj-1,jk) 
    513                zet = 1. / ( e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) ) 
     455               zet = 1. / ( e1e2t(ji,jj) * fse3t(ji,jj,jk) ) 
    514456               zhdiv(ji,jj,jk) = ( zu - zu1 + zv - zv1 ) * zet  
    515457            END DO 
    516458         END DO 
    517459      END DO 
     460      !                              !  update the horizontal divergence with the runoff inflow 
     461      IF( ln_dynrnf )  zhdiv(:,:,1) = zhdiv(:,:,1) - rnf(:,:) * r1_rau0 / fse3t(:,:,1) 
     462      ! 
    518463      CALL lbc_lnk( zhdiv, 'T', 1. )      ! Lateral boundary conditions on zhdiv 
    519       ! 
    520464      ! computation of vertical velocity from the bottom 
    521465      pw(:,:,jpk) = 0._wp 
     
    540484      REAL(wp), DIMENSION(jpi,jpj,jpk)     , INTENT(out) :: pwslpj   ! meridional diapycnal slopes 
    541485      !!--------------------------------------------------------------------- 
    542 #if defined key_ldfslp && ! defined key_c1d 
    543       CALL eos    ( pts, rhd, rhop, gdept_0(:,:,:) ) 
    544       CALL eos_rab( pts, rab_n )       ! now local thermal/haline expension ratio at T-points 
    545       CALL bn2    ( pts, rab_n, rn2  ) ! now    Brunt-Vaisala 
    546  
    547       ! Partial steps: before Horizontal DErivative 
    548       IF( ln_zps  .AND. .NOT. ln_isfcav)                            & 
    549          &            CALL zps_hde    ( kt, jpts, pts, gtsu, gtsv,  &  ! Partial steps: before horizontal gradient 
    550          &                                        rhd, gru , grv    )  ! of t, s, rd at the last ocean level 
    551       IF( ln_zps .AND.        ln_isfcav)                            & 
    552          &            CALL zps_hde_isf( kt, jpts, pts, gtsu, gtsv,  &    ! Partial steps for top cell (ISF) 
    553          &                                        rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
    554          &                                 gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi    ) ! of t, s, rd at the first ocean level 
    555  
    556       rn2b(:,:,:) = rn2(:,:,:)         ! need for zdfmxl 
    557       CALL zdf_mxl( kt )            ! mixed layer depth 
    558       CALL ldf_slp( kt, rhd, rn2 )  ! slopes 
    559       puslp (:,:,:) = uslp (:,:,:)  
    560       pvslp (:,:,:) = vslp (:,:,:)  
    561       pwslpi(:,:,:) = wslpi(:,:,:)  
    562       pwslpj(:,:,:) = wslpj(:,:,:)  
    563 #else 
    564       puslp (:,:,:) = 0.            ! to avoid warning when compiling 
    565       pvslp (:,:,:) = 0. 
    566       pwslpi(:,:,:) = 0. 
    567       pwslpj(:,:,:) = 0. 
    568 #endif 
     486      IF( l_ldfslp .AND. .NOT.lk_c1d ) THEN    ! Computes slopes (here avt is used as workspace)                        
     487         CALL eos    ( pts, rhd, rhop, gdept_0(:,:,:) ) 
     488         CALL eos_rab( pts, rab_n )       ! now local thermal/haline expension ratio at T-points 
     489         CALL bn2    ( pts, rab_n, rn2  ) ! now    Brunt-Vaisala 
     490 
     491         ! Partial steps: before Horizontal DErivative 
     492         IF( ln_zps  .AND. .NOT. ln_isfcav)                            & 
     493            &            CALL zps_hde    ( kt, jpts, pts, gtsu, gtsv,  &  ! Partial steps: before horizontal gradient 
     494            &                                        rhd, gru , grv    )  ! of t, s, rd at the last ocean level 
     495         IF( ln_zps .AND.        ln_isfcav)                            & 
     496            &            CALL zps_hde_isf( kt, jpts, pts, gtsu, gtsv,  &    ! Partial steps for top cell (ISF) 
     497            &                                        rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
     498            &                                 gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi    ) ! of t, s, rd at the first ocean level 
     499 
     500         rn2b(:,:,:) = rn2(:,:,:)         ! need for zdfmxl 
     501         CALL zdf_mxl( kt )            ! mixed layer depth 
     502         CALL ldf_slp( kt, rhd, rn2 )  ! slopes 
     503         puslp (:,:,:) = uslp (:,:,:)  
     504         pvslp (:,:,:) = vslp (:,:,:)  
     505         pwslpi(:,:,:) = wslpi(:,:,:)  
     506         pwslpj(:,:,:) = wslpj(:,:,:)  
     507     ELSE 
     508         puslp (:,:,:) = 0.            ! to avoid warning when compiling 
     509         pvslp (:,:,:) = 0. 
     510         pwslpi(:,:,:) = 0. 
     511         pwslpj(:,:,:) = 0. 
     512     ENDIF 
    569513      ! 
    570514   END SUBROUTINE dta_dyn_slp 
  • trunk/NEMOGCM/NEMO/OFF_SRC/nemogcm.F90

    r5504 r5836  
    2626   USE traqsr          ! solar radiation penetration    (tra_qsr_init routine) 
    2727   USE trabbl          ! bottom boundary layer          (tra_bbl_init routine) 
     28   USE traldf          ! lateral physics                (tra_ldf_init routine) 
    2829   USE zdfini          ! vertical physics: initialization 
    2930   USE sbcmod          ! surface boundary condition       (sbc_init     routine) 
     
    283284                            CALL     sbc_init   ! Forcings : surface module 
    284285 
    285 #if ! defined key_degrad 
    286286                            CALL ldf_tra_init   ! Lateral ocean tracer physics 
    287 #endif 
    288       IF( lk_ldfslp )       CALL ldf_slp_init   ! slope of lateral mixing 
     287                            CALL ldf_eiv_init   ! Eddy induced velocity param 
     288                            CALL tra_ldf_init   ! lateral mixing 
     289      IF( l_ldfslp )        CALL ldf_slp_init   ! slope of lateral mixing 
    289290 
    290291                            CALL tra_qsr_init   ! penetrative solar radiation qsr 
     
    444445      USE dom_oce,      ONLY: dom_oce_alloc 
    445446      USE zdf_oce,      ONLY: zdf_oce_alloc 
    446       USE ldftra_oce,   ONLY: ldftra_oce_alloc 
    447447      USE trc_oce,      ONLY: trc_oce_alloc 
    448448      ! 
     
    453453      ierr = ierr + dia_wri_alloc   () 
    454454      ierr = ierr + dom_oce_alloc   ()          ! ocean domain 
    455       ierr = ierr + ldftra_oce_alloc()          ! ocean lateral  physics : tracers 
    456455      ierr = ierr + zdf_oce_alloc   ()          ! ocean vertical physics 
    457456      ! 
  • trunk/NEMOGCM/NEMO/OPA_SRC/ASM/asmbkg.F90

    r5215 r5836  
    1818 
    1919   !!---------------------------------------------------------------------- 
    20    !!   'key_asminc' : Switch on the assimilation increment interface 
    21    !!---------------------------------------------------------------------- 
    2220   !!   asm_bkg_wri  : Write out the background state 
    2321   !!   asm_trj_wri  : Write out the model state trajectory (used with 4D-Var) 
     
    2725   USE zdf_oce            ! Vertical mixing variables 
    2826   USE zdfddm             ! Double diffusion mixing parameterization 
    29    USE ldftra_oce         ! Lateral tracer mixing coefficient defined in memory 
    30    USE ldfslp             ! Slopes of neutral surfaces 
     27   USE ldftra             ! Lateral diffusion: eddy diffusivity coefficients 
     28   USE ldfslp             ! Lateral diffusion: slopes of neutral surfaces 
    3129   USE tradmp             ! Tracer damping 
    3230#if defined key_zdftke 
     
    4139   USE asmpar             ! Parameters for the assmilation interface 
    4240   USE zdfmxl             ! mixed layer depth 
    43 #if defined key_traldf_c2d 
    44    USE ldfeiv             ! eddy induced velocity coef.      (ldf_eiv routine) 
    45 #endif 
    4641#if defined key_lim2 
    4742   USE ice_2 
     
    155150            CALL iom_rstput( kt, nitdin_r, inum, 'sshn'   , sshn              ) 
    156151#if defined key_lim2 || defined key_lim3 
    157             IF(( nn_ice == 2 ) .OR. ( nn_ice == 3 )) THEN 
    158           IF(ALLOCATED(frld)) THEN 
    159                   CALL iom_rstput( kt, nitdin_r, inum, 'iceconc', 1.0 - frld(:,:)   ) 
     152            IF( nn_ice == 2  .OR.  nn_ice == 3 ) THEN 
     153               IF( ALLOCATED(frld) ) THEN 
     154                  CALL iom_rstput( kt, nitdin_r, inum, 'iceconc', 1._wp - frld(:,:)   ) 
    160155               ELSE 
    161         CALL ctl_warn('Ice concentration not written to background as ice variable frld not allocated on this timestep') 
    162           ENDIF 
     156                  CALL ctl_warn('Ice concentration not written to background as ice variable frld not allocated on this timestep') 
     157               ENDIF 
    163158            ENDIF 
    164159#endif 
  • trunk/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

    r5541 r5836  
    1414 
    1515   !!---------------------------------------------------------------------- 
    16    !!   'key_asminc'   : Switch on the assimilation increment interface 
    17    !!---------------------------------------------------------------------- 
    1816   !!   asm_inc_init   : Initialize the increment arrays and IAU weights 
    1917   !!   calc_date      : Compute the calendar date YYYYMMDD on a given step 
     
    2826   USE domvvl           ! domain: variable volume level 
    2927   USE oce              ! Dynamics and active tracers defined in memory 
    30    USE ldfdyn_oce       ! ocean dynamics: lateral physics 
     28   USE ldfdyn           ! lateral diffusion: eddy viscosity coefficients 
    3129   USE eosbn2           ! Equation of state - in situ and potential density 
    3230   USE zpshde           ! Partial step : Horizontal Derivative 
     
    5654    LOGICAL, PUBLIC, PARAMETER :: lk_asminc = .FALSE.  !: No assimilation increments 
    5755#endif 
    58    LOGICAL, PUBLIC :: ln_bkgwri = .FALSE.      !: No output of the background state fields 
    59    LOGICAL, PUBLIC :: ln_asmiau = .FALSE.      !: No applying forcing with an assimilation increment 
    60    LOGICAL, PUBLIC :: ln_asmdin = .FALSE.      !: No direct initialization 
    61    LOGICAL, PUBLIC :: ln_trainc = .FALSE.      !: No tracer (T and S) assimilation increments 
    62    LOGICAL, PUBLIC :: ln_dyninc = .FALSE.      !: No dynamics (u and v) assimilation increments 
    63    LOGICAL, PUBLIC :: ln_sshinc = .FALSE.      !: No sea surface height assimilation increment 
    64    LOGICAL, PUBLIC :: ln_seaiceinc             !: No sea ice concentration increment 
    65    LOGICAL, PUBLIC :: ln_salfix = .FALSE.      !: Apply minimum salinity check 
     56   LOGICAL, PUBLIC :: ln_bkgwri     !: No output of the background state fields 
     57   LOGICAL, PUBLIC :: ln_asmiau     !: No applying forcing with an assimilation increment 
     58   LOGICAL, PUBLIC :: ln_asmdin     !: No direct initialization 
     59   LOGICAL, PUBLIC :: ln_trainc     !: No tracer (T and S) assimilation increments 
     60   LOGICAL, PUBLIC :: ln_dyninc     !: No dynamics (u and v) assimilation increments 
     61   LOGICAL, PUBLIC :: ln_sshinc     !: No sea surface height assimilation increment 
     62   LOGICAL, PUBLIC :: ln_seaiceinc  !: No sea ice concentration increment 
     63   LOGICAL, PUBLIC :: ln_salfix     !: Apply minimum salinity check 
    6664   LOGICAL, PUBLIC :: ln_temnofreeze = .FALSE. !: Don't allow the temperature to drop below freezing 
    67    INTEGER, PUBLIC :: nn_divdmp                !: Apply divergence damping filter nn_divdmp times 
     65   INTEGER, PUBLIC :: nn_divdmp     !: Apply divergence damping filter nn_divdmp times 
    6866 
    6967   REAL(wp), PUBLIC, DIMENSION(:,:,:), ALLOCATABLE ::   t_bkg   , s_bkg      !: Background temperature and salinity 
     
    9088   !! * Substitutions 
    9189#  include "domzgr_substitute.h90" 
    92 #  include "ldfdyn_substitute.h90" 
    9390#  include "vectopt_loop_substitute.h90" 
    9491   !!---------------------------------------------------------------------- 
     
    139136      ! Read Namelist nam_asminc : assimilation increment interface 
    140137      !----------------------------------------------------------------------- 
    141       ln_seaiceinc = .FALSE. 
     138      ln_seaiceinc   = .FALSE. 
    142139      ln_temnofreeze = .FALSE. 
    143140 
     
    428425 
    429426      IF ( ln_dyninc .AND. nn_divdmp > 0 ) THEN 
    430  
    431          CALL wrk_alloc(jpi,jpj,hdiv)  
    432  
    433          DO  jt = 1, nn_divdmp 
    434  
     427         ! 
     428         CALL wrk_alloc( jpi,jpj,   hdiv )  
     429         ! 
     430         DO jt = 1, nn_divdmp 
     431            ! 
    435432            DO jk = 1, jpkm1 
    436  
    437433               hdiv(:,:) = 0._wp 
    438  
    439434               DO jj = 2, jpjm1 
    440435                  DO ji = fs_2, fs_jpim1   ! vector opt. 
     
    444439                         + e1v(ji  ,jj  ) * fse3v(ji  ,jj  ,jk) * v_bkginc(ji  ,jj  ,jk)     & 
    445440                         - e1v(ji  ,jj-1) * fse3v(ji  ,jj-1,jk) * v_bkginc(ji  ,jj-1,jk)  )  & 
    446                          / ( e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) ) 
     441                         / ( e1e2t(ji,jj) * fse3t(ji,jj,jk) ) 
    447442                  END DO 
    448443               END DO 
    449  
    450444               CALL lbc_lnk( hdiv, 'T', 1. )   ! lateral boundary cond. (no sign change) 
    451  
     445               ! 
    452446               DO jj = 2, jpjm1 
    453447                  DO ji = fs_2, fs_jpim1   ! vector opt. 
    454                      u_bkginc(ji,jj,jk) = u_bkginc(ji,jj,jk) + 0.2_wp * ( e1t(ji+1,jj)*e2t(ji+1,jj) * hdiv(ji+1,jj)   & 
    455                                                                         - e1t(ji  ,jj)*e2t(ji  ,jj) * hdiv(ji  ,jj) ) & 
    456                                                                       / e1u(ji,jj) * umask(ji,jj,jk)  
    457                      v_bkginc(ji,jj,jk) = v_bkginc(ji,jj,jk) + 0.2_wp * ( e1t(ji,jj+1)*e2t(ji,jj+1) * hdiv(ji,jj+1)   & 
    458                                                                         - e1t(ji,jj  )*e2t(ji,jj  ) * hdiv(ji,jj  ) ) & 
    459                                                                       / e2v(ji,jj) * vmask(ji,jj,jk)  
     448                     u_bkginc(ji,jj,jk) = u_bkginc(ji,jj,jk) + 0.2_wp * ( e1e2t(ji+1,jj) * hdiv(ji+1,jj)   & 
     449                        &                                               - e1e2t(ji  ,jj) * hdiv(ji  ,jj) ) & 
     450                        &                                             * r1_e1u(ji,jj) * umask(ji,jj,jk)  
     451                     v_bkginc(ji,jj,jk) = v_bkginc(ji,jj,jk) + 0.2_wp * ( e1e2t(ji,jj+1) * hdiv(ji,jj+1)   & 
     452                        &                                               - e1e2t(ji,jj  ) * hdiv(ji,jj  ) ) & 
     453                        &                                             * r1_e2v(ji,jj) * vmask(ji,jj,jk)  
    460454                  END DO 
    461455               END DO 
    462  
    463456            END DO 
    464  
     457            ! 
    465458         END DO 
    466  
    467          CALL wrk_dealloc(jpi,jpj,hdiv)  
    468  
     459         ! 
     460         CALL wrk_dealloc( jpi,jpj,   hdiv )  
     461         ! 
    469462      ENDIF 
    470  
    471  
    472463 
    473464      !----------------------------------------------------------------------- 
     
    476467 
    477468      IF ( ln_asmdin ) THEN 
    478  
     469         ! 
    479470         ALLOCATE( t_bkg(jpi,jpj,jpk) ) 
    480471         ALLOCATE( s_bkg(jpi,jpj,jpk) ) 
     
    482473         ALLOCATE( v_bkg(jpi,jpj,jpk) ) 
    483474         ALLOCATE( ssh_bkg(jpi,jpj)   ) 
    484  
    485          t_bkg(:,:,:) = 0.0 
    486          s_bkg(:,:,:) = 0.0 
    487          u_bkg(:,:,:) = 0.0 
    488          v_bkg(:,:,:) = 0.0 
    489          ssh_bkg(:,:) = 0.0 
    490  
     475         ! 
     476         t_bkg(:,:,:) = 0._wp 
     477         s_bkg(:,:,:) = 0._wp 
     478         u_bkg(:,:,:) = 0._wp 
     479         v_bkg(:,:,:) = 0._wp 
     480         ssh_bkg(:,:) = 0._wp 
     481         ! 
    491482         !-------------------------------------------------------------------- 
    492483         ! Read from file the background state at analysis time 
    493484         !-------------------------------------------------------------------- 
    494  
     485         ! 
    495486         CALL iom_open( c_asmdin, inum ) 
    496  
     487         ! 
    497488         CALL iom_get( inum, 'rdastp', zdate_bkg )  
    498          
     489         ! 
    499490         IF(lwp) THEN 
    500491            WRITE(numout,*)  
    501             WRITE(numout,*) 'asm_inc_init : Assimilation background state valid at : ', & 
    502                &  NINT( zdate_bkg ) 
     492            WRITE(numout,*) 'asm_inc_init : Assimilation background state valid at : ', NINT( zdate_bkg ) 
    503493            WRITE(numout,*) '~~~~~~~~~~~~' 
    504494         ENDIF 
    505  
     495         ! 
    506496         IF ( NINT( zdate_bkg ) /= iitdin_date ) & 
    507497            & CALL ctl_warn( ' Validity time of assimilation background state does', & 
    508498            &                ' not agree with Direct Initialization time' ) 
    509  
     499         ! 
    510500         IF ( ln_trainc ) THEN    
    511501            CALL iom_get( inum, jpdom_autoglo, 'tn', t_bkg ) 
     
    514504            s_bkg(:,:,:) = s_bkg(:,:,:) * tmask(:,:,:) 
    515505         ENDIF 
    516  
     506         ! 
    517507         IF ( ln_dyninc ) THEN    
    518508            CALL iom_get( inum, jpdom_autoglo, 'un', u_bkg ) 
     
    521511            v_bkg(:,:,:) = v_bkg(:,:,:) * vmask(:,:,:) 
    522512         ENDIF 
    523          
     513         ! 
    524514         IF ( ln_sshinc ) THEN 
    525515            CALL iom_get( inum, jpdom_autoglo, 'sshn', ssh_bkg ) 
    526516            ssh_bkg(:,:) = ssh_bkg(:,:) * tmask(:,:,1) 
    527517         ENDIF 
    528  
     518         ! 
    529519         CALL iom_close( inum ) 
    530  
     520         ! 
    531521      ENDIF 
    532522      ! 
     
    574564      ! If kt = kit000 - 1 then set the date to the restart date 
    575565      IF ( kt == kit000 - 1 ) THEN 
    576  
    577566         kdate = ndastp 
    578567         RETURN 
    579  
    580568      ENDIF 
    581569 
     
    646634      !! ** Action  :  
    647635      !!---------------------------------------------------------------------- 
    648       INTEGER, INTENT(IN) :: kt               ! Current time step 
    649       ! 
    650       INTEGER :: ji,jj,jk 
    651       INTEGER :: it 
     636      INTEGER, INTENT(IN) ::   kt   ! Current time step 
     637      ! 
     638      INTEGER  :: ji, jj, jk 
     639      INTEGER  :: it 
    652640      REAL(wp) :: zincwgt  ! IAU weight for current time step 
    653641      REAL (wp), DIMENSION(jpi,jpj,jpk) :: fzptnz ! 3d freezing point values 
    654642      !!---------------------------------------------------------------------- 
    655  
     643      ! 
    656644      ! freezing point calculation taken from oc_fz_pt (but calculated for all depths)  
    657645      ! used to prevent the applied increments taking the temperature below the local freezing point  
    658  
    659646      DO jk = 1, jpkm1 
    660647        CALL eos_fzp( tsn(:,:,jk,jp_sal), fzptnz(:,:,jk), fsdept(:,:,jk) ) 
    661648      END DO 
    662  
    663       IF ( ln_asmiau ) THEN 
    664  
    665          !-------------------------------------------------------------------- 
    666          ! Incremental Analysis Updating 
    667          !-------------------------------------------------------------------- 
    668  
     649         ! 
     650         !                             !-------------------------------------- 
     651      IF ( ln_asmiau ) THEN            ! Incremental Analysis Updating 
     652         !                             !-------------------------------------- 
     653         ! 
    669654         IF ( ( kt >= nitiaustr_r ).AND.( kt <= nitiaufin_r ) ) THEN 
    670  
     655            ! 
    671656            it = kt - nit000 + 1 
    672657            zincwgt = wgtiau(it) / rdt   ! IAU weight for the current time step 
    673  
     658            ! 
    674659            IF(lwp) THEN 
    675660               WRITE(numout,*)  
     
    677662               WRITE(numout,*) '~~~~~~~~~~~~' 
    678663            ENDIF 
    679  
     664            ! 
    680665            ! Update the tracer tendencies 
    681666            DO jk = 1, jpkm1 
     
    700685               ENDIF 
    701686            END DO 
    702  
    703          ENDIF 
    704  
     687            ! 
     688         ENDIF 
     689         ! 
    705690         IF ( kt == nitiaufin_r + 1  ) THEN   ! For bias crcn to work 
    706691            DEALLOCATE( t_bkginc ) 
    707692            DEALLOCATE( s_bkginc ) 
    708693         ENDIF 
    709  
    710  
    711       ELSEIF ( ln_asmdin ) THEN 
    712  
    713          !-------------------------------------------------------------------- 
    714          ! Direct Initialization 
    715          !-------------------------------------------------------------------- 
    716              
     694         !                             !-------------------------------------- 
     695      ELSEIF ( ln_asmdin ) THEN        ! Direct Initialization 
     696         !                             !-------------------------------------- 
     697         !             
    717698         IF ( kt == nitdin_r ) THEN 
    718  
     699            ! 
    719700            neuler = 0  ! Force Euler forward step 
    720  
     701            ! 
    721702            ! Initialize the now fields with the background + increment 
    722703            IF (ln_temnofreeze) THEN 
     
    745726!!gm 
    746727 
    747  
    748             IF( ln_zps .AND. .NOT. lk_c1d .AND. .NOT. ln_isfcav)      & 
    749                &  CALL zps_hde    ( kt, jpts, tsb, gtsu, gtsv,        &  ! Partial steps: before horizontal gradient 
    750                &                              rhd, gru , grv          )  ! of t, s, rd at the last ocean level 
    751             IF( ln_zps .AND. .NOT. lk_c1d .AND.       ln_isfcav)      & 
    752                &  CALL zps_hde_isf( nit000, jpts, tsb, gtsu, gtsv,    &    ! Partial steps for top cell (ISF) 
    753                &                                  rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
    754                &                           gtui, gtvi, grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi    ) ! of t, s, rd at the last ocean level 
    755  
    756 #if defined key_zdfkpp 
    757             CALL eos( tsn, rhd, fsdept_n(:,:,:) )                      ! Compute rhd 
    758 !!gm fabien            CALL eos( tsn, rhd )                      ! Compute rhd 
    759 #endif 
    760  
     728            IF( ln_zps .AND. .NOT. lk_c1d ) THEN      ! Partial steps: before horizontal gradient 
     729               IF(ln_isfcav) THEN                        ! ocean cavities: top and bottom cells (ISF) 
     730                  CALL zps_hde_isf( nit000, jpts, tsb, gtsu, gtsv, gtui, gtvi,     & 
     731                     &                            rhd, gru , grv , aru , arv , gzu , gzv , ge3ru , ge3rv ,   & 
     732                     &                     grui, grvi, arui, arvi, gzui, gzvi, ge3rui, ge3rvi    ) 
     733               ELSE                                      ! no ocean cavities: bottom cells 
     734                  CALL zps_hde    ( kt, jpts, tsb, gtsu, gtsv,        &  !  
     735                     &                        rhd, gru , grv          )  ! of t, s, rd at the last ocean level 
     736               ENDIF 
     737            ENDIF 
     738            ! 
    761739            DEALLOCATE( t_bkginc ) 
    762740            DEALLOCATE( s_bkginc ) 
     
    767745      ENDIF 
    768746      ! Perhaps the following call should be in step 
    769       IF   ( ln_seaiceinc  )   CALL seaice_asm_inc ( kt )   ! apply sea ice concentration increment 
     747      IF ( ln_seaiceinc  )   CALL seaice_asm_inc ( kt )   ! apply sea ice concentration increment 
    770748      ! 
    771749   END SUBROUTINE tra_asm_inc 
     
    788766      REAL(wp) :: zincwgt  ! IAU weight for current time step 
    789767      !!---------------------------------------------------------------------- 
    790  
    791       IF ( ln_asmiau ) THEN 
    792  
    793          !-------------------------------------------------------------------- 
    794          ! Incremental Analysis Updating 
    795          !-------------------------------------------------------------------- 
    796  
     768      ! 
     769      !                          !-------------------------------------------- 
     770      IF ( ln_asmiau ) THEN      ! Incremental Analysis Updating 
     771         !                       !-------------------------------------------- 
     772         ! 
    797773         IF ( ( kt >= nitiaustr_r ).AND.( kt <= nitiaufin_r ) ) THEN 
    798  
     774            ! 
    799775            it = kt - nit000 + 1 
    800776            zincwgt = wgtiau(it) / rdt   ! IAU weight for the current time step 
    801  
     777            ! 
    802778            IF(lwp) THEN 
    803779               WRITE(numout,*)  
    804                WRITE(numout,*) 'dyn_asm_inc : Dynamics IAU at time step = ', & 
    805                   &  kt,' with IAU weight = ', wgtiau(it) 
     780               WRITE(numout,*) 'dyn_asm_inc : Dynamics IAU at time step = ', kt,' with IAU weight = ', wgtiau(it) 
    806781               WRITE(numout,*) '~~~~~~~~~~~~' 
    807782            ENDIF 
    808  
     783            ! 
    809784            ! Update the dynamic tendencies 
    810785            DO jk = 1, jpkm1 
     
    812787               va(:,:,jk) = va(:,:,jk) + v_bkginc(:,:,jk) * zincwgt 
    813788            END DO 
    814             
     789            ! 
    815790            IF ( kt == nitiaufin_r ) THEN 
    816791               DEALLOCATE( u_bkginc ) 
    817792               DEALLOCATE( v_bkginc ) 
    818793            ENDIF 
    819  
    820          ENDIF 
    821  
    822       ELSEIF ( ln_asmdin ) THEN  
    823  
    824          !-------------------------------------------------------------------- 
    825          ! Direct Initialization 
    826          !-------------------------------------------------------------------- 
    827           
     794            ! 
     795         ENDIF 
     796         !                          !----------------------------------------- 
     797      ELSEIF ( ln_asmdin ) THEN     ! Direct Initialization 
     798         !                          !----------------------------------------- 
     799         !          
    828800         IF ( kt == nitdin_r ) THEN 
    829  
     801            ! 
    830802            neuler = 0                    ! Force Euler forward step 
    831  
     803            ! 
    832804            ! Initialize the now fields with the background + increment 
    833805            un(:,:,:) = u_bkg(:,:,:) + u_bkginc(:,:,:) 
    834806            vn(:,:,:) = v_bkg(:,:,:) + v_bkginc(:,:,:)   
    835  
     807            ! 
    836808            ub(:,:,:) = un(:,:,:)         ! Update before fields 
    837809            vb(:,:,:) = vn(:,:,:) 
    838   
     810            ! 
    839811            DEALLOCATE( u_bkg    ) 
    840812            DEALLOCATE( v_bkg    ) 
     
    864836      REAL(wp) :: zincwgt  ! IAU weight for current time step 
    865837      !!---------------------------------------------------------------------- 
    866  
    867       IF ( ln_asmiau ) THEN 
    868  
    869          !-------------------------------------------------------------------- 
    870          ! Incremental Analysis Updating 
    871          !-------------------------------------------------------------------- 
    872  
     838      ! 
     839      !                             !----------------------------------------- 
     840      IF ( ln_asmiau ) THEN         ! Incremental Analysis Updating 
     841         !                          !----------------------------------------- 
     842         ! 
    873843         IF ( ( kt >= nitiaustr_r ).AND.( kt <= nitiaufin_r ) ) THEN 
    874  
     844            ! 
    875845            it = kt - nit000 + 1 
    876846            zincwgt = wgtiau(it) / rdt   ! IAU weight for the current time step 
    877  
     847            ! 
    878848            IF(lwp) THEN 
    879849               WRITE(numout,*)  
     
    882852               WRITE(numout,*) '~~~~~~~~~~~~' 
    883853            ENDIF 
    884  
     854            ! 
    885855            ! Save the tendency associated with the IAU weighted SSH increment 
    886856            ! (applied in dynspg.*) 
     
    891861               DEALLOCATE( ssh_bkginc ) 
    892862            ENDIF 
    893  
    894          ENDIF 
    895  
    896       ELSEIF ( ln_asmdin ) THEN 
    897  
    898          !-------------------------------------------------------------------- 
    899          ! Direct Initialization 
    900          !-------------------------------------------------------------------- 
    901  
     863            ! 
     864         ENDIF 
     865         !                          !----------------------------------------- 
     866      ELSEIF ( ln_asmdin ) THEN     ! Direct Initialization 
     867         !                          !----------------------------------------- 
     868         ! 
    902869         IF ( kt == nitdin_r ) THEN 
    903  
    904             neuler = 0                    ! Force Euler forward step 
    905  
    906             ! Initialize the now fields the background + increment 
    907             sshn(:,:) = ssh_bkg(:,:) + ssh_bkginc(:,:)   
    908  
    909             ! Update before fields 
    910             sshb(:,:) = sshn(:,:)          
    911  
     870            ! 
     871            neuler = 0                                   ! Force Euler forward step 
     872            ! 
     873            sshn(:,:) = ssh_bkg(:,:) + ssh_bkginc(:,:)   ! Initialize the now fields the background + increment 
     874            ! 
     875            sshb(:,:) = sshn(:,:)                        ! Update before fields 
     876            ! 
    912877            IF( lk_vvl ) THEN 
    913878               DO jk = 1, jpk