args = dict(
rho_a=denss,
t_a=temps,
dsm=2.7,
csm=3.5e5,
dhvapm=3e10,
tboil=2.5e3,
ideg=1,
z=[13],
drek=[3.0],
fact=[100],
nsigma=[0],
)
zbar2 = pysnop.run(**args)['zbar']
#opp.plot_zbar(denss, temps, zbar, cond_data_i['zmax'], plt.figure())
#opp.plot_zbar(denss, temps, zbar-zbar2, cond_data_i['zmax'], plt.figure())
# <markdowncell>
# Ion EoS - Pressure
# ------------------
# <codecell>
# write ionmix file with
#RHO, TE = np.meshgrid(eos_data.origin['ioncc_dens'], eos_data.origin['ioncc_temps'])
#mask = 1e5*(np.abs(np.log10(RHO.T) - 1)**(-30))*np.exp(-np.abs(np.log(TE.T))**2)
#eos_data.origin["ioncc_pres"] = np.where(eos_data_i["ioncc_pres"] <= 0,
args= dict(
rho_a= denss,
t_a= temps,
dsm=2.7,
csm=3.5e5,
dhvapm=3e10,
tboil=2.5e3,
ideg=1,
z=[13],
drek=[3.0],
fact=[100],
nsigma=[0],
)
zbar2 = pysnop.run(**args)['zbar']
#opp.plot_zbar(denss, temps, zbar, cond_data_i['zmax'], plt.figure())
#opp.plot_zbar(denss, temps, zbar-zbar2, cond_data_i['zmax'], plt.figure())
# <markdowncell>
# Ion EoS - Pressure
# ------------------
# <codecell>
# write ionmix file with
#RHO, TE = np.meshgrid(eos_data.origin['ioncc_dens'], eos_data.origin['ioncc_temps'])
#mask = 1e5*(np.abs(np.log10(RHO.T) - 1)**(-30))*np.exp(-np.abs(np.log(TE.T))**2)
#eos_data.origin["ioncc_pres"] = np.where(eos_data_i["ioncc_pres"] <= 0,
# compute Zbar
#eos_w['ele_dens'] *= 0.5
# merge Ion and Ele grids
eos_o = opp.adapt.EosMergeGrids(eos_w,
filter_dens=lambda x: (x>0),
filter_temps=lambda x: (x>0.025),
thresh=[])
dens_arr, temp_arr = np.meshgrid(eos_o['ele_dens'], eos_o['ele_temps'] )
zbar = pysnop.run( dens_arr, temp_arr, np.logspace(0,4,100),
lte=True,
**mat_dict.snop)['zbar']
zbar_tf = thomas_fermi_ionization(dens_arr, temp_arr, mat_dict.Z, mat_dict.A).T
# ok there is a messup somewhere!
zbar = (zbar.T).reshape(eos_o['ele_pres'].shape)
#zbar = np.zeros(eos_w['ele_pres'].shape)
#zbar[1:,1:] = zbar2.T
#zbar[0,:] = eos_w['zbar'][1,:]
#zbar[:,0] = eos_w['zbar'][:,1]
CheckEosConsistency(eos_o)
# Lookig what is happenning at room temperature and solid density
rho0 = hedp.matdb(mat).rho0
rho0_idx = np.argmin(np.abs(eos_o['ele_dens']-0.92))
