Changeset 2772

Timestamp:

2011-06-01T08:31:07+02:00 (13 years ago)

Author:

sga

Message:

NEMO trunk, bugfix to ldfslp.F90: correct jk indices

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/LDF/ldfslp.F90 (modified) (1 diff)

trunk/NEMOGCM/NEMO/OPA_SRC/LDF/ldfslp.F90

@@ -691 +691 @@
             !                          !- bound the slopes: abs(zw.)<= 1/100 and zb..<0
             !                                kxz max= ah slope max =< e1 e3 /(pi**2 2 dt)
-            zbu = MIN(  zbu , -100._wp* ABS( zau ) , -7.e+3_wp/fse3u(ji,jj,ik)* ABS( zau )  )
-            zbv = MIN(  zbv , -100._wp* ABS( zav ) , -7.e+3_wp/fse3v(ji,jj,ik)* ABS( zav )  )
+            zbu = MIN(  zbu , -100._wp* ABS( zau ) , -7.e+3_wp/fse3u(ji,jj,iku)* ABS( zau )  )
+            zbv = MIN(  zbv , -100._wp* ABS( zav ) , -7.e+3_wp/fse3v(ji,jj,ikv)* ABS( zav )  )
             !                          !- Slope at u- & v-points (uslpml, vslpml)
-            uslpml(ji,jj) = zau / ( zbu - zeps ) * umask(ji,jj,ik)
-            vslpml(ji,jj) = zav / ( zbv - zeps ) * vmask(ji,jj,ik)
+            uslpml(ji,jj) = zau / ( zbu - zeps ) * umask(ji,jj,iku)
+            vslpml(ji,jj) = zav / ( zbv - zeps ) * vmask(ji,jj,ikv)
             !
             !                    !==   i- & j-slopes at w-points just below the Mixed Layer   ==!