导入表面张力系数的udf文件后报错

py19870531

请问有没有大侠遇到过这种问题
在导入表面张力系数的udf文件后
初始化会有这样的错误
见下图
请问是什么地方出了问题啊
ps：后来发现是这行的问题 real MF=C_YI(c,t,0);//质量浓度
     如果将C_YI(c,t,0)换成其他的就可行了
是不是这个宏里面不能使用C_YI(c,t,0)呢？但是其他的比如密度就可以用啊
或者是因为表面张力hook的地方和一般物性不一样？
有哪位大侠帮忙解释下？谢谢

附上程序
#include "udf.h"
#include <math.h>

#define licl 0
#define water_liquid 1

#define air 0
#define solution 1

DEFINE_PROPERTY(sigma_licl,c,t)
{
     real sigma_licl;//表面张力，单位mN/m，达因/厘米
         real sigma_water;//水的表面张力
         Thread *thread_l = THREAD_SUB_THREAD(t, solution); //液相的指针
         real T=C_T(c,t);//温度，k
     real MF=C_YI(c,t,0);//质量浓度
         
        real a1=2.757115;
     real a2=-12.011299;
         real a3=14.751818;
         real a4=2.443204;
         real a5=-3.147739;//计算licl表面张力的系数

        real a6=235.8;
         real a7=0.625;//公式有误
         real a8=1.256;

    sigma_water=a6*(1-a7*(1-T/647.3))*pow((1-T/647.3),a8)*0.001;
         sigma_licl=sigma_water*(1+a1*MF+a2*MF*T/647.3+a3*MF*pow(T/647.3,2)+a4*pow(MF,2)+a5*pow(MF,3));

        return sigma_water;
}