[vekh]

Hello. I have problem with my UDF for 2d case.
I want to get phase level at certain area/volume to define boundary between phases. I wrote simple UDF:

#include <udf.h>
# define r 0.0004
# define x0 0.002
# define y0 0.002
FILE *fout;

DEFINE_SOURCE(my_source, cell, thread, dS, eqn) 
{ 
  real x[ND_ND]; 
  real source, xc, yc, prom;

  Thread *tm = THREAD_SUPER_THREAD(thread); 
  Thread **pt = THREAD_SUB_THREADS(tm); 

  xc=x[0];
  yc=x[1]; 
  prom = sqrt(pow((xc-(x0)),2));

  if ((prom <= r)&&(C_VOF(cell, pt[1])<1))
  {
    fout=fopen("test.dat","a+");
    fprintf(fout,"%f\n",xc);
    fclose(fout);  
  }

source=1;
dS[eqn] =0; 
return source; 
}

and I've compiled it. But after first iteration Fluent returns:

Error: received a fatal signal (Segmentation fault).

Function my_source if hooked at Cell Zone/Source Terms/Energy.
I'm quite new to all VoF things so I need some help to start. Thanks in advance!

[`e`]

Move the line of code "FILE *fout;" to within the DEFINE_SOURCE macro. If the error persists, try simplifying your code until you determine the problem line (or print text to the screen using Message(); every few lines for debugging).