print 'no. of direntions ', ndir
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
