Changeset 1775 for trunk/NEMO

Timestamp:

2009-11-26T12:02:18+01:00 (14 years ago)

Author:

rblod

Message:

Correct arrays out of bounds in diaptr, see ticket #619

File:

• trunk/NEMO/OPA_SRC/DIA/diaptr.F90 (modified) (3 diffs)

trunk/NEMO/OPA_SRC/DIA/diaptr.F90

@@ -27 +27 @@
    USE in_out_manager
    USE lib_mpp
+   USE lbclnk
 
    IMPLICIT NONE
@@ -398 +399 @@
             ! T times V on T points (include bolus velocities)
 #if defined key_diaeiv 
-            DO jj = 1, jpj
+            DO jj = 2, jpj
                DO ji = 1, jpi
                   vt(ji,jj,:) = tn(ji,jj,:) * ( vn(ji,jj,:) + vn(ji,jj-1,:) + u_eiv(ji,jj,:) + u_eiv(ji,jj-1,:) )*0.5
@@ -405 +406 @@
             END DO
 #else
-            DO jj = 1, jpj
+            DO jj = 2, jpj
                DO ji = 1, jpi
                   vt(ji,jj,:) = tn(ji,jj,:) * ( vn(ji,jj,:) + vn(ji,jj-1,:) )*0.5
-                  vs(ji,jj,:) = sn(ji,jj,:) * ( vn(ji,jj,:) + vn(ji,jj-1,:)  )*0.5
+                  vs(ji,jj,:) = sn(ji,jj,:) * ( vn(ji,jj,:) + vn(ji,jj-1,:) )*0.5
                END DO
             END DO
 #endif 
+            CALL lbc_lnk( vs, 'V', -1. )   ;   CALL lbc_lnk( vt, 'V', -1. )
 
             ht_glo(:) = SUM( ptr_vjk( vt(:,:,:)), 2 )