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.