对流换热边界UDF问题

tianya

我的问题是一个壁面是对流换热边界条件
其中对流换热系数用下式计算 h=17.6*(0.024-0.86*y)^0.32；
Y是沿壁面的高度。程序如下，但是编译不通，哪位高手指点一下，感激不尽。
#include "udf.h"
#include <math.h>
DEFINE_PROFILE(Coefficient,t,i)
{
  real x[ND_ND];             /* this will hold the position vector */

  real y;
  face_t f;
  
  
  begin_f_loop(f,t)
    {
      F_CENTROID(x,f,t);
      y=x[1];
      F_PROFILE(f,t,i)=17.6*pow((0.024-0.86*y),0.32);
    }
  end_f_loop(f,t)
}