# 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")
