Directory:
    ${PENCIL_HOME}/samples/supernova-driven-turbulence/snowplough/
SVN Id:
    $Id: README 13067 2010-01-25 07:08:53Z f.a.gent@ncl.ac.uk $
Maintainer:
    Frederick Gent <fred.gent.ncl/gmail[dot]com>
Added:
    22-Jan-2010
Status:
    has been working in this form since 22-Jan-2010
Recommended resolution:
    at least 24x24x24 (A larger box will be required with resolution 4pc) 
    for 1 atom per cc a box radius > 120pc is required end time 5 Myr
    for 0.001 atom per cc a box radius > 500pc is required end time 2.5 Myr
Comments:
    designed to verify the code reproduces the Snowplough solution for a single
    supernova remnant in a monotomic uniform ISM at rest. 
    if using density < 0.5 atom per cc increase box to keep away from boundary
    injection of 10^51 ergs => 50% kinetic and 50% thermal energy.

    To calculate analytic profiles and plot 1D slices of SN radius and profile speed
    run 
    > python snowplough.py