#! /usr/bin/env python
import numpy as np
import h5py, os.path, sys
from Cosmology import *
import matplotlib
matplotlib.use("Agg")
from matplotlib import pyplot as plt
from distinct_colours import get_distinct, pault_cmap
import matplotlib.gridspec as gridspec
from scipy import ndimage
from stats import *
import mpl_style
plt.style.use(mpl_style.style1)
ptmap = pault_cmap(1)
path = '/gpfs/data/violeta/Galform_Out/v2.7.0/stable/MillGas/'
nvol = 64
sn = '41'
plotdir = '/gpfs/data/violeta/lines/cosmicweb/plots/modelplots/'
models = ['gp18', 'gp18.font', 'gp18.starvation', 'gp17']
inleg = ['This work', '10% stripping ', 'Starvation', 'GP18']
# Initialize GSMF
mmin = 8.5
mmax = 15.
dm = 0.1
mbins = np.arange(mmin, mmax, dm)
mhist = mbins + dm * 0.5
#! /usr/bin/env python
import numpy as np
import h5py, os.path, sys
from Cosmology import *
import matplotlib
matplotlib.use("Agg")
from matplotlib import pyplot as plt
from distinct_colours import get_distinct,pault_cmap
import matplotlib.gridspec as gridspec
from scipy import ndimage
from stats import *
import mpl_style
plt.style.use(mpl_style.style1) ; ptmap=pault_cmap(1)
path = '/gpfs/data/violeta/Galform_Out/v2.7.0/stable/MillGas/'
nvol = 64
sn = '61'
plotdir = '/gpfs/data/violeta/lines/cosmicweb/plots/modelplots/gaea.'
models = ['gp19']
gaea = '/gpfs/data/violeta/gaea/'
inleg = ['This work','GAEA']#,'10% stripping ','Starvation','GP18']
# Initialize GSMF
mmin = 8.5
mmax = 15.
dm = 0.1
mbins = np.arange(mmin,mmax,dm)
mhist = mbins + dm*0.5