1	!
2	!	HEADS.B4S - calculate head loss etc in pipe
3	!
4	!					TRUTH.B4S
5	.DEFINE .TRUE%=-1%
6	.DEFINE .FALSE%=0%
7	!
8	more.data% = .true%
9	!
10	print "RICHARD WITTENOOM AND ASSOCIATES PTY LTD"
11	print
12	print "Pipe Flow Calculations"
13	print \ print
14	input "C (e.g. 120): "; pipe.coeff
15	!
16	WHILE more.data%
17	input "flow (l/sec) "; q
18	input "diam (mm) "; diam
19	input "length (m) "; length
20	!
21	q = q/1000	! to cu. m.
22	!
23	veloc = q/(PI * (diam^2)/4.0E6)
24	K = ((1/(PI * (diam^2)/4.0E6)) / (0.849 * pipe.coeff * ((diam/4.0E3)^0.63)))^1.85
25	hy.grad = K * (q ^ 1.85)
26	head = hy.grad * length
27	!
28	print
29	print "velocity:                "; veloc; " m/sec"
30	print "hydraulic gradient:      "; (hy.grad * 1000); " m/1000 m"
31	print "head loss:               "; head; " m"
32	print
33	input "more "; ans$
34	!
35	more.data% = fn.ans%(ans$)
36	!
37	NEXT
38	!
39	GOTO 32767
40	!	
41	!	FNANSCHK.B4S
42	!
43	DEF fn.ans$ (text$) = edit$ (left$ (text$, 1%), 32%)
44	DEF fn.upper.case$ (text$) = edit$ (text$, 32%)
45	!
46	DEF fn.ans% (text$)
47	IF	((edit$ (left$ (text$, 1%), 32%) = "Y") OR (text$ = "1")) &
	THEN	fn.ans% = -1%						&
	ELSE	fn.ans% =  0%
48	FNEND
49	!
32767	END