Halo radial profiles
Hello, I have two question: 1. I am trying to get radial profiles for my halos and trying to use the following instructions: https://yt-project.org/doc/cookbook/halo_analysis_example.html#halo-analysis... I am directly copy-pasting the code but once I try to create the catalog with hc.create() it gives me error, saying: "bins must be monotonically increasing or decreasing". I do not understand what can be the problem, since as I can see from the output it creates the catalog object for many halos but at some point it gives the error. 2. I wonder if yt features can be used for hdf5 file (output from ENZO, written in hdf5)? I.e if I can analyze hdf5 data with yt. Thanks, Salome
Hi Salome, All of the yt features you discuss should work just fine with Enzo data. I can't seem to find that exact error message, so it's possible you're working with an outdated version of yt, at least with respect to this functionality. If you run "yt instinfo", you can find out what version you're on. If you're able to update to a later version and are still getting such an error, I would guess it's because the halos are too small to be profiled. Sorry it took so long for you to get a response to this message. Please do follow up if you have additional questions and I should be able to get back to you sooner. Britton On Sat, May 4, 2019 at 3:14 PM Salome Mtchedlidze <salomchedlidze@gmail.com> wrote:
Hello,
I have two question:
1. I am trying to get radial profiles for my halos and trying to use the following instructions:
https://yt-project.org/doc/cookbook/halo_analysis_example.html#halo-analysis...
I am directly copy-pasting the code but once I try to create the catalog with
hc.create()
it gives me error, saying: "bins must be monotonically increasing or decreasing". I do not understand what can be the problem, since as I can see from the output it creates the catalog object for many halos but at some point it gives the error.
2. I wonder if yt features can be used for hdf5 file (output from ENZO, written in hdf5)? I.e if I can analyze hdf5 data with yt.
Thanks, Salome _______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org
Hi Britton, Thank you very much, turns out that it was the problem. I have one more question - I wanted to create not only dm halos but visible mass halos as well. I was following these instructions: http://yt-project.org/doc/reference/api/yt.analysis_modules.halo_finding.hal... where for HOP halo finder it was said that dm_only is depreciated and I should use ptype =all ifor including all particles. I tried this with the code: hc = HaloCatalog(data_ds=data_ds, finder_method='hop', ptype=all, output_dir="halo_catalogs/..../") but it was giving me the error: __init__() got an unexpected keyword argument 'ptype' I tried this once I updated the yt version but the result was same.
Hi Salome, Great to hear that's solved. When running the halo finders through the HaloCatalog, keyword arguments meant for the halo finder are passed via the "finder_kwargs" keyword. This should be a dictionary of all keywords meant to passed on. In your case, this would look something like: finder_kwargs={"ptype": "all"} There is more documentation here: https://yt-astro-analysis.readthedocs.io/en/latest/halo_catalogs.html#halo-f... Note, this is for using the yt_astro_analysis, which is an external package that has replaced the internal yt analysis modules. You can find more information on installation and modifying your imports in those docs. Britton On Thu, May 23, 2019 at 1:07 PM Salome Mtchedlidze <salomchedlidze@gmail.com> wrote:
Hi Britton,
Thank you very much, turns out that it was the problem. I have one more question - I wanted to create not only dm halos but visible mass halos as well. I was following these instructions:
http://yt-project.org/doc/reference/api/yt.analysis_modules.halo_finding.hal...
where for HOP halo finder it was said that dm_only is depreciated and I should use ptype =all ifor including all particles. I tried this with the code:
hc = HaloCatalog(data_ds=data_ds, finder_method='hop', ptype=all, output_dir="halo_catalogs/..../")
but it was giving me the error:
__init__() got an unexpected keyword argument 'ptype'
I tried this once I updated the yt version but the result was same. _______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org
participants (2)
-
Britton Smith -
Salome Mtchedlidze