Source file:   README.txt
 Directory:   /Users/rjl/git/rjleveque/clawpack-4.6.3/book/chap23/acoustics
 Converted:   Mon Jan 21 2013 at 20:15:42   using clawcode2html
 This documentation file will not reflect any later changes in the source file.

$\phantom{******** If you see this on the webpage then the browser could not locate *********}$
$\phantom{******** the jsMath file load.js *********}$

$\newcommand{\vector}[1]{\left[\begin{array}{c} #1 \end{array}\right]}$ $\newenvironment{matrix}{\left[\begin{array}{cccccccccc}} {\end{array}\right]}$ $\newcommand{\A}{{\cal A}}$ $\newcommand{\W}{{\cal W}}$


CLAWPACK Sample Code

2d acoustics with a radially-symmetric solution.

Solved on a quadrilateral grid using the mapping specified in mapc2p.f. This is currently set to cover a chevron-shaped region with a grid that is not orthogonal and not smooth along one line.

See setaux.f for a description of how grid information is stored in the aux array. Variable coefficients for the acoustics equation are also allowed, though in this example the they are constant and set in setprob.f.

The data is a smooth radial hump specified in qinit.f The 1d equation with source term is solved in subdirectory 1drad to create an accurate reference solution.

Example [book/chap23/acoustics] to accompany the book
   Finite Volume Methods for Hyperbolic Problems by R. J. LeVeque.

Converted to Clawpack 4.4 form in 2009.

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 and specifies where output and plots should be directed. Type "make .help" at the Unix prompt for options.

driver.f [.html]
The driver routine allocates storage and then calls the main Clawpack routine.

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

rpn2acvq.f [.html]
Normal Riemann solver (normal to cell interface).

rpt2acvq.f [.html]
Transverse Riemann solver.

qinit.f [.html]
This subroutine sets the initial data at time $t=0$.

Python files [.html]
This file contains a function that specifies what run-time parameters will be used.

Some parameters that you might want to modify are described in the documentation . [.html]
This file contains a function that specifies what plots will be done and sets various plotting parameters.

Data files

Warning: These files are generally changed when setting up a run, usually in [.html]. [.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/2d/lib/claw2ez.f [.html]. [.html]
This file may contain various parameters used in setting the initial conditions or otherwise setting up the problem.

Library routines

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