Changeset 7991 for branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfric.F90

Timestamp:

2017-05-01T16:21:42+02:00 (7 years ago)

Author:

gm

Message:

#1880 (HPC-09) - step-5: OPA/ZDF shear production term computed outside closure schemes

File:

• branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfric.F90 (modified) (3 diffs)

branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfric.F90

@@ -21 +21 @@
    USE oce            ! ocean dynamics and tracers variables
    USE dom_oce        ! ocean space and time domain variables
-   USE zdf_oce        ! ocean vertical physics
+   USE zdf_oce        ! vertical physics: variables
+   USE zdfsh2         ! vertical physics: shear production term of TKE
    USE phycst         ! physical constants
    USE sbc_oce,  ONLY :   taum
@@ -104 +105 @@
       INTEGER  ::   ji, jj, jk                 ! dummy loop indices
       LOGICAL  ::   ll_av_weight = .TRUE.      ! local logical
-      REAL(wp) ::   zsh2, zcfRi, zav, zustar   ! local scalars
-      REAL(wp), DIMENSION(jpi,jpj) ::   zdu2, zdv2, zh_ekm   ! 2D workspace
+      REAL(wp) ::   zcfRi, zav, zustar   ! local scalars
+      REAL(wp), DIMENSION(jpi,jpj)     ::   zh_ekm   ! 2D workspace
+      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zsh2     ! 3D     -
       !!----------------------------------------------------------------------
       !
       IF( nn_timing == 1 )   CALL timing_start('zdf_ric')
       !
-      DO jk = 2, jpkm1        !==  avm and avt = F(Richardson number)  ==!
-         !
-      IF( ll_av_weight ) THEN    !== viscosity weighted shear & stratification terms
-         !
-         DO jj = 1, jpjm1           !* viscosity weighted shear term
+      !                       !==  avm and avt = F(Richardson number)  ==!
+      !
+      !                          !* Shear production at uw- and vw-points (energy conserving form)
+      CALL zdf_sh2( zsh2 )
+      !
+      DO jk = 2, jpkm1           !* Vertical eddy viscosity and diffusivity coefficients
+         DO jj = 1, jpjm1
             DO ji = 1, jpim1
-               zdu2(ji,jj) = 0.5 * ( avm(ji,jj,jk  ) + avm(ji+1,jj,jk) ) * wumask(ji,jj,jk)      &
-                  &              * (  un(ji,jj,jk-1) -  un(ji  ,jj,jk) )                         &
-                  &              * (  ub(ji,jj,jk-1) -  ub(ji  ,jj,jk) ) / (  e3uw_n(ji,jj,jk) * e3uw_b(ji,jj,jk) )
-               zdv2(ji,jj) = 0.5 * ( avm(ji,jj,jk  ) + avm(ji,jj+1,jk) ) * wumask(ji,jj,jk)      &
-                  &              * (  vn(ji,jj,jk-1) -  vn(ji,jj  ,jk) )                         &
-                  &              * (  vb(ji,jj,jk-1) -  vb(ji,jj  ,jk) ) / (  e3vw_n(ji,jj,jk) * e3vw_b(ji,jj,jk) )
-            END DO
-         END DO
-         DO jj = 2, jpjm1
-            DO ji = 2, jpim1        !* Richardson number dependent coefficient
-               !                          ! shear of horizontal velocity
-               zsh2  =  ( zdu2(ji-1,jj) + zdu2(ji,jj) ) / MAX( 1._wp , umask(ji-1,jj,jk) + umask(ji,jj,jk) )   &
-                  &   + ( zdv2(ji,jj-1) + zdv2(ji,jj) ) / MAX( 1._wp , vmask(ji,jj-1,jk) + vmask(ji,jj,jk) )    
-               !                          ! coefficient = F(richardson number)
-               zcfRi = 1._wp / (  1._wp + rn_alp * MAX(  0._wp , avt(ji,jj,jk) * rn2(ji,jj,jk) / ( zsh2 + 1.e-20 ) )  )
+               !                          ! coefficient = F(richardson number) (avm-weighted Ri)
+               zcfRi = 1._wp / (  1._wp + rn_alp * MAX(  0._wp , avm(ji,jj,jk) * rn2(ji,jj,jk) / ( zsh2(ji,jj,jk) + 1.e-20 ) )  )
+               zav   = rn_avmri * zcfRi**nn_ric
                !
-               !                    !* Vertical eddy viscosity and diffusivity coefficients
-               zav = rn_avmri * zcfRi**nn_ric
+               !                          ! avm and avt coefficients
                avm(ji,jj,jk) = MAX(  zav         , avmb(jk)  ) * wmask(ji,jj,jk)
                avt(ji,jj,jk) = MAX(  zav * zcfRi , avtb(jk)  ) * wmask(ji,jj,jk)
             END DO
          END DO
-      ELSE                    !==  classical Richardson number definition  ==!
-         DO jj = 1, jpjm1           !* shear term
-            DO ji = 1, jpim1
-               zdu2(ji,jj) = 0.5 * (  un(ji,jj,jk-1) -  un(ji  ,jj,jk) )                         &
-                  &              * (  ub(ji,jj,jk-1) -  ub(ji  ,jj,jk) ) / (  e3uw_n(ji,jj,jk) * e3uw_b(ji,jj,jk) )
-               zdv2(ji,jj) = 0.5 * (  vn(ji,jj,jk-1) -  vn(ji,jj  ,jk) )                         &
-                  &              * (  vb(ji,jj,jk-1) -  vb(ji,jj  ,jk) ) / (  e3vw_n(ji,jj,jk) * e3vw_b(ji,jj,jk) )
-            END DO
-         END DO
-         DO jj = 2, jpjm1
-            DO ji = 2, jpim1        !* Richardson number dependent coefficient
-               !                          ! shear of horizontal velocity
-               zsh2  =  ( zdu2(ji-1,jj) + zdu2(ji,jj) ) / MAX( 1._wp , umask(ji-1,jj,jk) + umask(ji,jj,jk) )   &
-                  &   + ( zdv2(ji,jj-1) + zdv2(ji,jj) ) / MAX( 1._wp , vmask(ji,jj-1,jk) + vmask(ji,jj,jk) )    
-               !                          ! coefficient = F(richardson number)
-               zcfRi = 1._wp / (  1._wp + rn_alp * MAX(  0._wp , rn2(ji,jj,jk) / ( zsh2 + 1.e-20 ) )  )
-               !
-               !                    !* Vertical eddy viscosity and diffusivity coefficients
-               zav = rn_avmri * zcfRi**nn_ric
-               avm(ji,jj,jk) = MAX(  zav         , avmb(jk)  ) * wmask(ji,jj,jk)
-               avt(ji,jj,jk) = MAX(  zav * zcfRi , avtb(jk)  ) * wmask(ji,jj,jk)
-            END DO
-         END DO
-      ENDIF
-         !
       END DO
+      !
+!!gm BUG <<<<====  This param can't work at low latitude 
+!!gm               it provides there much to thick mixed layer ( summer 150m in GYRE configuration !!! )
       !
       IF( ln_mldw ) THEN      !==  set a minimum value in the Ekman layer  ==!
@@ -184 +153 @@
             END DO
          END DO
-      END IF
+      ENDIF
       !
       IF( nn_timing == 1 )   CALL timing_stop('zdf_ric')