UDF求助

dgf

各位大神们帮我看看这个UDF有问题吗？问题是这样的：下游有1m水深，下游坝坡1:2，坝高6m（如图一）。大坝下游的压力边界条件用UDF设置，程序如下：
#include "udf.h"
#define C2 4905.0     /*C2为一个系数*/
real h2=1.0;           /*h2为下游水深*/
DEFINE_PROFILE(outlet_pressure,thread,position)
{
   real x[ND_ND];       /*this will hold the position vector*/
   face_t f;
   real y;
   begin_f_loop(f,thread)
{
   F_CENTROID(x,f,thread);    /*获取单元中心坐标存放在x数组中*/
   y=x[1];             /*把存在x[1]中的y坐标赋给变量y*/
   if(y<=1.0)
F_PROFILE(f,thread,position)=1000*9.81*(h2-y);
   else
     F_PROFILE(f,thread,position)=0.0;
}
end_f_loop(f,thread);
}
各位牛人帮忙看看，小弟初学fluent，还望各位赐教！