var/da/da_radar/da_radzicevar_calc_ice_abc.inc

   1   subroutine da_radzicevar_calc_ice_abc(phimean,sigma,ice_abc)
   2   !----------------------------------------------------------------------
   3   ! purpose: calculate the A B C in ice-phase equations
   4   !----------------------------------------------------------------------
   5   implicit none
   6   real :: phimean,sigma
   7   real :: ice_abc(3)  ! 1 for A, 2 for B, 3 for C
   8   real,parameter :: pi=3.1415926
   9   real :: a2a=pi/180.
  10
  11   ice_abc(1)=1./8.*(3.+4.*cos(2.*phimean*pi)*exp(-2.*sigma**2)  &
  12               +cos(4.*phimean)*exp(-8.*sigma**2))
  13   ice_abc(2)=1./8.*(3.-4.*cos(2.*phimean*pi)*exp(-2.*sigma**2) &
  14               +cos(4.*phimean)*exp(-8.*sigma**2))
  15   ice_abc(3)=1./8.*(1.                                          &
  16               -cos(4.*phimean)*exp(-8.*sigma**2))
  17
  18   end subroutine da_radzicevar_calc_ice_abc