Esempio n. 1
0
from pystellar.opacity import OpacityTable
from pystellar.threading import ObjectThread

from pyshell.config import DottedConfiguration

Config = DottedConfiguration()
Config.dn = DottedConfiguration
Config.load("Star.yml")
X = 1.000
Y = 0.000
rho = 1
Opacity = ObjectThread(OpacityTable,
                       ikwargs=dict(fkey='GN93hz', X=X, Y=Y, snap=True),
                       locking=True,
                       timeout=None)
Opacity.start()

T = np.logspace(6.5, 7.5, 100)

epsilon_pp = epp(T=T, rho=rho, X=X, c=Config["Data.Energy"])
epsilon_CNO = eCNO(T=T, rho=rho, X=X, XCNO=0.01, c=Config["Data.Energy"])

plt.loglog(T, epsilon_CNO + epsilon_pp, '-', label='Total')
plt.loglog(T, epsilon_pp, ':', label="PP")
plt.loglog(T, epsilon_CNO, '--', label="CNO")
plt.legend()
plt.show()

Opacity.stop()
Esempio n. 2
0
from pyshell.config import DottedConfiguration

Config = DottedConfiguration()
Config.dn = DottedConfiguration
Config.load("pystellar/Star.yml")
Config.load("Star.yml")
X = 0.700
Y = 0.280
XCNO = 0.7 * (1 - X - Y)
rho = 80

T = np.logspace(6.8, 7.8, 100)

epsilon_pp = epp(T=T, rho=rho, X=X, c=Config["Data.Energy"])
epsilon_CNO = eCNO(T=T, rho=rho, X=X, XCNO=XCNO, c=Config["Data.Energy"])

T1 = 18 * np.power(10, 6)
epsilon_pp1 = epp(T=T1, rho=rho, X=X, c=Config["Data.Energy"])
epsilon_CNO1 = eCNO(T=T1, rho=rho, X=X, XCNO=XCNO, c=Config["Data.Energy"])

print u"X        = %g" % X
print u"XCNO     = %g" % XCNO
print u"ρ        = %g" % rho
print u"T        = %g" % T1
print u"εCNO     = %g" % epsilon_CNO1
print u"εpp      = %g" % epsilon_pp1
print u"ε        = %g" % (epsilon_pp1 + epsilon_CNO1)
print u"εpp/εCNO = %g" % (epsilon_pp1 / epsilon_CNO1)

plt.loglog(T, epsilon_CNO + epsilon_pp, '-', label=r"$\epsilon$")
Esempio n. 3
0
import matplotlib.pyplot as plt

from pystellar.energy import epp, eCNO

from pystellar.opacity import OpacityTable
from pystellar.threading import ObjectThread

from pyshell.config import DottedConfiguration

Config = DottedConfiguration()
Config.dn = DottedConfiguration
Config.load("Star.yml")
X = 1.000
Y = 0.000
rho = 1
Opacity = ObjectThread(OpacityTable,ikwargs=dict(fkey='GN93hz',X=X,Y=Y,snap=True),locking=True,timeout=None)
Opacity.start()

T = np.logspace(6.5,7.5,100)

epsilon_pp = epp(T=T,rho=rho,X=X,c=Config["Data.Energy"])
epsilon_CNO = eCNO(T=T,rho=rho,X=X,XCNO=0.01,c=Config["Data.Energy"])

plt.loglog(T,epsilon_CNO+epsilon_pp,'-',label='Total')
plt.loglog(T,epsilon_pp,':',label="PP")
plt.loglog(T,epsilon_CNO,'--',label="CNO")
plt.legend()
plt.show()

Opacity.stop()
Esempio n. 4
0
from pyshell.config import DottedConfiguration

Config = DottedConfiguration()
Config.dn = DottedConfiguration
Config.load("pystellar/Star.yml")
Config.load("Star.yml")
X = 0.700
Y = 0.280
XCNO = 0.7 * (1 - X - Y)
rho = 80

T = np.logspace(6.8,7.8,100)

epsilon_pp = epp(T=T,rho=rho,X=X,c=Config["Data.Energy"])
epsilon_CNO = eCNO(T=T,rho=rho,X=X,XCNO=XCNO,c=Config["Data.Energy"])

T1 = 18 * np.power(10,6)
epsilon_pp1 = epp(T=T1,rho=rho,X=X,c=Config["Data.Energy"])
epsilon_CNO1 = eCNO(T=T1,rho=rho,X=X,XCNO=XCNO,c=Config["Data.Energy"])

print u"X        = %g" % X
print u"XCNO     = %g" % XCNO
print u"ρ        = %g" % rho
print u"T        = %g" % T1
print u"εCNO     = %g" % epsilon_CNO1
print u"εpp      = %g" % epsilon_pp1
print u"ε        = %g" % (epsilon_pp1 + epsilon_CNO1)
print u"εpp/εCNO = %g" % (epsilon_pp1/epsilon_CNO1)