did a little cleaning up of step.f95 and run.py but not much. More will probably happen over time in the process of doing other things.
altered formulation for gradients of slope in the anisodiffusion flag in diffusion.f95. However, I am not sure it is correct quite yet. I want to carefully evaluate the indices used in both h and the grid spacing.