da_radial_velocity_lin.inc

1 subroutine da_radial_velocity_lin(rv,p,u,v,w,qrn,ps,x,y,z,qrn9)
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(in)  :: qrn9
12    real, intent(out) :: rv
13 
14    real    :: r, alpha, vt
15 
16    if (trace_use) call da_trace_entry("da_radial_velocity_lin")
17 
18    r     = sqrt(x*x+y*y+z*z)
19    alpha = (ps/p)**0.4
20 
21    if (qrn9 <= 0.0) then
22       vt=0.0
23    end if
24 
25    if (qrn9 >  0.0) then
26       vt=0.675*alpha*qrn9**(-0.875)*qrn
27    end if
28 
29    rv = u*x+v*y+(w-vt)*z
30    rv = rv/r
31 
32    if (trace_use) call da_trace_exit("da_radial_velocity_lin")
33 
34 end subroutine da_radial_velocity_lin
35 
36