da_transform_xtoseasfcwind.inc

1 subroutine da_transform_xtoseasfcwind(U,V,Speed,zhmkz)
2 
3    !----------------------------------------------------------------------------
4    ! Purpose: Convert (U-V in m/s) components into wind speed (Speed in m/s)
5    !----------------------------------------------------------------------------
6 
7    implicit none
8 
9    real, intent (out) :: Speed
10    real, intent (in)  :: U, V, zhmkz
11 
12    real :: usfc, vsfc
13 
14    if (trace_use) call da_trace_entry("da_transform_xtoseasfcwind")
15 
16    usfc   = U*log(10./0.0001)/log(zhmkz/0.0001) ! roughness = 0.0001
17    vsfc   = V*log(10./0.0001)/log(zhmkz/0.0001) ! roughness = 0.0001
18    speed  = sqrt(usfc*usfc + vsfc*vsfc)
19 
20    if (trace_use) call da_trace_exit("da_transform_xtoseasfcwind")
21 
22 end subroutine da_transform_xtoseasfcwind
23 
24