c
c =====================================================
subroutine setaux(maxmx,maxmy,maxmz,meqn,mbc,ibx,iby,ibz,
& mx,my,mz,q,aux,maux,cornx,corny,cornz,dx,dy,dz,t,dt)
c =====================================================
c
c Copyright (C) 2003-2007 California Institute of Technology
c Ralf Deiterding, ralf@cacr.caltech.edu
c
implicit double precision (a-h, o-z)
dimension q(meqn, 1-ibx*mbc:maxmx+ibx*mbc,
& 1-iby*mbc:maxmy+iby*mbc, 1-ibz*mbc:maxmz+ibz*mbc)
dimension aux(maux, 1-ibx*mbc:maxmx+ibx*mbc,
& 1-iby*mbc:maxmy+iby*mbc, 1-ibz*mbc:maxmz+ibz*mbc)
c
do 10 j = 1-iby*mbc, my+iby*mbc
do 10 i = 1-ibx*mbc, mx+ibx*mbc
r = dsqrt((cornx+(i+ibx*mbc-.5)*dx)**2+
& (corny+(j+iby*mbc-.5)*dy)**2)
do 10 k = 1-ibz*mbc, mz+ibz*mbc
aux(1,i,j,k) = cornz+(k+ibz*mbc-.5)*dz
aux(2,i,j,k) = r
10 continue
return
end