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"
c
parameter( lin = 5, lmechout = 11 )
character *16 cwork
dimension cwork(3), Xkl(2), Xkr(2)
c
open(unit=lin,status='old',form='formatted',file='init.dat')
read (lin, *) xc0, yc0, rc0, wlev
read (lin, *) tnts, ts
read (lin, *) xm, ym, rm, hm
close (lin)
c
cwork(1)= 'Gas'
cwork(2)= 'Water'
Wk(1) = 29.08d-3
Wk(2) = 18.015d-3
g(1) = 1.4d0
g(2) = 7.415d0
pinf(1) = 0.d0
pinf(2) = 296.2d6
c
RU = 0.83140000E+01
PA = 0.10132500E+06
c
cwork(3)= 'Pinf'
Wkhelp = 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,3)
write (lmechout,403) (k,Wk(k),k=1,2)
write (lmechout,403) 3,Wkhelp
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
pi = 4.d0*atan(1.d0)
es = 4520d3*tnts
ps = (3.d0*g(1)-3.d0)/(4.d0*pi)*es/rc0**3
rhos = ps*Wk(1)/(RU*ts)
c
pai = 1.d5
c rhoai = 1.29d0
rhoai = 40.d0
c
pw = 1.d5
rhow = 1027.d0
grav = -9.81d0
c
Nsegs = 3
Narcs = 3
c
xs(1,1) = xm-rm
xs(2,1) = ym
vns(1,1) = 1.d0
vns(2,1) = 0.d0
vls(1) = hm
ifs(1) = -1
c
xs(1,2) = xm+rm
xs(2,2) = ym
vns(1,2) = 1.d0
vns(2,2) = 0.d0
vls(2) = hm
ifs(2) = 1
c
xs(1,3) = xm-rm
xs(2,3) = ym+hm
vns(1,3) = 0.d0
vns(2,3) = -1.d0
vls(3) = 2.d0*rm
ifs(3) = -1
c
xa(1,1) = xm
xa(2,1) = ym
vna(1,1,1) = 0.d0
vna(2,1,1) = -1.d0
vna(1,2,1) = 0.d0
vna(2,2,1) = -1.d0
vra(1) = rm
ifa(1) = 1
c
xa(1,2) = xm-rm
xa(2,2) = ym+hm
vna(1,1,2) = -1.d0
vna(2,1,2) = 0.d0
vna(1,2,2) = 0.d0
vna(2,2,2) = 1.d0
vra(2) = 0.d0
ifa(2) = 1
c
xa(1,3) = xm+rm
xa(2,3) = ym+hm
vna(1,1,3) = 1.d0
vna(2,1,3) = 0.d0
vna(1,2,3) = 0.d0
vna(2,2,3) = 1.d0
vra(3) = 0.d0
ifa(3) = 1
c
return
end