New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 8882 for branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90 – NEMO

Ignore:
Timestamp:
2017-12-01T18:44:09+01:00 (6 years ago)
Author:
flavoni
Message:

dev_CNRS_2017 branch: merged dev_r7881_ENHANCE09_RK3 with trunk r8864

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90

    r7753 r8882  
    2424   USE lib_mpp         ! distribued memory computing library 
    2525   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    26    USE wrk_nemo        ! Memory Allocation 
    2726 
    2827   IMPLICIT NONE 
     
    3332 
    3433   !                                                !!* Namelist namdyn_ldf : lateral mixing on momentum * 
     34   LOGICAL , PUBLIC ::   ln_dynldf_NONE  !: No operator (i.e. no explicit diffusion) 
    3535   LOGICAL , PUBLIC ::   ln_dynldf_lap   !: laplacian operator 
    3636   LOGICAL , PUBLIC ::   ln_dynldf_blp   !: bilaplacian operator 
     
    9696      REAL(wp) ::   zah0              ! local scalar 
    9797      ! 
    98       NAMELIST/namdyn_ldf/ ln_dynldf_lap, ln_dynldf_blp,                  & 
    99          &                 ln_dynldf_lev, ln_dynldf_hor, ln_dynldf_iso,   & 
    100          &                 nn_ahm_ijk_t , rn_ahm_0, rn_ahm_b, rn_bhm_0,   & 
    101          &                 rn_csmc      , rn_minfac, rn_maxfac 
     98      NAMELIST/namdyn_ldf/ ln_dynldf_NONE, ln_dynldf_lap, ln_dynldf_blp,   &   ! type of operator 
     99         &                 ln_dynldf_lev, ln_dynldf_hor, ln_dynldf_iso ,   &   ! acting direction of the operator 
     100         &                 nn_ahm_ijk_t , rn_ahm_0, rn_ahm_b, rn_bhm_0 ,   &   ! lateral eddy coefficient 
     101         &                 rn_csmc      , rn_minfac, rn_maxfac                 ! Smagorinsky settings 
    102102      !!---------------------------------------------------------------------- 
    103103      ! 
     
    118118         ! 
    119119         WRITE(numout,*) '      type :' 
     120         WRITE(numout,*) '         no explicit diffusion                ln_dynldf_NONE= ', ln_dynldf_NONE 
    120121         WRITE(numout,*) '         laplacian operator                   ln_dynldf_lap = ', ln_dynldf_lap 
    121122         WRITE(numout,*) '         bilaplacian operator                 ln_dynldf_blp = ', ln_dynldf_blp 
     
    131132         WRITE(numout,*) '         background viscosity (iso case)      rn_ahm_b      = ', rn_ahm_b, ' m2/s' 
    132133         WRITE(numout,*) '         lateral bilaplacian eddy viscosity   rn_bhm_0      = ', rn_bhm_0, ' m4/s' 
    133          WRITE(numout,*) '      smagorinsky settings (nn_ahm_ijk_t  = 32) :' 
     134         WRITE(numout,*) '      Smagorinsky settings (nn_ahm_ijk_t  = 32) :' 
    134135         WRITE(numout,*) '         Smagorinsky coefficient              rn_csmc       = ', rn_csmc 
    135136         WRITE(numout,*) '         factor multiplier for theorectical lower limit for ' 
     
    140141 
    141142      !                                ! Parameter control 
    142       IF( .NOT.ln_dynldf_lap .AND. .NOT.ln_dynldf_blp ) THEN 
     143      IF( ln_dynldf_NONE ) THEN 
    143144         IF(lwp) WRITE(numout,*) '   No viscous operator selected. ahmt and ahmf are not allocated' 
    144145         l_ldfdyn_time = .FALSE. 
     
    284285      !!---------------------------------------------------------------------- 
    285286      ! 
    286       IF( nn_timing == 1 )  CALL timing_start('ldf_dyn') 
     287      IF( ln_timing )   CALL timing_start('ldf_dyn') 
    287288      ! 
    288289      SELECT CASE(  nn_ahm_ijk_t  )       !== Eddy vicosity coefficients ==! 
     
    411412      CALL iom_put( "ahmf_3d", ahmf(:,:,:) )   ! 3D      v-eddy diffusivity coeff. 
    412413      ! 
    413       IF( nn_timing == 1 )  CALL timing_stop('ldf_dyn') 
     414      IF( ln_timing )   CALL timing_stop('ldf_dyn') 
    414415      ! 
    415416   END SUBROUTINE ldf_dyn 
Note: See TracChangeset for help on using the changeset viewer.