[lf.lopez18]

Hi!


I have a custom solver where I export the information of the Lagrangian parcels into a file. For this, I have the following code:

forAllIter(basicSprayCloud, parcels, pIter)
{
    sprayVarsFile << pIter().origId() << "\t";
    for (label j=0; j<3; j++)
    {
        sprayVarsFile << pIter().position().component(j) << "\t";
    }
    for (label j=0; j<3; j++)
    {
        sprayVarsFile << pIter().U().component(j) << "\t";
    }
    sprayVarsFile << pIter().d() << "\t" << pIter().T() << "\t" << pIter().rho() << "\t"
        << pIter().nParticle() << "\t" << pIter().Cp() << "\t" << pIter().mu() << "\t" 
        << pIter().sigma() << "\t" << pIter().age() << "\t" << pIter().Xi() << "\t" 
        << mag(pIter().gradXi()) << "\t" << pIter().elecCh() << "\t" 
        << "\n";
}

where sprayVarsFile is the output file. That code works well for cases with a single processor.



However, for parallelized cases (multiple processors), the program only exports the data of one of the processors, resulting in missing parcels in the output file.



How can I fix this? Is there an alternative that considers the parcels of all the processors?


Best regards,
Luis Lopez.