! write the file charges.bin for the test
program aer
  open(unit=125,file='charges.bin',form='unformatted')
  write(125) 3
  write(125) 0.D0,0.D0,0.D0,-0.8D0
  write(125) 0.0000000000D0,  -1.4774185392D0,   1.0830805621D0, 0.4D0
  write(125) 0.0000000000D0,   1.4774185392D0,   1.0830805621D0, 0.4D0
  close(125)
end program aer