Changeset 4427 for branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DYN/dynkeg.F90

Timestamp:

2014-02-04T13:14:00+01:00 (10 years ago)

Author:

trackstand2

Message:

First files changed on last FINISS work package. Stephen's work although
commited by Andy P.

File:

• branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DYN/dynkeg.F90 (modified) (2 diffs)

branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DYN/dynkeg.F90

@@ -87 +87 @@
       DO jj = 2, jpj            ! Horizontal kinetic energy at T-point
          DO ji = 2, jpi
-            DO jk = 1, jpkm1
+            DO jk = 1, mbkmax(ji,jj)-1
                zhke(ji,jj,jk) = 0.25 * (   un(ji-1,jj  ,jk) * un(ji-1,jj  ,jk)   &
                   &                      + un(ji  ,jj  ,jk) * un(ji  ,jj  ,jk)   &
-                                         + vn(ji  ,jj-1,jk) * vn(ji  ,jj-1,jk)   &
+                  &                      + vn(ji  ,jj-1,jk) * vn(ji  ,jj-1,jk)   &
                   &                      + vn(ji  ,jj  ,jk) * vn(ji  ,jj  ,jk)   )
             END DO  
@@ -97 +97 @@
       DO jj = 2, jpjm1          ! add the gradient of kinetic energy to the general momentum trends
          DO ji = 2, jpim1
-            DO jk = 1, jpkm1
+            DO jk = 1, mbkmax(ji,jj)-1
                ua(ji,jj,jk) = ua(ji,jj,jk) - ( zhke(ji+1,jj  ,jk) - zhke(ji,jj,jk) ) / e1u(ji,jj)
                va(ji,jj,jk) = va(ji,jj,jk) - ( zhke(ji  ,jj+1,jk) - zhke(ji,jj,jk) ) / e2v(ji,jj)