# abiding by the rules of distribution of free software. import numpy as np from scipy.constants import N_A try: from yt.data_objects.field_info_container import add_field except ImportError: def add_field(*cargs, **vargs): return None def _nion(field, data): return data['dens'] * data['sumy'] * N_A add_field('nion', function=_nion, take_log=True) def _nele(field, data): return data['dens'] * data['ye'] * N_A add_field('nele', function=_nele, take_log=True) def _abar(field, data): return 1.0 / data['sumy'] add_field('abar', function=_abar, take_log=False)
# This software is governed by the CeCILL-B license under French law and # abiding by the rules of distribution of free software. import numpy as np from scipy.constants import N_A try: from yt.data_objects.field_info_container import add_field except ImportError: def add_field(*cargs, **vargs): return None def _nion(field, data): return data['dens'] * data['sumy'] * N_A add_field ('nion', function=_nion, take_log=True) def _nele(field, data): return data['dens'] * data['ye'] * N_A add_field ('nele', function=_nele, take_log=True) def _abar(field, data): return 1.0 / data['sumy'] add_field ('abar', function=_abar, take_log=False) def _zbar(field, data): return data['ye'] / data['sumy']
"""By importing this module, common derivied fields are registered with yt.""" from yt.data_objects.field_info_container import add_field def _nele(field, data): return data['ye'] * data['dens'] * data['sumy'] * 6.022E23 add_field('nele', function=_nele, take_log=True) def _abar(field, data): return 1.0 / data['sumy'] add_field('abar', function=_abar, take_log=False)
"""By importing this module, common derivied fields are registered with yt.""" from yt.data_objects.field_info_container import add_field def _nion(field, data): return data['dens'] * data['sumy'] * 6.022E23 add_field('nion', function=_nion, take_log=True) def _abar(field, data): return 1.0 / data['sumy'] add_field('abar', function=_abar, take_log=False) def _velo(field, data): return (data['velx']**2 + data['vely']**2 + data['velz']**2)**0.5 add_field('velo', function=_velo, take_log=True) add_field('velr_rz', function=lambda field, data: data['velx'], take_log=True) add_field('velz_rz', function=lambda field, data: data['vely'], take_log=True) add_field('velt_rz', function=lambda field, data: data['velz'], take_log=True)