恳请各路高手帮我看看这个质量源项的定义有什么问题？

仁者无敌

DEFINE_SOURCE(ymass_source,c0,tf,dS,eqn)
{
            cell_t c0;    /*This is the cell in which the source term is to be applied*/
            real source;
            source=Density*r*dtime;
            dS[eqn]=0;
            return source;
}
其中C0是表面附近的一层网格单元，我想在表面附近这层网格中定义质量源项ρ*r*dtime

仁者无敌

高手帮忙啊！这个DEFINE_SOURCE宏能嵌套在DEFINE_GRID_MOTION中吗？先谢过了！！