F_FLUX是否是边界的质量流量

ltpreason

real NV_VEC(A);
thread_in=Lookup_Thread(domain,7);/*7是入口面的id*/
begin_f_loop(f,thread_in)
        {
               F_AREA(A,f,thread_in);
               sum_A+=NV_MAG(A);                 
               sum_M_A+=NV_MAG(A)*F_FLUX(f,thread_in);
        }
    end_f_loop(f,thread_in)
ave_M=sum_M_A/sum_A;
这样得出的ave_M和用fluent的surface report里得到的mass flow rate差的特别多，请各位帮忙看看哪里有问题，先谢过了！