Merge pull request #45 from ParticulateFlow/postProc

Post proc
This commit is contained in:
tlichtenegger
2018-01-30 08:38:43 +01:00
committed by GitHub
2 changed files with 11 additions and 1 deletions

View File

@ -62,7 +62,7 @@ int main(int argc, char *argv[])
Info<< "\nStarting time loop\n" << endl;
int count=0;
int DEM_dump_Interval=1000;
int DEM_dump_Interval(particleCloud.couplingProperties().lookupOrDefault<int>("dumpInterval",1000));
particleCloud.reAllocArrays();
double **positions_;
@ -109,6 +109,7 @@ int main(int argc, char *argv[])
particleCloud.voidFractionM().setvoidFraction(NULL,voidfractions_,particleWeights_,particleVolumes_,particleV_);
// make sure coupling interval = 1, otherwise update of voidfractionPrev and UsPrev necessary
voidfraction.ref() = particleCloud.voidFractionM().voidFractionInterp();
voidfraction.correctBoundaryConditions();
@ -121,6 +122,9 @@ int main(int argc, char *argv[])
NULL
);
Us = particleCloud.averagingM().UsInterp();
Us.correctBoundaryConditions();
runTime.write();
particleCloud.IOM().dumpDEMdata();

View File

@ -171,6 +171,12 @@ void oneWayVTK::getData
string just_read;
while (just_read.compare(name) != 0) input >> just_read; // read until we read "name"
input >> just_read; // skip text for dataType
// need to distinguish between file formats from dump custom/vtk and from LPP
if(just_read.compare("3") == 0)
{
input >> just_read;
input >> just_read;
}
for (int index = 0; index<particleCloud_.numberOfParticles(); ++index)
{
input >> field[index][0] >> field[index][1] >> field[index][2];