var/da/da_radar/da_radial_velocity_lin.inc

   1 subroutine da_radial_velocity_lin(rv,p,u,v,w,qrn,ps,x,y,z,qrn9,rho)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: Tangent linear of da_radial_velocity
   5    ! History:
   6    !    08/2017 - bug fix for Vt (Siou-Ying Jiang, CWB, Taiwan)
   7    !-----------------------------------------------------------------------
   8
   9    implicit none
  10
  11    real, intent(in)  :: x, y, z
  12    real, intent(in)  :: p, u, v, w, qrn, ps
  13    real, intent(in)  :: qrn9
  14    real, intent(in)  :: rho
  15    real, intent(out) :: rv
  16
  17    real    :: r, alpha, vt
  18    real    :: qrrc
  19    real    :: qrn_g, qrn9_g
  20
  21    qrn_g = qrn *1000. ! kg/kg -> g/kg
  22    qrn9_g= qrn9*1000. ! kg/kg -> g/kg
  23    qrrc = 0.01        ! g/kg
  24    vt = 0.0
  25
  26    if (trace_use) call da_trace_entry("da_radial_velocity_lin")
  27
  28    r     = sqrt(x*x+y*y+z*z)
  29    alpha = (ps/p)**0.4
  30
  31    if (qrn9_g <= qrrc)then
  32       vt=0.0
  33    else
  34       vt=0.675*alpha*qrn9_g**(-0.875)*qrn_g*rho**0.125
  35    end if
  36
  37    rv = u*x+v*y+(w-vt)*z
  38    rv = rv/r
  39
  40    if (trace_use) call da_trace_exit("da_radial_velocity_lin")
  41
  42 end subroutine da_radial_velocity_lin
  43
  44