CLAWPACK Code

These codes were developed to accompany the paper

Universality in the run-up of shock waves to the surface of a star by C. Gundlach and R.J. LeVeque.
Submitted to J. Fluid Mech., 2010.

Link to preprint  .

1d Euler equations with polytropic or ideal gas equation of state.

The EOS and parameters are specified in setrun.py [.html].

See Instructions

Plots of results

After running this code and creating plots via "make .plots", you should be able to view the plots in _plots/_PlotIndex.html .

Fortran files

Makefile [.html]

Determines which version of fortran files are used when compiling the code with make.

setprob.f [.html]

A routine by this name is called by the library routine clawpack/1d/lib/claw1ez.f [.html] and is generally used to set any values needed for the specific problem being solved.

rp1eu_geneos_wb.f [.html]

This is the Riemann solver, which takes the $q$ values stored in the arrays ql and qr and returns the waves in the array wave and speeds in the array s that result in solving the Riemann problem at each cell interface, and the fluctuations amdq and apdq. See $CLAW/doc/rp1.html for more information about 1d Riemann solvers.

pres_psi.f [.html]

Specifies the EOS.

qinit.f [.html]

This subroutine sets the initial data q(x,0) at time t=0.

b4step1.f [.html]

Handles the regridding when the shock approaches the right boundary.

Python files

setrun.py [.html]

This file contains a function that specifies what run-time parameters will be used.

setplot.py [.html]

This file contains a function that specifies what plots will be done and sets various plotting parameters.

make_plots.py [.html]

Runs the code and makes the plots for figures in the paper.

setrun2.py [.html]

Used by make_plots.py [.html].

setplot2.py [.html]

Used by make_plots.py [.html].

Data files

Warning: These files are generally changed when setting up a run and the versions here may not be the ones actually used.

claw.data [.html]

This file contains a number of parameter values that are used by CLAWPACK. The values in this file are read by the library routine clawpack/1d/lib/claw1ez.f [.html]. Each line contains one or more values to be read in, followed by comments that are ignored by the Fortran code but may be used by Pythons scripts.

Some parameters that you might want to modify are described in the documentation .

setprob.data [.html]

This file contains the advection velocity $u$ and various other parameters used in setting the initial conditions. Values in this file are read in by the subroutine setprob.f95 [.html].

Library routines

In addition to the Fortran routines in this library, several library routines from $CLAW/clawpack/1d/lib are used. See the Makefile [.html] to determine which ones are used.

Instructions

To run code, edit setrun.py [.html] and then:

    $ make .output

View plots interactively with Iplotclaw or use "make .plots" to create html files.

To run the code 4 times for the test cases used in the paper:

    $ python make_plots.py