c =====================================================
subroutine combl()
c =====================================================
c
c Create and initialize application specific common-blocks.
c
c Copyright (C) 2002 Ralf Deiterding
c Brandenburgische Universitaet Cottbus
c
implicit double precision (a-h,o-z)
character *16 cwork
parameter( lmechout = 10 )
c
include "cuser.i"
c
a2ref = 0.10132500E+06/1.18922d0
Wref = 0.83140000E+01/a2ref
c # Air
gamma = 1.4d0
gamma1 = 0.4d0
Wk = 29.08d-3/Wref
cwork= 'AIR'
RU = 0.83140000E+01/(Wref*a2ref)
PA = 0.10132500E+06/a2ref
c
rhoshk = 1.862068449d0
ushk = 0.821677818d0
vshk = 0.d0
pshk = 2.458332279d0
c
qin(1) = rhoshk
qin(2) = ushk*qin(1)
qin(3) = vshk*qin(1)
qin(4) = pshk/gamma1+0.5d0*(qin(2)**2+qin(3)**2)/qin(1)
c
rhoamb = 1.d0
uamb = 0.d0
vamb = 0.d0
pamb = 1.d0
c
alf = 1.d0
beta = 0.d0
idisc = 1
x0 = 1.25d0
y0 = 0.d0
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
return
end