二维进出口水蒸气体积分数传递问题

gjmfreeboy

我的模型是二维的，需要将出口的水蒸气体积分数传递回进口，下面是我的UDF，但是初始化后一求解就出错，不知道是什么问题。请问各位，知道是什么问题吗？
DEFINE_ADJUST(yangqi_avg_temp, domain)
{
        real NV_VEC(B);      
        real sum_T_B=0.0,sum_B=0.0;
        face_t f;
        int i=2;
        Thread*f_thread;
        f_thread = Lookup_Thread(domain, WALL1_ID);
        begin_f_loop(f, f_thread)
        {
             F_AREA(B,f,f_thread);
                sum_B+=NV_MAG(B);                  
                sum_T_B+=NV_MAG(B)*F_YI(f, f_thread,i);
              
        }
        end_f_loop(f, f_thread)
         avg_temp3=sum_T_B/sum_B;   
}
      

DEFINE_PROFILE(yangqi_Y,t,position)
{
        face_t f;
        begin_f_loop(f,t)
        {
                F_PROFILE(f,t,position) = avg_temp3;
  }
        end_f_loop(f,t)