da_radial_velocity_adj.inc

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