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
)
