[syavash]

Dear foamers,

I have compiled an external library for inflow turbulence generation b.c. on foam-extend 4.0. The code seems to be successfully compiled, no error in runtime. However, it does not work apparently.

There is an if-clause in the code telling to print some message if the master process is involved (decayingTurbulenceFvPatchVectorField.C):

    if (Pstream::master())
    {
        Pout<<"Starting inflgen, time = "<<this->db().time().elapsedClockTime()<<" s"<<endl;
        forAll(CF, I)
        {
            if (L[I] > 0.0)
            {
                scalar x    = CF[I].x();
                scalar ls   = lspot(L[I]);
                scalar ymin = CF[I].y() - L[I];
                scalar zmin = CF[I].z() - L[I];

                for (label k = 0; k < NUM_VORT; k++)
                {
                    vector v((direction_ > 0) ? x-ls : x+ls, ymin+ranGen.scalar01()*2*L[I], zmin+ranGen.scalar01()*2*L[I]);
                    bool allowed = true;
                    for (SLList<decayingVorton>::const_iterator it = vortons_.begin(); it != vortons_.end(); ++it)
                    {
                        if (mag(v - it().location()) < OVERLAP*ls)
                        {
                            allowed = false;
                            break;
                        }
                    }
                    if (!allowed)
                        continue;
                    else
                        vortons_.insert(decayingVorton(L[I], v, RF[I], (direction_ > 0) ? x+ls : x-ls));
                }
            }
        }

In other words, I shall be able to see the output message in either serial or parallel mode, but I don't see it in the solver log (pimpleFoam).

I have attached the code for people willing to help! I appreciate any hint

Kind regards,
syavash

P.S: The library originally developed for OpenFOAM 2.3.1, I tried to port it to foam-extend.