2楼: Originally posted by tumu_tu at 2017-09-02 22:05:35
F_PROFILE(f,thread,index)=sumh;没有在面上循环，这样只在一个网格面上定义了值，其他没定义。

3楼: Originally posted by summer同学 at 2017-09-03 08:34:13
感谢你的回复，F_PROFILE(f,thread,index)=sumh目的在于返回循环里面求和的最终值sumh，你的意思是把这一语句放入循环里面吗？...

#include "udf.h"
#define MW_W 0.018
#define pi 3.14159
#define a_c 0.04
#define R 8.314
#define N 185.371
#define L 22.6

DEFINE_PROFILE(hflux, thread, index)
{
        face_t f;
        real sumh = 0.0;
        real A[ND_ND];
        Thread *threadh;
        begin_f_loop(f, thread)
        {
                real Ti;
                Ti = F_T(f, thread);
                F_AREA(A, f, thread);
                sumh += (2 * a_c / (2 - a_c))*sqrt(MW_W / (2 * pi*R))*NV_MAG(A)*(611.2*(pow(2.718, (17.62*(Ti - 273.15) / (Ti - 30.03)))) / sqrt(Ti) - N)*L;
        }
        end_f_loop(f, thread)
        begin_f_loop(f, thread)
        {
                F_PROFILE(f, thread, index) = sumh;
        }
        end_f_loop(f, thread)
}