import pyMesaUtils as pym const_lib, const_def = pym.loadMod("const") crlibm_lib, _ = pym.loadMod("crlibm") crlibm_lib.crlibm_init() ierr = 0 const_lib.const_init(pym.MESA_DIR, ierr) chem_lib, chem_def = pym.loadMod("chem") chem_lib.chem_init('isotopes.data', ierr)
import pyMesaUtils as pym import sys utils_lib, utils_def = pym.loadMod("utils") utils_lib.utils_omp_get_max_threads() utils_lib.utils_omp_get_thread_num()
import pyMesaUtils as pym crlibm_lib, _ = pym.loadMod("crlibm") crlibm_lib.crlibm_init() crlibm_lib.ln10 crlibm_lib.exp_cr(100.0)
import pyMesaUtils as pym col_lib, col_def = pym.loadMod("colors")
from __future__ import print_function import pyMesaUtils as pym import numpy as np import os import matplotlib.pyplot as plt const_lib, const_def = pym.loadMod("const") crlibm_lib, _ = pym.loadMod("crlibm") chem_lib, chem_def = pym.loadMod("chem") rates_lib, rates_def = pym.loadMod("rates") ierr = 0 crlibm_lib.crlibm_init() const_lib.const_init(pym.MESA_DIR, ierr) chem_lib.chem_init('isotopes.data', ierr) if pym.MESA_VERSION >= 11213: rates_lib.rates_init('reactions.list', 'jina_reaclib_results_20171020_default', 'rate_tables', False, False, '', '', '', ierr) elif pym.MESA_VERSION >= 10000: #Function sig changed rates_lib.rates_init('reactions.list', 'jina_reaclib_results_20130213default2', 'rate_tables', False, False, '', '', '', ierr) else: rates_lib.rates_init('reactions.list', 'jina_reaclib_results_20130213default2', 'rate_tables', False, '', '', '', ierr)
import pyMesaUtils as pym import numpy as np pym.buildModule('kap') eos_lib, eos_def = pym.loadMod("eos") const_lib, const_def = pym.loadMod("const") crlibm_lib, _ = pym.loadMod("crlibm") chem_lib, chem_def = pym.loadMod("chem") kap_lib, kap_def = pym.loadMod("kap") ierr = 0 const_lib.const_init('', ierr) crlibm_lib.crlibm_init() chem_lib.chem_init('isotopes.data', ierr) if pym.MESA_VERSION >= 11354: kap_lib.kap_init('gs98', 'gs98_co', 'lowT_fa05_gs98', 3.88, 3.80, True, pym.KAP_CACHE, '', False, ierr) elif pym.MESA_VERSION >= 10398 and pym.MESA_VERSION < 11354: kap_lib.kap_init('gs98', 'gs98_co', 'lowT_fa05_gs98', 3.88, 3.80, True, pym.KAP_CACHE, '', ierr) else: kap_lib.kap_init('gs98', 'gs98_co', 'lowT_fa05_gs98', 3.88, 3.80, 3.80, True, pym.KAP_CACHE, '', ierr) kap_handle = kap_lib.alloc_kap_handle(ierr) if pym.MESA_VERSION >= 10398: kap_lib.kap_set_choices(kap_handle, False, False, True, True, True, 0.71,
import pyMesaUtils as pym import numpy as np eos_lib, eos_def = pym.loadMod("eos") const_lib, const_def = pym.loadMod("const") crlibm_lib, _ = pym.loadMod("crlibm") chem_lib, chem_def = pym.loadMod("chem") net_lib, net_def = pym.loadMod("net") rates_lib, rates_def = pym.loadMod("rates") ierr = 0 crlibm_lib.crlibm_init() const_lib.const_init(pym.MESA_DIR, ierr) chem_lib.chem_init('isotopes.data', ierr) ierr = 0 if pym.MESA_VERSION >= 10000: #Function sig changed rates_lib.rates_init('reactions.list', 'jina_reaclib_results_20130213default2', 'rate_tables', False, False, '', '', '', ierr) else: rates_lib.rates_init('reactions.list', 'jina_reaclib_results_20130213default2', 'rate_tables', False, '', '', '', ierr) net_lib.net_init(ierr) eos_lib.eos_init('mesa', '', '', '', False, ierr)
from __future__ import print_function import pyMesaUtils as pym import numpy as np import os eos_lib, eos_def = pym.loadMod("eos") const_lib, const_def = pym.loadMod("const") crlibm_lib, _ = pym.loadMod("crlibm") chem_lib, chem_def = pym.loadMod("chem") net_lib, net_def = pym.loadMod("net") rates_lib, rates_def = pym.loadMod("rates") kap_lib, kap_def = pym.loadMod("kap") ion_lib, ion_def = pym.loadMod("ionization") star_lib, star_def = pym.loadMod("star") pym.buildRunStarSupport() run_star, _ = pym.loadMod("run_star_support") ierr = 0 crlibm_lib.crlibm_init() const_lib.const_init(pym.MESA_DIR, ierr) chem_lib.chem_init('isotopes.data', ierr)
ax.set_xlabel(kwargs['arg2']) ax.set_ylabel(kwargs['arg']) else: ax.set_xlabel(kwargs['arg']) ax.set_ylabel(kwargs['arg2']) fig.savefig(name) # plt.show() plt.close(fig) mod="net" pym.buildModule(mod) eos_lib, eos_def = pym.loadMod("eos") const_lib, const_def = pym.loadMod("const") crlibm_lib, _ = pym.loadMod("crlibm") chem_lib, chem_def = pym.loadMod("chem") net_lib, net_def = pym.loadMod("net") rates_lib, rates_def = pym.loadMod("rates") kap_lib, kap_def = pym.loadMod("kap") ion_lib, ion_def = pym.loadMod("ionization") ierr=0 crlibm_lib.crlibm_init() const_lib.const_init(pym.MESA_DIR,ierr) chem_lib.chem_init('isotopes.data',ierr) if pym.MESA_VERSION >= 11213:
import pyMesaUtils as pym const_lib, const_def = pym.loadMod("const") print(const_def.a2rad) print(const_def.thermohaline_mixing) print(const_def.mev_to_ergs) ierr = 0 const_lib.const_init(pym.MESA_DIR, ierr) print(const_def.mev_to_ergs)
cmap='bwr') cb = fig.colorbar(cax) cb.set_label('log abs neu ') ax.set_title(title) ax.set_xlabel('logT') ax.set_ylabel('LogRho') fig.savefig(name) plt.close(fig) mod = "neu" pym.buildModule(mod) const_lib, const_def = pym.loadMod("const") crlibm_lib, _ = pym.loadMod("crlibm") ierr = 0 crlibm_lib.crlibm_init() const_lib.const_init(pym.MESA_DIR, ierr) neu_lib, neu_def = pym.loadMod(mod) prefix = sys.argv[1] steps = 200 tmin = 10**6.9 tmax = 10**10.0 rmin = 10**-1
import pyMesaUtils as pym import numpy as np neu_lib, neu_def = pym.loadMod("neu") T = 10**9.0 log10_T = np.log10(T) Rho = 10**9.0 log10_Rho = np.log10(T) abar = 1.0 zbar = 1.0 z2bar = 1.0 log10_Tlim = 7.5 flags = np.zeros(neu_def.num_neu_types.get()) flags[:] = True loss = np.zeros(neu_def.num_neu_rvs.get()) sources = np.zeros((neu_def.num_neu_types.get(), neu_def.num_neu_rvs.get())) info = 0 res = neu_lib.neu_get(T, log10_T, Rho, log10_Rho, abar, zbar, z2bar, log10_Tlim, flags, loss, sources, info)
import pyMesaUtils as pym import numpy as np const_lib, const_def = pym.loadMod("const") crlibm_lib, _ = pym.loadMod("crlibm") ion_lib, ion_def = pym.loadMod("ionization") ierr = 0 crlibm_lib.crlibm_init() const_lib.const_init(pym.MESA_DIR, ierr) res = ion_lib.ionization_init('ion', '', pym.ION_CACHE, False, ierr) Z = 0.02 X = 0.78 Rho = 10**9 log10Rho = np.log10(Rho) T = 10**9 log10T = np.log10(T) res = np.zeros(ion_def.num_ion_vals.get()) ierr = 0 res = ion_lib.eval_ionization(Z, X, Rho, log10Rho, T, log10T, res, ierr)
import pyMesaUtils as pym const_lib, const_def = pym.loadMod("const") crlibm_lib, _ = pym.loadMod("crlibm") crlibm_lib.crlibm_init() ierr = 0 const_lib.const_init(pym.MESA_DIR, ierr) atm_lib, atm_def = pym.loadMod("atm") atm_lib.eval_lrad(1.0, 1.0, 1.0, 1.0, 1.0, 1.0)