Implement ewald_box()
Add test_ewald3.f90
Add another test case
Use correct sign for forces