c
c =====================================================
subroutine ic(maxmx,meqn,mbc,mx,x,dx,q)
c =====================================================
implicit double precision (a-h,o-z)
c
include "cuser.i"
c
dimension q(meqn,1-mbc:maxmx+mbc)
dimension x(1-mbc:maxmx+mbc)
c
do 150 i=1,mx
if (x(i) .lt. sloc) then
q(1,i) = (1.d0-Yl)*rhol
q(2,i) = Yl*rhol
q(3,i) = ul*rhol
q(4,i) = pl/gamma1 + Yl*rhol*q0 + 0.5d0*rhol*ul**2
else
q(1,i) = (1.d0-Yr)*rhor
q(2,i) = Yr*rhor
q(3,i) = ur*rhor
q(4,i) = pr/gamma1 + Yr*rhor*q0 + 0.5d0*rhor*ur**2
endif
150 continue
c
return
end
c