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