Index of /clawpack/clawpack-4.6.3/book/chap16/vctraffic

Name	Last modified	Size

Parent Directory		-
Makefile	2013-01-21 18:05	1.7K
Makefile.html	2013-01-21 21:31	3.8K
_plots/	2013-01-21 19:56	-
b4step1.f	2013-01-21 18:05	615
b4step1.f.html	2013-01-21 21:31	2.4K
b4step1.o	2013-01-21 19:56	660
claw.data	2013-01-21 19:56	2.1K
claw.data.html	2013-01-21 21:31	4.4K
driver.f	2013-01-21 18:05	804
driver.f.html	2013-01-21 21:31	2.6K
driver.o	2013-01-21 19:56	1.1K
out1.f	2013-01-21 18:05	2.3K
out1.f.html	2013-01-21 21:31	4.5K
out1.o	2013-01-21 19:56	4.7K
pyclaw.log	2013-01-21 19:56	1.9K
qinit.f	2013-01-21 18:05	601
qinit.f.html	2013-01-21 21:31	2.4K
qinit.o	2013-01-21 19:56	1.0K
rp1trvfw.f	2013-01-21 18:05	3.2K
rp1trvfw.f.html	2013-01-21 21:31	5.7K
rp1trvfw.o	2013-01-21 19:56	3.8K
setaux.f	2013-01-21 18:05	716
setaux.f.html	2013-01-21 21:31	2.4K
setaux.o	2013-01-21 19:56	1.1K
setplot.py	2013-01-21 18:05	1.8K
setplot.py.html	2013-01-21 21:31	3.8K
setplot.pyc	2013-01-21 19:56	1.5K
setprob.f	2013-01-21 18:05	116
setprob.f.html	2013-01-21 21:31	1.6K
setprob.o	2013-01-21 19:56	564
setrun.py	2013-01-21 18:05	5.5K
setrun.py.html	2013-01-21 21:31	9.0K
step1fw.f	2013-01-21 18:05	4.3K
step1fw.f.html	2013-01-21 21:31	7.6K
step1fw.o	2013-01-21 19:56	6.2K
xclaw	2013-01-21 19:56	166K

README.html

README.html
Source file: README.txt
Directory: /Users/rjl/git/rjleveque/clawpack-4.6.3/book/chap16/vctraffic
Converted: Mon Jan 21 2013 at 20:15:35 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

Variable-coefficient traffic flow model. The flux function is f(q,x) = u(x)*q*(1-q) where the "speed limit" u(x) varies with x. It's value in cell i is stored in aux(i,1) (see setaux.f).

Here a Riemann problem is solved with u(x) = 2 for x<0 u(x) = 1 for x>0

The initial data is set in qinit.f. Set ql to 0.13 for Figure 16.9 Set ql to 0.2 for Figure 16.10

The method used is based on splitting the flux difference into f-waves Z as described in Section 15.5 and implemented in rp1trvfw.f.

Note that this requires a modified routine step1fw.f in place of the library routine step1.f in order to implement the modified wave-propagation algorithm.

The transonic rarefaction case is very sensitive to the flux splitting used. The approach taken in rp1trvfw.f gives the correct solution for this case but has not been extensively tested...

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

Converted to Clawpack 4.6 form in 2011.

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/1d/lib/claw1ez.f [.html] and is generally used to set any values needed for the specific problem being solved.

rp1.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.

qinit.f [.html]

This subroutine sets the initial data at time $t=0$.

Python files

setrun.py [.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 .

setplot.py [.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 setrun.py [.html].

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].

setprob.data [.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/1d/lib are used. See the Makefile [.html] to determine which ones are used.