navier1.f File Reference

Functions/Subroutines
subroutine	plan1 (igeom)
subroutine	crespuz (respr, g1, g2, g3, h1, h2, h2inv, intype)
subroutine	cresvuz (resv1, resv2, resv3)
subroutine	makeg (out1, out2, out3, h1, h2, intype)
subroutine	ctolspl (tolspl, respr)
subroutine	ortho (respr)
subroutine	cdabdtp (ap, wp, h1, h2, h2inv, intype)
subroutine	opgrad (out1, out2, out3, inp)
subroutine	cdtp (dtx, x, rm2, sm2, tm2, isd)
subroutine	multd (dx, x, rm2, sm2, tm2, isd, iflg)
subroutine	ophx (out1, out2, out3, inp1, inp2, inp3, h1, h2)
subroutine	opbinv (out1, out2, out3, inp1, inp2, inp3, h2inv)
subroutine	opbinv1 (out1, out2, out3, inp1, inp2, inp3, SCALE)
subroutine	uzprec (rpcg, rcg, h1m1, h2m1, intype, wp)
subroutine	eprec (z2, r2)
subroutine	convprn (iconv, rbnorm, rrpt, res, z, tol)
subroutine	convpr (res, tol, iconv, rbnorm)
subroutine	chktcg2 (tol, res, iconv)
subroutine	dudxyz (du, u, rm1, sm1, tm1, jm1, imsh, isd)
subroutine	convopo (conv, fi)
subroutine	conv2 (dtfi, fi)
subroutine	cmask (cmask1, cmask2)
subroutine	makef
subroutine	makeuf
subroutine	nekuf (f1, f2, f3)
subroutine	advab
subroutine	makebdf
subroutine	makeabf
subroutine	setab3 (ab0, ab1, ab2)
subroutine	setabbd (ab, dtlag, nab, nbd)
subroutine	setbd (bd, dtbd, nbd)
subroutine	bdsys (a, b, dt, nbd, ldim)
subroutine	tauinit (tau, ilag)
subroutine	velinit (vel1, vel2, vel3, ilag)
subroutine	velconv (vxn, vyn, vzn, tau)
subroutine	frkconv (y, x, mask)
subroutine	velchar (vel, vn, vlag, nbd, tau, dtbd)
subroutine	lagvel
subroutine	hypmsk3 (hv1msk, hv2msk, hv3msk)
subroutine	setordbd
subroutine	normsc (h1, semi, l2, linf, x, imesh)
subroutine	normvc (h1, semi, l2, linf, x1, x2, x3)
subroutine	genwp (wp, wm2, p)
subroutine	convuz (ifstuz)
subroutine	convsp (ifstsp)
subroutine	antimsk (y, x, xmask, n)
subroutine	opamask (vbdry1, vbdry2, vbdry3)
subroutine	opmask (res1, res2, res3)
subroutine	opadd2 (a1, a2, a3, b1, b2, b3)
subroutine	opsub2 (a1, a2, a3, b1, b2, b3)
subroutine	opsub3 (a1, a2, a3, b1, b2, b3, c1, c2, c3)
subroutine	opcolv3 (a1, a2, a3, b1, b2, b3, c)
subroutine	opcolv (a1, a2, a3, c)
subroutine	opcol2 (a1, a2, a3, b1, b2, b3)
subroutine	opchsgn (a, b, c)
subroutine	opcopy (a1, a2, a3, b1, b2, b3)
subroutine	rotate_cyc (r1, r2, r3, idir)
subroutine	opdssum (a, b, c)
subroutine	opdsop (a, b, c, op)
subroutine	opicol2 (a1, a2, a3, b1, b2, b3)
subroutine	oprzero (a, b, c)
subroutine	oprone (a, b, c)
subroutine	opcmult (a, b, c, const)
subroutine	opcolv2c (a1, a2, a3, b1, b2, c)
subroutine	opcolv2 (a1, a2, a3, b1, b2)
subroutine	opadd2col (a1, a2, a3, b1, b2, b3, c)
subroutine	opcolv3c (a1, a2, a3, b1, b2, b3, c, d)
subroutine	uzawa (rcg, h1, h2, h2inv, intype, iter)
subroutine	mapw (md, nd, m1, n1, mflg)
subroutine	mapwp (md, nd, m1, n1, mflg)
subroutine	specmp (b, nb, a, na, ba, ab, w)
subroutine	setmap (n1, nd)
subroutine	set_pnd (P, LkD, LkNt, N, D)
subroutine	convop (conv, fi)
subroutine	conv1d (dfi, fi)
subroutine	conv1 (du, u)
subroutine	conv1no (du, u)
subroutine	conv1rk (du, dv, dw, u, v, w)
subroutine	velconvv (vxn, vyn, vzn, tau)
subroutine	frkconvv (du, dv, dw, u, v, w, mu)
subroutine	conv1rk2 (du, dv, dw, u, v, w, cu, cv, cw, beta, wk)
subroutine	frkconvv2 (du, dv, dw, u, v, w, cu, cv, cw, beta, mu, wk)
subroutine	hypmsk3v (msk, mask)
subroutine	ophyprk (vel1, vel2, vel3, ilag)
subroutine	opdiv (outfld, inpx, inpy, inpz)
subroutine	opgradt (outx, outy, outz, inpfld)
subroutine	setproj (n1, nd)
subroutine	set_pndoi (Pt, P, LkNt, N, D)
subroutine	wgradm1 (ux, uy, uz, u, nel)
subroutine	makevis
subroutine	comp_siej (U1, U2, U3)
subroutine	wlaplacian (out, a, diff, ifld)
subroutine	explstrs
subroutine	expl_strs (w1, w2, w3, u1, u2, u3)
subroutine	expl_strs_e (w1, w2, w3, u1, u2, u3, e)
subroutine	expl_strs_e_3d (w1, w2, w3, u1, u2, u3, e)
subroutine	expl_strs_e_2d (w1, w2, u1, u2, e)
subroutine	gradl_rst_t (u, ur, us, ut, md, if3d)
subroutine	gradl_rst (ur, us, ut, u, md, if3d)
subroutine	local_grad3_t (u, ur, us, ut, N, e, D, Dt, w)
subroutine	local_grad2_t (u, ur, us, N, e, D, Dt, w)
subroutine	get_dgll_ptr (ip, mx, md)
subroutine	gen_dgll (dgl, dgt, mp, np, w)

Function/Subroutine Documentation

advab()

subroutine advab

antimsk()

subroutine antimsk	(	real, dimension(1)	y,
		real, dimension(1)	x,
		real, dimension(1)	xmask,
			n
	)

bdsys()

subroutine bdsys	(	real, dimension(ldim,9)	a,
		real, dimension(9)	b,
		real, dimension(9)	dt,
			nbd,
			ldim
	)

cdabdtp()

subroutine cdabdtp	(	real, dimension (lx2,ly2,lz2,1)	ap,
		real, dimension (lx2,ly2,lz2,1)	wp,
		real, dimension (lx1,ly1,lz1,1)	h1,
		real, dimension (lx1,ly1,lz1,1)	h2,
		real, dimension (lx1,ly1,lz1,1)	h2inv,
			intype
	)

cdtp()

subroutine cdtp	(	real, dimension (lx1*ly1*lz1,lelv)	dtx,
		real, dimension (lx2*ly2*lz2,lelv)	x,
		real, dimension (lx2*ly2*lz2,lelv)	rm2,
		real, dimension (lx2*ly2*lz2,lelv)	sm2,
		real, dimension (lx2*ly2*lz2,lelv)	tm2,
			isd
	)

chktcg2()

subroutine chktcg2	(		tol,
		real, dimension (lx2,ly2,lz2,lelv)	res,
			iconv
	)

cmask()

subroutine cmask	(	real, dimension (lx1,ly1,lz1,lelv)	cmask1,
		real, dimension (lx1,ly1,lz1,lelv)	cmask2
	)

comp_siej()

subroutine comp_siej	(	dimension(lx1,ly1,lz1,1)	U1,
		dimension(lx1,ly1,lz1,1)	U2,
		dimension(lx1,ly1,lz1,1)	U3
	)

conv1()

subroutine conv1	(	real, dimension (lx1*ly1*lz1,1)	du,
		real, dimension (lx1,ly1,lz1,1)	u
	)

conv1d()

subroutine conv1d	(	real, dimension (lx1,ly1,lz1,1)	dfi,
		real, dimension (lx1,ly1,lz1,1)	fi
	)

conv1no()

subroutine conv1no	(	real, dimension (lx1*ly1*lz1,1)	du,
		real, dimension (lx1,ly1,lz1,1)	u
	)

conv1rk()

subroutine conv1rk	(	real, dimension(lx1*ly1*lz1,1)	du,
		real, dimension(lx1*ly1*lz1,1)	dv,
		real, dimension(lx1*ly1*lz1,1)	dw,
		real, dimension (lx1,ly1,lz1,1)	u,
		real, dimension (lx1,ly1,lz1,1)	v,
		real, dimension (lx1,ly1,lz1,1)	w
	)

conv1rk2()

subroutine conv1rk2	(	real, dimension(lx1*ly1*lz1,1)	du,
		real, dimension(lx1*ly1*lz1,1)	dv,
		real, dimension(lx1*ly1*lz1,1)	dw,
		real, dimension (lx1,ly1,lz1,1)	u,
		real, dimension (lx1,ly1,lz1,1)	v,
		real, dimension (lx1,ly1,lz1,1)	w,
		real, dimension(lx1,ly1,lz1,1)	cu,
		real, dimension(lx1,ly1,lz1,1)	cv,
		real, dimension(lx1,ly1,lz1,1)	cw,
			beta,
		real, dimension(lx1,ly1,lz1,3)	wk
	)

conv2()

subroutine conv2	(	real, dimension (lx1,ly1,lz1,1)	dtfi,
		real, dimension (lx1,ly1,lz1,1)	fi
	)

convop()

subroutine convop	(	real, dimension (lx1,ly1,lz1,1)	conv,
		real, dimension (lx1,ly1,lz1,1)	fi
	)

convopo()

subroutine convopo	(	real, dimension (lx1,ly1,lz1,1)	conv,
		real, dimension (lx1,ly1,lz1,1)	fi
	)

convpr()

subroutine convpr	(	real, dimension (lx2,ly2,lz2,lelv)	res,
			tol,
			iconv,
			rbnorm
	)

convprn()

subroutine convprn	(		iconv,
			rbnorm,
			rrpt,
		real, dimension (1)	res,
		real, dimension (1)	z,
			tol
	)

convsp()

subroutine convsp

(

logical

ifstsp

)

convuz()

subroutine convuz

(

logical

ifstuz

)

crespuz()

subroutine crespuz	(	real, dimension (lx2,ly2,lz2,lelv)	respr,
		real, dimension (lx1,ly1,lz1,lelv)	g1,
		real, dimension (lx1,ly1,lz1,lelv)	g2,
		real, dimension (lx1,ly1,lz1,lelv)	g3,
		real, dimension (lx1,ly1,lz1,lelv)	h1,
		real, dimension (lx1,ly1,lz1,lelv)	h2,
		real, dimension (lx1,ly1,lz1,lelv)	h2inv,
			intype
	)

cresvuz()

subroutine cresvuz	(	real, dimension (lx1,ly1,lz1,1)	resv1,
		real, dimension (lx1,ly1,lz1,1)	resv2,
		real, dimension (lx1,ly1,lz1,1)	resv3
	)

ctolspl()

subroutine ctolspl	(		tolspl,
		real, dimension (lx2,ly2,lz2,lelv)	respr
	)

dudxyz()

subroutine dudxyz	(	real, dimension (lx1,ly1,lz1,1)	du,
		real, dimension (lx1,ly1,lz1,1)	u,
		real, dimension (lx1,ly1,lz1,1)	rm1,
		real, dimension (lx1,ly1,lz1,1)	sm1,
		real, dimension (lx1,ly1,lz1,1)	tm1,
		real, dimension (lx1,ly1,lz1,1)	jm1,
			imsh,
			isd
	)

eprec()

subroutine eprec	(	real, dimension (lx2,ly2,lz2,lelv)	z2,
		real, dimension (lx2,ly2,lz2,lelv)	r2
	)

expl_strs()

subroutine expl_strs	(	real, dimension(1)	w1,
		real, dimension(1)	w2,
		real, dimension(1)	w3,
		real, dimension(1)	u1,
		real, dimension(1)	u2,
		real, dimension(1)	u3
	)

expl_strs_e()

subroutine expl_strs_e	(	real, dimension(1)	w1,
		real, dimension(1)	w2,
		real, dimension(1)	w3,
		real, dimension(1)	u1,
		real, dimension(1)	u2,
		real, dimension(1)	u3,
		integer	e
	)

expl_strs_e_2d()

subroutine expl_strs_e_2d	(	real, dimension(1)	w1,
		real, dimension(1)	w2,
		real, dimension(1)	u1,
		real, dimension(1)	u2,
		integer	e
	)

expl_strs_e_3d()

subroutine expl_strs_e_3d	(	real, dimension(1)	w1,
		real, dimension(1)	w2,
		real, dimension(1)	w3,
		real, dimension(1)	u1,
		real, dimension(1)	u2,
		real, dimension(1)	u3,
		integer	e
	)

explstrs()

subroutine explstrs

frkconv()

subroutine frkconv	(	real, dimension (lx1,ly1,lz1,1)	y,
		real, dimension (lx1,ly1,lz1,1)	x,
		real, dimension (lx1,ly1,lz1,1)	mask
	)

frkconvv()

subroutine frkconvv	(	real, dimension(1)	du,
		real, dimension(1)	dv,
		real, dimension(1)	dw,
		real, dimension (1)	u,
		real, dimension (1)	v,
		real, dimension (1)	w,
		integer, dimension(0:1)	mu
	)

frkconvv2()

subroutine frkconvv2	(	real, dimension(1)	du,
		real, dimension(1)	dv,
		real, dimension(1)	dw,
		real, dimension (1)	u,
		real, dimension (1)	v,
		real, dimension (1)	w,
		real, dimension(1)	cu,
		real, dimension(1)	cv,
		real, dimension(1)	cw,
			beta,
		integer, dimension(0:1)	mu,
		real, dimension(lx1*ly1*lz1,3)	wk
	)

gen_dgll()

subroutine gen_dgll	(	real, dimension(mp,np)	dgl,
		real, dimension(np*mp)	dgt,
			mp,
			np,
		real, dimension(1)	w
	)

genwp()

subroutine genwp	(	real, dimension (lx2,ly2,lz2,lelv)	wp,
		real, dimension (lx2,ly2,lz2)	wm2,
		real, dimension (lx2,ly2,lz2,lelv)	p
	)

get_dgll_ptr()

subroutine get_dgll_ptr	(		ip,
			mx,
			md
	)

gradl_rst()

subroutine gradl_rst	(	real, dimension(1)	ur,
		real, dimension(1)	us,
		real, dimension(1)	ut,
		real, dimension(1)	u,
			md,
		logical	if3d
	)

gradl_rst_t()

subroutine gradl_rst_t	(	real, dimension(1)	u,
		real, dimension(1)	ur,
		real, dimension(1)	us,
		real, dimension(1)	ut,
			md,
		logical	if3d
	)

hypmsk3()

subroutine hypmsk3	(	real, dimension (lx1,ly1,lz1,1)	hv1msk,
		real, dimension (lx1,ly1,lz1,1)	hv2msk,
		real, dimension (lx1,ly1,lz1,1)	hv3msk
	)

hypmsk3v()

subroutine hypmsk3v	(	integer, dimension(0:1)	msk,
		real, dimension(lxyz1,lelt)	mask
	)

lagvel()

subroutine lagvel

local_grad2_t()

subroutine local_grad2_t	(	real, dimension (0:n,0:n,1)	u,
		real, dimension(0:n,0:n)	ur,
		real, dimension(0:n,0:n)	us,
			N,
		integer	e,
		real, dimension (0:n,0:n)	D,
		real, dimension(0:n,0:n)	Dt,
		real, dimension (0:n,0:n)	w
	)

local_grad3_t()

subroutine local_grad3_t	(	real, dimension (0:n,0:n,0:n,1)	u,
		real, dimension(0:n,0:n,0:n)	ur,
		real, dimension(0:n,0:n,0:n)	us,
		real, dimension(0:n,0:n,0:n)	ut,
			N,
		integer	e,
		real, dimension (0:n,0:n)	D,
		real, dimension(0:n,0:n)	Dt,
		real, dimension (0:n,0:n,0:n)	w
	)

makeabf()

subroutine makeabf

makebdf()

subroutine makebdf

makeg()

subroutine makeg	(	real, dimension (lx1,ly1,lz1,lelv)	out1,
		real, dimension (lx1,ly1,lz1,lelv)	out2,
		real, dimension (lx1,ly1,lz1,lelv)	out3,
		real, dimension (lx1,ly1,lz1,lelv)	h1,
		real, dimension (lx1,ly1,lz1,lelv)	h2,
			intype
	)

makeuf()

subroutine makeuf

makevis()

subroutine makevis

mapw()

subroutine mapw	(	real, dimension(lxd,lyd,lzd,lelv)	md,
			nd,
		real, dimension(lx1,ly1,lz1,lelv)	m1,
			n1,
			mflg
	)

mapwp()

subroutine mapwp	(	real, dimension(lxd,lyd,lzd,lelv)	md,
			nd,
		real, dimension(lx1,ly1,lz1,lelv)	m1,
			n1,
			mflg
	)

multd()

subroutine multd	(	real, dimension (lx2*ly2*lz2,lelv)	dx,
		real, dimension (lx1*ly1*lz1,lelv)	x,
		real, dimension (lx2*ly2*lz2,lelv)	rm2,
		real, dimension (lx2*ly2*lz2,lelv)	sm2,
		real, dimension (lx2*ly2*lz2,lelv)	tm2,
			isd,
			iflg
	)

nekuf()

subroutine nekuf	(	real, dimension (lx1,ly1,lz1,lelv)	f1,
		real, dimension (lx1,ly1,lz1,lelv)	f2,
		real, dimension (lx1,ly1,lz1,lelv)	f3
	)

normsc()

subroutine normsc	(	real	h1,
		real	semi,
		real	l2,
		real	linf,
		real, dimension (lx1,ly1,lz1,1)	x,
			imesh
	)

normvc()

subroutine normvc	(	real	h1,
		real	semi,
		real	l2,
		real	linf,
		real, dimension (lx1,ly1,lz1,1)	x1,
		real, dimension (lx1,ly1,lz1,1)	x2,
		real, dimension (lx1,ly1,lz1,1)	x3
	)

opadd2()

subroutine opadd2	(	real, dimension(1)	a1,
		real, dimension(1)	a2,
		real, dimension(1)	a3,
		real, dimension(1)	b1,
		real, dimension(1)	b2,
		real, dimension(1)	b3
	)

opadd2col()

subroutine opadd2col	(	real, dimension(1)	a1,
		real, dimension(1)	a2,
		real, dimension(1)	a3,
		real, dimension(1)	b1,
		real, dimension(1)	b2,
		real, dimension(1)	b3,
		real, dimension(1)	c
	)

opamask()

subroutine opamask	(	real, dimension (lx1,ly1,lz1,1)	vbdry1,
		real, dimension (lx1,ly1,lz1,1)	vbdry2,
		real, dimension (lx1,ly1,lz1,1)	vbdry3
	)

opbinv()

subroutine opbinv	(	real, dimension (1)	out1,
		real, dimension (1)	out2,
		real, dimension (1)	out3,
		real, dimension (1)	inp1,
		real, dimension (1)	inp2,
		real, dimension (1)	inp3,
		real, dimension (1)	h2inv
	)

opbinv1()

subroutine opbinv1	(	real, dimension (1)	out1,
		real, dimension (1)	out2,
		real, dimension (1)	out3,
		real, dimension (1)	inp1,
		real, dimension (1)	inp2,
		real, dimension (1)	inp3,
			SCALE
	)

opchsgn()

subroutine opchsgn	(	real, dimension(1)	a,
		real, dimension(1)	b,
		real, dimension(1)	c
	)

opcmult()

subroutine opcmult	(	real, dimension(1)	a,
		real, dimension(1)	b,
		real, dimension(1)	c,
			const
	)

opcol2()

subroutine opcol2	(	real, dimension(1)	a1,
		real, dimension(1)	a2,
		real, dimension(1)	a3,
		real, dimension(1)	b1,
		real, dimension(1)	b2,
		real, dimension(1)	b3
	)

opcolv()

subroutine opcolv	(	real, dimension(1)	a1,
		real, dimension(1)	a2,
		real, dimension(1)	a3,
		real, dimension(1)	c
	)

opcolv2()

subroutine opcolv2	(	real, dimension(1)	a1,
		real, dimension(1)	a2,
		real, dimension(1)	a3,
		real, dimension(1)	b1,
		real, dimension(1)	b2
	)

opcolv2c()

subroutine opcolv2c	(	real, dimension(1)	a1,
		real, dimension(1)	a2,
		real, dimension(1)	a3,
		real, dimension(1)	b1,
		real, dimension(1)	b2,
			c
	)

opcolv3()

subroutine opcolv3	(	real, dimension(1)	a1,
		real, dimension(1)	a2,
		real, dimension(1)	a3,
		real, dimension(1)	b1,
		real, dimension(1)	b2,
		real, dimension(1)	b3,
		real, dimension (1)	c
	)

opcolv3c()

subroutine opcolv3c	(	real, dimension(1)	a1,
		real, dimension(1)	a2,
		real, dimension(1)	a3,
		real, dimension(1)	b1,
		real, dimension(1)	b2,
		real, dimension(1)	b3,
		real, dimension (1)	c,
			d
	)

opcopy()

subroutine opcopy	(	real, dimension(1)	a1,
		real, dimension(1)	a2,
		real, dimension(1)	a3,
		real, dimension(1)	b1,
		real, dimension(1)	b2,
		real, dimension(1)	b3
	)

opdiv()

subroutine opdiv	(	real, dimension (lx2,ly2,lz2,1)	outfld,
		real, dimension (lx1,ly1,lz1,1)	inpx,
		real, dimension (lx1,ly1,lz1,1)	inpy,
		real, dimension (lx1,ly1,lz1,1)	inpz
	)

opdsop()

subroutine opdsop	(	real, dimension(1)	a,
		real, dimension(1)	b,
		real, dimension(1)	c,
		character*3	op
	)

opdssum()

subroutine opdssum	(	real, dimension(1)	a,
		real, dimension(1)	b,
		real, dimension(1)	c
	)

opgrad()

subroutine opgrad	(	real, dimension (lx2,ly2,lz2,1)	out1,
		real, dimension (lx2,ly2,lz2,1)	out2,
		real, dimension (lx2,ly2,lz2,1)	out3,
		real, dimension (lx1,ly1,lz1,1)	inp
	)

opgradt()

subroutine opgradt	(	real, dimension (lx1,ly1,lz1,1)	outx,
		real, dimension (lx1,ly1,lz1,1)	outy,
		real, dimension (lx1,ly1,lz1,1)	outz,
		real, dimension (lx2,ly2,lz2,1)	inpfld
	)

ophx()

subroutine ophx	(	real, dimension (lx1,ly1,lz1,1)	out1,
		real, dimension (lx1,ly1,lz1,1)	out2,
		real, dimension (lx1,ly1,lz1,1)	out3,
		real, dimension (lx1,ly1,lz1,1)	inp1,
		real, dimension (lx1,ly1,lz1,1)	inp2,
		real, dimension (lx1,ly1,lz1,1)	inp3,
		real, dimension (lx1,ly1,lz1,1)	h1,
		real, dimension (lx1,ly1,lz1,1)	h2
	)

ophyprk()

subroutine ophyprk	(	real, dimension (lx1,ly1,lz1,1)	vel1,
		real, dimension (lx1,ly1,lz1,1)	vel2,
		real, dimension (lx1,ly1,lz1,1)	vel3,
			ilag
	)

opicol2()

subroutine opicol2	(	real, dimension(1)	a1,
		real, dimension(1)	a2,
		real, dimension(1)	a3,
		real, dimension(1)	b1,
		real, dimension(1)	b2,
		real, dimension(1)	b3
	)

opmask()

subroutine opmask	(	real, dimension(1)	res1,
		real, dimension(1)	res2,
		real, dimension(1)	res3
	)

oprone()

subroutine oprone	(	real, dimension(1)	a,
		real, dimension(1)	b,
		real, dimension(1)	c
	)

oprzero()

subroutine oprzero	(	real, dimension(1)	a,
		real, dimension(1)	b,
		real, dimension(1)	c
	)

opsub2()

subroutine opsub2	(	real, dimension(1)	a1,
		real, dimension(1)	a2,
		real, dimension(1)	a3,
		real, dimension(1)	b1,
		real, dimension(1)	b2,
		real, dimension(1)	b3
	)

opsub3()

subroutine opsub3	(	real, dimension(1)	a1,
		real, dimension(1)	a2,
		real, dimension(1)	a3,
		real, dimension(1)	b1,
		real, dimension(1)	b2,
		real, dimension(1)	b3,
		real, dimension(1)	c1,
		real, dimension(1)	c2,
		real, dimension(1)	c3
	)

ortho()

subroutine ortho

(

real, dimension (lx2,ly2,lz2,lelv)

respr

)

plan1()

subroutine plan1

(

igeom

)

rotate_cyc()

subroutine rotate_cyc	(	real, dimension(lx1,ly1,lz1,1)	r1,
		real, dimension(lx1,ly1,lz1,1)	r2,
		real, dimension(lx1,ly1,lz1,1)	r3,
			idir
	)

set_pnd()

subroutine set_pnd	(	real, dimension(n,d)	P,
		real, dimension(0:n-1,d)	LkD,
		real, dimension(n,0:n-1)	LkNt,
		integer	N,
		integer	D
	)

set_pndoi()

subroutine set_pndoi	(	real, dimension(n,d)	Pt,
		real, dimension(d,n)	P,
		real, dimension(n,0:n-1)	LkNt,
		integer	N,
		integer	D
	)

setab3()

subroutine setab3	(		ab0,
			ab1,
			ab2
	)

setabbd()

subroutine setabbd	(	real, dimension(nab)	ab,
		real, dimension(nab)	dtlag,
			nab,
			nbd
	)

setbd()

subroutine setbd	(	real, dimension(1)	bd,
		real, dimension(1)	dtbd,
			nbd
	)

setmap()

subroutine setmap	(		n1,
			nd
	)

setordbd()

subroutine setordbd

setproj()

subroutine setproj	(		n1,
			nd
	)

specmp()

subroutine specmp	(	real, dimension(nb,nb,nb)	b,
			nb,
		real, dimension(na,na,na)	a,
			na,
			ba,
			ab,
		real, dimension(1)	w
	)

tauinit()

subroutine tauinit	(		tau,
			ilag
	)

uzawa()

subroutine uzawa	(	real, dimension (lx2,ly2,lz2,lelv)	rcg,
		real, dimension (lx1,ly1,lz1,lelv)	h1,
		real, dimension (lx1,ly1,lz1,lelv)	h2,
		real, dimension(lx1,ly1,lz1,lelv)	h2inv,
			intype,
			iter
	)

uzprec()

subroutine uzprec	(	real, dimension (lx2,ly2,lz2,lelv)	rpcg,
		real, dimension (lx2,ly2,lz2,lelv)	rcg,
		real, dimension (lx1,ly1,lz1,lelv)	h1m1,
		real, dimension (lx1,ly1,lz1,lelv)	h2m1,
			intype,
		real, dimension (lx2,ly2,lz2,lelv)	wp
	)

velchar()

subroutine velchar	(	real, dimension (lx1,ly1,lz1,lelv)	vel,
		real, dimension (lx1,ly1,lz1,lelv)	vn,
		real, dimension (lx1,ly1,lz1,lelv,9)	vlag,
			nbd,
			tau,
		real, dimension (nbd)	dtbd
	)

velconv()

subroutine velconv	(	real, dimension (lx1,ly1,lz1,lelv)	vxn,
		real, dimension (lx1,ly1,lz1,lelv)	vyn,
		real, dimension (lx1,ly1,lz1,lelv)	vzn,
			tau
	)

velconvv()

subroutine velconvv	(	real, dimension(1)	vxn,
		real, dimension(1)	vyn,
		real, dimension(1)	vzn,
			tau
	)

velinit()

subroutine velinit	(	real, dimension (lx1,ly1,lz1,lelv)	vel1,
		real, dimension (lx1,ly1,lz1,lelv)	vel2,
		real, dimension (lx1,ly1,lz1,lelv)	vel3,
			ilag
	)

wgradm1()

subroutine wgradm1	(	real, dimension(lxyz,1)	ux,
		real, dimension(lxyz,1)	uy,
		real, dimension(lxyz,1)	uz,
		real, dimension(lxyz,1)	u,
			nel
	)

wlaplacian()

subroutine wlaplacian	(	real, dimension(1)	out,
		real, dimension(1)	a,
		real, dimension(1)	diff,
			ifld
	)