SRJ works!  need to split this up into two notebooks, possibly make SRJ a numpy routine