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])