Directory:
    ${PENCIL_HOME}/samples/supernova-driven-turbulence/lyapunov_HDtest
SVN Id:
    $Id$
Maintainer:
    Fred A. Gent <fred.gent.ncl/gmail[dot]com>
Added:
    21-Jun-2016
Status:
    introduced 15-Aug-2017
Recommended resolution:
    at least 16x16x1000 (maintain grid resolution <= 0.004 for larger Lx,Ly)
Comments:
    Adapted from exisiting run to include Dhruba's tools to calculate the
    Lyapunov exponents for the SN turbulence models.
    Use the 1D model (an example is given in sample ../1D_ism_init) to obtain
    the file init_ism.dat containing density and temperatue profiles for the 
    initial condition in this 3D supernova driven turbulence model of the ISM.
    The 1D model has identical vertical grid and hydrodynamic run parameters.
    Magnetic field and cosmic rays can be added to this 3D model
    
    The averaged SN heating substitutes for the SN driven turbulence and will be
    smoothly superceded by the SN explosions as the turbulence evolves.
    
    SN are spatially distributed according to a uniform distribution in the 
    horizontal plane and gaussian distribution vertically centred about the 
    centre of mass, which fluctuates near to the midlplane. Default SN frequency
    rates are milky way estimates and default mass and gravitional acceleration
    are taken from the solar vicinity.
    
    To obtain the 1D data file run
    > python ism_save.dat
    from the command line then copy the init_ism.dat file into the new 3D run
    directory
References:
    1. Gent, F.A., Shukurov, A., Fletcher, A., Sarson, G.R.S., &
       Mantere, M.J.: ``Supernova Driven Turbulence. I. The Multiphase
       Structure'',
       MNRAS, 432, 1396-1423