rpn2from1d.f.html | ![]() |
Source file: rpn2from1d.f | |
Directory: /home/rjl/git/rjleveque/clawpack-4.x/amrclaw/2d/lib | |
Converted: Sun May 15 2011 at 19:16:15 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