mgrid = np.log10(M)
zgrid = old_div((z_edges[1:]+z_edges[:-1]),2.)
zz = np.arange(0.1,2.01,0.05)
MMexp = np.arange(13.5,15.71,0.1)
MM = 10**MMexp
hscgrid = np.loadtxt("data/HSC_DeltalnM_z0_z2_17_04_04.txt")
sngrid = old_div(1.,hscgrid)
pgrid = np.rot90(sngrid)
pl = Plotter(xlabel="$\\mathrm{log}_{10}(M)$",ylabel="$z$",ftsize=14)
pl.plot2d(pgrid,extent=[MMexp.min(),MMexp.max(),zz.min(),zz.max()],levels=[3.0,5.0],labsize=14,aspect="auto")
pl.done(outDir+"origHSCgrid.png")
print(hscgrid.shape)
#outmerr = interpolateGrid(hscgrid,MM,zz,M,zgrid,regular=False,kind="cubic",bounds_error=False,fill_value=np.inf)
outmerr = interpolateGrid(hscgrid,MM,zz,M,zgrid,regular=False,kind="cubic",bounds_error=False)
sngrid = old_div(1.,outmerr)
pgrid = np.rot90(sngrid)
pl = Plotter(xlabel="$\\mathrm{log}_{10}(M)$",ylabel="$z$",ftsize=14)
pl.plot2d(pgrid,extent=[mgrid.min(),mgrid.max(),zgrid.min(),zgrid.max()],levels=[3.0,5.0],labsize=14,aspect="auto")
pl.done(outDir+"interpHSCgrid.png")
pl._ax.set_ylim(1, 5.e4) # fsky
pl._ax.set_xlim(0., 3.)
pl.done(outDir + "clNofz.pdf")
fsense[fsense > 10.] = np.nan
from orphics.io import Plotter
import os
mmin = mgrid.min()
mmax = mgrid.max()
zmin = zgrid.min()
zmax = zgrid.max()
pgrid = np.rot90((fsense))
pl = Plotter(xlabel="$\\mathrm{log}_{10}(M)$", ylabel="$z$", ftsize=14)
pl.plot2d(pgrid,
extent=[mmin, mmax, zmin, zmax],
labsize=14,
aspect="auto",
lim=[0., 10.])
pl.done(outDir + "massSense.png")
rn[rn < 1] = np.nan
mmin = mgrid.min()
mmax = mgrid.max()
zmin = zgrid.min()
zmax = zgrid.max()
pgrid = np.rot90((rn))
pl = Plotter(labelX="$\\mathrm{log}_{10}(M)$", labelY="$z$", ftsize=14)
pl.plot2d(pgrid, extent=[mmin, mmax, zmin, zmax], labsize=14, aspect="auto")
pl.done(outDir + "rn.png")
#rn[rn<1]=np.nan
jointgrid = np.sqrt(old_div(1.,jointgridsqinv))
snjoint = old_div(1.,jointgrid)
for ind in mindicesList:
pl.add(outzgrid,snjoint[ind,:].ravel(),ls="-.")
print(mexpgrid[ind])
pl.legendOn(loc='upper right',labsize=10)
pl.done(outDir+"slice"+cmbtype+".pdf")
#pl.legendOn(loc='upper right',labsize=8)
#pl.done(outDir+"slice.pdf")
#from orphics.io import Plotter
pgrid = np.rot90(old_div(1.,jointgrid))
pl = Plotter(labelX="$\\mathrm{log}_{10}(M)$",labelY="$z$",ftsize=14)
pl.plot2d(pgrid,extent=[mmin,mmax,zmin,zmax],levels=[1.0,3.0,5.0],labsize=14)
pl.done(outDir+"joint"+cmbtype+".png")
#savefile = "/astro/astronfs01/workarea/msyriac/data/SZruns/v0.6/lensgrid_S4-1.0-0.4_grid-default_CMB_all_joint_v0.5.pkl"
#savefile = "/astro/astronfs01/workarea/msyriac/data/SZruns/v0.6/lensgridRayUp_S4-1.0-0.4_grid-default_CMB_all_joint_v0.5.pkl"
#savefile = "/astro/astronfs01/workarea/msyriac/data/SZruns/v0.6/lensgridRayDn_S4-1.0-0._grid-default_CMB_all_joint_v0.5.pkl"
#pickle.dump((outmedges,outzedges,jointgrid),open(savefile,'wb'))
save_file = save_func(bigDataDir,expName,cmbGrid,saveCalName,version)
pickle.dump((outmedges,outzedges,jointgrid),open(save_file,'wb'))
