!
	!	HEADS.B4S - calculate head loss etc in pipe
	!
	{BL:[101,5]TRUTH.B4S}
	!
	more.data% = .true%
	!
	print "RICHARD WITTENOOM AND ASSOCIATES PTY LTD"
	print
	print "Pipe Flow Calculations"
	print \ print
	input "C (e.g. 120): "; pipe.coeff
	!
	WHILE more.data%
	input "flow (l/sec) "; q
	input "diam (mm) "; diam
	input "length (m) "; length
	!
	q = q/1000	! to cu. m.
	!
	veloc = q/(PI * (diam^2)/4.0E6)
	K = ((1/(PI * (diam^2)/4.0E6)) / (0.849 * pipe.coeff * ((diam/4.0E3)^0.63)))^1.85
	hy.grad = K * (q ^ 1.85)
	head = hy.grad * length
	!
	print
	print "velocity:                "; veloc; " m/sec"
	print "hydraulic gradient:      "; (hy.grad * 1000); " m/1000 m"
	print "head loss:               "; head; " m"
	print
	input "more "; ans$
	!
	more.data% = fn.ans%(ans$)
	!
	NEXT
	!
	GOTO 32767
	!	
	{BL:[101,5]FNANSCHK.B4S}
	!
32767	END