looping over face threads in domain
 

Dear Users
I have a problem with my udf.
Could you help my to fix it?


/************************************************** *****
faces_1
01.07.2010
************************************************** *****/

#include "udf.h"


DEFINE_ADJUST(radiation_adjust, domain)
{

real x[ND_ND];
Thread *f_thread;
FILE *fout_face;
face_t f;

fout_face = fopen("data1.out", "w");

{

thread_loop_f(f_thread, domain)/* loops over all face threads in a domain*/
{
Message("The file was created: %d\n", fout_face);


begin_f_loop(f, f_thread) /* loops over faces in a face thread */
{
F_CENTROID(x,f,f_thread);
fprintf(fout_face, "%g %g %g %g\n", x[0], x[1], x[2], F_T(f, f_thread));
}
end_f_loop(f, f_thread)

}
fclose(fout_face); /* close the file pointed to by fout */
Message("The file was created: %d\n", fout_face);
}
}