Changeset 1230 for trunk/NEMO/OPA_SRC/SBC/sbc_oce.F90

Timestamp:

2008-12-01T08:10:47+01:00 (16 years ago)

Author:

smasson

Message:

new coupling inteface: compute wind speed from tau, see ticket:155

File:

• trunk/NEMO/OPA_SRC/SBC/sbc_oce.F90 (modified) (2 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
    
    !!----------------------------------------------------------------------
@@ -67 +69 @@
    !! 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