var/da/da_radar/da_radzicevar_waterfraction.inc

   1   subroutine da_radzicevar_waterfraction(qr,qice,waterfraction)
   2   !----------------------------------------------------------------------
   3   ! purpose: calculate the water fraction fw
   4   !----------------------------------------------------------------------
   5   implicit none
   6   real :: qr,qice
   7   real :: waterfraction
   8
   9   if(qr<1.e-8.and.qice<1.e-8) then
  10     waterfraction=1.e-8
  11   else
  12     waterfraction=qr/(qr+qice)
  13   endif
  14
  15   end subroutine da_radzicevar_waterfraction