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