最近在搞反应ODE看了一下代码,做个记录。
在chemistryModel.C中有定义:
forAll(RR_, fieldI)
{
RR_.set
(
fieldI,
new DimensionedField<scalar, volMesh>
(
IOobject
(
"RR." + Y_[fieldI].name(),
mesh.time().timeName(),
mesh,
IOobject::NO_READ,
IOobject::NO_WRITE
),
mesh,
dimensionedScalar("zero", dimMass/dimVolume/dimTime, 0.0)
)
);
}
RR并不是体场。所以在求解器层面,可以定义一个volScalarField,然后把RR的值复制给这个volScalarField然后进行输出,并在paraview查看。
也可以吧上面的NO_WRITE改成AUTO_WRITE,就是不能再paraview查看。