[ali-group]

Hi Foamers
I am a bit confused with using this codes for accessing Cp.
In createFields:
PtrList<volScalarField>& cpi = Y;
forAll( cpi, rowI){
cpi[rowI].setSize(mesh.cells().size(), 0);
}
In YEqn.H:
forAll(T, cellI)
{
forAll(Y, specieI)
{
cpi[specieI][cellI] = composition.Cp(specieI,p[cellI],T[cellI]);
}
}

They work well for accessing cp but the whole code crash.
I cant understand what is the problem. If i add these codes to my solver without any other changes for using Cp inside the code, my solver crash.

Is this function changing something else??

Thank you