Changeset 4389 for branches/2013

Timestamp:

2014-01-31T15:55:04+01:00 (10 years ago)

Author:

hliu

Message:

add a missed loop in wadlmt.F90

File:

• branches/2013/dev_r4050_NOC_WaD/NEMOGCM/NEMO/OPA_SRC/DYN/wadlmt.F90 (modified) (2 diffs)

branches/2013/dev_r4050_NOC_WaD/NEMOGCM/NEMO/OPA_SRC/DYN/wadlmt.F90

@@ -52 +52 @@
       REAL(wp) ::   zzflxp, zzflxn      ! local scalars
       REAL(wp) ::   ztmp                ! local scalars
-      INTEGER,  POINTER,  DIMENSION(:,:) ::   iwdflg                   !: W/D cell flag
+      INTEGER,  POINTER,  DIMENSION(:,:) ::   iwdflg                   !: W/D cell flag (should become global)
       REAL(wp), POINTER,  DIMENSION(:,:) ::   zwdlmtu, zwdlmtv         !: W/D flux limiters
       REAL(wp), POINTER,  DIMENSION(:,:) ::   zflxp,  zflxn            ! specific 2D workspace
@@ -158 +158 @@
           
            IF(jk1 == nn_waditr) THEN
-              
-              IF(zflxu1(ji,  jj  ) >= 0._wp) zwdlmtu(ji,  jj  ) = 0._wp
-              IF(zflxu1(ji-1,jj  ) <  0._wp) zwdlmtu(ji-1,jj  ) = 0._wp
-              IF(zflxv1(ji,  jj  ) >= 0._wp) zwdlmtv(ji,  jj  ) = 0._wp
-              IF(zflxv1(ji,  jj-1) <  0._wp) zwdlmtv(ji,  jj-1) = 0._wp
-
-              CALL lbc_lnk( zwdlmtu, 'U', 1. )
-              CALL lbc_lnk( zwdlmtv, 'V', 1. )
+              DO jj = 2, jpjm1
+                 DO ji = fs_2, fs_jpim1   ! vector opt.
+                 
+                 IF(zflxu1(ji,  jj  ) >= 0._wp) zwdlmtu(ji,  jj  ) = 0._wp
+                 IF(zflxu1(ji-1,jj  ) <  0._wp) zwdlmtu(ji-1,jj  ) = 0._wp
+                 IF(zflxv1(ji,  jj  ) >= 0._wp) zwdlmtv(ji,  jj  ) = 0._wp
+                 IF(zflxv1(ji,  jj-1) <  0._wp) zwdlmtv(ji,  jj-1) = 0._wp
+            
+                 CALL lbc_lnk( zwdlmtu, 'U', 1. )
+                 CALL lbc_lnk( zwdlmtv, 'V', 1. )
+                 END DO ! ji loop
+              END DO  ! jj loop
            END IF