Hello All, I'm trying to find halos in mutiple snapshot files from DM-only cosmological simulation generated by Gadget-2, using FOF halo finder in yt. The box size of the simulation is 13.8Mpc/h and contain 1024^3 particles. Here is my code, data_ds = yt.load('./snapshot_016.0') hc = HaloCatalog(data_ds=data_ds, finder_method='fof',finder_kwargs={"dm_only":False,"ptype":"Halo"}) hc.create() and I met the error like this, yt.utilities.exceptions.YTDomainOverflow: Particle bounds (-386.63174, -385.38947, -525.52057) and (6900.443, 6900.271, 13799.999) exceed domain bounds [0. 0. 0.] code_length and [13800. 13800. 13800.] code_length. I have already updated yt and it's dependecy to the version 3.6.1, but I still have this problem. It seems like yt treat the particle velocity as the bound, I have no idea. What make me more confusing is there is no such error when I run above code for another small scale simualtion (886^3 particles and L_box=12Mpc/h), So, does anyone know what happend here or has ever met the same problem?
Hi Zhang, I'm not sure, but maybe explicitly setting the bounding_box on the yt.load command would help, as in: https://yt-project.org/doc/cookbook/gadget_notebook.html -desika On Thu, May 20, 2021 at 5:44 AM Zhang Meng <zhangmeng182@mails.ucas.ac.cn> wrote:
Hello All,
I'm trying to find halos in mutiple snapshot files from DM-only cosmological simulation generated by Gadget-2, using FOF halo finder in yt. The box size of the simulation is 13.8Mpc/h and contain 1024^3 particles. Here is my code,
data_ds = yt.load('./snapshot_016.0') hc = HaloCatalog(data_ds=data_ds, finder_method='fof',finder_kwargs={"dm_only":False,"ptype":"Halo"}) hc.create()
and I met the error like this,
yt.utilities.exceptions.YTDomainOverflow: Particle bounds (-386.63174, -385.38947, -525.52057) and (6900.443, 6900.271, 13799.999) exceed domain bounds [0. 0. 0.] code_length and [13800. 13800. 13800.] code_length.
I have already updated yt and it's dependecy to the version 3.6.1, but I still have this problem. It seems like yt treat the particle velocity as the bound, I have no idea. What make me more confusing is there is no such error when I run above code for another small scale simualtion (886^3 particles and L_box=12Mpc/h), So, does anyone know what happend here or has ever met the same problem? _______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org https://mail.python.org/mailman3/lists/yt-users.python.org/ Member address: desika.narayanan@gmail.com
participants (2)
-
Desika Narayanan -
Zhang Meng