Script to plot the gas fraction, based on plot_sfr.py