UDF求解并在控制台打印每个迭代步的边界面压力问题

Firehole

自己编制了一个UDF来求解边界面上的压力，并在控制台打印出来。但是自己验算的时候发现是错误了，而将F_P改为F_W后计算速度时，速度的结果又是正确的。求高手指导下，UDF是哪里出的问题
#include "udf.h"
real NV_VEC(A);
real p;
real sum_P_A=0.0;
real sum_A=0.0;
Thread *thread_in;
face_t f;
Domain *domain;
DEFINE_ADJUST(inlet,d)
{
domain=Get_Domain(1);
thread_in=Lookup_Thread(domain,5);
begin_f_loop(f,thread_in)
{
  F_AREA(A,f,thread_in);
  sum_A+=NV_MAG(A);
  sum_P_A+=NV_MAG(A)*F_P(f,thread_in);
}
end_f_loop(f,thread_in);
p=sum_P_A/sum_A;
printf("P is:%g\n",p);
}