UDF求助

六月繁星

我依照例子编了一个密度随温度的变化，如下
/*********************************************************************/
/* UDF for specifying a temperature-dependent dencity property */
/*********************************************************************/
#include "udf.h"
DEFINE_PROPERTY(cell_dencity, cell, thread)
{
real dencity;
real temp = C_T(cell, thread);
dencity=1.225* pow(temp, 4.2559)/29369117821;
return dencity;
}
可以读入，但是计算的时候会出现错误：
Updating solution at time level N... done.
  iter continuity x-velocity y-velocity     energy          k    epsilon     time/iter

Error: Floating point error: invalid number

Error Object: ()

Interrupting...
Done.
请问这个问题怎么解决呀？是我的UDF编的不对么？可是我觉得不错呀？请高手指点！

六月繁星

没有人遇到这个问题么？我觉得这个UDF不错呀！为什么就是没法计算下去呢？请高手回答！