Changeset 5870


Ignore:
Timestamp:
2015-11-09T18:33:54+01:00 (5 years ago)
Author:
acc
Message:

Branch 2015/dev_r5803_NOC_WAD. Merge in trunk changes from 5803 to 5869 in preparation for merge. Also tidied and reorganised some wetting and drying code. Renamed wadlmt.F90 to wetdry.F90. Wetting drying code changes restricted to domzgr.F90, domvvl.F90 nemogcm.F90 sshwzv.F90, dynspg_ts.F90, wetdry.F90 and dynhpg.F90. Code passes full SETTE tests with ln_wd=.false.. Still awaiting test case for checking with ln_wd=.false.

Location:
branches/2015/dev_r5803_NOC_WAD/NEMOGCM
Files:
31 deleted
218 edited
8 copied
1 moved

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/ARCH/arch-X64_MOBILIS.fcm

    r5656 r5870  
    3838%NCDF_HOME           /home/acc/shared 
    3939%HDF5_HOME           /home/acc/shared 
    40 %XIOS_HOME           /home/acc/XIOS_1.0 
     40%XIOS_HOME           /home/acc/XIOS_1.0_r777 
    4141%OASIS_HOME           
    4242 
     
    6161%USER_LIB            %XIOS_LIB %OASIS_LIB %NCDF_LIB 
    6262 
    63 %CC                  cc 
     63%CC                  icc 
    6464%CFLAGS              -O0 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/AMM12/EXP00/namelist_cfg

    r5501 r5870  
    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") 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/AMM12/cpp_AMM12.fcm

    r5842 r5870  
    1  bld::tool::fppkeys  key_bdy key_tide key_dynspg_ts key_ldfslp  key_zdfgls  key_vvl key_diainstant key_mpp_mpi  
     1 bld::tool::fppkeys  key_bdy key_tide key_dynspg_ts  key_zdfgls  key_vvl key_diainstant key_mpp_mpi key_iomput 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/C1D_PAPA/EXP00/namelist_cfg

    r5407 r5870  
    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!----------------------------------------------------------------------- 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/GYRE/EXP00/namelist_cfg

    r5407 r5870  
    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!----------------------------------------------------------------------- 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/GYRE/cpp_GYRE.fcm

    r4990 r5870  
    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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist_cfg

    r5407 r5870  
    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!----------------------------------------------------------------------- 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/GYRE_BFM/EXP00/namelist_top_cfg

    r4152 r5870  
    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!----------------------------------------------------------------------- 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/GYRE_BFM/cpp_GYRE_BFM.fcm

    r4230 r5870  
    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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/GYRE_PISCES/EXP00/namelist_cfg

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

    r4340 r5870  
    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!----------------------------------------------------------------------- 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/GYRE_PISCES/cpp_GYRE_PISCES.fcm

    r4990 r5870  
    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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/GYRE_XIOS/EXP00/namelist_cfg

    r5407 r5870  
    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!----------------------------------------------------------------------- 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/GYRE_XIOS/cpp_GYRE_XIOS.fcm

    r4373 r5870  
    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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/1_namelist_cfg

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

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

    r5385 r5870  
    1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi key_diaobs key_asminc  
     1bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_zdftke key_zdfddm key_zdftmx  key_mpp_mpi key_iomput 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/1_namelist_cfg

    r5499 r5870  
    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!----------------------------------------------------------------------- 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_LIM3/EXP00/namelist_cfg

    r4995 r5870  
    9393/ 
    9494!----------------------------------------------------------------------- 
    95 &namcla        !   cross land advection 
    96 !----------------------------------------------------------------------- 
    97 / 
    98 !----------------------------------------------------------------------- 
    9995&nambfr        !   bottom friction 
    10096!----------------------------------------------------------------------- 
     
    115111&namtra_adv    !   advection scheme for tracer 
    116112!----------------------------------------------------------------------- 
     113   ln_traadv_fct =  .true.    !  FCT scheme 
     114      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     115      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     116      nn_fct_zts =  0               !  > 1 , 2nd order FCT scheme with vertical sub-timestepping 
     117      !                             !        (number of sub-timestep = nn_fct_zts) 
    117118/ 
    118119!----------------------------------------------------------------------- 
     
    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/ 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_LIM3/cpp_ORCA2_LIM3.fcm

    r5385 r5870  
    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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/1_namelist_cfg

    r5407 r5870  
    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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/1_namelist_top_cfg

    r4147 r5870  
    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!----------------------------------------------------------------------- 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/namelist_cfg

    r5407 r5870  
    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!----------------------------------------------------------------------- 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/EXP00/namelist_top_cfg

    r4340 r5870  
    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!----------------------------------------------------------------------- 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_LIM_CFC_C14b/cpp_ORCA2_LIM_CFC_C14b.fcm

    r4523 r5870  
    1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_top key_cfc key_c14b key_iomput key_mpp_mpi 
     1bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_zdftke key_zdfddm key_zdftmx key_top key_cfc key_c14b key_iomput key_mpp_mpi 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_cfg

    r4370 r5870  
    9393/ 
    9494!----------------------------------------------------------------------- 
    95 &namcla        !   cross land advection 
    96 !----------------------------------------------------------------------- 
    97 / 
    98 !----------------------------------------------------------------------- 
    9995&nambfr        !   bottom friction 
    10096!----------------------------------------------------------------------- 
     
    115111&namtra_adv    !   advection scheme for tracer 
    116112!----------------------------------------------------------------------- 
     113   ln_traadv_fct =  .true.    !  FCT scheme 
     114      nn_fct_h   =  2               !  =2/4, horizontal 2nd / 4th order  
     115      nn_fct_v   =  2               !  =2/4, vertical   2nd / COMPACT 4th order  
     116      nn_fct_zts =  0               !  > 1 , 2nd order FCT scheme with vertical sub-timestepping 
     117      !                             !        (number of sub-timestep = nn_fct_zts) 
    117118/ 
    118119!----------------------------------------------------------------------- 
     
    120121!----------------------------------------------------------------------- 
    121122/ 
    122 !----------------------------------------------------------------------- 
    123 &namtra_adv_mle !  mixed layer eddy parametrisation (Fox-Kemper param) 
    124 !----------------------------------------------------------------------- 
    125 / 
    126123!---------------------------------------------------------------------------------- 
    127124&namtra_ldf    !   lateral diffusion scheme for tracers 
    128125!---------------------------------------------------------------------------------- 
     126   !                       !  Operator type: 
     127   ln_traldf_lap   =  .true.   !    laplacian operator 
     128   ln_traldf_blp   =  .false.  !  bilaplacian operator 
     129   !                       !  Direction of action: 
     130   ln_traldf_lev   =  .false.  !  iso-level 
     131   ln_traldf_hor   =  .false.  !  horizontal (geopotential) 
     132   ln_traldf_iso   =  .true.   !  iso-neutral (standard operator) 
     133   ln_traldf_triad =  .false.  !  iso-neutral (triad    operator) 
     134   ! 
     135   !                       !  iso-neutral options:         
     136   ln_traldf_msc   =  .true.   !  Method of Stabilizing Correction (both operators) 
     137   rn_slpmax       =   0.01    !  slope limit                      (both operators) 
     138   ln_triad_iso    =  .false.  !  pure horizontal mixing in ML              (triad only) 
     139   rn_sw_triad     =  1        !  =1 switching triad ; =0 all 4 triads used (triad only) 
     140   ln_botmix_triad =  .false.  !  lateral mixing on bottom                  (triad only) 
     141   ! 
     142   !                       !  Coefficients: 
     143   nn_aht_ijk_t    = 20        !  space/time variation of eddy coef 
     144   !                                !   =-20 (=-30)    read in eddy_diffusivity_2D.nc (..._3D.nc) file 
     145   !                                !   =  0           constant 
     146   !                                !   = 10 F(k)      =ldf_c1d 
     147   !                                !   = 20 F(i,j)    =ldf_c2d 
     148   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     149   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
     150   !                                !   = 31 F(i,j,k,t)=F(local velocity) 
     151   rn_aht_0        = 2000.     !  lateral eddy diffusivity   (lap. operator) [m2/s] 
     152   rn_bht_0        = 1.e+12    !  lateral eddy diffusivity (bilap. operator) [m4/s] 
     153/ 
     154!---------------------------------------------------------------------------------- 
     155&namtra_ldfeiv !   eddy induced velocity param. 
     156!---------------------------------------------------------------------------------- 
     157   ln_ldfeiv     =.true.   ! use eddy induced velocity parameterization 
     158   ln_ldfeiv_dia =.true.   ! diagnose eiv stream function and velocities 
     159   rn_aeiv_0     = 2000.   ! eddy induced velocity coefficient   [m2/s] 
     160   nn_aei_ijk_t  = 21      ! space/time variation of the eiv coeficient 
     161   !                                !   =-20 (=-30)    read in eddy_induced_velocity_2D.nc (..._3D.nc) file 
     162   !                                !   =  0           constant 
     163   !                                !   = 10 F(k)      =ldf_c1d 
     164   !                                !   = 20 F(i,j)    =ldf_c2d 
     165   !                                !   = 21 F(i,j,t)  =Treguier et al. JPO 1997 formulation 
     166   !                                !   = 30 F(i,j,k)  =ldf_c2d + ldf_c1d 
    129167/ 
    130168!----------------------------------------------------------------------- 
     
    139177&namdyn_vor    !   option of physics/algorithm (not control by CPP keys) 
    140178!----------------------------------------------------------------------- 
     179   ln_dynvor_ene = .false. !  enstrophy conserving scheme 
     180   ln_dynvor_ens = .false. !  energy conserving scheme 
     181   ln_dynvor_mix = .false. !  mixed scheme 
     182   ln_dynvor_een = .true.  !  energy & enstrophy scheme 
     183      nn_een_e3f = 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!----------------------------------------------------------------------- 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/EXP00/namelist_top_cfg

    r5385 r5870  
    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!----------------------------------------------------------------------- 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_LIM_PISCES/cpp_ORCA2_LIM_PISCES.fcm

    r4523 r5870  
    1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_cfg

    r5385 r5870  
    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") 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/EXP00/namelist_top_cfg

    r5385 r5870  
    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!----------------------------------------------------------------------- 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_OFF_PISCES/cpp_ORCA2_OFF_PISCES.fcm

    r4523 r5870  
    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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/ORCA2_SAS_LIM/EXP00/namelist_cfg

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

    r4523 r5870  
    1  bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_diaeiv key_ldfslp key_traldf_c2d key_traldf_eiv key_dynldf_c3d key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi 
     1 bld::tool::fppkeys key_trabbl key_lim2 key_dynspg_flt key_zdftke key_zdfddm key_zdftmx key_iomput key_mpp_mpi 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/SHARED/domain_def.xml

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

    r5517 r5870  
    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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/SHARED/namelist_ref

    r5656 r5870  
     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!----------------------------------------------------------------------- 
     
    12931283   rn_htrmax         =  200.0   ! max. depth of transition range 
    12941284/ 
     1285!----------------------------------------------------------------------- 
     1286&namwad  !   Wetting and drying 
     1287!----------------------------------------------------------------------- 
     1288   ln_wd             = .false.  ! T/F activation of wetting and drying 
     1289   rn_wdmin1         =  0.1     ! Minimum wet depth on dried cells 
     1290   rn_wdmin2         =  0.01    ! Tolerance of min wet depth on dried cells 
     1291   rn_wdld           =  20.0    ! Land elevation below which wetting/dryins is allowed 
     1292   nn_wdit           =  10      ! Max iterations for W/D limiter 
     1293/ 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/CONFIG/SHARED/namelist_top_ref

    r5416 r5870  
    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!----------------------------------------------------------------------- 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_2/limdmp_2.F90

    r4624 r5870  
    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) )  )  ) 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_2/limrhg_2.F90

    r5123 r5870  
    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                  ! 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_2/limsbc_2.F90

    r5407 r5870  
    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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_2/limthd_2.F90

    r5407 r5870  
    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) 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_2/limthd_lac_2.F90

    r3625 r5870  
    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 ) ) ) 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_3/limadv.F90

    r5429 r5870  
    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               
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_3/limcons.F90

    r5183 r5870  
    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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_3/limctl.F90

    r5167 r5870  
    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)        
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_3/limdiahsb.F90

    r5215 r5870  
    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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_3/limdyn.F90

    r5123 r5870  
    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      :') 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_3/limhdf.F90

    r5429 r5870  
    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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90

    r5202 r5870  
    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      :') 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_3/limrhg.F90

    r5429 r5870  
    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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90

    r5407 r5870  
    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      :') 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_3/limtrp.F90

    r5202 r5870  
    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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_3/limupdate1.F90

    r5215 r5870  
    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        :') 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/LIM_SRC_3/limupdate2.F90

    r5410 r5870  
    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        :') 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/NST_SRC/agrif_opa_sponge.F90

    r5656 r5870  
    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) & 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/NST_SRC/agrif_top_sponge.F90

    r5656 r5870  
    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)  ) 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/OFF_SRC/domrea.F90

    r5504 r5870  
    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      ! 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/OFF_SRC/dtadyn.F90

    r5768 r5870  
    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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/OFF_SRC/nemogcm.F90

    r5504 r5870  
    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      ! 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/ASM/asmbkg.F90

    r5215 r5870  
    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 
  • branches/2015/dev_r5803_NOC_WAD/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

    r5541 r5870  
    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)