|
rpt2.f.html |
|
|
Source file: rpt2.f
|
|
Directory: /home/rjl/git/claworg/clawpack-4.x/amrclaw/2d/lib
|
|
Converted: Sat Aug 6 2011 at 21:59:31
using clawcode2html
|
|
This documentation file will
not reflect any later changes in the source file.
|
c
c
c =====================================================
subroutine rpt2(ixy,maxm,meqn,mwaves,mbc,mx,
& ql,qr,aux1,aux2,aux3,
& ilr,asdq,bmasdq,bpasdq)
c =====================================================
c
c # Riemann solver in the transverse direction.
c # This is a dummy routine that returns zeros and is only intended
c # to illustrate the format of this routine. See various example
c # directories for better examples.
c # This dummy routine can be used if transverse solves are not being
c # used, i.e. if method(3)=0.
c
c # Split asdq (= A^* \Delta q, where * = + or -)
c # into down-going flux difference bmasdq (= B^- A^* \Delta q)
c # and up-going flux difference bpasdq (= B^+ A^* \Delta q)
c
c
implicit double precision (a-h,o-z)
dimension ql(1-mbc:maxm+mbc, meqn)
dimension qr(1-mbc:maxm+mbc, meqn)
dimension asdq(1-mbc:maxm+mbc, meqn)
dimension bmasdq(1-mbc:maxm+mbc, meqn)
dimension bpasdq(1-mbc:maxm+mbc, meqn)
c
do i = 2-mbc, mx+mbc
do m=1,meqn
bmasdq(i,m) = 0.d0
bpasdq(i,m) = 0.d0
enddo
enddo
c
return
end