[求助] 帮忙修改一下圆柱绕流动网格udf，急 已有2人参与

DEFINE_ON_DEMAND(print_force)
{
   real A[ND_ND];
   face_t f;
   real dfx=0,dfy=0;
   int zone_ID=3;
   Thread*thread_name=Lookup_Thread(domain,zone_ID);
   begin_f_loop(f,thread)
   {
         F_AREA(A,f,thread)
         dfx+=F_P(f,thread)*A[1];
         dfy+=F_P(f,thread)*A[0];
     }
     end_f_loop(f,thread)
     printf("x-force=%f y-force=%f,dfx,dfy)