Fluent动网格燃面推移

需要编写，加质燃面退移代码，参考论文：基于动网格的装药燃烧的燃面退移仿真。
已经写了UDF还出现点错误。主要是DEFINE_GRID_MOTION宏的编写。
DEFINE_SOURCE(burn_rate, c, t, dS, eqn)
{
real p = C_P(c, t);
real r = a * pow(p, n);
dS[eqn] = 0.0;
return r;
}
DEFINE_SOURCE(mass_source, c, t, dS, eqn)
{
real p = C_P(c, t);
real rho_g = C_R(c, t);
real r = a * pow(p, n);
real v = -r * (rho_p / rho_g);
dS[eqn] = v * rho_g;
}

DEFINE_GRID_MOTION(dynamic_mesh, domain, dt, time, dtime)