Esempio n. 1
0
# 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)
Esempio n. 2
0
# 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']
Esempio n. 3
0
"""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)
Esempio n. 4
0
"""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)