Changeset 4427 for branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DYN/dynzad.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/dynzad.F90 (modified) (3 diffs)

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

@@ -99 +99 @@
             zwuw(ji,jj, 1 ) = 0.e0     ! Surface values set to zero
             zwvw(ji,jj, 1 ) = 0.e0
-            DO jk = 2, jpkm1
+            DO jk = 2, mbkmax(ji,jj)-1
                zwuw(ji,jj,jk) =   ( zww(ji+1,jj  )*wn(ji+1,jj  ,jk) + zww(ji,jj)*wn(ji,jj,jk) )   &
                   &             * ( un(ji,jj,jk-1)-un(ji,jj,jk) ) 
@@ -105 +105 @@
                   &             * ( vn(ji,jj,jk-1)-vn(ji,jj,jk) )
             END DO  
-            zwuw(ji,jj,jpk) = 0.e0     ! Bottom values set to zero
-            zwvw(ji,jj,jpk) = 0.e0
+            zwuw(ji,jj,mbkmax(ji,jj)) = 0.e0     ! Bottom values set to zero
+            zwvw(ji,jj,mbkmax(ji,jj)) = 0.e0
          END DO   
       END DO
@@ -136 +136 @@
       DO jj = 2, jpjm1              ! Vertical momentum advection at u- and v-points
          DO ji = 2, jpim1
-            DO jk = 1, jpkm1
+            DO jk = 1, mbkmax(ji,jj)-1
 #else
       DO jk = 1, jpkm1              ! Vertical momentum advection at u- and v-points