There was a bug when converting the indices back from fortran indexing to python indexing in the run loop. there was also a bug in readfields for a z depth isoslice in which the dzt calculation using isoslice had the input arguments switched in order such that the whole array ended up masked out.