Exemplo n.º 1
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")
Exemplo n.º 2
0
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")
Exemplo n.º 3
0
                      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")