Currently I'm using sonicTurbF
 

Currently I'm using sonicTurbFoam for my simulation. I'm trying to implement the heat capacity Cp which is temperature dependent as an input variable for the boundary condition. Now the case in sonicTurbFoam is that Cp is treated as a constant in thermophysicalProperties. Could anybody give some suggestion how to implement the temperature-dependent Cp? Thank you so much!

I guess the temperature-depend
 

I guess the temperature-dependent Cp should be written by myself. If I had to write the code for it, is there some code already existed so that I can reference? Thanks!