deltaC_x = np.zeros((n_c,))
deltaC_r = np.zeros((n_c,))
deltaC_r[-1] = 10 #second to last element, set to 10
deltaC = np.array(zip(deltaC_x,deltaC_r))
geom.deform(centerbody=deltaC)
print "Run Time: ", time.time()-start_time
start_time = time.time()
geom.writeSTL('new.stl', ascii=False)
print "STL Write Time: ", time.time()-start_time
start_time = time.time()
geom.writeFEPOINT('deform.dat')
print "FEPoint Write Time: ", time.time()-start_time
#geometry.writeSTL('new.stl',body.coords.cartesian,ascii=False)
#centerbody.writeFEPOINT('deform.dat',body.coords.cartesian)
deltaC_c = np.array(zip(deltaC_cx,deltaC_cr))
deltaC_tx = np.zeros((n_c,))
deltaC_tr = np.zeros((n_c,))
# deltaC_tr[:-1] = 1 #second to last element, set to 10
deltaC_t = np.array(zip(deltaC_tx,deltaC_tr))
geom.deform(cowl=(deltaC_c,deltaC_t))
print "Run Time: ", time.time()-start_time
start_time = time.time()
geom.writeSTL('new.stl', ascii=False)
print "STL Write Time: ", time.time()-start_time
start_time = time.time()
geom.writeFEPOINT('new.fepoint')
import pylab as p
profile = geom.project_profile()
for point_set in profile:
X = point_set[:,0]
Y = point_set[:,1]
p.plot(X,Y)
p.show()
