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)
}