<HTML> <BODY BGCOLOR=#ccccdd LINK=#0000aa VLINK=#0000ff ALINK=#ff0000 ><BASE TARGET="bottom_target"><PRE>
<A NAME='DA_RADIAL_VELOCITY'><A href='../../html_code/radar/da_radial_velocity.inc.html#DA_RADIAL_VELOCITY' TARGET='top_target'><IMG SRC="../../gif/bar_red.gif" border=0></A>
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