put in specific heat estimator into analysis scripts.  For this, the main driver now writes T, L to disk (and H, for use in future) along with E^2.  All analysis scripts have been changed accordingly, but now the spline fitting functionality of the plot_range is broken for some reason...