Hi Rob nThanks for letting me know. So how do I define a check for saturation temperature when I am determining it as a function of pressure? Also I tried initializing the domain with mass fraction value of h20 as 1 and absolute pressure as 10KPa Hence there is no possibility of sat_t becoming less than 0 anywhere in the domain. I dont understand from where the solver calculates this value.