예제 #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)
예제 #2
0
파일: flash_fields.py 프로젝트: fimay/hedp
# 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']
예제 #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)
예제 #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)