def readData(self):
if self.type != 'vti':
if self.xrange[0] <=0. :
self.data=read.load(self.offset,file=self.filenameout,type=self.type,mirrorPlane=0)
else:
self.data=read.load(self.offset,file=self.filenameout,type=self.type)
else:
self.data=read.loadvti(self.offset,file=self.filenameout)
self.nx=self.data.nx
self.ny=self.data.nx
if self.xrange == None:
self.xrange = np.array(self.data.getBounds()[0:2])
if self.yrange == None:
self.yrange = np.array(self.data.getBounds()[2:4])
self.Lx = self.xrange[1] - self.xrange[0]
self.Ly = self.yrange[1] - self.yrange[0]
def readData(self):
if self.type != 'vti':
if self.xrange[0] <= 0.:
self.data = read.load(self.offset,
file=self.filenameout,
type=self.type,
mirrorPlane=0)
else:
self.data = read.load(self.offset,
file=self.filenameout,
type=self.type)
else:
self.data = read.loadvti(self.offset, file=self.filenameout)
self.nx = self.data.nx
self.ny = self.data.nx
if self.xrange == None:
self.xrange = np.array(self.data.getBounds()[0:2])
if self.yrange == None:
self.yrange = np.array(self.data.getBounds()[2:4])
self.Lx = self.xrange[1] - self.xrange[0]
self.Ly = self.yrange[1] - self.yrange[0]
fontsize=10
dpi=300
filenames = ['output/data1D','output/data3D','output/data3DFCT']
labels = [r'1D', r'3D GLM', r'3D, FCT']
markers = ['ro','gs','bD']
filenameout = 'cuts_'
#slice 3D data at these indices:
ix=16
iy=16
simulations=[]
for filename in filenames:
simulations.append(read.loadvti(1,file=filename))
for (filename,sim) in zip(filenames,simulations):
outarr=np.empty((sim.nx,len(sim.getVarnames())+1))
varnames='x, '
outarr[:,0] = sim.x[:]
i=1
for varname in sim.getVarnames():
varnames=varnames+varname+', '
if sim.ndim==3:
# We have to extract one line in x-direction
exec "outarr[:,i]=sim.%s[:,ix,iy]" % (varname)
else:
exec "outarr[:,i]=sim.%s[:]" % (varname)
i=i+1
