How to set h as initial condition
 

Hi, everyone,

I have a small problem. As we know, when we give OF initial condition, we usually give temperature and pressure as initial condition, but this is not ok for two phase condition. Is there a way to give pressure and enthalpy as initial condition??

In hPsiThermo.C, the constructor function construct the h_ field, and calculate h at each cell and surface from pressure and temperature. My first thought is make the h_ field must_read and auto_wirte, and do not calculate h at the constructor function. but it doesn't work.

Does anyone have any idea about this problem??

thank you in advance.