import math 
	
#initialisation of variables
w= 78.9 	#lbf.ft**3 weight
d= 0.01 	#in bore
u= 8.67*10**-9 	#lbf/ hr ft**2 viscosity
h= 18.   	#ft
l= 10. 	    #ft length
	
#CALCULATIONS
Q= math.pi*w*(d/12)**4*(h+l)/(l*128*u)
	
#RESULTS
print  'Flow rate = %.2e ft**3/hr'%(Q)
	
#initialisation of variables
x= 0.1       	#ft
w= 62.4     	#lbf/ft**3 fluid density
v1= 10.      	#ft/sec velocity
u= 2.4*10**-5 	#lbf/ft viscosity
g= 32.2     	#ft/sec**2
k= 4.91
	
#CALCULATIONS
s= k*x*(w*v1*x/(u*g))**-0.5
Tw= 0.332*w*v1**2*(u*g/(w*x*v1))**0.5/g
R= 0.332*6*Tw
	
#RESULTS
print  'Thickness = %.2e*ft'%(s)
print  ' Shear stress = %.3f lbf/ft**2'%(Tw)
print  ' Shear stress = %.3f lbf/ft'%(R)
	
#initialisation of variables
r = 1
r1 = 1
	
#CALCULATIONS
e1= r+r1
e2= r-r1
	
#RESULTS
print  ' vorticity in a forced vortex = %.f*k'%(e1)
print  ' vorticity in a free vortex = %.f'%(e2)