入口速度满足1/7定律的UDF 求助

trackcom

我在算一个倾斜圆管的温度场跟速度场，其中入口速度想设成是满足1/7定律的情况，用UDF 来实现。
面心的中心坐标是（-0.185803，0.750705，0.816300）。在算的过程中，报错：

Error: > (greater-than): invalid argument [2]: wrong type [not a number]

UDF如下：
#include "udf.h"

#define R 0.15
#define UMEAN 0.02829
#define B 1/7


DEFINE_PROFILE(a_velocity,t,i)      
{
real r,x[ND_ND];
face_t f;



begin_f_loop(f,t)
{
F_CENTROID(x,f,t);
r=sqrt((x[0]+0.185803)*(x[0]+0.185803)+(x[1]-0.750705)*(x[1]-0.750705)+(x[2]-0.816300)*(x[2]-0.816300));  

F_PROFILE(f,t,i)=UMEAN*pow(1-r/R,B);
}
end_f_loop(f,t);
}

有高手能告诉我怎么修改这个UDF吗？
实在是改不了了。。。。
期待牛人。。。。。。。