clusterDict = dict_from_section(Config,clusterParams)
cc = ClusterCosmology(cosmoDict,constDict,pickling=True,clTTFixFile = "data/cltt_lensed_Feb18.txt")
# make an SZ profile example
SZProfExample = SZ_Cluster_Model(clusterCosmology=cc,clusterDict=clusterDict,rms_noises = noise,fwhms=beam,freqs=freq,lmax=lmax,lknee=lknee,alpha=alpha,dell=dell,pmaxN=pmaxN,numps=numps,v3mode=2,fsky=0.4)
#MM = 10**np.linspace(13.,14.,5)
#print SZProfExample.quickVar(MM,zz,tmaxN=tmaxN,numts=numts)
#sys.exit()
print(("quickvar " , np.sqrt(SZProfExample.quickVar(MM,zz,tmaxN=tmaxN,numts=numts))))
#print "filtvar " , np.sqrt(SZProfExample.filter_variance(MM,zz))
print(("y_m",SZProfExample.Y_M(MM,zz)))
R500 = cc.rdel_c(MM,zz,500.).flatten()[0]
print(R500)
print((cc.rhoc(0)))
print((cc.rhoc(zz)))
#sys.exit()
ells=ls,
Nls=Nls,
kellmax=kmax,
overdensity=overdensity,
critical=critical,
atClusterZ=atClusterZ,
saveId=None,
rayleighSigmaArcmin=rayDn,
arcStamp=arcStamp,
pxStamp=pxStamp)
else:
snRetDn = snRet
MerrGridDn[mindex, zindex] = old_div(1., snRetDn)
if doSZ:
var = SZCluster.quickVar(10**mass, z)
siggrid[mindex, zindex] = np.sqrt(var)
if rank != 0:
if doLens:
MerrGrid = MerrGrid.astype(np.float64)
comm.Send(MerrGrid, dest=0, tag=77)
if True: #doRayDeriv:
MerrGridUp = MerrGridUp.astype(np.float64)
MerrGridDn = MerrGridDn.astype(np.float64)
comm.Send(MerrGridUp, dest=0, tag=98)
comm.Send(MerrGridDn, dest=0, tag=99)
if doSZ:
siggrid = siggrid.astype(np.float64)
comm.Send(siggrid, dest=0, tag=78)