请问电场力作为动量源项进行UDF编程对吗？

Xiah

各位老师们好，我是用MHD模型加载了电场，只考虑电场力和界面张力的情况下，将麦克斯韦应力作为电场力加入到动量源项，这是我自己编写的程序，可以连接，但计算不起来，请问哪里需要修改呢？
#include "udf.h"
DEFINE_SOURCE(maxwell_source,c,t,dS,eqn)
{       
real x[ND_ND];
real E,a0,a1,a2,b1,b2,r;
real T1,T2,source;
C_CENTROID(x,c,t);
T1=(9*a0*a2*(b1/b2-1)*E*E)/(2*(a1/a2+2)*(a1/a2+2));
T2=9*a0*a2*((a1/a2)*(a1/a2)-2*(b1/b2)+1)*E*E*x[1]*x[1]/(r*r*2*(a1/a2+2)*(a1/a2+2));
source=T1+T2;
dS[eqn]=9*a0*a2*((a1/a2)*(a1/a2)-2*(b1/b2)+1)*E*E*x[1]/(r*r*(a1/a2+2)*(a1/a2+2));
return source;
}

jiangfan2008

计算时有什么提示没？