da_radial_velocity_lin.inc

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