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

trunk/NEMO/OPA_SRC/SBC/sbccpl.F90

@@ -545 +545 @@
       !!
       !! ** Action  :   update  utau, vtau   ocean stress at U,V grid 
+      !!                        taum, wndm   wind stres and wind speed module at T-point
       !!                        qns , qsr    non solar and solar ocean heat fluxes   ('ocean only case)
       !!                        emp = emps   evap. - precip. (- runoffs) (- calving) ('ocean only case)
@@ -557 +558 @@
       REAL(wp) ::   zcumulneg, zcumulpos   ! temporary scalars     
       REAL(wp) ::   zcoef                  ! temporary scalar
+      REAL(wp) ::   zrhoa  = 1.22          ! Air density kg/m3
+      REAL(wp) ::   zcdrag = 1.5e-3        ! drag coefficient
+      REAL(wp) ::   zzx, zzy               ! temporary variables
       REAL(wp), DIMENSION(jpi,jpj) ::   ztx, zty   ! 2D workspace 
       !!----------------------------------------------------------------------
@@ -625 +629 @@
       ! u(v)tau will be modified by ice model -> need to be reset before each call of the ice/fsbc      
       IF( MOD( kt-1, k_fsbc ) == 0 ) THEN
+         !
          utau(:,:) = frcv(:,:,jpr_otx1)                   
          vtau(:,:) = frcv(:,:,jpr_oty1)
-         IF( .NOT. srcv(jpr_w10m)%laction )   CALL sbc_tau2wnd
+         !
+         ! wind stress module at T-point
+!CDIR NOVERRCHK
+         DO jj = 2, jpjm1
+!CDIR NOVERRCHK
+            DO ji = fs_2, fs_jpim1   ! vect. opt.
+               zzx = utau(ji-1,jj  ) + utau(ji,jj) 
+               zzy = vtau(ji  ,jj-1) + vtau(ji,jj) 
+               taum(ji,jj) = 0.5 * SQRT( zzx * zzx + zzy * zzy )
+            END DO
+         END DO
+         CALL lbc_lnk( taum(:,:), 'T', 1. )
+         !
+         ! wind speed module at T-point
+         IF( .NOT. srcv(jpr_w10m)%laction ) THEN 
+            zcoef = 1. / ( zrhoa * zcdrag ) 
+            !CDIR NOVERRCHK
+            DO jj = 1, jpj
+               !CDIR NOVERRCHK
+               DO ji = 1, jpi 
+                  wndm(ji,jj) = SQRT( taum(ji,jj) * zcoef )
+               END DO
+            END DO
+         ENDIF
+         !  
       ENDIF
       !                                                      ! ========================= !
@@ -678 +707 @@
          !                                                           ! 10 m wind speed
          IF( srcv(jpr_w10m)%laction )   wndm(:,:) = frcv(:,:,jpr_w10m)
-         ! it not, we call sbc_tau2wnd in sbc_cpl_rcv (or later, after the ice???)
          !
 #if defined  key_cpl_carbon_cycle