在做降膜蒸發器的過程用，應用了網上的一個水蒸發的UDF，具體程序如下
#include "udf.h"          /*包括常規宏*/
#include "sg_mphase.h"    /*包括體積分數宏*/
#define T_SAT 373         /*定義飽和溫度*/
#define LAT_HT 1.e3       /*定義水蒸汽潛熱*/

DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn)  /*定義液相源項*/
{
  Thread *mix_th, *sec_th;   /*混合相、第二相 定義計算區域指針*/
  real m_dot_l;             /*蒸發冷凝速率  定義液相質量轉移*/
  mix_th = THREAD_SUPER_THREAD(pri_th);   /*指向混合區的主相即液相的指針*/
  sec_th = THREAD_SUB_THREAD(mix_th, 1); /* 指向單相控制區的氣相的指針*/
  if(C_T(cell, mix_th)>=T_SAT){
    m_dot_l = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
                fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;  /*如果液相單元的溫度高于蒸發溫度，液相向氣相的質量轉移*/
   dS[eqn] = -0.1*C_R(cell, pri_th)*
      fabs(C_T(cell, pri_th) - T_SAT)/T_SAT; /*定義源項對質量轉移偏導*/
                               }
  else {
    m_dot_l = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
                fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    dS[eqn] = 0.; /*于是氣相向液相轉移，所以液相的質量源項對質量轉移的偏導為零*/
       }
   
  return m_dot_l;
}

DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) /*定義氣相源項*/
{
  Thread * mix_th, *pri_th;    /*混合相、第一相*/
  real m_dot_v;
  mix_th = THREAD_SUPER_THREAD(sec_th);   /*指向混合區的第二相即氣相的指針*/
  pri_th = THREAD_SUB_THREAD(mix_th, 0);  /*指向單相控制區的液相的指針，液相為主相*/
  if(C_T(cell, mix_th)>=T_SAT){
    m_dot_v = 0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
                fabs(C_T(cell, mix_th) - T_SAT)/T_SAT; /*如果混合區單元的溫度高于蒸發溫度，液相向氣相的質量質量轉移*/
    dS[eqn] = 0.;    /*由于是液相向氣相轉移，所以氣相的質量源項對來自液相的質量轉移的偏導為零*/
                               }
  else {
    m_dot_v = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
                fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    dS[eqn] = -0.1*C_R(cell, sec_th)*
      fabs(C_T(cell, sec_th) - T_SAT)/T_SAT;  /*由于是氣相向液相轉移，所以氣相的質量源項對自身的質量轉移的偏導不為零*/
       }
   
  return m_dot_v;
}


DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn)   /*混合模型能量源項*/
{
  Thread *pri_th, *sec_th;
  real m_dot;
  pri_th = THREAD_SUB_THREAD(mix_th, 0);   /*指向混合區的液相的指針*/
  sec_th = THREAD_SUB_THREAD(mix_th, 1);   /*指向混合區的氣相的指針*/
  if(C_T(cell, mix_th)>=T_SAT){
    m_dot = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)*
                fabs(C_T(cell, pri_th) - T_SAT)/T_SAT;  /*如果混合區的單元溫度高于蒸發溫度。質量轉移由液相向氣相轉移，吸熱，質量轉移量前有負號*/
    dS[eqn] = -0.1*C_VOF(cell, pri_th)*C_R(cell, pri_th)/T_SAT;
                               }
  else {
    m_dot = 0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)*
                fabs(T_SAT-C_T(cell,mix_th))/T_SAT;
    dS[eqn] = -0.1*C_VOF(cell, sec_th)*C_R(cell, sec_th)/T_SAT;}  /*相反，氣相向液相轉移則放熱*/
  return LAT_HT*m_dot; 
   }                            /*能量大小氣化潛熱與質量轉移率相乘得能量源項*/


在FLUENT中解釋時報錯line 66: label "store_vof_norm" not found (pc=56).請問這個問題應該怎么解決呢???