求大神指导物性随两个变量改变的UDF怎么写啊

#include "udf.h"

DEFINE_PROPERTY(cell_viscosity ,cell ,thread)

{

real mu_lam;

real temp = C_T(cell,thread);

real density = 1011.63-0.2205194*temp-1.92249*pow(10,-5)*pow(temp,2)+5.63788*pow(10,-9)*pow(temp,3);

if (temp > 773)

mu_lam = 0.0851*pow(10,-4)*pow(density,1/3)*exp((1.04*density)/temp);

else

mu_lam = 0.1235*pow(10,-4)*pow(density,1/3)*exp((1.04*density)/temp);

return mu_lam;

}

麻烦大神帮忙看下哪里不对？图是物性的公式，密度本身就是温度的函数，粘度又是密度和温度的函数。下面是我写的UDF。