Code from publication  http://dx.doi.org/10.1086/589889