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)
c
include 'air.i'
include 'cuser.i'
c
call stdcdat()
c
idisc = 1
alf = 1.d0
beta = 0.d0
x0 = 0.1d0
y0 = 0.d0
c
c
Minf = 10.d0
rhoamb = 1.4d0
uamb = 0.d0
vamb = 0.d0
wamb = 0.d0
pamb = 1.d0
eamb = pamb/gamma1 + .5d0*rhoamb*(uamb**2+vamb**2+wamb**2)
s1amb = 0.d0
s2amb = 0.d0
dcamb = 0.d0
t1 = 0.05d0
c
rhoshk = rhoamb
ushk = 10.d0
vshk = 0.d0
wshk = 0.d0
pshk = pamb
eshk = pshk/gamma1 + .5d0*rhoshk*(ushk**2+vshk**2+wshk**2)
s1shk = 0.d0
s2shk = 0.d0
dcshk = 0.d0
c
c
c Two spheres
NSph = 1
c
c The number of approximation points for
c integration on the surfaces is 8*Np()**2.
c
c Sphere 1
c
xm(1,1) = 0.32d0
xm(2,1) = 0.4d0
xm(3,1) = 0.4d0
rd(1) = 0.16d0
Np(1) = 180
c
c
return
end
c =====================================================
subroutine cblread(nst,xmt,rdt,npt)
c =====================================================
c
c Export values midpoint coordinates and radiuses of
c spheres into Problem.h.
c
implicit double precision (a-h,o-z)
c
include "cuser.i"
dimension xmt(15),rdt(5),npt(5)
c
nst = NSph
do n = 1, NSph
xmt(3*n-2) = xm(1,n)
xmt(3*n-1) = xm(2,n)
xmt(3*n) = xm(3,n)
rdt(n) = rd(n)
npt(n) = Np(n)
enddo
c
return
end