Hi, dear all,

Somebody help me to check the following UDF code and Error message.

/************************************/

/*UDF for fixed flux on the boundary*/

#include "udf.h"

#include "sg.h"

DEFINE_PROFILE(secretion, thread, position)

{

/* real y;*/

real k=1.8;

face_t f;

real dz;

real A[ND_ND], At;

real dr0[ND_ND], es[ND_ND], ds, A_by_es;

/* Thread *t0=thread->t0;*/

begin_f_loop(f, thread)

{

cell_t c0 = F_C0(f,thread);

BOUNDARY_FACE_GEOMETRY(f,thread,A,ds,es,A_by_es,dr 0);

dz=ds;

/*dz=0.00001;*/

y=C_YI(c0, thread, position);

F_PROFILE(f, thread, position)=k*dz+y;

end_f_loop(f, thread)

}

ERROR message is:

Error: udfbottom.c: line 21(which is the line: BOUNDARY_FACE_GEOMETRY(f,thread,A,ds,es,A_by_es,dr 0): structure reference not implemented

Thank you very much for your attention.

Regards

Steven

can you tell us what you are trying to do with this UDF ?

I am trying to apply a Neumann boundary condition (the value of variable's gradient at the boudary is knwon) to one face of the domain. BOUNDARY_FACE_GEOMETRY was used to get the distance beween the cell centroid and cell boudary face.