源项DEFINE_SOURCE怎样加入速度和面积的内积 MHD

ajene

我这段时间想做磁场的UDF，源项中矢量的内积，即在源项中加入速度和表面积矢量的点积值。我将UDF写好了，能通过编译，但是一旦加载就出现读取错误（ACCESS VIOLATION），不知是什么原因，我现在被卡了很长时间了，麻烦您帮忙看看，希望得到你的指点！谢谢了！ UDF如下：
DEFINE_SOURCE(C_a, cell, t, dS, eqn)
{
real source;
real rho=6000;
int n;
face_t f;
Thread *tf;
real V[ND_ND],A[ND_ND];

     c_face_loop(cell,t,n)

     {
                 f=C_FACE(cell,t,n);
                 tf=C_FACE_THREAD(cell,t,n);           
                 NV_DS(V, =, C_U(cell,t), C_V(cell,t),C_W(cell,t), *,rho);     
               
                 F_AREA(A,f,tf);
                 source=NV_DOT(V,A);
           }
    dS[eqn]=0;
   return source;

}