c =====================================================
subroutine combl()
c =====================================================
c
c Create and initialize application specific common-blocks.
c
c Copyright (C) 2003-2007 California Institute of Technology
c Ralf Deiterding, ralf@amroc.net
c
implicit double precision (a-h,o-z)
include "cuser.i"
parameter( lin = 5 )
c
call stdcdat()
c
open(unit=lin,status='old',form='formatted',file='init.dat')
read (lin, *) rhol, ul, pl
read (lin, *) rhor, ur, pr
read (lin, *) sloc
close (lin)
c
gamma = 1.4d0
gamma1 = gamma-1.d0
c # data in left state:
rhoul = rhol * ul
el = pl/gamma1 + 0.5d0*rhoul**2/rhol
c
c # data in right state:
rhour = rhor * ur
er = pr/gamma1 + 0.5d0*rhour**2/rhor
c
return
end