JAYMES- added 22 Apr 2015 This function calculates saturation vapor pressure. Separate ice and liquid functions are used (identical to those in module_mp_thompson.F, v3.6). Then, the final returned value is a temperature-dependant "blend". Because the final value is "phase-aware", this formulation may be preferred for use throughout the module (replacing "svp").