Changeset 5324 for trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynkeg.F90

Timestamp:

2015-06-01T10:54:30+02:00 (9 years ago)

Author:

pabouttier

Message:

Fix a mistake in jpkm1 loop in dynkeg.F90

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynkeg.F90 (modified) (4 diffs)

trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynkeg.F90

@@ -98 +98 @@
       zhke(:,:) = 0._wp
 
-      SELECT CASE ( kscheme )             !== Horizontal kinetic energy at T-point  ==!
-      !
-      CASE ( nkeg_C2 )                          !--  Standard scheme  --!
-         DO jk = 1, jpkm1
+      DO jk = 1, jpkm1
+
+         SELECT CASE ( kscheme )             !== Horizontal kinetic energy at T-point  ==!
+         !
+         CASE ( nkeg_C2 )                          !--  Standard scheme  --!
+         
             DO jj = 2, jpj
                DO ji = fs_2, jpi   ! vector opt.
@@ -111 +113 @@
                END DO
             END DO
-         END DO
          !
-      CASE ( nkeg_HW )                          !--  Hollingsworth scheme  --!
-         DO jk = 1, jpkm1
+         CASE ( nkeg_HW )                          !--  Hollingsworth scheme  --!
+
             DO jj = 2, jpjm1
                DO ji = fs_2, jpim1   ! vector opt.
@@ -129 +130 @@
                END DO
             END DO
-         END DO
-         CALL lbc_lnk( zhke, 'T', 1. )
+            CALL lbc_lnk( zhke, 'T', 1. )
+            !
+         END SELECT
          !
-      END SELECT
-      !
-      DO jk = 1, jpkm1                    !==  grad( KE ) added to the general momentum trends  ==!
          DO jj = 2, jpjm1
             DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -141 +140 @@
             END DO
          END DO
+
       END DO
       !