Changeset 1025 for trunk/NEMO/OPA_SRC/SBC/sbcana.F90

Timestamp:

2008-05-30T13:26:09+02:00 (16 years ago)

Author:

cetlod

Message:

adding wind speed module variable, see ticket 172

File:

• trunk/NEMO/OPA_SRC/SBC/sbcana.F90 (modified) (7 diffs)

trunk/NEMO/OPA_SRC/SBC/sbcana.F90

@@ -66 +66 @@
       !!----------------------------------------------------------------------
       INTEGER, INTENT(in) ::   kt       ! ocean time step
-      !
+      !!
       INTEGER             ::   ji, jj          ! dummy loop indices
       REAL(wp)            ::   zfacto          ! local scalar
       REAL(wp)            ::   ztx, zty, ztau  ! local scalar
+      !!
+      NAMELIST/namsbc_ana/ nn_tau000, rn_utau0, rn_vtau0, rn_qns0, rn_qsr0, rn_emp0
       !!---------------------------------------------------------------------
-      NAMELIST/namsbc_ana/ nn_tau000, rn_utau0, rn_vtau0, rn_qns0, rn_qsr0, rn_emp0
       !
       IF( kt == nit000 ) THEN
@@ -77 +78 @@
          REWIND ( numnam )                   ! Read Namelist namsbc : surface fluxes
          READ   ( numnam, namsbc_ana )
-
+         !
          IF(lwp) WRITE(numout,*)' '
          IF(lwp) WRITE(numout,*)' sbc_ana : Constant surface fluxes read in namsbc_ana namelist'
@@ -87 +88 @@
          IF(lwp) WRITE(numout,*)'              solar heat flux        rn_qsr0   = ', rn_qsr0  , ' W/m2'
          IF(lwp) WRITE(numout,*)'              net heat flux          rn_emp0   = ', rn_emp0  , ' Kg/m2/s'
-
+         !
          nn_tau000 = MAX( nn_tau000, 1 )   ! must be >= 1
          qns   (:,:) = rn_qns0
@@ -97 +98 @@
    
       ! Increase the surface stress to its nominal value during the first nn_tau000 time-steps
-         
       IF( kt <= nn_tau000 ) THEN
          zfacto = 0.5 * (  1. - COS( rpi * FLOAT( kt ) / FLOAT( nn_tau000 ) )  )
@@ -103 +103 @@
          vtau(:,:) = zfacto * rn_vtau0
       ENDIF
+
+      ! Estimation of wind speed as a function of wind stress ( |tau|=rhoa*Cd*|U|^2 )
+      zfacto = 0.5 / ( rhoa * cdrag )
+!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 * zfacto ) * tmask(ji,jj,1)
+         END DO
+      END DO
+      CALL lbc_lnk( wndm(:,:) , 'T', 1. )
       !
-   
-      ! Estimation of wind speed as a function of wind stress
-!CDIR NOVERRCHK
-      DO jj = 1, jpj
-!CDIR NOVERRCHK
-         DO ji = 1, jpi
-            ztx  = utau(ji,jj) * umask(ji,jj,1)
-            zty  = vtau(ji,jj) * vmask(ji,jj,1)
-            ztau = SQRT( ztx * ztx + zty * zty )
-            wndm(ji,jj) = SQRT ( ztau / ( rhoa * cdrag ) )
-         ENDDO
-      ENDDO
-
    END SUBROUTINE sbc_ana
 
@@ -151 +152 @@
       REAL(wp) ::   zcos_sais1, zcos_sais2, ztrp, zconv, t_star
       REAL(wp) ::   zsumemp, zsurf
-      REAL(wp) ::   ztx, zty
+      REAL(wp) ::   ztx, zty, zfacto
       !!---------------------------------------------------------------------
          
@@ -268 +269 @@
       END DO
 
-      ! Estimation of wind speed as a function of wind stress
+      ! Estimation of wind speed as a function of wind stress ( |tau|=rhoa*Cd*|U|^2 )
+      zfacto = 0.5 / ( rhoa * cdrag )
 !CDIR NOVERRCHK
-      DO jj = 1, jpj
+      DO jj = 2, jpjm1
 !CDIR NOVERRCHK
-         DO ji = 1, jpi
-            ztx  = utau(ji,jj) * umask(ji,jj,1)
-            zty  = vtau(ji,jj) * vmask(ji,jj,1)
+         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 / ( rhoa * cdrag ) )
-         ENDDO
-      ENDDO
+            wndm(ji,jj) = SQRT ( ztau * zfacto ) * tmask(ji,jj,1)
+         END DO
+      END DO
+      CALL lbc_lnk( wndm(:,:) , 'T', 1. )
 
       ! ---------------------------------- !