This is a more complex example that demonstrates quality of the shock wave
reproduction in the code. Complex problem was experimentally studied by Takayama.
(Cited from E.F. Toro, M. Spruce, W. Speares, Restoration of the Contact Surface
in the Harten-Lax-van Leer Riemann solver, J. Shock Waves 4 (1994) 25-34.).

We have in this folder pictures showing results of our calculation and the
experimental data. You can visually compare these results. Also can compare
values for the characteristic parameters such as velocities, pressure, density
in triple point.

To perform the calculations please follow same procedure as for the first example:


1. Prepare data file com00000.nc for the first step. Run the command:

comfile -k test.geom

   It will ask you some question (see User's Guide for details). To simplify
   the process of answering the questions we provide the test.ini file. It
   contains all the answer in correct order. You can simply pass content of
   this file as input for the comfile i.e.

comfile -k test.geom < test.ini

2. Inflow conditions are contained in vi-01.dat. Edit this file if you want
   to alter them.

3. Start the calculation. In case if you have three processors you can start
   it with the command:

mpirun -np 3 com -s com000000.nc

3. While the calculation runs start the vizir with command:

vizir

4. Connect vizir with the calculation. (see User Guide)
5. From Vizir:
	Select KYLCOM combustion method from Vizir's menu.
	Select 100,000 steps and pres play button.
	Press the "show" button
	Watch the evolution of the calculation model by selecting "continuous play".