[yt-users] help: yt Halo finder for Nyx data

[Dr. Surajit Paul]

Hi,
   I was just trying to compute the halos from a Nyx data. While running the halo finder it stops saying
........
yt.data_objects.field_info_container.NeedsDataField: (['particle_mass'])

   Below is the code, results and the list of fields available. In the list I have "particle_mass_density" instead of
"particle_mass". Please suggest me how I can make it running.

Thanks,
surajit

from yt.mods import *
import fields

# Load the dataset.
pf = load("/gfs1/work/nipiagws/MUSIC/102-26468_512zoom2a_amr2_sgs2_si_600_05_flat/plt02593_z01")
halos = HaloFinder(pf)

halos.write_out("%s_halos.txt")



Traceback (most recent call last):
  File "halos.py", line 14, in <module>
    halos = HaloFinder(pf)
  File "/home/h/nipiagws/yt-x86_64/src/yt-hg/yt/analysis_modules/halo_finding/halo_objects.py", line 2484, in __init__
    self.comm.mpi_allreduce((self._data_source["ParticleMassMsun"][select]).sum(dtype='float64'), op='sum')
  File "/home/h/nipiagws/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 318, in __getitem__
    self.get_data(key)
  File "/home/h/nipiagws/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 2619, in get_data
    if self._generate_field(field):
  File "/home/h/nipiagws/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 345, in _generate_field
    self._generate_field_in_grids(field, ngt_exception.ghost_zones)
  File "/home/h/nipiagws/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 355, in _generate_field_in_grids
    grid[field] = self.__touch_grid_field(grid, field)
  File "/home/h/nipiagws/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 81, in save_state
    tr = func(self, grid, field, *args, **kwargs)
  File "/home/h/nipiagws/yt-x86_64/src/yt-hg/yt/data_objects/data_containers.py", line 359, in __touch_grid_field
    return grid[field]
  File "/home/h/nipiagws/yt-x86_64/src/yt-hg/yt/data_objects/grid_patch.py", line 147, in __getitem__
    self.get_data(key)
  File "/home/h/nipiagws/yt-x86_64/src/yt-hg/yt/data_objects/grid_patch.py", line 190, in get_data
    self._generate_field(field)
  File "/home/h/nipiagws/yt-x86_64/src/yt-hg/yt/data_objects/grid_patch.py", line 122, in _generate_field
    self.pf.field_info[field].check_available(self)
  File "/home/h/nipiagws/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py", line 356, in check_available
    validator(data)
  File "/home/h/nipiagws/yt-x86_64/src/yt-hg/yt/data_objects/field_info_container.py", line 450, in __call__
    raise NeedsDataField(doesnt_have)
yt.data_objects.field_info_container.NeedsDataField: (['particle_mass'])


filed_list

AveMomX
AveMomY
AveMomZ
ErrCovarX
ErrCovarY
ErrCovarZ
Ne
Temp
density
diss_sgs
divu
grav_x
grav_y
grav_z
magvort
particle_mass_density
phi_grav
prod_sgs
rho_E
rho_H
rho_He
rho_K
rho_e
support_sgs
support_therm
support_turb
turb_src
xmom
ymom
zmom






_______________________________________________
yt-users mailing list
yt-users at lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org