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(3,i) = el else q(1,i) = rhor q(2,i) = rhour q(3,i) = er endif 150 continue c return end