Changeset 11653 for NEMO/trunk/src/OCE/LDF/ldfdyn.F90

Timestamp:

2019-10-04T14:34:02+02:00 (5 years ago)

Author:

davestorkey

Message:

Fix for bug in Smagorinsky scheme #2318.

File:

• NEMO/trunk/src/OCE/LDF/ldfdyn.F90 (modified) (1 diff)

NEMO/trunk/src/OCE/LDF/ldfdyn.F90

@@ -315 +315 @@
             DO jj = 1, jpj             ! Set local gridscale values
                DO ji = 1, jpi
-                  esqt(ji,jj) = ( e1e2t(ji,jj) / ( e1t(ji,jj) + e2t(ji,jj) ) )**2 
-                  esqf(ji,jj) = ( e1e2f(ji,jj) / ( e1f(ji,jj) + e2f(ji,jj) ) )**2 
+                  esqt(ji,jj) = ( 2._wp * e1e2t(ji,jj) / ( e1t(ji,jj) + e2t(ji,jj) ) )**2 
+                  esqf(ji,jj) = ( 2._wp * e1e2f(ji,jj) / ( e1f(ji,jj) + e2f(ji,jj) ) )**2 
                END DO
             END DO