Add a script to plot column density