请教大虾：二维边界条件上怎样计算y方向边界上的温度梯度积分！！！

lgdyw

请教大虾：二维边界条件上怎样计算y方向边界上的温度梯度积分！！！

lgdyw

#include "udf.h"
DEFINE_PROFILE(th_heatlines,thread,position)
{
cell_t c;
face_t f;
real V;
Thread *tc;      /* Do nothing if areas aren't computed yet or not next to fluid. */
   begin_f_loop(f,thread)
{
c=F_C0(f,thread);/*获得边界处相邻网格C0的cell index*/
   tc=THREAD_T0(thread);/*获得C0的cell thread*/
  VC_T_RG(c,tc)[0]*C_VOLUME(c,tc）;
  F_PROFILE(f,thread,position)=V; /*13*/
}
end_f_loop(f,thread)
}/