之前在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"))
{}