Split interp_zeta_level.py into two files. One does the interpolation to a set height and saves the result. The other plots the data at the saved height. A better more general height interpolation scheme would be good, maybe ones is built in?