该udf文件出错在哪里

fastbb

/*******************************************************/
/* inlet velocity.c                                                              */
/* UDF for specifying inlet velocity profile condition */
/*******************************************************/

#include "udf.h"

#define aa1 -2395.91088
#define bb1 75.804571
#define cc1 14.76751972
#define dd1 4.0085397
#define ee1 0.0122174
#define ff1 0.00159555


DEFINE_PROFILE(inlet_velocity,thread,np)
{
  real x[ND_ND];
  real y;
  cell_t c;

  begin_c_loop(c,thread)
    {
          C_CENTROID(x,c,thread);
          y=x[1];
          F_PROFILE(c,thread,np) =
aa1*pow(y,5.0)+bb1*pow(y,4.0)+cc1*pow(y,3.0)+dd1*pow(y,2.0)+ee1*y+ff1;
    }
  end_c_loop(c,thread)
}



DEFINE_PROPERTY(cell_density1,cell,thread)
{
  real den_fluid3;

den_fluid3 =1022.0;

return den_fluid3;


有没有高手帮忙指点一下







}

topshen

DEFINE_PROFILE(inlet_velocity,thread,np)
{
  real x[ND_ND];
  real y;
  cell_t c;

  begin_c_loop(c,thread)
    {
          C_CENTROID(x,c,thread);
          y=x[1];
          F_PROFILE(c,thread,np) =
aa1*pow(y,5.0)+bb1*pow(y,4.0)+cc1*pow(y,3.0)+dd1*pow(y,2.0)+ee1*y+ff1;
    }
  end_c_loop(c,thread)
}

改成：

DEFINE_PROFILE(inlet_velocity,thread,np)
{
  real x[ND_ND];
  real y;
  face_t c;

  begin_f_loop(c,thread)
    {
          F_CENTROID(x,c,thread);
          y=x[1];
          F_PROFILE(c,thread,np) =
aa1*pow(y,5.0)+bb1*pow(y,4.0)+cc1*pow(y,3.0)+dd1*pow(y,2.0)+ee1*y+ff1;
    }
  end_f_loop(c,thread)
}

fastbb

谢谢楼上的回答，我试过了，好像行不通，总是提示line1：syntax error

gbd20022002

请仔细检查你的 程序，你可以选择将以下代码直接复制，进行编译，我编译过，没问题。
#include "udf.h"

#define aa1 -2395.91088
#define bb1 75.804571
#define cc1 14.76751972
#define dd1 4.0085397
#define ee1 0.0122174
#define ff1 0.00159555

DEFINE_PROFILE(inlet_velocity,thread,np)
{
  real x[ND_ND];
  real y;
  face_t f;

  begin_f_loop(f,thread)
    {
          F_CENTROID(x,f,thread);
          y=x[1];
          F_PROFILE(f,thread,np)=aa1*pow(y,5.0)+bb1*pow(y,4.0)+cc1*pow(y,3.0)+dd1*pow(y,2.0)+ee1*y+ff1;
    }
  end_f_loop(f,thread)
}

topshen

最近line1：syntax error问题好多喔。呵呵呵