applications/euler/1d/BlastWaves/src/init1.f

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. sloc1) then
             q(1,i) = rhol
             q(2,i) = rhoul
             q(3,i) = el
          else if (x(i) .ge. sloc1 .and. x(i) .lt. sloc2) then
             q(1,i) = rhom
             q(2,i) = rhoum
             q(3,i) = em
          else
             q(1,i) = rhor
             q(2,i) = rhour
             q(3,i) = er
          endif
 150   continue
c     
       return
       end