How do I get density
 

Could some one please help me out.
I need to get a density in a cell. I am using rhoPimpleFoam and the code is

forAll(mesh.C(), celli)
{
scalar rro=rho[celli];
.....
}

I am getting an error. "error: invalid types ‘float[Foam::label {aka int}]’ for array subscript"

I spend two day tryng to understand why I can not get density.

Please help me out.

Hi,

have a look at hPsiThermo.C

Kind Regards,
Christian

Thank you Christian for your reply, but it did not help me. I looked to this file and it did not tell me how to get density.

Could you please describe what I should check in this file.

Hi,

get the pressure and psi of the cell (as seen in class function calculate) and multiply them (have a look at hPsiThermo.H)

Best Regards,
Christian

One more question.

My rho[celli]=psi[celli]*p[celli]

I can get p[celli] but when I try to get psi[celli] I am getting error.

How can I get a p[celli] ???

Thank you one more time.

QUOTE=Chris Lucas;365934]Hi,

get the pressure and psi of the cell (as seen in class function calculate) and multiply them (have a look at hPsiThermo.H)

Best Regards,
Christian[/QUOTE]

One more question.

My rho[celli]=psi[celli]*p[celli]

I can get p[celli] but when I try to get psi[celli] I am getting error.

How can I get a p[celli] ???

Thank you one more time.

QUOTE=Chris Lucas;365934]Hi,

get the pressure and psi of the cell (as seen in class function calculate) and multiply them (have a look at hPsiThermo.H)

Best Regards,
Christian[/QUOTE]