Simulations in three spatial dimensions

It would be nice to extend the code to square lattices in three spatial dimensions--this would add relevance for ALF to the nuclear physics community, for example.

I would be happy to take on this extension, after discussions on how general we would want to make it.