[Foam_User5]

Hi,
Im running a case in rhoPimpleFoam. I have set the initial conditions in 0/file sub directory as follows. Pressure
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 101325;
and then, I have set the boundary field types to each boundary.
For velocity, I have given the inlet, absolute values which is something between 10-20 m/s.
My problem is this, after running the case for several iterations, the pressure values reach to abt 1.9e07, which is huge. the thermal field also seems not to be calculated correctly, however the velocities seems to have captured the flow filed. Does anyone know, what could be the problem.

Another question is this, I have calculated the case in incompressible solvers correctly (since it is a M<0.1case) but cannot find a suitable way to add the thermal filed to it. I have tried thermoFoam application but was unable to run it. Is there any tutorial to help running with thermoFoam?
Thanks