Updated H_SMOOTHING routine so it is more similar to the FARGO version.  That is, smoothing is done by calculating a global average of H/R and then setting the local smoothing to R*H_frac*<H/R>