c =====================================================
subroutine rpn3(ixyz,maxm,meqn,mwaves,mbc,mx,ql,qr,maux,
& auxl,auxr,wave,s,amdq,apdq)
c =====================================================
c
c # Dummy function for 3D Riemann-solver in normal direction.
c # A problem-dependent rpn3()-function is always necessary!
c
c # Solve Riemann problem in x-direction if ixy=1
c # or in y-direction if ixy=2
c # or in z-direction if ixyz=3.
c
c # auxl(i,ma,2) contains auxiliary data for cells along this slice,
c # where ma=1,maux in the case where maux=method(7) > 0.
c # auxl(i,ma,1) and auxl(i,ma,3) contain auxiliary data along
c # neighboring slices that generally aren't needed in the rpn3 routine.
c
c # Note that the i'th Riemann problem has left state qr(i-1,:)
c # and right state ql(i,:)
c # From the basic routine step1, rp is called with ql = qr = q.
c
implicit double precision (a-h,o-z)
dimension ql(1-mbc:maxm+mbc, meqn)
dimension qr(1-mbc:maxm+mbc, meqn)
dimension auxl(1-mbc:maxm+mbc, maux, 3)
dimension auxr(1-mbc:maxm+mbc, maux, 3)
dimension s(1-mbc:maxm+mbc, mwaves)
dimension wave(1-mbc:maxm+mbc, meqn, mwaves)
dimension amdq(1-mbc:maxm+mbc, meqn)
dimension apdq(1-mbc:maxm+mbc, meqn)
c
return
end