== Single Halo ==

For this test, we take the grid from the [wiki:VOBOZ/NullResult null result test], and add a blob around the center of one of the tiles. The code used to create the halo files is in [attachment:single_halo.c single_halo.c]. Again, you'll need to modify {{{readfiles.c}}} to read this and the velocity field file in.

I'm going to skip the diatribe about what the programs do, go back to the [wiki:VOBOZ/NullResult null result test] if things are unclear.

=== Create Grid ===

{{{
cable:~/packages/voboz1.1.1 rpwagner$ gcc single_halo.c
cable:~/packages/voboz1.1.1 rpwagner$ a.out
cable:~/packages/voboz1.1.1 rpwagner$ ls -l single_halo*
-rw-r--r--   1 rpwagner  rpwagner  1132 Jan 14 13:52 single_halo.raw
-rw-r--r--   1 rpwagner  rpwagner  1132 Jan 14 13:52 single_halo_vel.raw
cable:~/packages/voboz1.1.1 rpwagner$ 
}}}

=== {{{vozinit}}} ===

The blob (halo) is inside of the {{{8}}} particles in the tile, so it should be well blocked from the guard particles.

{{{
cable:~/packages/voboz1.1.1 rpwagner$ ./vozinit single_halo.raw 0.3 1 2 goo
}}}

{{{
cable:~/packages/voboz1.1.1 rpwagner$ cat scrgoo 
#!/bin/bash -f
./voz1b1 single_halo.raw 0.300000 1.000000 goo 2 0 0 0
./voz1b1 single_halo.raw 0.300000 1.000000 goo 2 0 0 1
./voz1b1 single_halo.raw 0.300000 1.000000 goo 2 0 1 0
./voz1b1 single_halo.raw 0.300000 1.000000 goo 2 0 1 1
./voz1b1 single_halo.raw 0.300000 1.000000 goo 2 1 0 0
./voz1b1 single_halo.raw 0.300000 1.000000 goo 2 1 0 1
./voz1b1 single_halo.raw 0.300000 1.000000 goo 2 1 1 0
./voz1b1 single_halo.raw 0.300000 1.000000 goo 2 1 1 1
./voztie 2 goo
cable:~/packages/voboz1.1.1 rpwagner$ 
}}}


=== {{{voz1b1}}} and {{{voztie}}} ===

{{{
cable:~/packages/voboz1.1.1 rpwagner$ ./scrgoo 
<< blah, blah, blah >>
}}}

We've got the results from each tile,

{{{
cable:~/packages/voboz1.1.1 rpwagner$ ls -1 part.goo.0*
part.goo.00.00.00
part.goo.00.00.01
part.goo.00.01.00
part.goo.00.01.01
part.goo.01.00.00
part.goo.01.00.01
part.goo.01.01.00
part.goo.01.01.01
cable:~/packages/voboz1.1.1 rpwagner$ 
}}}

and the data files from {{{voztie}}}.

{{{
cable:~/packages/voboz1.1.1 rpwagner$ ls -1 *goo.dat
adjgoo.dat
volgoo.dat
cable:~/packages/voboz1.1.1 rpwagner$ 
}}}

=== {{{jovoz}}} ===

{{{
cable:~/packages/voboz1.1.1 rpwagner$ ./jovoz adjgoo.dat volgoo.dat goozone.dat goozone.txt 0.01
About to jump ...
Post-jump ...
1 initial haloes found
Found zone adjacencies
Volumes range from 1.401526e-02 to 1.468750e+00.
Maxprob = 104.796516.
about to sort ...
cable:~/packages/voboz1.1.1 rpwagner$ 
}}}

Again, some files.

{{{
cable:~/packages/voboz1.1.1 rpwagner$ ls -1 goozone.*
goozone.dat
goozone.txt
}}}

And now for something new. The new clump creates a zone with a different volume.

{{{
cable:~/packages/voboz1.1.1 rpwagner$ head goozone.txt 
1
0       94      1.047965e+02    1.468750e+00
cable:~/packages/voboz1.1.1 rpwagner$}}}

=== {{{boz}}} ===

{{{boz}}} gives us even more new stuff, like finding a single halo.

{{{
cable:~/packages/voboz1.1.1 rpwagner$ ./boz 1.0 0.3 1.0 single_halo.raw single_halo_vel.raw goozone.dat goohalos.dat goohalos.txt 1.2
a = 1.000000 => Multiplying H0 by 1.000000e+00
cell, potconst = 0.007812,1.707547e-04
np = 94
np: 94, x: 0.000000,0.875000; y: 0.000000,0.875000; z: 0.000000,0.875000
vx: 0.000000,0.000000; vy: 0.000000,0.000000; vz: 0.000000,0.000000
np = 94, nhaloes = 1
Done! Outputting ...
Nhunb = 1
cable:~/packages/voboz1.1.1 rpwagner$
}}}

Here's the info on the halo found. Interestingly, all of the particles are included.

{{{
cable:~/packages/voboz1.1.1 rpwagner$ cat goohalos.txt 
1       1
0       94      94      1.047965e+02 1.401526e-02 -inf -inf     70      0       0       1 1 0.000000 0.000000 0.000000 0.000000e+00 0.000000e+00 0.000000e+00
cable:~/packages/voboz1.1.1 rpwagner$ 
}}}