UDF:DEFINE_DPM_BC(AD_SCAF,p,t,f,f_normal,dim)

Adriano

hI, with this little UDF I have to impose this simple BC to a Porous Jump interior surface: 1 - if particle identity is equal to selected number, trap particle; else 2 - apply the 'interior' attribute to my surface, so continue with tracking particle.

the problem is: in case 2 - particles stop their run! It is wrong!!!!! Sigh please any suggestions?!? I'll be crazy very soon!

Many thanx to helper!!

Cheers Adriano

ps: this is the code.

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;

}