Add `interp_hybrid_to_height` functionality?

[AndrewILWilliams]

Firstly, I've spent the day marvelling at how well interp_hybrid_to_pressure() works for interpolating CMIP6 model data on hybrid sigma pressure coordinates! So thank you!

Now, my request...

There are a couple of CMIP6 models (primarily the HadGEM model), which use a hybrid height coordinate instead of hybrid pressure. The formula for height is pretty simple z(level)=a(level)+b(level)*orog(lat,lon), so would it be possible to adapt the existing hybrid pressures function to interpolate a hybrid-height coordinate to a defined set of height levels?

Not sure how easy this would be?

Cheers! Andrew