Fixes to interpolation to pressure levels

[djgagne]

I made some significant revisions to the pressure level interpolation module to accommodate the issues @yingkaisha 's PR raised and to make the pressure interpolation process for geopotential and temperature consistent with CESM/ECMWF. New features:

• surface geopotential is now passed into credit.interp.full_state_pressure_interpolation as a numpy array rather than embedded in xarray Dataset with other state variables.
• output.py updated to reflect changed API.
• temperature and geopotential are extrapolated below ground based on Trenberth et al. 1993.
• mean sea level pressure reduction implemented based on Trenberth et al. 1993 and CAM fortran code. The formula adjusts the temperature and lapse rate correction based on the surface height. According to the paper you might still see some wonkiness in the Himalayas/Tibetan plateau, which would require some west-east interpolation over mountains.
• pressure to hybrid coordinate interpolation implemented. Interpolation routine ignores levels below ground but also requires surface pressure field to do so.