diff --git a/applications/utilities/cfdemPostproc/cfdemPostproc.C b/applications/utilities/cfdemPostproc/cfdemPostproc.C index 90faf24e..a832bacc 100644 --- a/applications/utilities/cfdemPostproc/cfdemPostproc.C +++ b/applications/utilities/cfdemPostproc/cfdemPostproc.C @@ -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("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(); diff --git a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/oneWayVTK/oneWayVTK.C b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/oneWayVTK/oneWayVTK.C index 4a2774df..af104f38 100644 --- a/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/oneWayVTK/oneWayVTK.C +++ b/src/lagrangian/cfdemParticle/subModels/dataExchangeModel/oneWayVTK/oneWayVTK.C @@ -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> field[index][0] >> field[index][1] >> field[index][2];