[cpplabs]

Hi,
In simpleFoam.C, runTime.Write() is used to write results, and simpleFoam p.storePrevIter() is called to store the pressure field of the previous iteration. Because the WOpt_ member in the copy constructor of GeometricField is always set to be NO_Write, and it seems that the write function always export the pressure field for the previous iteration (not p itself), pressue is always missing from the export.
Can someone tell me how to fix the problem or what part of my understanding is wrong?

Thanks.

============================================
// construct as copy
template<class> class PatchField, class GeoMesh>
Foam::GeometricField<type,>::GeometricField
(
const GeometricField<type,>& gf
)
:
DimensionedField<type,>(gf),
timeIndex_(gf.timeIndex()),
field0Ptr_(NULL),
fieldPrevIterPtr_(NULL),
boundaryField_(*this, gf.boundaryField_)
{
if (debug)
{
Info<<>::GeometricField : "
"constructing as copy"
<<>info() << endl;
}

if (gf.field0Ptr_)
{
field0Ptr_ = new GeometricField<type,>
(
*gf.field0Ptr_
);
}

this->writeOpt() = IOobject::NO_WRITE;
}