Should ma.log10() be used at the end of the interp.hght() function instead of np.log10()?
return generic_interp_pres(ma.log10(p), prof.logp[::-1], prof.hght[::-1])
I would like to avoid the log10 of zero warning, but don't want to break anything. I tried a pull request, but couldn't figure out how to isolate this change from the rest of my changes.