Boundary_face_geometry
 

its a udf to set a boundary condition for UDS0
whan i choose the UDF to the boundary condition pannel , we see fettal error , please advise me
my UDF is :

DEFINE_PROFILE(wire_profile,t,i)
{
real A[ND_ND];
real ds ;/*distance between the cell centroid and the face centroid */
real es[ND_ND]; /* unit normal vector in the direction from centroid of cell c0 to the face centroid */
real A_by_es;
real dr0[ND_ND];
real y;
face_t f;
cell_t c;
Domain *d;
c = F_C0(f,t);
thread_loop_f(t,d)
{
begin_f_loop(f,t)
{
BOUNDARY_FACE_GEOMETRY(f,t,A,ds,es,A_by_es,dr0);
y = dr0[1];
F_PROFILE(f,t,i) = -epsilon_0*( NV_MAG(C_UDSI_G(c,t,1))-(43.25e+6))/y ;
}
end_f_loop(f,t)
}
}


thanks
IVI