fluent用户自定义函数求助！！！

雨季不再来

我现在在正在做一个模型，需要对一个二维圆筒面的外边界条件进行设置，外边界条件热流密度是一个极坐标下关于方位的函数q=600*cosα，边界面任一点的位置是（r，α）请高手指导如何用udf编译

雨季不再来

这是我做的一个udf编译，但总是提醒出现错误，请高手指导一下
#include "udf.h"
DEFINE_PROFILE

(WALL_cemian,thread,position)
{
face_t f;
real x[ND_ND];
real y;
begin_f_loop(f,thread)
{
F_CENTROID(x,f,thread);
y=x[1];
F_PROFILE(f,thread,position)=600*sqrt

(double(1-(y/x)*(y/x)));
}
end_f_loop(f,thread);
}

[ 本帖最后由 雨季不再来 于 2014-1-10 15:04 编辑 ]

aztpb

F_PROFILE(f,thread,position)=600*sqrt(double(1-(y/x)*(y/x)));

Y是实数，x是实数数组，怎么相除？？？