Changeset 1695 for trunk/NEMO/OPA_SRC/ZDF/zdfkpp.F90

Timestamp:

2009-10-30T17:42:38+01:00 (14 years ago)

Author:

smasson

Message:

wind stress module directly at T-point, see ticket:577

File:

• trunk/NEMO/OPA_SRC/ZDF/zdfkpp.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/ZDF/zdfkpp.F90

@@ -200 +200 @@
                       !
       REAL(wp) ::                       &  !
-         ztx, zty, ztau, zflageos,      &  !
+         ztx, zty, zflageos,            &  !
          zstabl, zbuofdep,zucube,       &  !
          zrhos, zalbet, zbeta,          &  !
@@ -469 +469 @@
             zrhos         = rhop(ji,jj,1) + zflageos * rau0 * ( 1. - tmask(ji,jj,1) )  
             ! Friction velocity (zustar), at T-point : LMD94 eq. 2
-            ztx           = 0.5 * ( utau(ji,jj) + utau(ji - 1, jj    ) )
-            zty           = 0.5 * ( vtau(ji,jj) + vtau(ji    , jj - 1) )
-            ztau          = SQRT( ztx * ztx + zty * zty )
-            zustar(ji,jj) = SQRT( ztau / ( zrhos +  epsln ) )
+            zustar(ji,jj) = SQRT( taum(ji,jj) / ( zrhos +  epsln ) )
          ENDDO
       ENDDO