初学，请高手指点关于UDF的

spring025

<P>fluent中一个圆平面的动载荷,圆周固定,平面做余弦曲面周期振动（圆心为顶点），下面这个C程序不能通过编译，请高手指点一下。<BR>#include "udf.h"<BR>#include "math.h"<BR>DEFINE_GRID_MOTION(membrane_moving, domain, dt, time, dtime)<BR>{<BR>  Thread *tf=DT_THREAD(dt);<BR>  face_t f;<BR>  Node *v;<BR>  real NV_VEC(axis);<BR>  real NV_VEC(origin);/* NV_VEC(rvec);*/<BR>  real sign, x[3], a, y, c,signx; <BR>  int n;<BR>  float d=0.0008, fre=50;                 /* maxim deformation through all timesteps and frequency */<BR>  </P>
<P>     /*Message("magnitude is: %g\n", sign);<BR>   Message("magnitude of center is: %g\n",r);<BR>  Message("radius is: %g\n",a); */<BR>  /*NV_D(axis, =, 1.0, 1.0, 0.0);   */   <BR>  /*NV_D(origin, =, 0.0, 0.0,0.0);*/   <BR> <BR>  begin_f_loop(f,tf)<BR>   {<BR>     f_node_loop(f, tf, n)<BR>       {<BR>         v=F_NODE(f, tf, n);        /* check nodes */<BR>         c=NODE_X(v); <BR>   y=NODE_Y(v); <BR>          if(NODE_POS_NEED_UPDATE(v))  /* ensure each node only be used once */<BR>           {   <BR>            NODE_POS_UPDATED(v); <BR>   signx=d*cos(sqrt(c*c+y*y));<BR>            x[0]=NODE_X (v); <BR>     x[1]=NODE_Y (v);<BR>            x[2]=signx*sin(2*3.1415926*fre*time);<BR>           /* NV_VV(rvec, =, x, -, origin);     */        <BR>             <BR>            NV_V(NODE_COORD(v), =, x);     /* return new position to node */<BR>           /* Message("Node-coordinate x: %g\n", c);<BR>            Message("Node-coordinate y: %g\n", x[1]);*/<BR>           } <BR>       }<BR>    }<BR>  end_f_loop(f,tf);<BR>}</P>