rpt2.f.html | clawcode2html |
Source file: rpt2.f | |
Directory: /home/rjl/www/pubs/cise08/cise08levequeV2/lib | |
Converted: Wed Jul 2 2008 at 13:39:55 | |
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