Following issue #30 (closed), I have added a new script to generate arbitrarily large initial conditions.
The scripts downloads an initial tile and then replicates it along each axis multiple times and/or downsamples the initial particle load. The script takes two arguments:
python makeIC.py factor copies
where factor is the downsampling factor [0, 1] and copies is an integer giving the number of copies along each axis.