SlicePlot of dark matter density fields from an enzo sim with ParallelParticleIO=1.
Hi, Want to touch base in analyzing cosmological simulation results with ParallelParticleIO on. I'm in the middle of testing my own initial condition generator (named BCCOMICS). Starting from a unigrid IC from BCCOMICS, I ran two cosmological (with cooling on) simulations, one with AMR on and one with AMR off. For the AMR-on one, I turned on ParallelParticleIO and this figure ( https://www.chosun.ac.kr/user/kjahn/tmpFile/RD0001_Slice_y_Dark_Matter_Densi... ) is the plot I got by yting as follows: ------- import yt fname="RD0001/RD0001" ds=yt.load(fname) pc=yt.SlicePlot(ds,"y","Dark_Matter_Density") pc.annotate_quiver("all_cic_velocity_z","all_cic_velocity_x") pc.save() -------- The ring command I used is: ring -np 20 pv ParticlePositions ParticleVelocities So I see those rectangular artifacts in the plot. I don't see any artifact in the baryon density field (ParallelRootGridIO=1 too). Q1: I guess that enzo is internally getting the dark matter particle information and gravity correctly? In the AMR-off simulation I did not turn on ParallelParticleIO, and surely the DM field plotted by the same yt script looks fine. Q2: So, is this the right way you analyze the DM density field when ParallelParticleIO was on? Many thanks in advance. Kyungjin
This is probably a better question for the Enzo-users or Enzo-dev mailing lists. It may also be worth making a dark matter particle filter and then plotting the deposited dark matter density in yt using that filter. The Dark_Matter_Density field is computed internally in Enzo and yt is just displaying what Enzo has written to disk. -Nathan On Fri, Mar 9, 2018 at 2:03 AM, arrau08 . <kjyaya@gmail.com> wrote:
Hi,
Want to touch base in analyzing cosmological simulation results with ParallelParticleIO on.
I'm in the middle of testing my own initial condition generator (named BCCOMICS). Starting from a unigrid IC from BCCOMICS, I ran two cosmological (with cooling on) simulations, one with AMR on and one with AMR off. For the AMR-on one, I turned on ParallelParticleIO and this figure ( https://www.chosun.ac.kr/user/kjahn/tmpFile/RD0001_Slice_y_D ark_Matter_Density.png ) is the plot I got by yting as follows: ------- import yt
fname="RD0001/RD0001" ds=yt.load(fname)
pc=yt.SlicePlot(ds,"y","Dark_Matter_Density") pc.annotate_quiver("all_cic_velocity_z","all_cic_velocity_x") pc.save() --------
The ring command I used is: ring -np 20 pv ParticlePositions ParticleVelocities
So I see those rectangular artifacts in the plot. I don't see any artifact in the baryon density field (ParallelRootGridIO=1 too).
Q1: I guess that enzo is internally getting the dark matter particle information and gravity correctly? In the AMR-off simulation I did not turn on ParallelParticleIO, and surely the DM field plotted by the same yt script looks fine.
Q2: So, is this the right way you analyze the DM density field when ParallelParticleIO was on?
Many thanks in advance.
Kyungjin _______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org _______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org
Thank you Nathan, I thought it was more of a yt related question. Anyway your answer leads to the DM filter, which I'll try to learn and apply. Many thanks! Kyungjin On Sat, Mar 10, 2018 at 12:29 AM, Nathan Goldbaum <nathan12343@gmail.com> wrote:
This is probably a better question for the Enzo-users or Enzo-dev mailing lists.
It may also be worth making a dark matter particle filter and then plotting the deposited dark matter density in yt using that filter. The Dark_Matter_Density field is computed internally in Enzo and yt is just displaying what Enzo has written to disk.
-Nathan
On Fri, Mar 9, 2018 at 2:03 AM, arrau08 . <kjyaya@gmail.com> wrote:
Hi,
Want to touch base in analyzing cosmological simulation results with ParallelParticleIO on.
I'm in the middle of testing my own initial condition generator (named BCCOMICS). Starting from a unigrid IC from BCCOMICS, I ran two cosmological (with cooling on) simulations, one with AMR on and one with AMR off. For the AMR-on one, I turned on ParallelParticleIO and this figure ( https://www.chosun.ac.kr/user/kjahn/tmpFile/RD0001_Slice_y_D ark_Matter_Density.png ) is the plot I got by yting as follows: ------- import yt
fname="RD0001/RD0001" ds=yt.load(fname)
pc=yt.SlicePlot(ds,"y","Dark_Matter_Density") pc.annotate_quiver("all_cic_velocity_z","all_cic_velocity_x") pc.save() --------
The ring command I used is: ring -np 20 pv ParticlePositions ParticleVelocities
So I see those rectangular artifacts in the plot. I don't see any artifact in the baryon density field (ParallelRootGridIO=1 too).
Q1: I guess that enzo is internally getting the dark matter particle information and gravity correctly? In the AMR-off simulation I did not turn on ParallelParticleIO, and surely the DM field plotted by the same yt script looks fine.
Q2: So, is this the right way you analyze the DM density field when ParallelParticleIO was on?
Many thanks in advance.
Kyungjin _______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org _______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org
_______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org
_______________________________________________ yt-users mailing list -- yt-users@python.org To unsubscribe send an email to yt-users-leave@python.org
-- Kyungjin Ahn Associate Professor Department of Earth Sciences, Astronomy Division Chosun University Tel: 062-230-7340 +82-62-230-7340 (from abroad) homepage: www.chosun.ac.kr/kjahn
participants (3)
-
arrau08 . -
Kyungjin Ahn -
Nathan Goldbaum