[duaiduaihu]

the boundary condition is dT/dn=f(t) .
I think we should use DEFINE_PROFILE to difine the boundary condition,
#include "udf.h"
DEFINE_PROFILE(currentdensity,thread,position)
{
begin_f_loop(f,thread)
{
a=f(t)/////suppose we have known f(t);
F_PROFILE(f,thread,position)=a;
}
end_f_loop(f,thread)
}
i think the one with red color defines T=f(t), not dT/dn=f(t);
how to modify F_PROFILE(f,thread,position)=a; to define dT/dn=f(t);

thank every dear friends!