# Chris Fuller, April 2014
#import mods
from atpy import Table
import numpy as np
from os.path import join as pj
import matplotlib.pyplot as plt
from lmfit import minimize, Parameters, report_fit
from matplotlib.ticker import MaxNLocator
#Inputs
folder = "/Users/chrisfuller/Dropbox/phd/herchel/coma/final_outputs/" # input/output folder
fname = 'coma_supercluster_cal12.fits' #input name
cat = Table(pj(folder,fname))
cat.VELOCITY_1 = cat.VELOCITY_1/1000.0
if False:
cat = cat.where((cat.RADIUS_VIR > 1.0))
types = ['late', 'inter', 'early']
colours = ['b', 'g', 'r']
sigmas = ['SIGMA1', 'SIGMA5', 'SIGMA10']
sig_labs = ['$\log_{10}$($\Sigma_{1}$) (Mpc$^{-2}$)', '$\log_{10}$($\Sigma_{5}$) (Mpc$^{-2}$)', '$\log_{10}$($\Sigma_{10}$) (Mpc$^{-2}$)']
fig, subs = plt.subplots(nrows=3, ncols=1, sharex=True, sharey=True, squeeze=True, figsize = (8,9), facecolor='w',edgecolor='w')
for j in range(3):
sub = subs[j]
