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()
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()
