def saveRestartFile(self, outDir, n): saveFileName = outDir + 'checkpoint_' + str(n) + '.hdf5' f = Persistence(saveFileName,'w') f.saveAttribute('globalCount', self.globalCount) f.saveAttribute('globalTime', self.globalTime) ts = self.timeStep tsN1 = self.poptions.timeStepN1 tsN2 = self.poptions.timeStepN2 f.saveAttribute('timeStep', ts) f.saveAttribute('timeStepN1', tsN1) f.saveAttribute('timeStepN2', tsN2) f.saveFluidMeshes(self.fluidMeshes) f.saveSolidMeshes(self.solidMeshes) f.savePlateModel(self.plateFields,self.pmodel,self.solidMeshes) f.saveElectricModel(self.elecFields,self.emodel,self.fluidMeshes) f.close()
DistFunc=esbgk1.getdsf() for i in range(0,ndir): dsfList.append( DistFunc.getField(i) ) if esbgk1options.transient: DistFunc1=esbgk1.getdsf1() for i in range(0,ndir): dsfList.append( DistFunc1.getField(i) ) if esbgk1options.timeDiscretizationOrder > 1: DistFunc2=esbgk1.getdsf2() for i in range(0,ndir): dsfList.append( DistFunc2.getField(i) ) if restartFile is not None: restartFile.readKineticModel(macroFields,esbgk1,meshes,dsfList,ndir) restartFile.close() print "read f,fgamma,macropr from restart file" #esbgk1.initializeMaxwellian() tecplotESBGKEntireDomain.esbgkTecplotEntireDomain(1,meshes,meshes_fluent,options.type,macroFields,"dump.dat") esbgk1.ComputeMacroparameters() if (fgamma==0): esbgk1.initializeMaxwellianEq() else: esbgk1.EquilibriumDistributionBGK() if(fgamma==2): esbgk1.EquilibriumDistributionESBGK() #collision frequency based on Prandlt esbgk1.ComputeCollisionfrequency()