UDF定义阻力方向

lwj120

请教用UDF定义阻力方向，需要用到什么命令？下面这段程序是一个定义通过多孔介质阻力方向的程序，看起来像是定义重心的坐标，不是很明白，求看明白下列程序的高手给予指点？
程序1：
#include "udf.h"
#include "unsteady.h"
#include "stdio.h"
#include "stdlib.h"
#include "threads.h"
/*********************************************************/
DEFINE_PROFILE(dir1, t, nv)
{
  real xc[ND_ND],y[ND_ND],X;
   cell_t c;
begin_c_loop(c,t)
{
     C_CENTROID(xc,c,t)
             X=xc[0];
        y[0]=-X/5;
   C_PROFILE(c,t,nv)=y[0];
}
  end_c_loop(c,t)     
}
程序2：
#include "udf.h"
#include "unsteady.h"
#include "stdio.h"
#include "stdlib.h"
#include "threads.h"
/*********************************************************/
DEFINE_PROFILE(dir2, t, nv)
{
  real xc[ND_ND],M[ND_ND],X;
  cell_t c;
begin_c_loop(c,t)
{
     C_CENTROID(xc,c,t)
            X=xc[1];
           M[1]=-X/5;
   C_PROFILE(c,t,nv)=M[1];
}
  end_c_loop(c,t)     
}

fty0083

这个程序定义了阻力方向为 (-x/5, -y/5)，其中点(x,y)表示控制体中心点的坐标。仅此而已。