

This directory contains an example driver routine for solving the 2D
Euler equations of gas dynamics  with piecewise
constant initial data having 4 different values in 4 quadrants.
The particular problem is specified in the routine ic2rp2.f 
and corresponds to the problem shown in Figure 4 of

@article{csr-col-glaz,
  author="C. W. Schulz-Rinne and J. P. Collins and H. M. Glaz",
  title="Numerical Solution of the {R}iemann Problem for Two-Dimensional Gas
   Dynamics",
  journal="SIAM J. Sci. Comput.",
  volume="14",
  year="1993",
  pages="1394-1414"
}


To run, type
   a.out < data
This is currently set up to use a 200 x 200 grid, but this can be changed in
the data files.  With this grid, the program runs for 416 time steps and
takes about 20 minutes on a DEC 3000/600 Alpha workstation.  The output file
fort.101 is 2.6 MB and fort.dump is 4MB (creating the latter can be
suppressed by deleting the call to dump2 in driver.f).

matlab command:   plotq
