66       subroutine initSponge( y,x)  argument
71       real*8   y(nshg,nflow), x(numnp,3)  local
102          radsqr=x(id,2)**2+x(id,1)**2
106             zval=max(x(id,3),zstart)
128       subroutine initBCt( x, iBC, BC )  argument
133       real*8   x(numnp,nsd), BC(nshg,ndofBC), rj1,rj2,rj3,rj4,distd,epsd  local
196                   dd= distds(x1,x2,x3,x(i,1),x(i,2),x(i,3))
297       subroutine BCint(timel,shp,shgl,shpb,shglb,x,BC,iBC)  argument
304       real*8   x(numnp,nsd),     local
339       real*8 x1,y1,z1,x2,y2,z2,x,y,z  local
340       x=x1-x2
343       distds=x*x+y*y+z*z