Have fixed up z grid calculations to be more modern using octant code from Rob. Surface case is now explicitly built into run.py. Currently working on cleaning up the code with better headers.