Add analogue code to simplify a couple of code chunks

[gavinsimpson]

This addresses issue #3 plus a minor tweak to the Stratiplot() call as the y axis limits get taken care of internally.