纠结啊！！！——C_H_G问题！！！

chuanyue

我在能量源项中包含了C_H_G这个宏，结果fluent提示
Error:
FLUENT received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
3. Exit program and restart to continue.
4. Report error to your distributor.
Error Object: ()

我的源项是：
DEFINE_SOURCE(energy_source,c,t,dS,eqn)
{
real x[ND_ND];
real source;
real jr,jx;
C_CENTROID(x,c,t);
jr=-b*C_UDSI_G(c,t,0)[1];
jx=-b*C_UDSI_G(c,t,0)[0];
source=(pow(jx,2)+pow(jr,2))/b
        +(2.5*kB/e)*((jx*C_H_G(c,t)[0]+jr*C_H_G(c,t)[1])/C_CP(c,t))
        -sr;
dS[eqn]=0;
return source;
}

源项编译能通过，但初始化后一计算就出错。源项中各变量我都用常数一一替换运行过，就只有C_H_G这里不行。

大家有谁知道我错在哪里，或者C_H_G这个宏有什么用法要点，能帮小弟的忙，不胜感激涕零！！！