da_transform_xtospeed_adj.inc

1 subroutine da_transform_xtospeed_adj(ADJ_speed,U,V,ADJ_u,ADJ_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 (in  ) :: ADJ_speed
10    real,   intent (in  ) :: U, V
11    real,   intent (inout) :: ADJ_u, ADJ_v
12    real                   :: speed
13 
14    if (trace_use) call da_trace_entry("da_transform_xtospeed_adj")
15 
16    speed = sqrt(U*U+V*V+ 1.E-6)
17 
18    ADJ_u = U*ADJ_speed/speed + ADJ_u
19    ADJ_v = V*ADJ_speed/speed + ADJ_v
20 
21    if (trace_use) call da_trace_exit("da_transform_xtospeed_adj")
22 
23 end subroutine da_transform_xtospeed_adj
24 
25