Changeset 1694 for trunk/NEMO/OPA_SRC/DOM

Timestamp:

2009-10-30T17:06:23+01:00 (15 years ago)

Author:

smasson

Message:

avoid out of bounds access, see ticket:576

Location:

trunk/NEMO/OPA_SRC/DOM

Files:

• dommsk.F90 (modified) (1 diff)
• domvvl.F90 (modified) (1 diff)
• domzgr.F90 (modified) (1 diff)

trunk/NEMO/OPA_SRC/DOM/dommsk.F90

@@ -213 +213 @@
                vmask(ji,jj,jk) = tmask(ji,jj  ,jk) * tmask(ji  ,jj+1,jk)
             END DO
-            DO ji = 1, jpim1 
+            DO ji = 1, jpim1      ! NO vector opt.
                fmask(ji,jj,jk) = tmask(ji,jj  ,jk) * tmask(ji+1,jj  ,jk)   &
                   &            * tmask(ji,jj+1,jk) * tmask(ji+1,jj+1,jk)

trunk/NEMO/OPA_SRC/DOM/domvvl.F90

@@ -117 +117 @@
 
       DO jj = 1, jpjm1                          ! initialise before and now Sea Surface Height at u-, v-, f-points
-         DO ji = 1, jpim1
+         DO ji = 1, jpim1   ! NO vector opt.
             zcoefu = 0.5  * umask(ji,jj,1) / ( e1u(ji,jj) * e2u(ji,jj) )
             zcoefv = 0.5  * vmask(ji,jj,1) / ( e1v(ji,jj) * e2v(ji,jj) )

trunk/NEMO/OPA_SRC/DOM/domzgr.F90

@@ -1253 +1253 @@
       hbatf(:,:) = rn_sbot_min
       DO jj = 1, jpjm1
-        DO ji = 1, jpim1
+        DO ji = 1, jpim1   ! NO vector opt.
            hbatu(ji,jj) = 0.5 * ( hbatt(ji  ,jj) + hbatt(ji+1,jj  ) )
            hbatv(ji,jj) = 0.5 * ( hbatt(ji  ,jj) + hbatt(ji  ,jj+1) )