Changeset 14072 for NEMO/trunk/src/OCE/ZDF/zdfsh2.F90

Timestamp:

2020-12-04T08:48:38+01:00 (3 years ago)

Author:

laurent

Message:

Merging branch "2020/dev_r13648_ASINTER-04_laurent_bulk_ice", ticket #2369

File:

• NEMO/trunk/src/OCE/ZDF/zdfsh2.F90 (modified) (4 diffs)

NEMO/trunk/src/OCE/ZDF/zdfsh2.F90

@@ -2 +2 @@
    !!======================================================================
    !!                       ***  MODULE  zdfsh2  ***
-   !! Ocean physics:  shear production term of TKE 
+   !! Ocean physics:  shear production term of TKE
    !!=====================================================================
    !! History :   -   !  2014-10  (A. Barthelemy, G. Madec)  original code
@@ -36 +36 @@
 CONTAINS
 
-   SUBROUTINE zdf_sh2( Kbb, Kmm, p_avm, p_sh2  ) 
+   SUBROUTINE zdf_sh2( Kbb, Kmm, p_avm, p_sh2  )
       !!----------------------------------------------------------------------
       !!                   ***  ROUTINE zdf_sh2  ***
@@ -44 +44 @@
       !! ** Method  : - a stable discretization of this term is linked to the
       !!                time-space discretization of the vertical diffusion
-      !!                of the OGCM. NEMO uses C-grid, a leap-frog environment 
+      !!                of the OGCM. NEMO uses C-grid, a leap-frog environment
       !!                and an implicit computation of vertical mixing term,
       !!                so the shear production at w-point is given by:
-      !!                   sh2 = mi[   mi(avm) * dk[ub]/e3ub * dk[un]/e3un   ] 
-      !!                       + mj[   mj(avm) * dk[vb]/e3vb * dk[vn]/e3vn   ] 
+      !!                   sh2 = mi[   mi(avm) * dk[ub]/e3ub * dk[un]/e3un   ]
+      !!                       + mj[   mj(avm) * dk[vb]/e3vb * dk[vn]/e3vn   ]
       !!                NB: wet-point only horizontal averaging of shear
       !!
@@ -81 +81 @@
                zsh2u(ji,jj) = ( p_avm(ji+1,jj,jk) + p_avm(ji,jj,jk) ) &
                   &         * (   uu(ji,jj,jk-1,Kmm) -   uu(ji,jj,jk,Kmm) ) &
-                  &         * (   uu(ji,jj,jk-1,Kbb) -   uu(ji,jj,jk,Kbb) ) & 
+                  &         * (   uu(ji,jj,jk-1,Kbb) -   uu(ji,jj,jk,Kbb) ) &
                   &         / ( e3uw(ji,jj,jk  ,Kmm) * e3uw(ji,jj,jk,Kbb) ) &
                   &         * wumask(ji,jj,jk)