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