|
rpn2from1d.f.html |
|
|
Source file: rpn2from1d.f
|
|
Directory: /home/rjl/git/claworg/clawpack-4.x/amrclaw/2d/lib
|
|
Converted: Sat Aug 6 2011 at 21:53:21
using clawcode2html
|
|
This documentation file will
not reflect any later changes in the source file.
|
c
c
c =====================================================
subroutine rpn2(ixy,maxm,meqn,mwaves,mbc,mx,ql,qr,auxl,auxr,
& wave,s,amdq,apdq)
c =====================================================
c
c # Solve Riemann problems for the 2D hyperbolic problem.
c
c # This is a wrapper for a 1d Riemann solver rp1 so that it
c # can be used with AMR by using the 2d amrclaw routines with my=1.
c
implicit double precision (a-h,o-z)
c
dimension wave(1-mbc:maxm+mbc, meqn, mwaves)
dimension s(1-mbc:maxm+mbc, mwaves)
dimension ql(1-mbc:maxm+mbc, meqn)
dimension qr(1-mbc:maxm+mbc, meqn)
dimension apdq(1-mbc:maxm+mbc, meqn)
dimension amdq(1-mbc:maxm+mbc, meqn)
c
if (ixy.eq.2) then
c write(6,*) '*** Error, this Riemann solver should only be used'
c write(6,*) '*** for a 1d problem with my=1'
c write(6,*) '*** ixy=2 and my = ',mx
c stop
endif
call rp1(maxm,meqn,mwaves,mbc,mx,ql,qr,auxl,auxr,
& wave,s,amdq,apdq)
return
end