def get_atomic_repo_config(): """ Get the repo configuration dictionary for the atomic data Returns ------- : dict """ atomic_repo_fname = get_internal_data_path("atomic_data_repo.yml") return yaml.load(open(atomic_repo_fname), Loader=yaml.CLoader)
import tardis from tardis.io.util import get_internal_data_path from IPython import get_ipython k_B_cgs = constants.k_B.cgs.value c_cgs = constants.c.cgs.value h_cgs = constants.h.cgs.value m_e_cgs = constants.m_e.cgs.value e_charge_gauss = constants.e.gauss.value logger = logging.getLogger(__name__) tardis_dir = os.path.realpath(tardis.__path__[0]) ATOMIC_SYMBOLS_DATA = (pd.read_csv( get_internal_data_path("atomic_symbols.dat"), delim_whitespace=True, names=["atomic_number", "symbol"], ).set_index("atomic_number").squeeze()) ATOMIC_NUMBER2SYMBOL = OrderedDict(ATOMIC_SYMBOLS_DATA.to_dict()) SYMBOL2ATOMIC_NUMBER = OrderedDict( (y, x) for x, y in ATOMIC_NUMBER2SYMBOL.items()) synpp_default_yaml_fname = get_internal_data_path("synpp_default.yaml") NUMERAL_MAP = tuple( zip( (1000, 900, 500, 400, 100, 90, 50, 40, 10, 9, 5, 4, 1), ("M", "CM", "D", "CD", "C", "XC", "L", "XL", "X", "IX", "V", "IV", "I"),
def get_atomic_repo_config(): atomic_repo_fname = get_internal_data_path('atomic_data_repo.yml') return yaml.load(open(atomic_repo_fname))
from astropy import units as u from pyne import nucname import tardis from tardis.io.util import get_internal_data_path k_B_cgs = constants.k_B.cgs.value c_cgs = constants.c.cgs.value h_cgs = constants.h.cgs.value m_e_cgs = constants.m_e.cgs.value e_charge_gauss = constants.e.gauss.value logger = logging.getLogger(__name__) tardis_dir = os.path.realpath(tardis.__path__[0]) ATOMIC_SYMBOLS_DATA = pd.read_csv(get_internal_data_path('atomic_symbols.dat'), delim_whitespace=True, names=['atomic_number', 'symbol']).set_index('atomic_number').squeeze() ATOMIC_NUMBER2SYMBOL = OrderedDict(ATOMIC_SYMBOLS_DATA.to_dict()) SYMBOL2ATOMIC_NUMBER = OrderedDict((y, x) for x, y in ATOMIC_NUMBER2SYMBOL.items()) synpp_default_yaml_fname = get_internal_data_path('synpp_default.yaml') NUMERAL_MAP = tuple(zip( (1000, 900, 500, 400, 100, 90, 50, 40, 10, 9, 5, 4, 1), ('M', 'CM', 'D', 'CD', 'C', 'XC', 'L', 'XL', 'X', 'IX', 'V', 'IV', 'I') )) class MalformedError(Exception): pass
return 'Expecting an atomic symbol (e.g. Fe) - supplied %s' % self.malformed_element_symbol class MalformedQuantityError(MalformedError): def __init__(self, malformed_quantity_string): self.malformed_quantity_string = malformed_quantity_string def __str__(self): return 'Expecting a quantity string(e.g. "5 km/s") for keyword - supplied %s' % self.malformed_quantity_string logger = logging.getLogger(__name__) tardis_dir = os.path.realpath(tardis.__path__[0]) ATOMIC_SYMBOLS_DATA = pd.read_csv( get_internal_data_path('atomic_symbols.dat'), delim_whitespace=True, names=['atomic_number', 'symbol']).set_index('atomic_number').squeeze() ATOMIC_NUMBER2SYMBOL = OrderedDict(ATOMIC_SYMBOLS_DATA.to_dict()) SYMBOL2ATOMIC_NUMBER = OrderedDict( (y, x) for x, y in ATOMIC_NUMBER2SYMBOL.items()) synpp_default_yaml_fname = get_internal_data_path('synpp_default.yaml') NUMERAL_MAP = tuple( zip((1000, 900, 500, 400, 100, 90, 50, 40, 10, 9, 5, 4, 1), ('M', 'CM', 'D', 'CD', 'C', 'XC', 'L', 'XL', 'X', 'IX', 'V', 'IV', 'I')))
class MalformedQuantityError(MalformedError): def __init__(self, malformed_quantity_string): self.malformed_quantity_string = malformed_quantity_string def __str__(self): return 'Expecting a quantity string(e.g. "5 km/s") for keyword - supplied %s' % self.malformed_quantity_string logger = logging.getLogger(__name__) tardis_dir = os.path.realpath(tardis.__path__[0]) ATOMIC_SYMBOLS_DATA = pd.read_csv(get_internal_data_path('atomic_symbols.dat'), delim_whitespace=True, names=['atomic_number', 'symbol']).set_index('atomic_number').squeeze() ATOMIC_NUMBER2SYMBOL = OrderedDict(ATOMIC_SYMBOLS_DATA.to_dict()) SYMBOL2ATOMIC_NUMBER = OrderedDict((y, x) for x, y in ATOMIC_NUMBER2SYMBOL.items()) synpp_default_yaml_fname = get_internal_data_path('synpp_default.yaml') NUMERAL_MAP = tuple(zip( (1000, 900, 500, 400, 100, 90, 50, 40, 10, 9, 5, 4, 1), ('M', 'CM', 'D', 'CD', 'C', 'XC', 'L', 'XL', 'X', 'IX', 'V', 'IV', 'I') )) def int_to_roman(i): result = []