body = Body(P,Cc,cartesian=False) #move the control points for centerline deltaC_r = np.array([0,2,0,0]) deltaC_x = np.array([0,0,0,1]) deltaCc = np.array(zip(deltaC_x,deltaC_r)) #calculate new P's body.deform(deltaCc) fig = pylab.figure() ax = fig.add_subplot(3,1,1) ax.set_title('Centerline b-sPiine Interpolant') body.plot_spline(ax) ax.legend(loc=2) ax = fig.add_subplot(3,1,3) ax.set_title('Initial and Final Geometries') body.plot_geom(ax) ax.legend() fig = pylab.figure() ax = fig.add_subplot(111, projection='3d') ax.plot_surface(body.Xo,body.Yo,body.Zo,rstride=1, cstride=1, alpha=0.2) #ax.plot_surface(body.Xi,body.Yi,body.Zi,rstride=1, cstride=1, alpha=1.0,color='r') pylab.show()