subroutine da_radial_velocity_adj(rv,p,u,v,w,qrn,ps,x,y,z,qrn9) !----------------------------------------------------------------------- ! Purpose: TBD !----------------------------------------------------------------------- implicit none real, intent(in) :: x, y, z real, intent(in) :: p real, intent(in) :: qrn9 real, intent(in) :: ps real, intent(inout) :: rv real, intent(inout) :: u, v, w, qrn real :: r, alpha, vt if (trace_use) call da_trace_entry("da_radial_velocity_adj") r = sqrt(x*x+y*y+z*z) alpha = (ps/p)**0.4 rv = rv/r u = u + rv*x v = v + rv*y w = w + rv*z vt = -rv*z if (qrn9 > 0.0) then qrn = qrn + vt*0.675*alpha*qrn9**(-0.875) end if if (trace_use) call da_trace_exit("da_radial_velocity_adj") end subroutine da_radial_velocity_adj