Hi,
I've been trying to construct a CiC mesh but couldn't get an answer comparable with theory for the mean density of the mesh yt gave me.
I've put all the relevant code and output in a pastebin here: http://pastebin.com/hbVM0rYN
I've spoken with Britton Smith and he thinks the physical size of the mesh might be wrong by a factor of (1+z) which translates to the densities all being off by (1+z)^3. There also seems to be a missing h**3 somewhere which brings everything in "rough" alignment.
I've run it for two redshifts z = 46 and z = 0 and put the output in the above pastebin. I've added a correction where I forcefully put /(1+z)^3 term at the end and it gives reasonable results. I've also added another two values including the missing /h**3 term.
Can this be corrected or can someone tell me where I went wrong in what I'm doing? Thanks!
Brendan
Hi all,
I am trying to debug this, but it would help if I had a small cosmological Gadget dataset at a redshift greater (significantly if possible) than 0. Does anyone have anything like this they could share?
Britton
On Mon, Aug 4, 2014 at 4:17 PM, Brendan Griffen <brendan.f.griffen@gmail.com
wrote:
Hi,
I've been trying to construct a CiC mesh but couldn't get an answer comparable with theory for the mean density of the mesh yt gave me.
I've put all the relevant code and output in a pastebin here: http://pastebin.com/hbVM0rYN
I've spoken with Britton Smith and he thinks the physical size of the mesh might be wrong by a factor of (1+z) which translates to the densities all being off by (1+z)^3. There also seems to be a missing h**3 somewhere which brings everything in "rough" alignment.
I've run it for two redshifts z = 46 and z = 0 and put the output in the above pastebin. I've added a correction where I forcefully put /(1+z)^3 term at the end and it gives reasonable results. I've also added another two values including the missing /h**3 term.
Can this be corrected or can someone tell me where I went wrong in what I'm doing? Thanks!
Brendan
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Hi Britton,
Yes, I have some. Look on galactica.pa.msu.edu in /data/gomez/halo1_128 - there are data outputs at many different redshifts. Each output is 67 MB, 128^3 particles (or thereabouts).
--Brian
On Mon, Aug 4, 2014 at 12:07 PM, Britton Smith brittonsmith@gmail.com wrote:
Hi all,
I am trying to debug this, but it would help if I had a small cosmological Gadget dataset at a redshift greater (significantly if possible) than 0. Does anyone have anything like this they could share?
Britton
On Mon, Aug 4, 2014 at 4:17 PM, Brendan Griffen < brendan.f.griffen@gmail.com> wrote:
Hi,
I've been trying to construct a CiC mesh but couldn't get an answer comparable with theory for the mean density of the mesh yt gave me.
I've put all the relevant code and output in a pastebin here: http://pastebin.com/hbVM0rYN
I've spoken with Britton Smith and he thinks the physical size of the mesh might be wrong by a factor of (1+z) which translates to the densities all being off by (1+z)^3. There also seems to be a missing h**3 somewhere which brings everything in "rough" alignment.
I've run it for two redshifts z = 46 and z = 0 and put the output in the above pastebin. I've added a correction where I forcefully put /(1+z)^3 term at the end and it gives reasonable results. I've also added another two values including the missing /h**3 term.
Can this be corrected or can someone tell me where I went wrong in what I'm doing? Thanks!
Brendan
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
-
Brendan Griffen -
Brian O'Shea -
Britton Smith