Description
===========
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.

Usage
=====

$ python enzo-mrp-music.py <config-file> <level>

<config-file> is the configuration file. See sample.conf for an
example.

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

This script can be run in parallel that will use multiple cores for
the particle identification and searching in yt.

$ mpirun -n <number-of-cores> python enzo-mrp.music.py <level>

Procedure
=========

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. Create a set of initial conditions for a unigrid simulation with
   MUSIC. Its parameter file needs to be references in this script's
   configuration file (e.g. sample.conf) as the variable
   template_config.  See template.conf for an example.
   
2. Run unigrid simulation (If an iteration after step #7, run zoom-in
   simulation.)

3. Find halo of interest in the last output

4. Modify parameters in the [region] section in the configuration file
   (e.g. sample.conf). The parameters are documented above their
   declaration.
   
5. Run "python enzo-mrp-music.py sample.conf 1"

6. 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.
   
7. If another nested grid is needed, goto step #2 and increase the
   level argument by 1.  If not, you're finished!