Changeset 1699 for trunk

Timestamp:

2009-11-02T11:40:30+01:00 (14 years ago)

Author:

cetlod

Message:

Correction in estimation of wind module with masked wind stress, see changeset:1643

File:

• trunk/NEMO/OFF_SRC/dtadyn.F90 (modified) (1 diff)

trunk/NEMO/OFF_SRC/dtadyn.F90

@@ -675 +675 @@
       zcoef = 1. / ( zrhoa * zcdrag )
 !CDIR NOVERRCHK
-      DO jj = 1, jpj
+      DO jj = 2, jpjm1
 !CDIR NOVERRCHK
-         DO ji = 1, jpi   ! vect. opt.
-            ztx = ptaux(ji,jj) * umask(ji,jj,1)
-            zty = ptauy(ji,jj) * vmask(ji,jj,1)
-            ztau = SQRT( ztx * ztx + zty * zty )
+         DO ji = fs_2, fs_jpim1   ! vector opt.
+            ztx = ptaux(ji,jj) * umask(ji,jj,1) + ptaux(ji-1,jj  ) * umask(ji-1,jj  ,1)
+            zty = ptauy(ji,jj) * vmask(ji,jj,1) + ptauy(ji  ,jj-1) * vmask(ji  ,jj-1,1)
+            ztau = 0.5 * SQRT( ztx * ztx + zty * zty )
             pwspd(ji,jj) = SQRT ( ztau * zcoef ) * tmask(ji,jj,1)
          END DO
       END DO
+      CALL lbc_lnk( pwspd(:,:), 'T', 1. )
 
    END SUBROUTINE tau2wnd