PROGRAM test_mie

  real :: Qe, Qs, Qb, g
  integer, parameter:: NANG=101
  complex, dimension(2*NANG-1) :: S1, S2
  
  complex:: m
  
  real  :: x
  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., m)
  
  print *, 'm = ', m

  
 END