输出数据 用udf怎么写 用哪个宏？附程序

moplye

我现在想用udf输出两相流界面处的速度和压强，也就是体积分数为0和1之间的界面处的速度和压强。每个时刻都要输出来，不知道该采用哪个宏？输出的结果不要存在case、data里，单独放一个文本中。请高手或者版主指点。
下面是我写的部分程序：由于宏定义错了 导进去没效果。
   #include "udf.h"
DEFINE_PROFILE(r_w_para,t,i)
{
FILE *fd;
//Thread *t;
int curr_ts;
cell_t c;
real p,a,b,x1,y1,ft,ti;
real x[ND_ND];
fd=fopen("data.txt","a+");
curr_ts = N_TIME;

//if (curr_ts%10==0.0)
//{

ft=RP_Get_Real("flow-time");
fprintf(fd,"%f\n",ft);
begin_c_loop(c,t)
  {
  
    ti=C_VOF(c,t);
    if(ti>0. &&ti<1.0)
        { p=C_P(c,t);
          a=C_U(c,t);
          b=C_V(c,t);
      C_CENTROID(x,c,t);
             x1=x[0];
          y1=x[1];

          fprintf(fd,"%f,%f,%f,%f,%f,%f,%f\n",x1,y1,a,b,p,ti);
        }

  }
end_c_loop(c,t)

//}
fclose(fd);
}

[ 本帖最后由 moplye 于 2009-9-10 15:49 编辑 ]