[求助]一个密度随压力变化的udf。总是出错请指教

zxkzxk

#include "udf.h"
DEFINE_PROPERTY(cell_p,cell,thread)
{
real midu;
real p=C_P(cell,thread);
midu=1.225+(1.+p)/p;
return midu;
}

运行之后只显示在多少个面上有回流，然后就不动了。很久以后显示错误：
Error: Floating point error: invalid number
Error Object: ()
可是如果不用这个udf，密度直接用常数就没有问题

hp270388026

#include "udf.h"
DEFINE_PROPERTY(cell_p,cell,thread)
{
real midu;
real p=C_P(cell,thread);
if(p>1e5) midu=1.225+(1.+p)/p;
else midu=1.225;
return midu;
}