Changeset 1695 for trunk/NEMO/OPA_SRC/SBC/sbc_oce.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/SBC/sbc_oce.F90 (modified) (3 diffs)

trunk/NEMO/OPA_SRC/SBC/sbc_oce.F90

@@ -11 +11 @@
    IMPLICIT NONE
    PRIVATE
-   
-   PUBLIC sbc_tau2wnd   ! compute wind speed based on the wind stress
    
    !!----------------------------------------------------------------------
@@ -41 +39 @@
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   utau      !: sea surface i-stress (ocean referential)     [N/m2]
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   vtau      !: sea surface j-stress (ocean referential)     [N/m2]
-   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   wndm      !: wind speed module at T-point (=|U10m-Uoce|)  [m/s]
+   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   taum      !: module of sea surface stress (at T-point)    [N/m2] 
+   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   wndm      !: wind speed module at T-point (=|U10m-Uoce|)  [m/s] Used only in PISCES
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   qsr       !: sea heat flux:     solar                     [W/m2]
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   qns       !: sea heat flux: non solar                     [W/m2]
@@ -73 +72 @@
    !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
    !!======================================================================
-CONTAINS
-
-   SUBROUTINE sbc_tau2wnd
-      !!---------------------------------------------------------------------
-      !!                    ***  ROUTINE sbc_tau2wnd  ***
-      !!                   
-      !! ** Purpose : Estimation of wind speed as a function of wind stress   
-      !!
-      !! ** Method  : |tau|=rhoa*Cd*|U|^2
-      !!---------------------------------------------------------------------
-      USE dom_oce         ! ocean space and time domain
-      USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
-      REAL(wp) ::   zrhoa  = 1.22         ! Air density kg/m3
-      REAL(wp) ::   zcdrag = 1.5e-3       ! drag coefficient
-      REAL(wp) ::   ztx, zty, ztau, zcoef ! temporary variables
-      INTEGER  ::   ji, jj                ! dummy indices
-      !! * Substitutions
-#  include "vectopt_loop_substitute.h90"
-      !!---------------------------------------------------------------------
-      zcoef = 0.5 / ( zrhoa * zcdrag ) 
-!CDIR NOVERRCHK
-      DO jj = 2, jpjm1
-!CDIR NOVERRCHK
-         DO ji = fs_2, fs_jpim1   ! vect. opt.
-            ztx = utau(ji-1,jj  ) + utau(ji,jj) 
-            zty = vtau(ji  ,jj-1) + vtau(ji,jj) 
-            ztau = SQRT( ztx * ztx + zty * zty )
-            wndm(ji,jj) = SQRT ( ztau * zcoef ) * tmask(ji,jj,1)
-         END DO
-      END DO
-      CALL lbc_lnk( wndm(:,:) , 'T', 1. )
-
-   END SUBROUTINE sbc_tau2wnd
 
 END MODULE sbc_oce