[petr.f.]

Dear Foamers,

I've programmed a custom boundary condition which prescribes velocity according to pressure over patch. In the code I'm accessing all the fields in the following manner :

const fvPatchField<scalar> & pp =
patch().lookupPatchField<volScalarField, scalar>(pName_);

scalar patchArea = gSum(patch().magSf());

avgP = gSum(patch().magSf()*pp)/patchArea;

However, if I run the computation in parallel, it blows up in the first iteration (pressure values are nan).

Do you have any idea, what might be the cause of these problems?

Thanks...