var/da/da_radar/da_radzicevar_rain_tl.inc

   1   subroutine da_radzicevar_rain_tl(zrhv,para_rain,dqra,dqsn,dqgr,dtmk,dqvp,rhoair,  &
   2                           prain_coef,zmm_ref,pdfrrs,pdfsrs,pdfgrg,         &
   3                           pdfrhot,pdfrhoq,pdfrrg,qra,prain)
   4
   5    !----------------------------------------------------------------------
   6    ! purpose: calculate the tl linear for rainwater
   7    !----------------------------------------------------------------------
   8     implicit none
   9     real :: zrhv
  10     real :: para_rain
  11     real :: dqra,dqsn,dqgr,dtmk,dqvp
  12     real :: rhoair,prain_coef,zmm_ref
  13     real :: pdfrrs,pdfsrs,pdfgrg,pdfrhot,pdfrhoq,pdfrrg
  14     real :: qra,prain
  15
  16     zrhv=para_rain*rhoair*prain_coef*dqra         &
  17         +para_rain*prain*pdfrhot*dtmk             &
  18         +para_rain*prain*pdfrhoq*dqvp             &
  19         -para_rain*(pdfrrs+pdfrrg)*qra*dqra       &
  20         -para_rain*pdfsrs*qra*dqsn                &
  21         -para_rain*pdfgrg*qra*dqgr
  22
  23   end subroutine da_radzicevar_rain_tl