c c c ===================================================== subroutine bc2wall(maxmx,maxmy,meqn,mbc,mx,my,q,aux,t) c ===================================================== c c # solid wall boundary conditions: c ------------------------------------------------ c c # Extend the data from the interior cells (1:mx, 1:my) c # to the ghost cells outside the region: c # (i, 1-jbc) for jbc = 1,mbc, i = 1-mbc, mx+mbc c # (i, my+jbc) for jbc = 1,mbc, i = 1-mbc, mx+mbc c # (1-ibc, j) for ibc = 1,mbc, j = 1-mbc, my+mbc c # (mx+ibc, j) for ibc = 1,mbc, j = 1-mbc, my+mbc c implicit double precision (a-h,o-z) dimension q(1-mbc:maxmx+mbc, 1-mbc:maxmy+mbc, meqn) dimension aux(1-mbc:maxmx+mbc, 1-mbc:maxmy+mbc, *) dimension isx(3),isy(3) data isx /1,-1,1/ data isy /1,1,-1/ c c do 112 m=1,meqn do 111 jbc=1,mbc do 110 i = 1-mbc, mx+mbc q(i,1-jbc,m) = q(i,jbc,m)*isy(m) q(i,my+jbc,m) = q(i,my+1-jbc,m)*isy(m) 110 continue 111 continue 112 continue c do 122 m=1,meqn do 121 j = 1-mbc, my+mbc do 120 ibc=1,mbc q(1-ibc,j,m) = q(ibc,j,m)*isx(m) q(mx+ibc,j,m) = q(mx+1-ibc,j,m)*isx(m) 120 continue 121 continue 122 continue return end