[求助]udf问题

wulong2477

我做的是用VOF计算两相流问题，出口设为压力出口，并用udf来定义，我想根据水面波动来得到压力值，但是在中间出现了错误，/************************************************************************/
    /* udfexample.c                                           ***************/
    /*UDF for specifying a steady-state velocity profile boundary condition */
    /*************************************************************************/
     #include "udf.h"
     DEFINE_PROFILE(unsteady_pressure, thread, position)
     {
      real x[ND_ND];
      real y;
      real vof;
      cell_t c;
      Thread **pt;      
      face_t f;
      begin_f_loop(f, thread)
     {
      vof=C_VOF(c,pt[0]);
      if (vof==1)
      {
      F_CENTROID(x,f,thread);
      y = x[1];      
      F_PROFILE(f, thread, position) =9.8*1000*y;
      }
      }
      end_f_loop(f, thread)
      }
编译的时候没有出现错误，我用的是interpreted,但是在我设置我边界条件时，进行all-zones initilize 时，出现了这样的错误
Error:
FLUENT received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
3. Exit program and restart to continue.
4. Report error to your distributor.
Error Object: ()
是不是不能采用获得体积分数来判断？我想在出口处有一定的水位，而且水面是自由波动的，有没有好的建议来一下以得到所期望的的结果？
请高手指点！！谢谢！