da_transform_xtospeed_lin.inc

1 subroutine da_transform_xtospeed_lin(TGL_speed,U,V,TGL_u,TGL_v)
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) :: TGL_speed
10    real, intent(in)  :: U,V
11    real, intent(in)  :: TGL_u,TGL_v
12 
13    real :: speed
14 
15    if (trace_use) call da_trace_entry("da_transform_xtospeed_lin")
16 
17    speed = sqrt(U*U + V*V + 1.0e-6)
18 
19    TGL_speed = (U*TGL_u + V*TGL_v)/speed
20 
21    if (trace_use) call da_trace_exit("da_transform_xtospeed_lin")
22 
23 end subroutine da_transform_xtospeed_lin
24 
25