- Mail actions: [ respond to this message ] [ mail a new topic ]
- Contemporary messages sorted: [ by date ] [ by thread ] [ by subject ] [ by author ] [ by messages with attachments ]

From: Jared Coughlin <Jared.W.Coughlin.29_at_nd.edu>

Date: Tue, 7 Apr 2015 16:42:55 -0400

Hello! I recently downloaded the sph version of the Aquila project initial

conditions from here:

http://www.aip.de/People/cscannapieco/aquila/ICs/Gadget/

The note says that the ics are in gadget format with units of 10^10 M_sun

/h, Mpc/h, and km/s. I obtained the cosmological parameters from their

paper:

http://arxiv.org/pdf/1112.0315v3.pdf

I tried running the Aq-C-6 ics but encountered several problems. First,

gadget told me that it detected incorrect block sizes. I have since

written a code that fixed this and also checked to make sure that the data

itself was untouched by this procedure.

The second problem, however, is that gadget, in the function check_omega in

init.c, is obtaining a value of omega = 2.6e13, when it is supposed to be

.25. I'm fairly certain that this is a units problem, and I was just

wondering if anyone had worked with these ics before?

For cosmology I have:

Omega0 = .25

OmegaLambda = .75

OmegaBaryone = 0.04

HubbleParam = 0.73

BoxSize = 100.0 (This should be in Mpc / h, which is the

unit of length that the Aquila project uses. I got this from Springel et

al. 2008)

For the system of units I have:

UnitLength_in_cm = 3.085678e24

UnitMass_in_g = 1.989e43

UnitVelocity_in_cm_per_s = 1e5

I'm pretty sure this system matches the note on the Aquila project site,

but I could be wrong. Anyways, any help would be greatly appreciated! Thank

you very much!

Sincerely,

-Jared

Received on 2015-04-07 22:43:00

Date: Tue, 7 Apr 2015 16:42:55 -0400

Hello! I recently downloaded the sph version of the Aquila project initial

conditions from here:

http://www.aip.de/People/cscannapieco/aquila/ICs/Gadget/

The note says that the ics are in gadget format with units of 10^10 M_sun

/h, Mpc/h, and km/s. I obtained the cosmological parameters from their

paper:

http://arxiv.org/pdf/1112.0315v3.pdf

I tried running the Aq-C-6 ics but encountered several problems. First,

gadget told me that it detected incorrect block sizes. I have since

written a code that fixed this and also checked to make sure that the data

itself was untouched by this procedure.

The second problem, however, is that gadget, in the function check_omega in

init.c, is obtaining a value of omega = 2.6e13, when it is supposed to be

.25. I'm fairly certain that this is a units problem, and I was just

wondering if anyone had worked with these ics before?

For cosmology I have:

Omega0 = .25

OmegaLambda = .75

OmegaBaryone = 0.04

HubbleParam = 0.73

BoxSize = 100.0 (This should be in Mpc / h, which is the

unit of length that the Aquila project uses. I got this from Springel et

al. 2008)

For the system of units I have:

UnitLength_in_cm = 3.085678e24

UnitMass_in_g = 1.989e43

UnitVelocity_in_cm_per_s = 1e5

I'm pretty sure this system matches the note on the Aquila project site,

but I could be wrong. Anyways, any help would be greatly appreciated! Thank

you very much!

Sincerely,

-Jared

Received on 2015-04-07 22:43:00

*
This archive was generated by hypermail 2.3.0
: 2022-09-01 14:03:42 CEST
*