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)
include "cuser.i"
c
parameter( lin = 5, lmechout = 11 )
character *16 cwork
dimension cwork(2)
c
cwork(1)= 'Product'
cwork(2)= 'Reactant'
c
Wk(1) = 1.d0
Wk(2) = 1.d0
RU = 1.d0
PA = 1.d0
c
open(unit=lmechout, status='unknown', form='formatted',
& file='chem.dat')
write (lmechout,400) RU
write (lmechout,401) PA
write (lmechout,402) (k,cwork(k),k=1,2)
write (lmechout,403) (k,Wk(k),k=1,2)
close (lmechout)
c
400 format('RU ',e16.8)
401 format('PA ',e16.8)
402 format('Sp(',i2.2,') ',a16)
403 format('W(',i2.2,') ',e16.8)
c
open(unit=lin,status='old',form='formatted',file='init.dat')
read (lin, *) gamma, q0, E0, f
read (lin, *) sloc, moving
close (lin)
c
gamma1 = gamma-1.d0
c
return
end