[ynsudo]

Dear all,
I am struggling with the following.
I am using interFoam (version 1.71) and want to obtain the each force-term appered in the code UEqn.H for a certain cell to know which term is large.

Those force terms are

        fvm::ddt(rho, U)
      + fvm::div(rhoPhi, U)
      - fvm::laplacian(muEff, U)
      - (fvc::grad(U) & fvc::grad(muEff))

                (
                    fvc::interpolate(interface.sigmaK())*fvc::snGrad(alpha1)
                  - ghf*fvc::snGrad(rho)
                  - fvc::snGrad(p_rgh)
                ) * mesh.magSf()

What should I add in the file UEqn.H to get those force values?
Any hints are appreciated.
Best regards,
Yu