banner

The site's hierarchy



Home page :: Private area :: Site Map :: Credits
Science / Conditions Initiales / Conditions initiales cosmologiques / Codes / Constrfield / index

Navigation

Mar 2009
MTWTFSS
2324252627281
2345678
9101112131415
16171819202122
23242526272829
303112345
News
  • Les premiers zoom de galaxies sont accessibles
  • Le halo FOF 6133 de la boite Horizon L et le halo FOF 544 de Horizon S ont été resimulés avec plusieurs techniques de zoom. Il est important que chaque "zoomer" valide sa méthode, avant de lancer une campagne de zoom sur un plus grand cataloque.
  • Les données de la simulation Mare Nostrum sont disponibles
  • 34 snapshots jusqu’à z=4 sont accessibles sur horizon3 et sur le serveur de fichiers de l’IDRIS à la collaboration Horizon.
  • Méso Machine HPC1 opérationnelle
  • Depuis le 23 Octobre 2005, la Méso machine du site horizon est operationnelle. Elle correspond à 3 quadriprocesseurs avec chacun 64 Giga de RAM reliés par infiniband, ainsi qu’un access conséquent (sur une base de projet dédié) au reste de la ferme). Son acces est ouvert a toute personne de la collaboration ayant acces à la minigrille et qui en fait la demande a admin-minigrille
  • http://
  • Workshop Horizon le 14 et 15 novembre 2005
  • Il aura lieu à Paris les 14 et 15 novembre 2005 (prévoir une nuit sur place). L’enregistrement est ouvert dans la rubrique "meeting!"
  • http://

index


by Pichon Christophe (Saturday 30 December 2006)

Constrfield

written by D. Pogosyan: implementation of baryon wiggles and F90 for HORIZON.

JPG - 1.9 Mb
constrained realisation
produced with a constrfield and mpgrafic

JPG - 52.6 kb
Power spectrum of CGRF
using mpgrafic and constfield. Here the constraints correspond to a regular grid.

An example of script to run constrfield is


#!/bin/csh -f
#Output size, smoothing, file
set Nf=g
set Rf=0.
set boxsize=142.857
set denfile=out.dat

# do prewhiting ?

set prewhite=y

#Number of constraint points
set Npoints=50

# cosmology & powerspectrum (Ed Bertshinger style)
set cosmo='b'
set powerspectype=4
set Omegam=0.3
set Omegav=0.7
set Omegab=0.05
set H0=70
set primordialspectrum=1
set sigma8=-0.92

# If random field is added, one need to include the following
# two parameters after mean_no_random
set random_seed=12345

set exec_dir="$HOME/soft/shear/F90/"

$exec_dir/constr_field << EOF
$random_seed
$boxsize
$Nf
$Rf
$prewhite
$cosmo
$powerspectype
$Omegam,$Omegav,$H0
$Omegab
$primordialspectrum
$sigma8
0,0
$denfile
$Npoints
n
g    # point number 01
5.000000
9.000000
9.000000
9.000000
y
25.000000
y
0.864533
0.700956
0.163927
n
n
n
y
g    # point number 02
5.000000
45.000000
45.000000
45.000000
y
25.000000
y
0.514600
0.490192
0.429319
n
n
n
y
EOF

An alternative type of constraint, corresponding to top hat with velocities would be


t    # point number 01
6.832125
96.876959
109.499716
106.630902
y
1.686000
y
30.031500
-159.169000
253.191000
n
e
n
y
t    # point number 02
6.832125
73.352816
116.384816
78.229673
y
1.686000
y
-428.348000
446.312000
35.606500
n
e
n
y