c c c ===================================================== subroutine setaux(maxmx,maxmy,meqn,mbc,ibx,iby,mx,my,q, & aux,maux,xc,yc,dx,dy,t,dt) c ===================================================== implicit double precision (a-h, o-z) dimension q(meqn, 1-ibx*mbc:maxmx+ibx*mbc, & 1-iby*mbc:maxmy+iby*mbc) dimension aux(maux, 1-ibx*mbc:maxmx+ibx*mbc, & 1-iby*mbc:maxmy+iby*mbc) c do 10 j = 1-iby*mbc, my+iby*mbc aux(1,1-ibx*mbc,j) = yc + (j+iby*mbc-.5)*dy do 10 i = 2-ibx*mbc, mx+ibx*mbc aux(1,i,j) = aux(1,1-ibx*mbc,j) 10 continue return end