# Copyright 2012 Alexander Rudy. All rights reserved.
#
import numpy as np
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')
import numpy as np
from astropysics.constants import Rs, Ms, Lsun
from pystellar.initial import inner_boundary, outer_boundary
from pystellar.density import mmw
from pystellar.opacity import OpacityTable
from pystellar.stellar import derivatives
from pystellar.threading import ObjectsManager, ObjectPassthrough, EngineManager, ObjectThread
from pystellar.star import Star
from pyshell.config import DottedConfiguration
Config = DottedConfiguration()
Config.dn = DottedConfiguration
Config.load("pystellar/Star.yml")
Config.load("Star.yml")
print "Stellar Structure Problem Set #3"
X = Config["Star.Composition.X"]
Y = Config["Star.Composition.Y"]
# P_Guess_Iteration = 2
P_Guess_Iteration = 0
Opacity = ObjectThread(OpacityTable,
ikwargs=dict(fkey='GN93hz', X=X, Y=Y),
locking=True)
Opacity.start()
star = Star(config=Config, optable_args=Opacity.duplicator)
print ""
print "Problem #4&5:"
# Copyright 2012 Alexander Rudy. All rights reserved.
#
import numpy as np
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()
# Copyright 2012 Alexander Rudy. All rights reserved.
#
import numpy as np
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("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"])
