The obove was extracted from hThermo.C which I think is responsible for writing the value into the Cp file after a writeInterval. So by referencing other variable like p or T. I changed IOobject::NO_WRITE into IOobject::AUTO_WRITE. But still not saw Cp. How do I change the code correctly? Thank you so much!

The current situation is that after running a case, I intend to use postprocessing method like calcMassFlow or something like that to output the Cp value as well as h. My idea is that to read the Cp value from for example "0.003" directory. But the prerequisite must be that I write the Cp data into a file called Cp the same as p, T, U and so on. I also took a look at the basicThermo.C and .H, but I cannot move forward now, can someone give me a hand? Thanks!