M | T | W | T | F | S | S |
---|---|---|---|---|---|---|
23 | 24 | 25 | 26 | 27 | 28 | 1 |
2 | 3 | 4 | 5 | 6 | 7 | 8 |
9 | 10 | 11 | 12 | 13 | 14 | 15 |
16 | 17 | 18 | 19 | 20 | 21 | 22 |
23 | 24 | 25 | 26 | 27 | 28 | 29 |
30 | 31 | 1 | 2 | 3 | 4 | 5 |
written by D. Pogosyan: implementation of baryon wiggles and F90 for HORIZON.
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