关于dpm边界的udf ,请教.

小旭

通过这个udf在多孔板(过滤板)上给颗粒加简单边界条件:1-捕获颗粒.2-颗粒继续跟踪.问题是2中颗粒也停止了.帮别人问.谢谢关注.
DEFINE_DPM_BC(ADESIONE_SCAF,p,t,f,f_normal,dim)
{
int select,l;
FILE *fp;
/*TO call function that makes the comparison */
select=confronta(p->part_id, a);
if (select==1)
{
Message("select = %d part_id ferma = %d\n",select,p->part_id);
fp=fopen("ADESIONE_SCAF.txt","w");
fprintf(fp,"((%10.6g %10.6g %10.6g %10d %10d %10.6g %10.6g %10.6g) %d)\n",p->state.pos[0],p->state.pos[1],p->state.pos[2],t,f,f_normal[0],f_normal[1],f_normal[2],p->part_id);
fclose(fp);
Trap_Particle(p);
l=PATH_ABORT;
}
if (select==0)
{
l=PATH_ACTIVE;
}
return l;
}