var/da/da_radar/da_radzicevar_wetice_tl.inc

   1   subroutine da_radzicevar_wetice_tl(zshv,para_wetice_dqwx,para_wetice_dfwx,pxkhv,pxkhv_tlr,pxkhv_tlx, &
   2                              dqsn,dqra,dtmk,dqvp,qsn,wsnow,rhoair,dsnow_coef,zmm_ref,pdfsrs, &
   3                              pdfrrs,pdfrhot,pdfrhoq)
   4    !----------------------------------------------------------------------
   5    ! purpose: calculate the tl linear for wet snow
   6    !----------------------------------------------------------------------
   7     implicit none
   8     real :: zshv
   9     real :: para_wetice_dqwx,para_wetice_dfwx
  10     real :: pxkhv,pxkhv_tlr,pxkhv_tlx
  11     real :: dqsn,dqra,dtmk,dqvp,qsn,wsnow
  12     real :: rhoair,dsnow_coef,zmm_ref
  13     real :: pdfsrs,pdfrrs,pdfrhot,pdfrhoq
  14
  15     zshv= para_wetice_dqwx*pxkhv*rhoair*(1-dsnow_coef)*dqsn   &
  16          +para_wetice_dqwx*pxkhv*rhoair*(1-dsnow_coef)*dqra   &
  17          +para_wetice_dqwx*pxkhv*wsnow*pdfrhot*dtmk           &
  18          +para_wetice_dqwx*pxkhv*wsnow*pdfrhoq*dqvp           &
  19          +para_wetice_dqwx*pxkhv*pdfrrs*qsn*dqra              &
  20          +para_wetice_dqwx*pxkhv*pdfsrs*qsn*dqsn              &
  21          +para_wetice_dfwx*pxkhv_tlr*dqra                     &
  22          +para_wetice_dfwx*pxkhv_tlx*dqsn
  23
  24   end subroutine da_radzicevar_wetice_tl