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 
13    real :: speed
14 
15    if (trace_use) call da_trace_entry("da_transform_xtospeed_adj")
16 
17    speed = sqrt(U*U+V*V+ 1.0e-6)
18 
19    ADJ_u = U*ADJ_speed/speed + ADJ_u
20    ADJ_v = V*ADJ_speed/speed + ADJ_v
21 
22    if (trace_use) call da_trace_exit("da_transform_xtospeed_adj")
23 
24 end subroutine da_transform_xtospeed_adj
25 
26