各位大佬大家好，我最近准备在rhoSimpleFoam的基础上加入自己的算法，算法中比较关键的点是有另一组速度、压力和温度场（暂且叫它们Uf，pf，Tf）。
之前在simpleFoam中顺利地植入了，但是在可压缩这里我发现：thermo.rho()这些函数在调用时，读取的物理量是求解器原来的U、p，如果我再另外定义Uf，pf的话，好像不能直接使用thermo.rho()函数进行另一组速度压力场的密度修正？
附上一段src/thermophysicalModels/basic/fluidThermo/fluidThermo.C中的代码，麻烦大佬们看下是不是这样的
Foam::fluidThermo::implementation::implementation
(
const fvMesh& mesh,
const word& phaseName
)
:
p_(lookupOrConstruct(mesh, "p"))
{}
Foam::fluidThermo::implementation::implementation
(
const fvMesh& mesh,
const dictionary& dict,
const word& phaseName
)
:
p_(lookupOrConstruct(mesh, "p"))
{}
