Changeset 12788 for NEMO/branches/UKMO/NEMO_4.0.1_biharmonic_GM/src/OCE/DYN/dynldf_lap_blp.F90

Timestamp:

2020-04-21T13:21:43+02:00 (4 years ago)

Author:

davestorkey

Message:

UKMO/NEMO_4.0.1_biharmonic_GM : New calculation for zmu with new rn_bgm_msc namelist parameter.

File:

• NEMO/branches/UKMO/NEMO_4.0.1_biharmonic_GM/src/OCE/DYN/dynldf_lap_blp.F90 (modified) (2 diffs)

NEMO/branches/UKMO/NEMO_4.0.1_biharmonic_GM/src/OCE/DYN/dynldf_lap_blp.F90

@@ -146 +146 @@
       CALL dyn_ldf_lap( kt, tmask, fmask, pub, pvb, zulap, zvlap, 1 )   
       !
-      ! Calculate ratio of bhm and avm for stabilising correction terms
+      ! Calculate zmu for stabilising correction terms
       ! and add to avm to be included in the vertical diffusion calculation later.
       zmu(:,:,:) = 0.0_wp
@@ -152 +152 @@
          DO jj = 2, jpjm1
             DO ji = fs_2, jpim1   ! vector opt.
-               zmu(ji,jj,jk) = ( bhm(ji,jj,jk) / (avm(ji,jj,jk) + rsmall) ) * wmask(ji,jj,jk)
+               zmu(ji,jj,jk) = ( rn_bgm_msc * bhm(ji,jj,jk) / e1e2t(ji,jj) ) * wmask(ji,jj,jk)
                avm(ji,jj,jk) = avm(ji,jj,jk) + zmu(ji,jj,jk)
             ENDDO