[求助]那位大侠帮俺看看我的udf为何不能用，我是在例子的基础上改的，急，不胜感激

suzhaohe

#include "udf.h"
DEFINE_ON_DEMAND(on_demand_calc)
{
  Domain *d;
  Thread *t;
  /* Integrate dissipation. */
  real sum_dot=0.;
  real x=0.0;
  real y=0.0;
  real z=0.0;
  cell_t c;
  real u=0.0;
  real v=0.0;
  real w=0.0;
d = Get_Domain(1);   
thread_loop_c (t,d)
  {
    begin_c_loop (c,t)
        x=C_T_G(c,t)[0];
        y=C_T_G(c,t)[1];
        z=C_T_G(c,t)[2];
        ND_SET(u, v, w, C_U(c, t), C_V(c, t), C_W(c, t));
        
        sum_dot += ND_DOT(x, y, z, u, v, w)*
       C_VOLUME(c,t);
    end_c_loop (c,t)
   /* y=C_CENTROID(x,c,t)*/
  }
printf("x %g\n", y);
  printf("Volume integral of temperature: %g\n", sum_dot);
}

suzhaohe

说明：我想计算速度与温度梯度的的点积值，那位帮忙非常感激，我闷了一个月了