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

c     =====================================================
      subroutine combl()
c     =====================================================
c     
c     Create and initialize application specific common-blocks.
c
      implicit double precision (a-h,o-z)
c
      include  "cuser.i"
      parameter( lin = 8, lmechout = 11 )
      character *16 cwork
c
      Wref    = 10.d6
      a2ref   = 10.d3
c     Air
      gamma = 1.4d0  
      Wk   = 29.08d0*10.d-3/Wref
      cwork= 'Air'
      RU = 0.83140000E+01/(Wref*a2ref)
      PA = 0.10132500E+06/a2ref
c
      open(unit=lmechout, status='unknown', form='formatted', 
     &     file='chem.dat')
      write (lmechout,400) RU
      write (lmechout,401) PA
      write (lmechout,402) gamma
      write (lmechout,403) cwork     
      write (lmechout,404) Wk        
      close (lmechout)
c
 400  format('RU       ',e16.8)
 401  format('PA       ',e16.8)
 402  format('Gamma    ',e16.8)
 403  format('Sp       ',a16)
 404  format('W        ',e16.8)
c
      gamma1 = gamma-1.d0
c
c     # data in left state:
      rhol = 1.d0
      ul   = 0.d0
      pl   = 1000.d0
      rhoul = rhol * ul
      el = pl/gamma1 + 0.5d0*rhoul**2/rhol
      sloc1 = 0.1d0
c
c     # data in middle state:
      rhom = 1.d0
      um   = 0.d0
      pm   = 0.01d0
      rhoum = rhom * um
      em = pm/gamma1 + 0.5d0*rhoum**2/rhom
      sloc2 = 0.9d0
c
c     # data in right state:
      rhor = 1.d0
      ur   = 0.d0
      pr   = 100.d0
      rhour = rhor * ur
      er = pr/gamma1 + 0.5d0*rhour**2/rhor
c
      return
      end