# 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)
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)