PROGRAM test_mie

  implicit none

  real :: Qe, Qs, Qb, g
  integer, parameter:: NANG=101
  complex, dimension(2*NANG-1) :: S1, S2
  
  complex:: m
  
  real  :: x, mr, mi
  real, parameter  :: wl = 1.0 ! in micron
  real, parameter  :: r = 1.0 ! in micron
  real, parameter  :: pi = 3.14159265

  x = 2 * PI * r / wl
  CALL index_water(0, wl, 273., mr, mi)
  
  m = CMPLX(mr,mi)
  
  CALL bhmie(x, m, nang, s1, s2, qe, qs, qb, g)
  
  print *, 'm = ', m
  print *, 'x = ', x
  print *, 'Qe = ', qe
  print *, 'Qs = ', qs
  print *, 'ssa = ', qs/qe
  print *, 'g = ', g
  
 END