Posted a version of this somewhere can't remember where, but anyway this is an updated version using Alma as one of the smoothing ma's and optimized the code somewhat.