Changeset 1643 for trunk/NEMO

Timestamp:

2009-10-07T15:44:29+02:00 (15 years ago)

Author:

cetlod

Message:

Correction in estimation of wind module with masked wind stress

File:

• trunk/NEMO/OFF_SRC/dtadyn.F90 (modified) (2 diffs)

trunk/NEMO/OFF_SRC/dtadyn.F90

@@ -20 +20 @@
    USE ldftra_oce      ! ocean tracer   lateral physics
    USE zdfmxl
-   USE trabbl          ! tracers: bottom boundary layer
+   USE trabbl
    USE eosbn2
    USE zdfddm          ! vertical  physics: double diffusion
@@ -673 +673 @@
       INTEGER  ::   ji, jj                ! dummy indices
       !!---------------------------------------------------------------------
-      zcoef = 0.5 / ( zrhoa * zcdrag )
+      zcoef = 1. / ( zrhoa * zcdrag )
 !CDIR NOVERRCHK
-      DO jj = 2, jpjm1
+      DO jj = 1, jpj
 !CDIR NOVERRCHK
-         DO ji = fs_2, fs_jpim1   ! vect. opt.
-            ztx = ptaux(ji-1,jj  ) + ptaux(ji,jj)
-            zty = ptauy(ji  ,jj-1) + ptauy(ji,jj)
+         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 )
             pwspd(ji,jj) = SQRT ( ztau * zcoef ) * tmask(ji,jj,1)
          END DO
       END DO
-      CALL lbc_lnk( pwspd(:,:) , 'T', 1. )
 
    END SUBROUTINE tau2wnd