Seem to have interpolation and plotting working. Now need to do the whole thing for all Zs we are interested in.