var/da/da_ssmi/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