c
c =====================================================
subroutine setaux(maxmx,meqn,mbc,ibx,mx,q,
& aux,maux,xc,dx,t,dt)
c =====================================================
c
c Copyright (C) 2002 Ralf Deiterding
c Brandenburgische Universitaet Cottbus
c
implicit double precision (a-h, o-z)
c
dimension q(meqn, 1-ibx*mbc:maxmx+ibx*mbc)
dimension aux(maux, 1-ibx*mbc:maxmx+ibx*mbc)
c
c # define the velocity field
u(x) = 2.d0
do 10 i=1-ibx*mbc,mx+ibx*mbc
aux(1,i) = u(xc+(i+ibx*mbc-1.d0)*dx)
10 continue
c
return
end