def ps():
for i in range(len(redshifts)):
data = IO.readmap("dbt_" + str("%.3f" % redshifts[i]))
ps = power_spectrum_1d(data)
ps = np.asarray(ps)
print np.shape(ps), type(ps)
IO.writedata(ps[0, :], "data/powerSpectra_dbt_" + str("%.3f" % redshifts[i]) + ".dat")
IO.writedata(ps[1, :], "data/powerSpectraFrequencies_dbt_" + str("%.3f" % redshifts[i]) + ".dat")
def plotdbt():
mini = -457
maxi=83
print mini,maxi
for i in range(start,len(redshifts)):
print "Doing redshift " + str(redshifts[i])+"..."
dbt = IO.readmap("dbt_"+str('%.3f' % redshifts[i]))
print np.min(dbt)
plot(dbt[mesh/2,:,:],i,"$\delta T_b$","dbt",mini,maxi,'seismic')
print "Complete"
def compare_dbtmaps():
cmapdbt = {'red': ((0.0, 0.0, 0.0),
(0.85, 0.0, 1.0),
(1.0, 1.0, 1.0)),
'green': ((0.0, 0.0, 0.0),
(0.85,0.0, 0.0),
(1.0, 1.0, 1.0)),
'blue': ((0.0, 0.0, 0.0),
(0.85, 1.0, 0.0),
(1.0, 1.0, 1.0))
}
plt.register_cmap(name='dbtmap',data=cmapdbt)
mini = -457
maxi=83
redshifts = [14.912,14.294,13.733]
fig, axes = plt.subplots(2,3,figsize=(12,7))
fig.text(0.5,0.02,"Distance (Mpc)",ha='center',va='center',fontsize=fontsize)
fig.text(0.06,0.55,"Distance (Mpc)",ha='center',va='center',rotation='vertical',fontsize=fontsize)
plt.gcf().subplots_adjust(right=0.15)
plt.gcf().subplots_adjust(bottom=0.15)
plt.subplots_adjust(left=None,bottom=None,right=None,top=None,wspace=0.02,hspace=0.02)
for i in range(len(redshifts)):#start,len(redshifts)):
print "Doing redshift " + str(redshifts[i])+"..."
dbt_wpls = IO.readmap("dbt_"+str('%.3f' % redshifts[i]),'data_wpls/')
im = axes[0,i].imshow(dbt_wpls[250/2,:,:],cmap='dbtmap',vmin=mini,vmax=maxi,origin='lower')
axes[0,i].text(20,220,"z = "+str(redshifts[i]),fontsize=fontsize,color='white')
dbt_wstars= IO.readmap("dbt_"+str('%.3f' % redshifts[i]),'data_wstars/')
axes[1,i].imshow(dbt_wstars[250/2,:,:],cmap='dbtmap',vmin=mini,vmax=maxi,origin='lower')
fig.subplots_adjust(right=0.85)
cbar_ax = fig.add_axes([0.87, 0.17, 0.03, 0.7])
cbar = fig.colorbar(im, cax=cbar_ax,ticks=[-450,-400,-350,-300,-250,-200,-150,-100,-50,0,50,100,150,200,250,300])
cbar.ax.set_yticklabels
cbar.set_label("$\delta T_b \ (mK)$",size=fontsize)
#cbar_ax = fig.add_axes([0.87,0.54,0.03,0.35])
#cbar = fig.colorbar(im,cax=cbar_ax,ticks=[-300,-150,0, 150, 300])
#cbar.ax.set_yticklabels
#cbar.set_label("$\delta T_b \ (mK)$",size=fontsize)
cnv=250.0/244.0
for i in range(len(axes[:,0])):
for j in range(len(axes[0,:])):
axes[i,j].tick_params(axis='both', which='major', labelsize=numberfontsize, width = tickwidth, length = 11, direction='out',pad=14.0,top='off',right='off')
axes[i,j].set_xlim([0,250])
axes[i,j].set_ylim([0,250])
m0 = MultipleLocator(50.0*cnv)
axes[i,j].yaxis.set_major_locator(m0)
axes[i,j].xaxis.set_major_locator(m0)
axes[i,j].yaxis.set_major_formatter(mpl.ticker.FuncFormatter(mjfFormatter))
axes[i,j].xaxis.set_major_formatter(mpl.ticker.FuncFormatter(mjfFormatter))
m0 = MultipleLocator(10.0*cnv)
axes[i,j].yaxis.set_minor_locator(m0)
axes[i,j].xaxis.set_minor_locator(m0)
axes[0,0].tick_params(axis='both', which='minor', width = tickwidth, length = 5.5,direction='out',top='off',right='off',bottom='off')
axes[0,1].tick_params(axis='both', which='minor', width = tickwidth, length = 5.5,direction='out',top='off',right='off',bottom='off',left='off')
axes[0,2].tick_params(axis='both', which='minor', width = tickwidth, length = 5.5,direction='out',top='off',right='off',bottom='off',left='off')
axes[1,0].tick_params(axis='both', which='minor', width = tickwidth, length = 5.5,direction='out',top='off',right='off')
axes[1,1].tick_params(axis='both', which='minor', width = tickwidth, length = 5.5,direction='out',top='off',right='off',left='off')
axes[1,2].tick_params(axis='both', which='minor', width = tickwidth, length = 5.5,direction='out',top='off',right='off',left='off')
plt.setp(axes[0,:],xticks=[])
plt.setp(axes[:,1:3],yticks=[])
plt.savefig("compare_wstars_wpls/map_dbt.png")
plt.close()
