[Kummi]

Hello Foamers,
I'm currently working on a custom solver and there is something I can't figure out. The solver calculates thermal conductivity (K) based on various scalar variables (pi, pR, pT, px) as shown below.

template<class thermo>
inline Foam::scalar Foam::conductivitySolidTransport<thermo>::K //thermal conductivity
(
const scalar T
) const
{
return
((1-pR)*(0.012848335*pow(T,0.5))) + (pT*(0.000496*T)) + (pi*(0.0000000000342*pow(T,3))) + (px*(0.00000000228*pow(T,3)));
}

Since, the thermal conductivity (K) is a field variable (volScalarField), the variable gets extracted and printed in every time folder while simulation runs and could able to plot easily in the post-processing. However, the scalar variables (pi, pR, pT, px) mentioned above are not extracted, since its not a field variable.
How to extract such scalar variables (pi, pR, pT, px) in OpenFOAM ?
My idea was to use an IOobject to write this variable in each time directory, but I can't figure out how.
Kindly someone suggest me some ideas or references.
Thank You