var/da/da_radar/da_radzicevar_upper_f.inc

   1   subroutine da_radzicevar_upper_f(upper_f,qr,qice,qthres,flg)
   2   !----------------------------------------------------------------------
   3   ! purpose: calculate the fraction F
   4   !----------------------------------------------------------------------
   5   implicit none
   6   real :: qr, qice,qthres
   7   real :: upper_f
   8   real :: fmax=0.5
   9   integer :: flg
  10
  11   if(flg==1) fmax=0.5 ! for snow
  12   if(flg==2) fmax=0.3 ! for hail/graupel
  13
  14   upper_f=0
  15   if(qr>1.0e-12 .and. qice>1.0e-12) then
  16     upper_f=fmax*min(qice/qr,qr/qice)**0.3
  17   endif
  18
  19   end subroutine da_radzicevar_upper_f