New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 1000 for trunk/NEMO/OPA_SRC/SBC/sbcflx.F90 – NEMO

Ignore:
Timestamp:
2008-05-30T11:22:27+02:00 (16 years ago)
Author:
cetlod
Message:

computation of wind speed module, see ticket:164

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/SBC/sbcflx.F90

    r888 r1000  
    5050   TYPE(FLD), ALLOCATABLE, DIMENSION(:) ::   sf    ! structure of input fields (file informations, fields read) 
    5151 
     52   REAL(wp) ::   rhoa  = 1.22         ! Air density kg/m3 
     53   REAL(wp) ::   cdrag = 1.5e-3       ! drag coefficient 
     54 
    5255   !! * Substitutions 
    5356#  include "domzgr_substitute.h90" 
     
    9699      INTEGER  ::   ierror       ! return error code 
    97100      REAL(wp) ::   zfact        ! temporary scalar 
     101      REAL(wp) ::   ztx, zty, ztau 
    98102      !! 
    99103      CHARACTER(len=100) ::  cn_dir                               ! Root directory for location of flx files 
     
    179183      END DO 
    180184 
     185      ! Estimation of wind speed as a function of wind stress 
     186!CDIR NOVERRCHK 
     187      DO jj = 1, jpj 
     188!CDIR NOVERRCHK 
     189         DO ji = 1, jpi 
     190            ztx  = utau(ji,jj) * umask(ji,jj,1) 
     191            zty  = vtau(ji,jj) * vmask(ji,jj,1) 
     192            ztau = SQRT( ztx * ztx + zty * zty ) 
     193            wndm(ji,jj) = SQRT ( ztau / ( rhoa * cdrag ) ) 
     194         ENDDO 
     195      ENDDO 
     196 
    181197      ! control print (if less than 100 time-step asked) 
    182198      IF( nitend-nit000 <= 100 .AND. lwp ) THEN 
Note: See TracChangeset for help on using the changeset viewer.