from yt.mods import * from matplotlib import use; use('Agg') import matplotlib.pyplot as plt pf = load("RedshiftOutput0002") vl=pf["RefineRegionLeftEdge"] vr=pf["RefineRegionRightEdge"] width=(vr-vl).max() width=0.1650390625/2 c=(vl+vr)/2 sv = pf.h.region((vl+vr)/2, c-width, c+width) #pc = PlotCollection(pf, center=c) #pc.add_projection('Density',2,fig_size=[5,5],use_colorbar=False,weight_field='Density',data_source=sv) #pc.set_width(2*width,'1') #pc.set_zlim(1e-27,2e-25) #pc.save("refined_region") #del pc #pc = PlotCollection(pf, center=c) #pc.add_projection('Density',2,fig_size=[5,5],use_colorbar=True,weight_field='Density',data_source=sv) #pc.set_width(2*width,'1') #pc.set_zlim(5e-28,1e-25) #pc.save("refined_region_colorbar") #del pc #pc = PlotCollection(pf, center=c) #pc.set_cmap('hot') #pc.add_projection('Temperature',2,fig_size=[5,5],use_colorbar=False,weight_field='Density',data_source=sv) #pc.set_width(2*width,'1') #pc.set_zlim(10,5e4) #pc.set_cmap('hot') #pc.save("refined_region") #del pc #pc = PlotCollection(pf, center=c) #pc.set_cmap('gist_stern') #pc.add_projection('TotalMetallicity',2,fig_size=[5,5],use_colorbar=False,weight_field='Density',data_source=sv) #pc.set_width(2*width,'1') #pc.set_zlim(1e-6,1e-1) #pc.set_cmap('gist_stern') #pc.save("refined_region") #del pc c=na.array([0.4882009076, 0.4725569205, 0.4654072348]) sv = pf.h.region(c,c-5.0/pf['kpc'],c+5.0/pf['kpc']) pc = PlotCollection(pf, center=c) pc.add_projection('Density',2,fig_size=[5,5],use_colorbar=False,weight_field='Density',data_source=sv) pc.set_width(10,'kpc') pc.set_zlim(1e-26,5e-24) pc.save("most_massive_halo") del pc #pc = PlotCollection(pf, center=c) #pc.set_cmap('hot') #pc.add_projection('Temperature',2,fig_size=[5,5],use_colorbar=False,weight_field='Density',data_source=sv) #pc.set_width(10,'kpc') #pc.set_zlim(1e3,1e6) #pc.set_cmap('hot') #pc.save("most_massive_halo") #del pc #pc = PlotCollection(pf, center=c) #pc.set_cmap('gist_stern') #pc.add_projection('TotalMetallicity',2,fig_size=[5,5],use_colorbar=False,weight_field='Density',data_source=sv) #pc.set_width(10,'kpc') #pc.set_zlim(1e-6,1e0) #pc.set_cmap('gist_stern') #pc.save("most_massive_halo") #del pc