c
c
c
c     =====================================================
       subroutine qinit(maxmx,maxmy,maxmz,meqn,mbc,mx,my,mz,
     &                   xlower,ylower,zlower,dx,dy,dz,q,maux,aux)
c     =====================================================
c
c     # Set initial conditions for q.
c
       implicit double precision (a-h,o-z)
c
       dimension q(1-mbc:maxmx+mbc, 1-mbc:maxmy+mbc, 
     &             1-mbc:maxmz+mbc, meqn)
c
c
      do 215 i=1,mx
         xcell = xlower + (i-0.5d0)*dx
	 do 210 j=1,my
            ycell = ylower + (j-0.5d0)*dy
	    do 205 k=1,mz
               zcell = zlower + (k-0.5d0)*dz
	       r = dsqrt(xcell**2 + ycell**2 + zcell**2)
	       rho = 1.d0 + 10.0d0*dexp(-20*(r-0.0d0)**2)
	       q(i,j,k,1) = rho
	       q(i,j,k,2) = 0.d0
	       q(i,j,k,3) = 0.d0
	       q(i,j,k,4) = 0.d0
	       q(i,j,k,5) = rho
  205           continue
  210        continue
  215     continue
       return
       end