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