#
from __future__ import division
from pystellar.density import density
from pystellar.opacity import OpacityTable
from pystellar.threading import ObjectsManager, EngineManager
import numpy as np
import time
import logging
X = 0.70
Y = 0.28
log = logging.getLogger('pystellar.opacity')
log.setLevel(logging.DEBUG)
console = logging.StreamHandler()
console.setLevel(logging.DEBUG)
log.addHandler(console)
OP = OpacityTable(fkey='OP17', X=0.70, Y=0.28, efkey='cunha06')
P, T = [8.726086186677213013e+07, 4.576702504411734481e+03]
rho = density(P=P, T=T, X=X, Y=Y)
print OP.kappa(rho=rho, T=T)
print OP.make_points(logT=np.log10(T), logrho=np.log10(rho))
T = 11264.2381423
rho = 0.0044603404639
print OP.kappa(T=T, rho=rho)
Opacity.start() print "" print "Settings for the Sun" print "M = %10.6g g" % Ms print "X = %10.3f" % X print "Y = %10.3f" % Y print u"μ = %10.3f" % mmw(X=X,Y=Y) print "" m = 1.9890999983295044e+33 print "Inner Boundary Conditions" print " Initial Step:" print " m = %10.6e g" % m (r,l,P,T) = [ 6.96000000e+10, 2.13492744e+33, 3.64246107e+27, 5.77489402e+03] print "At m=m:" print " r = %10.6e cm" % r print " l = %10.6e erg/s" % l print " P = %10.6e Dyne/cm^2" % P print " T = %10.6e K" % T print u" = %10.6e " % density(P=P,T=T,X=X,Y=Y) print "Derivatives" (dr,dl,dP,dT) = derivatives(np.array([m]),np.atleast_2d([r,l,P,T]),mu=mmw(X=X,Y=Y),optable=Opacity,X=X,XCNO=(1-X-Y),cfg=Config["Data.Energy"]) print " dr = %10.6e" % dr print " dl = %10.6e" % dl print " dP = %10.6e" % dP print " dT = %10.6e" % dT print "" Opacity.stop()
from __future__ import division
from pystellar.density import density
from pystellar.opacity import OpacityTable
from pystellar.threading import ObjectsManager, EngineManager
import numpy as np
import time
import logging
X = 0.70
Y = 0.28
log = logging.getLogger('pystellar.opacity')
log.setLevel(logging.DEBUG)
console = logging.StreamHandler()
console.setLevel(logging.DEBUG)
log.addHandler(console)
OP = OpacityTable(fkey='OP17',X=0.70,Y=0.28,efkey='cunha06')
P, T = [8.726086186677213013e+07, 4.576702504411734481e+03]
rho = density(P=P,T=T,X=X,Y=Y)
print OP.kappa(rho=rho,T=T)
print OP.make_points(logT=np.log10(T),logrho=np.log10(rho))
T= 11264.2381423
rho=0.0044603404639
print OP.kappa(T=T,rho=rho)