Large parts of the package re-written. This version gives exactly the same results for the scenarios that could be run with the original script of P. Burton.