can anyone help me with this udf.i get errors when i compile it.

#include "udf.h" #define hp 2.0

DEFINE_SOURCE (cell_source,cell,thread,dS,eqn) {

real x[ND_ND], x_coord, y_coord,r,p,source;

cell_t c;

begin_c_loop(c, thread) { C_CENTROID(x,c,thread);

x_coord = x[0];

y_coord = x[1];

p = SQR(x_coord) + SQR(y_coord); r = SQRT(p); source = hp/r; } end_c_loop(c,thread) }

add 'return source;' at last row

Your problem is 2D or 3D?

In DEFINE_SOURCE you can't use "begin_c_loop" and you must specify dS[eqn].

Good luck!!!

Giordano Bruno

its a 2D problem

Why did not ou try the function? face_t f thread *thread begin_f_loop(f, thread) {

}end_f_loop Good lucky

i want the centroid of the cell to find the radial distance of each cell from the center of the model.

thank you

Why do you want to define it as a source then? Just use a define on demand and save the values in User Defined Memory.

You can always access the distance with over the memory then C_UDMI(c,t,0).

DEFINE_ON_DEMAND(init_UDS_distance) { Domain *d; Thread *t; cell_t c; real x[ND_ND], x_coord, y_coord,r,p,source; d = Get_Domain(1); t = Lookup_Thread(d, thread_id); /* Loop over all the cells in the given thread */

begin_c_loop (c,t)

{

C_CENTROID(x,c,t);

x_coord = x[0];

y_coord = x[1];

p = SQR(x_coord) + SQR(y_coord);

r = SQRT(p);

C_UDMI(c,t,0) = hp/r;

}

end_c_loop(c,t) }

Where does hp (hp/r) come from ?