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