Directory:
    ${PENCIL_HOME}/samples/0d-tests/heating_ionize
SVN Id:
    $Id$
Maintainer:
    Axel Brandenburg <brandenb/nordita[dot]org>
Added:
    23-Apr-2004
Status:
    Works since 2003/08/12 (the date when doc/figs/pTTss.ps was made)
Recommended resolution:
    1x1x1 (zero-dimensional run)
Comments:
    In this test we solve an ODE with the Pencil Code (hence 0-D).
        rho T ds/dt = heat_uniform
    We have variable degree of hydrogen ionization, but Helium stays neutral.
    Our units are [rho]=1e-6 g/cm^3, [u]=km/s, [x]=Mm, [T]=K.
    idl.png shows the resulting temperature as a function of thermal energy.
    This result is plotted using plot_result.csh
    In order for the run not to last too long on slow machines, we have
    changed the file run.in to run longer and fewer time steps. The original
    file used for making the figure is still available as "run.in-fig.pTTss".

    To reproduce the plot pTTss.ps in the manual, one must adjust the timestep
    by hand, i.e. it longer as time goes on. The density value can be adjusted
    by changing simply unit_density in start.in.

    In print.in we should add one more digit for better accuracy.