applications/euler/2d/Shockbubble/src/sauxrx2.f

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) = corny + (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