Changeset 4959 for branches

Timestamp:

2014-12-02T15:57:40+01:00 (9 years ago)

Author:

cetlod

Message:

2014/dev_MERGE_2014 : fix out of bounds in ldfslp

File:

• branches/2014/dev_MERGE_2014/NEMOGCM/NEMO/OPA_SRC/LDF/ldfslp.F90 (modified) (1 diff)

branches/2014/dev_MERGE_2014/NEMOGCM/NEMO/OPA_SRC/LDF/ldfslp.F90

@@ -408 +408 @@
          ! set the slope of diffusion to the slope of s-surfaces 
          !      ( c a u t i o n : minus sign as fsdep has positive value ) 
-         DO jk = 1, jpk 
+         DO jj = 2, jpjm1 
+            DO ji = fs_2, fs_jpim1   ! vector opt. 
+               uslp(ji,jj,1) = -1./e1u(ji,jj) * ( fsdept_b(ji+1,jj,1) - fsdept_b(ji ,jj ,1) )  * umask(ji,jj,1) 
+               vslp(ji,jj,1) = -1./e2v(ji,jj) * ( fsdept_b(ji,jj+1,1) - fsdept_b(ji ,jj ,1) )  * vmask(ji,jj,1) 
+               wslpi(ji,jj,1) = -1./e1t(ji,jj) * ( fsdepw_b(ji+1,jj,1) - fsdepw_b(ji-1,jj,1) ) * tmask(ji,jj,1) * 0.5 
+               wslpj(ji,jj,1) = -1./e2t(ji,jj) * ( fsdepw_b(ji,jj+1,1) - fsdepw_b(ji,jj-1,1) ) * tmask(ji,jj,1) * 0.5 
+            END DO 
+         END DO 
+
+         DO jk = 2, jpk 
             DO jj = 2, jpjm1 
                DO ji = fs_2, fs_jpim1   ! vector opt.