2D AMRCLAW
rpn2from1d.f
Go to the documentation of this file.
1 
2 c
3 c
4 c =====================================================
5  subroutine rpn2(ixy,maxm,meqn,mwaves,maux,mbc,mx,ql,qr,auxl,auxr,
6  & wave,s,amdq,apdq)
7 c =====================================================
8 c
9 c # Solve Riemann problems for the 2D hyperbolic problem.
10 c
11 c # This is a wrapper for a 1d Riemann solver rp1 so that it
12 c # can be used with AMR by using the 2d amrclaw routines with my=1.
13 c
14  implicit double precision (a-h,o-z)
15 c
16  dimension wave(meqn, mwaves, 1-mbc:maxm+mbc)
17  dimension s(mwaves, 1-mbc:maxm+mbc)
18  dimension ql(meqn, 1-mbc:maxm+mbc)
19  dimension qr(meqn, 1-mbc:maxm+mbc)
20  dimension apdq(meqn, 1-mbc:maxm+mbc)
21  dimension amdq(meqn, 1-mbc:maxm+mbc)
22 c
23  if (ixy.eq.2) then
24 c write(6,*) '*** Error, this Riemann solver should only be used'
25 c write(6,*) '*** for a 1d problem with my=1'
26 c write(6,*) '*** ixy=2 and my = ',mx
27 c stop
28  endif
29 
30  call rp1(maxm,meqn,mwaves,maux,mbc,mx,ql,qr,auxl,auxr,
31  & wave,s,amdq,apdq)
32 
33  return
34  end
subroutine rpn2(ixy, maxm, meqn, mwaves, maux, mbc, mx, ql, qr, auxl, auxr, wave, s, amdq, apdq)
Definition: rpn2from1d.f:7