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

Nathan Goldbaum nathan12343 at gmail.com
Wed Dec 23 14:27:27 PST 2015


I think the nyx frontend probably needs to create a particle_mass field.
This likely needs to be patched in yt itself.

Unfortunately I don't think we have any nyx cosmological datasets suitable
for testing with for this purpose. Would you be ok with sharing a test
dataset?

-Nathan

On Mon, Dec 21, 2015 at 9:25 AM, Dr. Surajit Paul <
surajit at physics.unipune.ac.in> wrote:

> 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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20151223/cd46be30/attachment.html>
-------------- next part --------------
_______________________________________________
yt-users mailing list
yt-users at lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org


More information about the yt-users mailing list