Changeset 1000

Timestamp:

2008-05-30T11:22:27+02:00 (16 years ago)

Author:

cetlod

Message:

computation of wind speed module, see ticket:164

Location:

trunk/NEMO/OPA_SRC/SBC

Files:

• sbc_oce.F90 (modified) (1 diff)
• sbcana.F90 (modified) (5 diffs)
• sbcblk_clio.F90 (modified) (1 diff)
• sbcblk_core.F90 (modified) (1 diff)
• sbcflx.F90 (modified) (3 diffs)

trunk/NEMO/OPA_SRC/SBC/sbc_oce.F90

@@ -16 +16 @@
    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      !: 10m wind speed                               [m/s]
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   qns       !: sea heat flux: non solar                     [W/m2]
    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   qsr       !: sea heat flux:     solar                     [W/m2]

trunk/NEMO/OPA_SRC/SBC/sbcana.F90

@@ -35 +35 @@
    REAL(wp) ::   rn_qsr0   = 0.e0   !     solar heat flux
    REAL(wp) ::   rn_emp0   = 0.e0   ! net freshwater flux
-   
+ 
+   REAL(wp) ::   rhoa  = 1.22         ! Air density kg/m3
+   REAL(wp) ::   cdrag    = 1.5e-3       ! drag coefficient
+  
    !! * Substitutions
 #  include "domzgr_substitute.h90"
@@ -64 +67 @@
       INTEGER, INTENT(in) ::   kt       ! ocean time step
       !
-      REAL(wp)            ::   zfacto   ! local scalar
+      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
@@ -98 +103 @@
          vtau(:,:) = zfacto * rn_vtau0
       ENDIF
-      !   
+      !
+   
+      ! 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
 
@@ -133 +151 @@
       REAL(wp) ::   zcos_sais1, zcos_sais2, ztrp, zconv, t_star
       REAL(wp) ::   zsumemp, zsurf
+      REAL(wp) ::   ztx, zty
       !!---------------------------------------------------------------------
          
@@ -248 +267 @@
          END DO
       END DO
+
+      ! 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
 
       ! ---------------------------------- !

trunk/NEMO/OPA_SRC/SBC/sbcblk_clio.F90

@@ -294 +294 @@
       END DO
 
+      !------------------------------------!
+      !   store the wind speed  (wndm )    !
+      !------------------------------------!
+!CDIR COLLAPSE
+      DO jj = 1 , jpj
+         DO ji = 1, jpi
+            wndm(ji,jj) = sf(jp_wndm)%fnow(ji,jj)
+         END DO
+      END DO
+
+
       !------------------------------------------------!
       !   Shortwave radiation for ocean and snow/ice   !

trunk/NEMO/OPA_SRC/SBC/sbcblk_core.F90

@@ -247 +247 @@
       !      0   Wind components and module at T-point relative to the moving ocean   !
       ! ----------------------------------------------------------------------------- !
+
+      ! Store the wind speed U10m
+!CDIR NOVERRCHK
+!CDIR COLLAPSE
+      DO jj = 1, jpj
+         DO ji = 1, jpi   ! vect. opt.
+            wndm(ji,jj) = SQRT(  sf(jp_wndi)%fnow(ji,jj) * sf(jp_wndi)%fnow(ji,jj)   &
+         &                     + sf(jp_wndj)%fnow(ji,jj) * sf(jp_wndj)%fnow(ji,jj) )
+         ENDDO
+      ENDDO
+
       ! ... components ( U10m - U_oce ) at T-point (unmasked)
       zwnd_i(:,:) = 0.e0  

trunk/NEMO/OPA_SRC/SBC/sbcflx.F90

@@ -50 +50 @@
    TYPE(FLD), ALLOCATABLE, DIMENSION(:) ::   sf    ! structure of input fields (file informations, fields read)
 
+   REAL(wp) ::   rhoa  = 1.22         ! Air density kg/m3
+   REAL(wp) ::   cdrag = 1.5e-3       ! drag coefficient
+
    !! * Substitutions
 #  include "domzgr_substitute.h90"
@@ -96 +99 @@
       INTEGER  ::   ierror       ! return error code
       REAL(wp) ::   zfact        ! temporary scalar
+      REAL(wp) ::   ztx, zty, ztau
       !!
       CHARACTER(len=100) ::  cn_dir                               ! Root directory for location of flx files
@@ -179 +183 @@
       END DO
 
+      ! 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
+
       ! control print (if less than 100 time-step asked)
       IF( nitend-nit000 <= 100 .AND. lwp ) THEN