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