Address #180.

"Post-compute" basal melt rate (i.e. use enthalpy at the *end* of the
time-step to compute the upward flux at the base of the column).

Ed: please take a look at this.