var/da/da_radar/da_radzicevar_wetice_adj.inc

   1   subroutine da_radzicevar_wetice_adj(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    !----------------------------------------------------------------------
   6    ! purpose: calculate the adj associated with the contribution from wet snow
   7    !----------------------------------------------------------------------
   8     implicit none
   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     dqsn=dqsn+para_wetice_dqwx*pxkhv*rhoair*(1-dsnow_coef)*zmm_ref &
  16              +para_wetice_dqwx*pxkhv*pdfsrs*qsn*zmm_ref
  17     dqra=dqra+para_wetice_dfwx*pxkhv_tlr*zmm_ref &
  18              +para_wetice_dqwx*pxkhv*pdfrrs*qsn*zmm_ref &
  19              +para_wetice_dqwx*pxkhv*rhoair*(1-dsnow_coef)*zmm_ref
  20     dqsn=dqsn+para_wetice_dfwx*pxkhv_tlx*zmm_ref
  21     dtmk=dtmk+para_wetice_dqwx*pxkhv*wsnow*pdfrhot*zmm_ref
  22     dqvp=dqvp+para_wetice_dqwx*pxkhv*wsnow*pdfrhoq*zmm_ref
  23
  24   end subroutine da_radzicevar_wetice_adj