The scripts included in this package run MUSIC after setting the
necessary parameters to their correct values to focus on a Lagrangian
volume of some halo.


$ python <level>

<level> is the finest level in the initial conditions to be generated.


1. Collect parameters from the original set of initial conditions
   (ICs) to obtain the shift of the box.

2. Given a halo position and mass, it will calculate the Lagrangian
   volume (LV) at the initial redshift.

3. With this LV, it will modify the appropriate parameters and supply
   a particle position list to MUSIC.  Afterwards, MUSIC is
   automatically executed.

   * Four options for the LV are given: box, ellipsoid, convex hull,
     or exact.

4. The initial conditions directory is then moved to some
   user-specified location.

Example workflow

1. Run a uniform grid simulation
2. Find halo of interest in the last output
3. Modify halo_info and other parameters in
   The parameters are documented above their declaration.
4. Run "python 1"
5. Run simulation with 1 nested grid with the appropriate Enzo
   parameter file.  MUSIC will only give a skeleton parameter file,
   but the must-refine and nested grid parameters are given in it.
6. If another nested grid is needed, goto step #2 and increase the
   level argument by 1.  If not, you're finished!