帮忙看一下我的用来定义inlet velocity 的c文件

lxblxb

这下面是我的c文件。请看一下那里有错。谢谢！
这个c文件是定义一个2维风洞的inlet velocity随高度变化：Vx=20*((y+50)/396.24)^0.4
#include "udf.h"
DEFINE_PROFILE(inlet_x_velocity, thread, index)
{
real x[ND_ND]; /* this will hold the position vector */
real y;
face_t f;
begin_f_loop(f, thread) /* loops over all faces in the thread
passed in the DEFINE macro argument */
{
F_CENTROID(x,f,thread);
y = x[1];
F_PROFILE(f, thread, index) =20.*sqrt((y+50)/396.24);
}
end_f_loop(f, thread)
}[br][br][以下内容由 lxblxb 在 2006年11月29日 06:30am 时添加] [br]
上面c文件中的
F_PROFILE(f, thread, index) =20.*sqrt((y+50)/396.24);
应该改成
F_PROFILE(f, thread, index) =20.*pow((y+50)/396.24,0.4);

bppz

你编译后哪里报错？