The site's hierarchy
Navigation
Mar 2009
| 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 |
- 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://
Grafic2
Multi zoom GRF
by Pichon Christophe
(Thursday 21 December 2006)
Grafic2
Written by Ed Bertschinger; modified by S. Prunet and C. Pichon
 grafic2 exemple
nested (four levels) GRF obeying a given powerspectrum.
|  zoom of grafic2 inner region
|
What it does ? Grafic-2 produces Gaussian random fields with multiple levels of resolution. It implements the theory of adaptive mesh refinement of Gaussian random fields.
How it does it do it? It relies on contrained realisations of the inner cubes given the tidal field imposed by the outer cube. See the paper
How to make it work ?
Install grafic2 from here
Compile grafic1 grafic2 (with ifort, not g77) at the desired resolutions: see script below as an example. Note that in order to produce a 
multi nested zoom, one need to compile
grafic1 with np1,2,3=128 and grafic2 with np1,2,3=64; see the README.
Save this script, say as four.csh and type
csh four.csh
Find below an example of script which produces a nested zoom with four levels of zoom which was used to produce the figure above.
In order to zoom on a specific region (rather than
the center of the box) one needs to load a shifted whitenoise
file.
The shift can be done by a simple roll or through a sub pixel
phaser (ie a filter in Fourier space 
).
An example of 4 level of zoom script which uses
baryon wiggles and WMAP 3 cosmology together with a
reference whitenoise is given at the bottom.
Find below (at the bottom) the yorick code which produces automatically such scripts, compiles grafic1,2 and produces the plot.
#! /bin/csh -f
set echo
# Script for grafic1 and grafic2
# WARNING !: make sure that grafic{1,2} have been *compiled* with np1=128, np1=64
cd test_res
cp $HOME/soft/grafic-2/linger.lcdm .
$HOME/soft/grafic-2/grafic1.128 >& test.out <<!
1 # uses linger.dat file for transfer functions
linger.lcdm
1 # spectral index n
-0.9 # normalization (sigma_8=-0.9)
0 0 # don't output power.dat
0.5 # Grid spacing (Mpc)
16 # Ultimate refinement factor (used only to set astart)
0 # further refinement (1 or 0 for yes or no)
64 64 64 # next grid size (used only if further refinement)
32 32 32 # next grid offsets (used only if further refinement)
64 64 64 # final subvolume size (used only if further refinement)
32 32 32 # final subvolume offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
39156 # iseed
wnsub.lev0
!
date >>& test.out
mkdir level_0
\mv -f ic_* level_0
#---------------------------------#
echo ' loop vertically 1'
#---------------------------------#
# Now run grafic1 again to produce intermediate level_1 output.
$HOME/soft/grafic-2/grafic1.128 >>& test.out <<!
1 # use linger.dat file for transfer functions
linger.lcdm
1 # spectral index n
-0.9 # normalization (sigma_8=-0.9)
0 0 # don't output power.dat
0.5 # Grid spacing (Mpc)
16 # Ultimate refinement factor (used only to set astart)
1 # further refinement (1 or 0 for yes or no)
64 64 64 # next grid size (used only if further refinement)
32 32 32 # next grid offsets (used only if further refinement)
64 64 64 # final subvolume size (used only if further refinement)
32 32 32 # final subvolume offsets (used only if further refinement)
2 # irand =1 to generate new nos, =2 use old ones
39156 # iseed
wnsub.lev0
!
date >>& test.out
\mv -f grafic2.top grafic2.lev1
#----
#---------------------------------#
echo ' loop horizontally 1 1'
#---------------------------------#
# Now run grafic2 to produce intermediate level_1-1 output.
$HOME/soft/grafic-2/grafic2.64 >>& test.out <<!
grafic2.lev1
1 # use linger.dat file for transfer functions
linger.lcdm
1 # spectral index n
-0.9 # normalization (sigma_8=-0.9)
0 0 # don't output power.dat
0 # further refinement (1 or 0 for yes or no)
64 64 64 # next grid size (used only if further refinement)
32 32 32 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
2079 # iseed
wnsub.lev1
!
date >>& test.out
# Save output in level_1 subdirectory
mkdir level_1
mv -f ic_* level_1
#---------------------------------#
echo ' loop vertically 2'
#---------------------------------#
# Now run grafic1 again to produce intermediate level_2 output.
$HOME/soft/grafic-2/grafic1.128 >>& test.out <<!
1 # use linger.dat file for transfer functions
linger.lcdm
1 # spectral index n
-0.9 # normalization (sigma_8=-0.9)
0 0 # don't output power.dat
0.5 # Grid spacing (Mpc)
16 # Ultimate refinement factor (used only to set astart)
1 # further refinement (1 or 0 for yes or no)
64 64 64 # next grid size (used only if further refinement)
32 32 32 # next grid offsets (used only if further refinement)
32 32 32 # final subvolume size (used only if further refinement)
48 48 48 # final subvolume offsets (used only if further refinement)
2 # irand =1 to generate new nos, =2 use old ones
39156 # iseed
wnsub.lev0
!
date >>& test.out
\mv -f grafic2.top grafic2.lev1
#----
#---------------------------------#
echo ' loop horizontally 2 1'
#---------------------------------#
# Now run grafic2 to produce intermediate level_2-1 output.
$HOME/soft/grafic-2/grafic2.64 >>& test.out <<!
grafic2.lev1
1 # use linger.dat file for transfer functions
linger.lcdm
1 # spectral index n
-0.9 # normalization (sigma_8=-0.9)
0 0 # don't output power.dat
1 # further refinement (1 or 0 for yes or no)
64 64 64 # next grid size (used only if further refinement)
32 32 32 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
2079 # iseed
wnsub.lev1
!
date >>& test.out
\mv -f grafic2.dat grafic2.lev2
#---------------------------------#
echo ' loop horizontally 2 2'
#---------------------------------#
# Now run grafic2 to produce intermediate level_2-2 output.
$HOME/soft/grafic-2/grafic2.64 >>& test.out <<!
grafic2.lev2
1 # use linger.dat file for transfer functions
linger.lcdm
1 # spectral index n
-0.9 # normalization (sigma_8=-0.9)
0 0 # don't output power.dat
0 # further refinement (1 or 0 for yes or no)
64 64 64 # next grid size (used only if further refinement)
32 32 32 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
35080 # iseed
wnsub.lev2
!
date >>& test.out
# Save output in level_2 subdirectory
mkdir level_2
mv -f ic_* level_2
#---------------------------------#
echo ' loop vertically 3'
#---------------------------------#
# Now run grafic1 again to produce intermediate level_3 output.
$HOME/soft/grafic-2/grafic1.128 >>& test.out <<!
1 # use linger.dat file for transfer functions
linger.lcdm
1 # spectral index n
-0.9 # normalization (sigma_8=-0.9)
0 0 # don't output power.dat
0.5 # Grid spacing (Mpc)
16 # Ultimate refinement factor (used only to set astart)
1 # further refinement (1 or 0 for yes or no)
64 64 64 # next grid size (used only if further refinement)
32 32 32 # next grid offsets (used only if further refinement)
16 16 16 # final subvolume size (used only if further refinement)
56 56 56 # final subvolume offsets (used only if further refinement)
2 # irand =1 to generate new nos, =2 use old ones
39156 # iseed
wnsub.lev0
!
date >>& test.out
\mv -f grafic2.top grafic2.lev1
#----
#---------------------------------#
echo ' loop horizontally 3 1'
#---------------------------------#
# Now run grafic2 to produce intermediate level_3-1 output.
$HOME/soft/grafic-2/grafic2.64 >>& test.out <<!
grafic2.lev1
1 # use linger.dat file for transfer functions
linger.lcdm
1 # spectral index n
-0.9 # normalization (sigma_8=-0.9)
0 0 # don't output power.dat
1 # further refinement (1 or 0 for yes or no)
64 64 64 # next grid size (used only if further refinement)
32 32 32 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
2079 # iseed
wnsub.lev1
!
date >>& test.out
\mv -f grafic2.dat grafic2.lev2
#---------------------------------#
echo ' loop horizontally 3 2'
#---------------------------------#
# Now run grafic2 to produce intermediate level_3-2 output.
$HOME/soft/grafic-2/grafic2.64 >>& test.out <<!
grafic2.lev2
1 # use linger.dat file for transfer functions
linger.lcdm
1 # spectral index n
-0.9 # normalization (sigma_8=-0.9)
0 0 # don't output power.dat
1 # further refinement (1 or 0 for yes or no)
64 64 64 # next grid size (used only if further refinement)
32 32 32 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
35080 # iseed
wnsub.lev2
!
date >>& test.out
\mv -f grafic2.dat grafic2.lev3
#---------------------------------#
echo ' loop horizontally 3 3'
#---------------------------------#
# Now run grafic2 to produce intermediate level_3-3 output.
$HOME/soft/grafic-2/grafic2.64 >>& test.out <<!
grafic2.lev3
1 # use linger.dat file for transfer functions
linger.lcdm
1 # spectral index n
-0.9 # normalization (sigma_8=-0.9)
0 0 # don't output power.dat
0 # further refinement (1 or 0 for yes or no)
64 64 64 # next grid size (used only if further refinement)
32 32 32 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
76812 # iseed
wnsub.lev3
!
date >>& test.out
# Save output in level_3 subdirectory
mkdir level_3
mv -f ic_* level_3
#---------------------------------#
echo ' loop vertically 4'
#---------------------------------#
# Now run grafic1 again to produce intermediate level_4 output.
$HOME/soft/grafic-2/grafic1.128 >>& test.out <<!
1 # use linger.dat file for transfer functions
linger.lcdm
1 # spectral index n
-0.9 # normalization (sigma_8=-0.9)
0 0 # don't output power.dat
0.5 # Grid spacing (Mpc)
16 # Ultimate refinement factor (used only to set astart)
1 # further refinement (1 or 0 for yes or no)
64 64 64 # next grid size (used only if further refinement)
32 32 32 # next grid offsets (used only if further refinement)
8 8 8 # final subvolume size (used only if further refinement)
60 60 60 # final subvolume offsets (used only if further refinement)
2 # irand =1 to generate new nos, =2 use old ones
39156 # iseed
wnsub.lev0
!
date >>& test.out
\mv -f grafic2.top grafic2.lev1
#----
#---------------------------------#
echo ' loop horizontally 4 1'
#---------------------------------#
# Now run grafic2 to produce intermediate level_4-1 output.
$HOME/soft/grafic-2/grafic2.64 >>& test.out <<!
grafic2.lev1
1 # use linger.dat file for transfer functions
linger.lcdm
1 # spectral index n
-0.9 # normalization (sigma_8=-0.9)
0 0 # don't output power.dat
1 # further refinement (1 or 0 for yes or no)
64 64 64 # next grid size (used only if further refinement)
32 32 32 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
2079 # iseed
wnsub.lev1
!
date >>& test.out
\mv -f grafic2.dat grafic2.lev2
#---------------------------------#
echo ' loop horizontally 4 2'
#---------------------------------#
# Now run grafic2 to produce intermediate level_4-2 output.
$HOME/soft/grafic-2/grafic2.64 >>& test.out <<!
grafic2.lev2
1 # use linger.dat file for transfer functions
linger.lcdm
1 # spectral index n
-0.9 # normalization (sigma_8=-0.9)
0 0 # don't output power.dat
1 # further refinement (1 or 0 for yes or no)
64 64 64 # next grid size (used only if further refinement)
32 32 32 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
35080 # iseed
wnsub.lev2
!
date >>& test.out
\mv -f grafic2.dat grafic2.lev3
#---------------------------------#
echo ' loop horizontally 4 3'
#---------------------------------#
# Now run grafic2 to produce intermediate level_4-3 output.
$HOME/soft/grafic-2/grafic2.64 >>& test.out <<!
grafic2.lev3
1 # use linger.dat file for transfer functions
linger.lcdm
1 # spectral index n
-0.9 # normalization (sigma_8=-0.9)
0 0 # don't output power.dat
1 # further refinement (1 or 0 for yes or no)
64 64 64 # next grid size (used only if further refinement)
32 32 32 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
76812 # iseed
wnsub.lev3
!
date >>& test.out
\mv -f grafic2.dat grafic2.lev4
#---------------------------------#
echo ' loop horizontally 4 4'
#---------------------------------#
# Now run grafic2 to produce intermediate level_4-4 output.
$HOME/soft/grafic-2/grafic2.64 >>& test.out <<!
grafic2.lev4
1 # use linger.dat file for transfer functions
linger.lcdm
1 # spectral index n
-0.9 # normalization (sigma_8=-0.9)
0 0 # don't output power.dat
0 # further refinement (1 or 0 for yes or no)
64 64 64 # next grid size (used only if further refinement)
32 32 32 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
48903 # iseed
wnsub.lev4
!
date >>& test.out
# Save output in level_4 subdirectory
mkdir level_4
mv -f ic_* level_4
Here is an example of 8 level of zoom which correspond to an
effective resolution of 
:
8 level grafic2 zoom
Here is an example of 4 level of zoom script which uses baryon wiggles and WMAP 3 cosmology together with a reference whitenoise (which may have been shifted).
#! /bin/csh -f
set echo
# Script for grafic1 and grafic2
# WARNING !: make sure that grafic2 has been *compiled* with np1=256
cd test_res
$HOME/soft/grafic-2/grafic1.256 >& test.out <<!
4 # uses baryon wiggles
0.042 # omegab
0.24 0.76 73 # cosmology Om Ov H0
0.958 # spectral index n
-0.77 # sigma8
0.01 1000 # output power.dat
-100 # boxsize
16 # Ultimate refinement factor (used only to set astart)
0 # further refinement (1 or 0 for yes or no)
128 128 128 # next grid size (used only if further refinement)
64 64 64 # next grid offsets (used only if further refinement)
128 128 128 # final subvolume size (used only if further refinement)
64 64 64 # final subvolume offsets (used only if further refinement)
2 # irand =1 to generate new nos, =2 use old ones
0 # iseed
../white.dat
!
date >>& test.out
mkdir level_0
\mv -f ic_* level_0
#---------------------------------#
echo ' loop vertically 1'
#---------------------------------#
# Now run grafic1 again to produce intermediate level_1 output.
$HOME/soft/grafic-2/grafic1.256 >>& test.out <<!
4 # uses baryon wiggles
0.042 # omegab
0.24 0.76 73 # cosmology Om Ov H0
0.958 # spectral index n
-0.77 # sigma8
0,0 # output power.dat
-100 # boxsize
16 # Ultimate refinement factor (used only to set astart)
1 # further refinement (1 or 0 for yes or no)
128 128 128 # next grid size (used only if further refinement)
64 64 64 # next grid offsets (used only if further refinement)
128 128 128 # final subvolume size (used only if further refinement)
64 64 64 # final subvolume offsets (used only if further refinement)
2 # irand =1 to generate new nos, =2 use old ones
91632 # iseed
../white.dat
!
date >>& test.out
\mv -f grafic2.top grafic2.lev1
#----
#---------------------------------#
echo ' loop horizontally 1 1'
#---------------------------------#
# Now run grafic2 to produce intermediate level_1-1 output.
$HOME/soft/grafic-2/grafic2.128 >>& test.out <<!
grafic2.lev1
4 # uses baryon wiggles
0.042 # omegab
0.958 # spectral index n
-0.77 # sigma8
0.01,1000 # output power.dat
0 # further refinement (1 or 0 for yes or no)
128 128 128 # next grid size (used only if further refinement)
64 64 64 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
9785 # iseed
wnsub.lev1
!
date >>& test.out
# Save output in level_1 subdirectory
mkdir level_1
mv -f ic_* level_1
#---------------------------------#
echo ' loop vertically 2'
#---------------------------------#
# Now run grafic1 again to produce intermediate level_2 output.
$HOME/soft/grafic-2/grafic1.256 >>& test.out <<!
4 # uses baryon wiggles
0.042 # omegab
0.24 0.76 73 # cosmology Om Ov H0
0.958 # spectral index n
-0.77 # sigma8
0,0 # output power.dat
-100 # boxsize
16 # Ultimate refinement factor (used only to set astart)
1 # further refinement (1 or 0 for yes or no)
128 128 128 # next grid size (used only if further refinement)
64 64 64 # next grid offsets (used only if further refinement)
64 64 64 # final subvolume size (used only if further refinement)
96 96 96 # final subvolume offsets (used only if further refinement)
2 # irand =1 to generate new nos, =2 use old ones
91632 # iseed
../white.dat
!
date >>& test.out
\mv -f grafic2.top grafic2.lev1
#----
#---------------------------------#
echo ' loop horizontally 2 1'
#---------------------------------#
# Now run grafic2 to produce intermediate level_2-1 output.
$HOME/soft/grafic-2/grafic2.128 >>& test.out <<!
grafic2.lev1
4 # uses baryon wiggles
0.042 # omegab
0.958 # spectral index n
-0.77 # sigma8
0.01,1000 # output power.dat
1 # further refinement (1 or 0 for yes or no)
128 128 128 # next grid size (used only if further refinement)
64 64 64 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
9785 # iseed
wnsub.lev1
!
date >>& test.out
\mv -f grafic2.dat grafic2.lev2
#---------------------------------#
echo ' loop horizontally 2 2'
#---------------------------------#
# Now run grafic2 to produce intermediate level_2-2 output.
$HOME/soft/grafic-2/grafic2.128 >>& test.out <<!
grafic2.lev2
4 # uses baryon wiggles
0.042 # omegab
0.958 # spectral index n
-0.77 # sigma8
0.01,1000 # output power.dat
0 # further refinement (1 or 0 for yes or no)
128 128 128 # next grid size (used only if further refinement)
64 64 64 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
71371 # iseed
wnsub.lev2
!
date >>& test.out
# Save output in level_2 subdirectory
mkdir level_2
mv -f ic_* level_2
#---------------------------------#
echo ' loop vertically 3'
#---------------------------------#
# Now run grafic1 again to produce intermediate level_3 output.
$HOME/soft/grafic-2/grafic1.256 >>& test.out <<!
4 # uses baryon wiggles
0.042 # omegab
0.24 0.76 73 # cosmology Om Ov H0
0.958 # spectral index n
-0.77 # sigma8
0,0 # output power.dat
-100 # boxsize
16 # Ultimate refinement factor (used only to set astart)
1 # further refinement (1 or 0 for yes or no)
128 128 128 # next grid size (used only if further refinement)
64 64 64 # next grid offsets (used only if further refinement)
32 32 32 # final subvolume size (used only if further refinement)
112 112 112 # final subvolume offsets (used only if further refinement)
2 # irand =1 to generate new nos, =2 use old ones
91632 # iseed
../white.dat
!
date >>& test.out
\mv -f grafic2.top grafic2.lev1
#----
#---------------------------------#
echo ' loop horizontally 3 1'
#---------------------------------#
# Now run grafic2 to produce intermediate level_3-1 output.
$HOME/soft/grafic-2/grafic2.128 >>& test.out <<!
grafic2.lev1
4 # uses baryon wiggles
0.042 # omegab
0.958 # spectral index n
-0.77 # sigma8
0.01,1000 # output power.dat
1 # further refinement (1 or 0 for yes or no)
128 128 128 # next grid size (used only if further refinement)
64 64 64 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
9785 # iseed
wnsub.lev1
!
date >>& test.out
\mv -f grafic2.dat grafic2.lev2
#---------------------------------#
echo ' loop horizontally 3 2'
#---------------------------------#
# Now run grafic2 to produce intermediate level_3-2 output.
$HOME/soft/grafic-2/grafic2.128 >>& test.out <<!
grafic2.lev2
4 # uses baryon wiggles
0.042 # omegab
0.958 # spectral index n
-0.77 # sigma8
0.01,1000 # output power.dat
1 # further refinement (1 or 0 for yes or no)
128 128 128 # next grid size (used only if further refinement)
64 64 64 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
71371 # iseed
wnsub.lev2
!
date >>& test.out
\mv -f grafic2.dat grafic2.lev3
#---------------------------------#
echo ' loop horizontally 3 3'
#---------------------------------#
# Now run grafic2 to produce intermediate level_3-3 output.
$HOME/soft/grafic-2/grafic2.128 >>& test.out <<!
grafic2.lev3
4 # uses baryon wiggles
0.042 # omegab
0.958 # spectral index n
-0.77 # sigma8
0.01,1000 # output power.dat
0 # further refinement (1 or 0 for yes or no)
128 128 128 # next grid size (used only if further refinement)
64 64 64 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
46979 # iseed
wnsub.lev3
!
date >>& test.out
# Save output in level_3 subdirectory
mkdir level_3
mv -f ic_* level_3
#---------------------------------#
echo ' loop vertically 4'
#---------------------------------#
# Now run grafic1 again to produce intermediate level_4 output.
$HOME/soft/grafic-2/grafic1.256 >>& test.out <<!
4 # uses baryon wiggles
0.042 # omegab
0.24 0.76 73 # cosmology Om Ov H0
0.958 # spectral index n
-0.77 # sigma8
0,0 # output power.dat
-100 # boxsize
16 # Ultimate refinement factor (used only to set astart)
1 # further refinement (1 or 0 for yes or no)
128 128 128 # next grid size (used only if further refinement)
64 64 64 # next grid offsets (used only if further refinement)
16 16 16 # final subvolume size (used only if further refinement)
120 120 120 # final subvolume offsets (used only if further refinement)
2 # irand =1 to generate new nos, =2 use old ones
91632 # iseed
../white.dat
!
date >>& test.out
\mv -f grafic2.top grafic2.lev1
#----
#---------------------------------#
echo ' loop horizontally 4 1'
#---------------------------------#
# Now run grafic2 to produce intermediate level_4-1 output.
$HOME/soft/grafic-2/grafic2.128 >>& test.out <<!
grafic2.lev1
4 # uses baryon wiggles
0.042 # omegab
0.958 # spectral index n
-0.77 # sigma8
0.01,1000 # output power.dat
1 # further refinement (1 or 0 for yes or no)
128 128 128 # next grid size (used only if further refinement)
64 64 64 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
9785 # iseed
wnsub.lev1
!
date >>& test.out
\mv -f grafic2.dat grafic2.lev2
#---------------------------------#
echo ' loop horizontally 4 2'
#---------------------------------#
# Now run grafic2 to produce intermediate level_4-2 output.
$HOME/soft/grafic-2/grafic2.128 >>& test.out <<!
grafic2.lev2
4 # uses baryon wiggles
0.042 # omegab
0.958 # spectral index n
-0.77 # sigma8
0.01,1000 # output power.dat
1 # further refinement (1 or 0 for yes or no)
128 128 128 # next grid size (used only if further refinement)
64 64 64 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
71371 # iseed
wnsub.lev2
!
date >>& test.out
\mv -f grafic2.dat grafic2.lev3
#---------------------------------#
echo ' loop horizontally 4 3'
#---------------------------------#
# Now run grafic2 to produce intermediate level_4-3 output.
$HOME/soft/grafic-2/grafic2.128 >>& test.out <<!
grafic2.lev3
4 # uses baryon wiggles
0.042 # omegab
0.958 # spectral index n
-0.77 # sigma8
0.01,1000 # output power.dat
1 # further refinement (1 or 0 for yes or no)
128 128 128 # next grid size (used only if further refinement)
64 64 64 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
46979 # iseed
wnsub.lev3
!
date >>& test.out
\mv -f grafic2.dat grafic2.lev4
#---------------------------------#
echo ' loop horizontally 4 4'
#---------------------------------#
# Now run grafic2 to produce intermediate level_4-4 output.
$HOME/soft/grafic-2/grafic2.128 >>& test.out <<!
grafic2.lev4
4 # uses baryon wiggles
0.042 # omegab
0.958 # spectral index n
-0.77 # sigma8
0.01,1000 # output power.dat
0 # further refinement (1 or 0 for yes or no)
128 128 128 # next grid size (used only if further refinement)
64 64 64 # next grid offsets (used only if further refinement)
1 # irand =1 to generate new nos, =2 use old ones
81864 # iseed
wnsub.lev4
!
date >>& test.out
# Save output in level_4 subdirectory
mkdir level_4
mv -f ic_* level_4
