[求助]大家帮我看一下UDF错在哪

songxy

&#35;include <stdio.h>
&#35;include "udf.h"
DEFINE_CG_MOTION(vel, dt, cg_vel, cg_omega, time, dtime)
{
Thread *t;
face_t f;
real x[ND_ND];
  real y;
  real y1;  
  real y2;           
  face_t f;           
  begin_f_loop(f, thread)
    {
      F_CENTROID(x,f,thread);
      y = x[1];      
      y1 = y-0.2;
      y2 = y1/0.21;
      F_PROFILE(f, thread, position) =1.49*y2*sqrt(1-y2);}
  end_f_loop(f, thread)
cg_vel[1] = 1.49*y2*sqrt(1-y2);
}

yxg

看不懂想算什么，把要算的物理模型说一下，还有计算过程中出的错误也说一下！