var/da/da_radar/da_radzicevar_sigma_in_abc.inc

   1   subroutine da_radzicevar_sigma_in_abc(qg,fw,gsflag,sigma_in_abc)
   2   !----------------------------------------------------------------------
   3   ! purpose: calculate the standard deviation of canting angle
   4   !----------------------------------------------------------------------
   5   implicit none
   6   real :: qg,fw
   7   integer :: gsflag
   8   real :: sigma_in_abc
   9   real,parameter :: pi=3.1415926
  10   real :: a2a=pi/180.
  11
  12   if(gsflag==1) then ! for snow
  13     sigma_in_abc=20.*a2a
  14   elseif(gsflag==2) then ! for graupel
  15     if(qg>0.2*1e-3) then
  16       sigma_in_abc=60.*a2a*(1-0.8*fw)
  17     else
  18       sigma_in_abc=60.*a2a*(1-4.*qg*1e3*fw)
  19     endif
  20   endif
  21
  22   end subroutine da_radzicevar_sigma_in_abc