c 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) = rhol q(2,i) = rhoul q(5,i) = el else q(1,i) = rhor q(2,i) = rhour q(5,i) = er endif q(3,i) = 0.d0 q(4,i) = 0.d0 q(6,i) = 0.d0 q(7,i) = 0.d0 q(8,i) = 0.d0 150 continue c return end