我用loop循环获取单元内的压力和温度。
但是标准初始化会卡死
DEFINE_SPECIFIC_HEAT(UDF_cp, T, Tref, h, yi)
{
real Tz, P, cp;
Domain *domain = Get_Domain(1);
cell_t c;
Thread *cthread;
thread_loop_c(cthread, domain)
{
begin_c_loop(c, cthread)
{
Tz = C_T(c, cthread);
P = C_P(c, cthread);
cp = /* 定义比热容与压力温度的函数 */
*h = cp * (T - Tref);
}
end_c_loop(c, cthread)
}
return cp;
}