|
1 | 1 | #!/usr/bin/env python |
2 | 2 |
|
| 3 | +''' |
| 4 | +This script creates a single PNG from a pltfile/ |
| 5 | +coloring regions by RGB code for composition. |
| 6 | +Red is proportional to iron group elements, |
| 7 | +blue is proportional to silicon group, |
| 8 | +green is proportional to Mg24 and lighter elements. |
| 9 | +''' |
| 10 | + |
3 | 11 | import matplotlib |
4 | 12 | matplotlib.use('agg') |
5 | 13 |
|
@@ -84,17 +92,17 @@ def _light_nuclei(field, data): |
84 | 92 |
|
85 | 93 |
|
86 | 94 | fig = plt.figure() |
87 | | -fig.set_size_inches(12.0, 9.0) |
| 95 | +fig.set_size_inches(12.0, 6.0) |
88 | 96 |
|
89 | | -width_frac = 0.1 |
| 97 | +# set this to 0.1 instead to zoom in to the core |
| 98 | +width_frac = 1 |
90 | 99 |
|
91 | | -center=[xmin + 0.25*width_frac*L_x, yctr, 0.0*cm] |
92 | | -width=[0.5*width_frac*L_x, width_frac*L_y] |
| 100 | +center=[xmin + 0.5*width_frac*L_x, yctr, 0.0*cm] |
| 101 | +width=[width_frac*L_x, width_frac*L_y, 0.0*cm] |
93 | 102 |
|
94 | 103 | slc = yt.SlicePlot(ds, "theta", |
95 | 104 | fields=[("gas", "iron_group"), ("gas", "si_group"), ("gas", "light_nuclei")], |
96 | | - center=[xmin + 0.25*width_frac*L_x, yctr, 0.0*cm], |
97 | | - width=[0.5*width_frac*L_x, width_frac*L_y, 0.0*cm], fontsize="12") |
| 105 | + center=center, width=width, fontsize="12") |
98 | 106 |
|
99 | 107 | res = (1024, 512) |
100 | 108 | frb = slc.data_source.to_frb(width[0], res, height=width[1]) #width, res)#, center=center) |
|
0 commit comments