1d/equations/euler/rp/flx1eu.f

c
c =========================================================
      subroutine flx1(maxmx,meqn,mbc,mx,q,maux,aux,f)
c =========================================================
c
c     # Flux-function for the one-dimensional Euler equations.
c
      implicit double precision (a-h,o-z)
      dimension q(1-mbc:maxmx+mbc, meqn)
      dimension aux(1-mbc:maxmx+mbc, maux)
      dimension f(1-mbc:maxmx+mbc, meqn)
      common /param/  gamma,gamma1
c
      do i=1-mbc,mx+mbc
         u = q(i,2)/q(i,1)
         p = gamma1*(q(i,3) - 0.5d0*q(i,1)*u**2)
         f(i,1) =   q(i,2)
         f(i,2) = u*q(i,2)+p
         f(i,3) = u*(q(i,3)+p)
      enddo
c
      return
      end