Changes

I am committing a bug fix to the code that calculates mixed Brunt
Vaisala frequency. The variable ice_supersat_frac was used directly
in calculations that are supposed to take place on a momentum level.
The fix is to interpolated ice_supersat_frac to momentum levels before
use. This commit is BIT_CHANGING.