lupton_rgb has magic values for initial colour intensities

[mwalmsley]

These were chosen to set the relative intensities of each band to be approximately equal.
I should find a more rigorous way to do this - perhaps an average over the whole set of galaxies

grzscales = dict(g=(2, 0.00526), r=(1, 0.008), z=(0, 0.0135) )