Spacetime Diagrams

In [1]:

    

%pylab nbagg
from ppmpy import ppm


    

In [ ]:

    

FV_cmap_str = '''
Cnot:   0 0.000000000 0.000000000 0.000000000
Cnot:   24 0.000000000 0.000000000 0.392156900
Cnot:   41 0.000000000 0.501960800 0.501960800
Cnot:   81 0.000000000 1.000000000 1.000000000
Cnot:   115 1.000000000 1.000000000 1.000000000
Cnot:   153 1.000000000 1.000000000 0.000000000
Cnot:   197 1.000000000 0.000000000 0.000000000
Cnot:   238 0.250980400 0.000000000 0.000000000
Cnot:   255 0.000000000 0.000000000 0.000000000
Anot:   0 0.000000000
Anot:   8 0.031372550
Anot:   20 0.062745100
Anot:   32 0.784313700
Anot:   133 1.000000000
Anot:   153 0.894117700
Anot:   197 0.243137300
Anot:   241 0.258823500
Anot:   255 0.000000000
'''


    

In [24]:

    

prof= ppm.yprofile('/data/ppm_rpod2/YProfiles/AGBTP_M2.0Z1.e-5/F5/')
cmap = ppm.colourmap_from_str(FV_cmap_str, segment=(5, 251))
ifig = 1; pl.close(ifig); fig = pl.figure(ifig, figsize = (7., 2.25), dpi=150)
prof.spacetime_diagram('FV H+He', 1170, fig, rlim=(8., 33.), vlim=(1e-8, 1e0), cmap=cmap)
fig.tight_layout()


    

    















    












    




Processing will be done in 258 s.
Processing will be done in 238 s.
Processing will be done in 215 s.
Processing will be done in 192 s.
Processing will be done in 147 s.
Processing will be done in 112 s.
Processing will be done in 80 s.
Processing will be done in 55 s.
Processing will be done in 33 s.
Processing will be done in 10 s.






    




/usr/local/lib/python2.7/dist-packages/matplotlib/axes/_base.py:1298: UnicodeWarning: Unicode equal comparison failed to convert both arguments to Unicode - interpreting them as being unequal
  if aspect == 'normal':
/usr/local/lib/python2.7/dist-packages/matplotlib/axes/_base.py:1303: UnicodeWarning: Unicode equal comparison failed to convert both arguments to Unicode - interpreting them as being unequal
  elif aspect in ('equal', 'auto'):

In [26]:

    

airmu = 1.39165
cldmu = 0.725
fkair = 0.203606102635
fkcld = 0.885906040268
AtomicNoair = 6.65742024965
AtomicNocld = 1.34228187919
zero_intervals = None
ifig = 30; pl.close(ifig); fig = pl.figure(ifig, figsize = (7., 2.25), dpi=150)
prof.spacetime_diagram('enuc_C12pg', 1170, fig, rlim=(8., 33.), vlim=(1e9, 1e13), \
                       zero_intervals=zero_intervals, cmap='inferno', airmu=airmu, \
                       cldmu=cldmu, fkair=fkair, fkcld=fkcld, AtomicNoair=AtomicNoair, \
                       AtomicNocld=AtomicNocld, corr_fact=1.5)
fig.tight_layout()


    

    















    












    




Processing will be done in 1648 s.
Processing will be done in 1764 s.
Processing will be done in 1727 s.
Processing will be done in 1720 s.
Processing will be done in 1723 s.
Processing will be done in 1696 s.
Processing will be done in 1684 s.
Processing will be done in 1656 s.
Processing will be done in 1626 s.
Processing will be done in 1607 s.
Processing will be done in 1580 s.
Processing will be done in 1550 s.
Processing will be done in 1526 s.
Processing will be done in 1503 s.
Processing will be done in 1479 s.
Processing will be done in 1459 s.
Processing will be done in 1435 s.
Processing will be done in 1412 s.
Processing will be done in 1390 s.
Processing will be done in 1367 s.
Processing will be done in 1342 s.
Processing will be done in 1319 s.
Processing will be done in 1284 s.
Processing will be done in 1233 s.
Processing will be done in 1186 s.
Processing will be done in 1137 s.
Processing will be done in 1093 s.
Processing will be done in 1054 s.
Processing will be done in 1010 s.
Processing will be done in 969 s.
Processing will be done in 929 s.
Processing will be done in 891 s.
Processing will be done in 855 s.
Processing will be done in 820 s.
Processing will be done in 785 s.
Processing will be done in 750 s.
Processing will be done in 716 s.
Processing will be done in 683 s.
Processing will be done in 649 s.
Processing will be done in 615 s.
Processing will be done in 583 s.
Processing will be done in 552 s.
Processing will be done in 521 s.
Processing will be done in 490 s.
Processing will be done in 459 s.
Processing will be done in 428 s.
Processing will be done in 397 s.
Processing will be done in 365 s.
Processing will be done in 335 s.
Processing will be done in 304 s.
Processing will be done in 276 s.
Processing will be done in 251 s.
Processing will be done in 225 s.
Processing will be done in 195 s.
Processing will be done in 166 s.
Processing will be done in 138 s.
Processing will be done in 110 s.
Processing will be done in 82 s.
Processing will be done in 53 s.
Processing will be done in 24 s.

In [28]:

    

ifig = 3; pl.close(ifig); fig = pl.figure(ifig, figsize = (7., 2.25), dpi=150)
prof.spacetime_diagram('Ek', 1170, fig, rlim=(8., 33.), vlim=(1e10, 3e12))
fig.tight_layout()


    

    















    












    




Processing will be done in 233 s.
Processing will be done in 221 s.
Processing will be done in 202 s.
Processing will be done in 179 s.
Processing will be done in 138 s.
Processing will be done in 103 s.
Processing will be done in 73 s.
Processing will be done in 43 s.
Processing will be done in 14 s.