I do not remember what this branch was originally for, but I did some changes in readfields to improve the speed. The interpolation, flipping of arrays, and depth calculations should all be improved for speed if possible.
switched so that choices in run (3d/2d and turbulence/diffusion type) are chosen with flags in run.py instead of through compiler flags. Also worked on a problem in the anisodiffusion scheme.