da_radial_velocity.inc

1 subroutine da_radial_velocity(rv,p,u,v,w,qrn,ps,x,y,z)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: TBD
5    !-----------------------------------------------------------------------
6 
7    implicit none
8 
9    real, intent(in)  :: x, y, z
10    real, intent(in)  :: p, u, v, w, qrn, ps
11    real, intent(out) :: rv
12 
13    real :: r, alpha, vt
14 
15    if (trace_use) call da_trace_entry("da_radial_velocity")
16 
17    r=sqrt(x*x+y*y+z*z)
18    alpha=(ps/p)**0.4
19 
20    if (qrn <= 0.0) vt=0.0
21    if (qrn >  0.0) vt=5.4*alpha*qrn**0.125
22    rv=u*x+v*y+(w-vt)*z
23    rv=rv/r
24 
25    if (trace_use) call da_trace_exit("da_radial_velocity")
26 
27 end subroutine da_radial_velocity
28 
29