def deidti(opimx, rho, tion, log=False): return opp.interpDT(opimx.eion, opimx.dens, opimx.temps, rho, tion, log=False, lookup=opp.INTERP_DFDT, bctmin=opp.BC_EXTRAP_ZERO)
def deedte(opimx, rho, tele, log=False): return opp.interpDT(opimx.eele, opimx.dens, opimx.temps, rho, tele, log=False, lookup=opp.INTERP_DFDT, bctmin=opp.BC_EXTRAP_ZERO)
def dpidd(opimx, rho, tion, log=False): return opp.interpDT(opimx.pion, opimx.dens, opimx.temps, rho, tion, log=False, lookup=opp.INTERP_DFDD, bctmin=opp.BC_EXTRAP_ZERO, bcdmin=opp.BC_EXTRAP_ZERO)
def pele(opimx, rho, tele, log=False): return opp.interpDT(opimx.pele, opimx.dens, opimx.temps, rho, tele, log=False, lookup=opp.INTERP_FUNC, bctmin=opp.BC_EXTRAP_ZERO, bcdmin=opp.BC_EXTRAP_ZERO)
def dpidd(opimx, rho, tion, log=False): return opp.interpDT( opimx.pion, opimx.dens, opimx.temps, rho, tion, log=False, lookup=opp.INTERP_DFDD, bctmin=opp.BC_EXTRAP_ZERO, bcdmin=opp.BC_EXTRAP_ZERO, )
def dpedte(opimx, rho, tele, log=False): return opp.interpDT( opimx.pele, opimx.dens, opimx.temps, rho, tele, log=False, lookup=opp.INTERP_DFDT, bctmin=opp.BC_EXTRAP_ZERO, bcdmin=opp.BC_EXTRAP_ZERO, )
import opacplot2 as opp from opacplot2 import RADCONST as AR, KB import numpy as np from matplotlib import pyplot as plt mpi = 4.002602 / opp.NA opimx = opp.OpacIonmix("examples/data/al-imx-002.cn4", mpi, twot=True, man=True, verbose=True) t1 = 0.0 t2 = 20.0 tvals = np.linspace(t1, t2, 1001) dens = 2.16 # PLOT EELE: vals = np.array([opp.interpDT(opimx.eele, opimx.dens, opimx.temps, dens, t, bctmin=opp.BC_EXTRAP_ZERO, log=False) for t in tvals]) plt.plot(tvals, vals, "r", lw=2) plt.grid(True) plt.ylabel("Electron Specific Internal Energy (ergs/g)") plt.xlabel("Electron Temperature (eV)") # PLOT CVEELE: plt.figure() vals = np.array([opp.interpDT(opimx.eele, opimx.dens, opimx.temps, dens, t, bctmin=opp.BC_EXTRAP_ZERO, log=False, lookup=opp.INTERP_DFDT) for t in tvals]) plt.plot(tvals, vals, "r", lw=2) plt.grid(True) plt.ylabel("Electron Specific Heat (ergs/g/eV)")
opp.plot_zbar(cond_data_i['zbar_dens'], cond_data_i['zbar_temps'], cond_data_i['zbar'], cond_data_i['zmax'], plt.figure()) # <codecell> ntemp = eos_data["ele_ntemp"] ndens = eos_data["ele_ndens"] temps = eos_data["ele_temps"] denss = eos_data["ele_dens"] # this part is not working as it is expected to... zbar = np.empty((ndens, ntemp)) for jd in xrange(ndens): for jt in xrange(ntemp): zbar[jd,jt] = opp.interpDT(cond_data_i["zbar"], cond_data_i["zbar_dens"], cond_data_i["zbar_temps"], denss[jd], temps[jt], bctmin=opp.BC_EXTRAP_ZERO) # <codecell> import pysnop args= dict( rho_a= denss, t_a= temps, dsm=2.7, csm=3.5e5, dhvapm=3e10, tboil=2.5e3,
cond_data_i['zbar'], cond_data_i['zmax'], plt.figure()) # <codecell> ntemp = eos_data["ele_ntemp"] ndens = eos_data["ele_ndens"] temps = eos_data["ele_temps"] denss = eos_data["ele_dens"] # this part is not working as it is expected to... zbar = np.empty((ndens, ntemp)) for jd in xrange(ndens): for jt in xrange(ntemp): zbar[jd, jt] = opp.interpDT(cond_data_i["zbar"], cond_data_i["zbar_dens"], cond_data_i["zbar_temps"], denss[jd], temps[jt], bctmin=opp.BC_EXTRAP_ZERO) # <codecell> import pysnop args = dict( rho_a=denss, t_a=temps, dsm=2.7, csm=3.5e5, dhvapm=3e10, tboil=2.5e3, ideg=1,
def deidti(opimx, rho, tion, log=False): return opp.interpDT( opimx.eion, opimx.dens, opimx.temps, rho, tion, log=False, lookup=opp.INTERP_DFDT, bctmin=opp.BC_EXTRAP_ZERO )
def eele(opimx, rho, tele, log=False): return opp.interpDT( opimx.eele, opimx.dens, opimx.temps, rho, tele, log=False, lookup=opp.INTERP_FUNC, bctmin=opp.BC_EXTRAP_ZERO )