Python loadMod Beispiele

Beispiel #1

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)

Beispiel #2

Datei: utils.py Projekt: darthoctopus/pyMesa

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

Beispiel #3

import pyMesaUtils as pym

crlibm_lib, _ = pym.loadMod("crlibm")

crlibm_lib.crlibm_init()

crlibm_lib.ln10
crlibm_lib.exp_cr(100.0)

Beispiel #4

Datei: colors.py Projekt: darthoctopus/pyMesa

import pyMesaUtils as pym

col_lib, col_def = pym.loadMod("colors")

Beispiel #5

Datei: rates.py Projekt: darthoctopus/pyMesa

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)

Beispiel #6

Datei: kap.py Projekt: darthoctopus/pyMesa

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,

Beispiel #7

Datei: eos.py Projekt: darthoctopus/pyMesa

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)

Beispiel #8

Datei: star.py Projekt: darthoctopus/pyMesa

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)

Beispiel #9

Datei: pydfriddr_rates.py Projekt: darthoctopus/pyMesa

        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:

Beispiel #10

Datei: const.py Projekt: darthoctopus/pyMesa

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)

Beispiel #11

                    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

Beispiel #12

Datei: neu.py Projekt: darthoctopus/pyMesa

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)

Beispiel #13

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)

Beispiel #14

Datei: atm.py Projekt: darthoctopus/pyMesa

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)