added new script for plotting stage shift data

Also added functions for easy plotting of line densities in sf.py