var/da/da_radar/da_radzicevar_rhoair_tl.inc

   1   subroutine da_radzicevar_rhoair_tl(tlopt,rhoair,prs,rgas,tmk,qvp,pdfrhot,pdfrhoq)
   2    !----------------------------------------------------------------
   3    ! purpose: calculate the air density and its tl linear
   4    !----------------------------------------------------------------
   5    implicit none
   6    integer   :: tlopt
   7    real      :: rhoair,prs,rgas,tmk,qvp
   8    real      :: virtual
   9    real      :: pdfrhot,pdfrhoq
  10
  11     call da_radzicevar_virtual(virtual,tmk,qvp)
  12     rhoair=prs*1./(rgas*virtual)      ! air density
  13     if(tlopt>=1) then
  14       pdfrhot=-prs*1./(rgas*virtual)*(1./tmk)
  15       pdfrhoq=-prs*1./(rgas*virtual)*(0.61/(1+0.61*qvp))
  16     endif
  17     pdfrhot=0 ! not used in current version
  18     pdfrhoq=0 ! not used in current version
  19
  20   end subroutine da_radzicevar_rhoair_tl