Bandeau-IAP IAP Ohio CRAL
accueil 

Mocks

Ramses simulation ("Horizon-MareNostrum"):
WMAP1z=2.51 Lbox=50 Mpc/h10243 DM particles+hydro

Post-treatment: hydro spectra at the BOSS resolution:







Download binary or fits files:
GridsSlices
Real Space1024x256x256 (~256 M)       .bin 1024x1024 (~4 M)      .bin      .fits
Redshift Space1024x256x256 (~256 M)       .bin1024x1024 (~4 M)      .bin      .fits
The code read_grid.c can be used to read the binary files.

Structure of binary files:

The binary files are organized as follows (blocks are delimited by dummy variables indicating the size of the blocks for FORTRAN compatibility)

fieldtypesizecomment
dummyint(4B) 1= 28
dim1int(4B)1 resolution of spectra (direction of the line of sight)
dim2int(4B)1 number of spectra per slice
dim3int(4B)1 number of slices
Lboxfloat(4B)1 simulation Box side
zfloat(4B)1 redshift
z-spaceint(4B)10: real-space    1:redshift-space
F_smoothingfloat(4B)1 Flux smoothing scale in Mpc/h (here BOSS scale)
dummyint(4B) 1= 28
dummyint(4B) 1= dim1*dim2*4
Slice 1float(4B)dim1*dim21...dim1:                  trans. flux 1
1+dim1...2dim1:     trans. flux 2
                        .
                        .
dummyint(4B) 1= dim1*dim2*4
dummyint(4B) 1= dim1*dim2*4
Slice 2float(4B)dim1*dim21...dim1:                  trans. flux dim1*dim2+1
1+dim1...2dim1:     trans. flux dim1*dim2+2
                        .
                        .
dummyint(4B) 1= dim1*dim2*4
etc...
dummyint(4B) 1= dim1*dim2*4
Slice dim3float(4B)dim1*dim21...dim1:              trans. flux dim1*dim2*(dim3-1)+1
1+dim1...2dim1: trans. flux dim1*dim2*(dim3-1)+2
                        .
                        .
dummyint(4B) 1= dim1*dim2*4

Structure of fits files (for slice only):
SIMPLE= T / true FITS file
BITPIX= -32 / bits per pixel
NAXIS= 2 / number of dimensions
NAXIS1= 1024 / length of 1st dimension: Fs
NAXIS2= 1024 / length of 2nd dimension
LBOX= 50 / Simulation box size (Mpc/h)
z= 2.51 / Redshift
F_SMOOTH= 0.696 / Flux smoothing scale in Mpc/h (here BOSS scale)

TAB = array(double, NAXIS1, NAXIS2)

ex: TAB(456,) = trans.flux 456

Institut d'Astrophysique de Paris - 98 bis boulevard Arago - 75014 Paris  |  Contact: peirani AT iap.fr