subroutine da_radial_velocity(rv,p,u,v,w,qrn,ps,x,y,z) 1,2

   !-----------------------------------------------------------------------
   ! Purpose: TBD
   !-----------------------------------------------------------------------

   implicit none

   real, intent(in)  :: x, y, z
   real, intent(in)  :: p, u, v, w, qrn, ps
   real, intent(out) :: rv

   real :: r, alpha, vt
   real    :: qrrc

   qrrc = 1.0e-3
   vt = 0.0

   if (trace_use) call da_trace_entry("da_radial_velocity")

   r=sqrt(x*x+y*y+z*z)
   alpha=(ps/p)**0.4

!   if (qrn <= 0.0) vt=0.0
!   if (qrn >  0.0) vt=5.4*alpha*qrn**0.125

   if (use_radar_rf .or. use_radar_rhv)then
      if (qrn <= qrrc)then
         vt=0.0
      else
         vt=5.4*alpha*qrn**0.125
      end if
   end if
   rv=u*x+v*y+(w-vt)*z
   rv=rv/r

   if (trace_use) call da_trace_exit("da_radial_velocity")

end subroutine da_radial_velocity