da_radial_velocity_adj.inc

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