# H2 profile sphere = ds.sphere( keskusta, leveys) # keskitetään H2 fraction -kuvaaja säteilylähteeseen H2_prof = yt.ProfilePlot(sphere, "radius", ["H2_fraction"]) H2_prof.set_unit("radius", "kpc") H2_prof.set_xlim(prof_xrange[0], prof_xrange[1]) H2_prof.set_ylim("H2_fraction", 1e-12, 2e-6) H2_prof.set_line_property("linewidth", linewidth) H_prof.x_log = False plots.append(H2_prof) colorbar_flags.append(False) xaxis_flags.append(0) yaxis_flags.append(1) #H2_prof.save("kuvat/H2prof.png") multi = eps.multiplot(sarakkeet, rivit, plots, bare_axes=False, cb_labels=colorbar_label, cb_flags=colorbar_flags, xaxis_flags=xaxis_flags, yaxis_flags=yaxis_flags) multi.scale_line(label="%.1f kpc" % (0.5 * leveys), size=0.5, loc=(0.05, 1.08)) multi.title_box("z = %2.6f" % (ds.current_redshift), loc=(0.1, 1.95), color=white, bgcolor=None, text_opts=[pyx.text.size.large]) multi.save_fig("kuvat/nelikko", format="png")
H_prof.set_unit("radius", "kpc") H_prof.set_xlim(prof_xrange[0], prof_xrange[1]) H_prof.set_ylim("H_number_density", 0.015, 0.13) H_prof.set_line_property("linewidth", linewidth) H_prof.x_log = False #H_prof.save("kuvat/Hprof.png") plots.append(H_prof) colorbar_flags.append(False) xaxis_flags.append(0) yaxis_flags.append(0) # H2 profile sphere = ds.sphere(keskusta, leveys) # keskitetään H2 fraction -kuvaaja säteilylähteeseen H2_prof = yt.ProfilePlot(sphere, "radius", ["H2_fraction"]) H2_prof.set_unit("radius", "kpc") H2_prof.set_xlim(prof_xrange[0], prof_xrange[1]) H2_prof.set_ylim("H2_fraction", 1e-12, 2e-6) H2_prof.set_line_property("linewidth", linewidth) H_prof.x_log = False plots.append(H2_prof) colorbar_flags.append(False) xaxis_flags.append(0) yaxis_flags.append(1) #H2_prof.save("kuvat/H2prof.png") multi = eps.multiplot(sarakkeet, rivit, plots, bare_axes=False, cb_labels=colorbar_label, cb_flags = colorbar_flags, xaxis_flags = xaxis_flags, yaxis_flags = yaxis_flags) multi.scale_line(label="%.1f kpc" % (0.5*leveys),size=0.5, loc=(0.05,1.08)) multi.title_box("z = %2.6f" % (ds.current_redshift), loc=(0.1,1.95), color=white, bgcolor=None, text_opts=[pyx.text.size.large]) multi.save_fig("kuvat/nelikko", format="png")
width=(leveys, 'kpc')) tiheys.set_unit("Density", "g/cm**3") tiheys.set_cmap(field="Density", cmap='hot') tiheys.save("kuvat/tiheys.png") plotit.append(tiheys) xaxis_flags.append(0) # x-akseli alapuolelle yaxis_flags.append(0) # y-akseli vasemmalle cb_location.append("top") # väripalkki yläpuolelle #toinen kuva lampo = yt.SlicePlot(ds, 1, "Temperature", center=tiheysmaksimi, width=(leveys, 'kpc')) lampo.set_cmap(field="Temperature", cmap='hot') lampo.save("kuvat/lampo.png") plotit.append(lampo) xaxis_flags.append(0) # x-akseli alapuolelle yaxis_flags.append(1) # y-akseli oikealle cb_location.append("top") # väripalkki yläpuolelle multi = eps.multiplot(sarakkeet, rivit, plotit, bare_axes=False, xaxis_flags=xaxis_flags, yaxis_flags=yaxis_flags, cb_location=cb_location) multi.save_fig("kuvat/EPSMultiPlot", format="png")