Don't use meshgrid_nd anymore
make median line share the shading of the whole plot
Make io.fits the default over pyfits
rename Jy2Mag to Jy2mag
tweak Mstar value