from __future__ import division
from tecplot import Tecplot
from GeoMACH.PAM.configurations.conventional import Conventional
if __name__ == '__main__':
name = 'conventional'
aircraft = Conventional()
aircraft.oml0.addVars(['der'])
t = Tecplot()
counter = 100
for c in aircraft.comps.keys():
comp = aircraft.comps[c]
for p in comp.params.keys():
par = comp.params[p]
for j in range(par.P.shape[1]):
for i in range(par.P.shape[0]):
print c, p, i, j
der = aircraft.getDerivatives(c, p, (i,j), clean=False, FD=False)
aircraft.oml0.P0[:,6] = aircraft.oml0.exportPjtn(der)
aircraft.oml0.write2Tec(name+str(counter))
t.importDataSet(name+str(counter),aircraft.oml0.var,True)
t.setTransparency(False)
t.setCamera(66.0214, 133.832, -140.403, -251, 268, 170.3)
t.plotContours(7)
t.writeImage('temp%03d'%(counter))
counter += 1
from __future__ import division
from tecplot import Tecplot
import numpy, time
import sys
from GeoMACH.PAM.configurations.conventional import Conventional
name = 'conventional'
aircraft = Conventional()
aircraft.oml0.addVars(['der'])
t = Tecplot()
counter = 100
for c in aircraft.comps.keys():
comp = aircraft.comps[c]
for p in comp.params.keys():
par = comp.params[p]
for j in range(par.P.shape[1]):
for i in range(par.P.shape[0]):
print c, p, i, j
der = aircraft.getDerivatives(c, p, (i,j), clean=False, FD=False)
aircraft.oml0.P0[:,6] = aircraft.oml0.exportPjtn(der)
aircraft.oml0.write2Tec(name+str(counter))
t.importDataSet(name+str(counter),aircraft.oml0.var,True)
t.setTransparency(False)
t.setCamera(66.0214, 133.832, -140.403, -251, 268, 170.3)
t.plotContours(7)
t.writeImage('temp%03d'%(counter))
counter += 1
aircraft.writeStructure('temp'+str(counter))
counter += 1
n = 50
for i in range(n):
v = i/(n-1)
c['fuse'].setSections(sections=[2,3,4,5], t1L=0.35*v, t2L=1.0-0.35*v)
c['fuse'].props['posy'].set([0.5-0.2*v,0.5,0.5],[0,0.15,1],w=[1.0,0,0],d=[1,0,0])
c['fuse'].props['ry'].set([0.3-0.2*v,0.5,0.5,0.1],[0,0.15,0.75,1.0],w=[0.9985,0,0,0],d=[1,0,0,0])
c['fuse'].props['rz'].set([0.3-0.2*v,0.5,0.5,0.1],[0,0.15,0.75,1.0],w=[0.9985,0,0,0],d=[1,0,0,0])
c['fuse'].props['noseL'] = 0.3-0.2*v
aircraft.computePoints()
aircraft.export.write2Tec('temp'+str(counter))
aircraft.writeStructure('temp'+str(counter))
counter += 1
t = Tecplot()
for i in range(counter):
t.importDataSet('temp'+str(i),True)
t.createMirror(1,165,3)
t.importDataSet('temp'+str(i)+'_str')
t.setTranslucency(1,2617,50)
t.setTranslucency(166,330,1)
# t.setTransparency(False)
t.setRelCameraPosition(-20, 10, 20, -140)
t.writeImage('temp%03d'%(i))
t.runTecplot()
for i in range(counter,counter+30):
os.system('cp temp%03d.png temp%03d.png'%(counter-1,i))
t.makeVideo()
