Changeset 12788 for NEMO/branches
- Timestamp:
- 2020-04-21T13:21:43+02:00 (4 years ago)
- Location:
- NEMO/branches/UKMO/NEMO_4.0.1_biharmonic_GM
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0.1_biharmonic_GM/cfgs/SHARED/namelist_ref
r12535 r12788 949 949 ! ! iso-neutral laplacian operator (ln_dynldf_iso=T) : 950 950 rn_ahm_b = 0.0 ! background eddy viscosity [m2/s] 951 nn_bhm_ijk_t = 0 ! space/time variation of bilaplacian GM coefficient :951 nn_bhm_ijk_t = 0 ! space/time variation of bilaplacian GM coefficient : 952 952 ! ! = 11 uniform coefficient in interior, zero top and bottom 953 953 ! ! = 12 linear tapering with dz to surface, zero at top and bottom … … 955 955 rn_bhm_b = 0.0 ! bilaplacian GM coefficient 956 956 rn_bgmzcrit = 0.0 ! critical depth for depth variation of bilaplacian GM coefficient 957 rn_bgm_msc = 16.0 ! multiplicative parameter for Method of Stabilising Correction (BGM) 957 958 / 958 959 !----------------------------------------------------------------------- -
NEMO/branches/UKMO/NEMO_4.0.1_biharmonic_GM/src/OCE/DYN/dynldf_lap_blp.F90
r12777 r12788 146 146 CALL dyn_ldf_lap( kt, tmask, fmask, pub, pvb, zulap, zvlap, 1 ) 147 147 ! 148 ! Calculate ratio of bhm and avmfor stabilising correction terms148 ! Calculate zmu for stabilising correction terms 149 149 ! and add to avm to be included in the vertical diffusion calculation later. 150 150 zmu(:,:,:) = 0.0_wp … … 152 152 DO jj = 2, jpjm1 153 153 DO ji = fs_2, jpim1 ! vector opt. 154 zmu(ji,jj,jk) = ( bhm(ji,jj,jk) / (avm(ji,jj,jk) + rsmall) ) * wmask(ji,jj,jk)154 zmu(ji,jj,jk) = ( rn_bgm_msc * bhm(ji,jj,jk) / e1e2t(ji,jj) ) * wmask(ji,jj,jk) 155 155 avm(ji,jj,jk) = avm(ji,jj,jk) + zmu(ji,jj,jk) 156 156 ENDDO -
NEMO/branches/UKMO/NEMO_4.0.1_biharmonic_GM/src/OCE/LDF/ldfdyn.F90
r12777 r12788 55 55 REAL(wp), PUBLIC :: rn_bhm_b !: lateral bi-Laplacian GM background eddy viscosity [m4/s] 56 56 REAL(wp), PUBLIC :: rn_bgmzcrit !: critical depth (>=0) for linear tapering of bi-Lap. GM coef. to zero at surface [m] 57 REAL(wp), PUBLIC :: rn_bgm_msc !: multiplicative parameter for Method of Stabilising Correction: zero = no correction term. 57 58 ! !!* Parameter to control the type of lateral viscous operator 58 59 INTEGER, PARAMETER, PUBLIC :: np_ERROR =-10 !: error in setting the operator … … 125 126 & ln_dynldf_lev, ln_dynldf_hor, ln_dynldf_iso, & ! acting direction of the operator 126 127 & nn_ahm_ijk_t , rn_Uv , rn_Lv, rn_ahm_b, & ! lateral eddy coefficient 127 & nn_bhm_ijk_t, rn_bhm_b, rn_bgmzcrit, &! lateral bi-Lap. GM coefficient128 & nn_bhm_ijk_t, rn_bhm_b, rn_bgmzcrit, rn_bgm_msc, & ! lateral bi-Lap. GM coefficient 128 129 & rn_csmc , rn_minfac , rn_maxfac ! Smagorinsky settings 129 130 !!---------------------------------------------------------------------- … … 150 151 WRITE(numout,*) ' bi-Laplacian GM operator ln_dynldf_bgm = ', ln_dynldf_bgm 151 152 WRITE(numout,*) ' critical depth for taper (if used) rn_bgmzcrit = ', rn_bgmzcrit 153 WRITE(numout,*) ' MSC parameter rn_bgm_msc = ', rn_bgm_msc 152 154 ! 153 155 WRITE(numout,*) ' direction of action :'
Note: See TracChangeset
for help on using the changeset viewer.