急急急！UDF编程出错 请高手指点

巧遇角落2012

用户自定义函数定义边界条件，fluent仿真出的各点压力值控制各点的多相流气体体积分数，请各位大侠指点，万分感谢！！！
#include "udf.h"
#define A 1.5729e-11
#define B 1.7930e-8
#define C 7.6574e-6
#define D 1.4822e-3
#define E 0.1133

DEFINE_PROFILE(v_fration, thread, position)   
{
  cell_t c;
  face_t f;
  begin_f_loop(f, thread)
    {
      Thread *mix_th;
      real p=C_P(c,mix_th);     /*提取各点仿真压力赋值给p*/
      real w=0.549;
      real v;
      if(p<400000)
        {v = A*pow(p,4)-B*pow(p,3)+C*pow(p,2)-D*p+E;
         if(v>w)
            v=w;
         else  if(v<w&&v>0)
                v=v;
               else v=0;
       return v;
        }
         
    }  

  end_f_loop(f, thread)
}

巧遇角落2012

UDF编程出错 请给位高手百忙之中指点一下 看一下哪里有问题
用户自定义函数定义边界条件，fluent仿真出的各点压力值控制各点的多相流气体体积分数，请各位大侠指点，万分感谢！！！
#include "udf.h"
#define A 1.5729e-11
#define B 1.7930e-8
#define C 7.6574e-6
#define D 1.4822e-3
#define E 0.1133

DEFINE_PROFILE(v_fration, thread, position)   
{
  cell_t c;
  face_t f;
  begin_f_loop(f, thread)
    {
      Thread *mix_th;
      real p=C_P(c,mix_th);     /*提取各点仿真压力赋值给p*/
      real w=0.549;
      real v;
      if(p<400000)
        {v = A*pow(p,4)-B*pow(p,3)+C*pow(p,2)-D*p+E;
         if(v>w)
            v=w;
         else  if(v<w&&v>0)
                v=v;
               else v=0;
       return v;
        }
         
    }  

  end_f_loop(f, thread)
}