[yt-users] units error in cloud-in-cell construction in yt

[Britton Smith]

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 at 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 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/20140804/0ae84686/attachment.html>