print s1.radius
s2 = stadium.stadium('thickradius',0.5,0.2)
print s2.thick
print s2.radius
height = 5
#s1.plot(ax,'r')
#s2.plot(ax,'g')
pos = np.array([0,0,0])
direct = np.array([0,0,0])
theta = pos
S0 = stadium.stadiumsolid(densities.Ds[0],pos,direct,'r',s1,s2,height)
S0.draw(ax)
S0.calcLocalProperties()
pos = np.array([0,0,0])
direct = np.array([np.pi/10,0,0])
theta = pos
print np.arange(2)
S1 = stadium.stadiumsolid(densities.Ds[0],pos,direct,'g',s1,s2,height)
S1.draw(ax)
S1.calcLocalProperties()
pos = np.array([0,0,0])
