- Timestamp:
- 2017-12-01T18:44:09+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_CNRS_2017/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90
r7753 r8882 24 24 USE lib_mpp ! distribued memory computing library 25 25 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 26 USE wrk_nemo ! Memory Allocation27 26 28 27 IMPLICIT NONE … … 33 32 34 33 ! !!* Namelist namdyn_ldf : lateral mixing on momentum * 34 LOGICAL , PUBLIC :: ln_dynldf_NONE !: No operator (i.e. no explicit diffusion) 35 35 LOGICAL , PUBLIC :: ln_dynldf_lap !: laplacian operator 36 36 LOGICAL , PUBLIC :: ln_dynldf_blp !: bilaplacian operator … … 96 96 REAL(wp) :: zah0 ! local scalar 97 97 ! 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 102 102 !!---------------------------------------------------------------------- 103 103 ! … … 118 118 ! 119 119 WRITE(numout,*) ' type :' 120 WRITE(numout,*) ' no explicit diffusion ln_dynldf_NONE= ', ln_dynldf_NONE 120 121 WRITE(numout,*) ' laplacian operator ln_dynldf_lap = ', ln_dynldf_lap 121 122 WRITE(numout,*) ' bilaplacian operator ln_dynldf_blp = ', ln_dynldf_blp … … 131 132 WRITE(numout,*) ' background viscosity (iso case) rn_ahm_b = ', rn_ahm_b, ' m2/s' 132 133 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) :' 134 135 WRITE(numout,*) ' Smagorinsky coefficient rn_csmc = ', rn_csmc 135 136 WRITE(numout,*) ' factor multiplier for theorectical lower limit for ' … … 140 141 141 142 ! ! Parameter control 142 IF( .NOT.ln_dynldf_lap .AND. .NOT.ln_dynldf_blp) THEN143 IF( ln_dynldf_NONE ) THEN 143 144 IF(lwp) WRITE(numout,*) ' No viscous operator selected. ahmt and ahmf are not allocated' 144 145 l_ldfdyn_time = .FALSE. … … 284 285 !!---------------------------------------------------------------------- 285 286 ! 286 IF( nn_timing == 1 )CALL timing_start('ldf_dyn')287 IF( ln_timing ) CALL timing_start('ldf_dyn') 287 288 ! 288 289 SELECT CASE( nn_ahm_ijk_t ) !== Eddy vicosity coefficients ==! … … 411 412 CALL iom_put( "ahmf_3d", ahmf(:,:,:) ) ! 3D v-eddy diffusivity coeff. 412 413 ! 413 IF( nn_timing == 1 )CALL timing_stop('ldf_dyn')414 IF( ln_timing ) CALL timing_stop('ldf_dyn') 414 415 ! 415 416 END SUBROUTINE ldf_dyn
Note: See TracChangeset
for help on using the changeset viewer.