DEFINE_DPM_BC(capture_bounce_bc, p, t, f, f_normal, dim)
{
real deposit_mass;
Thread *t, *t0;
cell_t c,c0;
c0=F_C0(f,t);
t0=THREAD_T0(t);
deposit_mass=0.0;
begin_f_loop(c,t)
{
deposit_mass=deposit_mass+P_MASS0(p)*TIMESTEP_SIZE;
C_UDMI(c0,t0,0)=deposit_mass;
printf("Deposit_mass=%18.16lf\n",C_UDMI(c0,t0,0));
return PATH_END;
}
end_f_loop(c,t)
} |