求助：fluent 子程序编写-热源和温度耦合问题！

lxfsure

一个传热问题，两种材料组成的一个物体，其中一种材料是一个体热源，体热源随；另一种材料的某个边界的平均温度变化，函数是
source=1000000*(1-0.15*(1-0.0045*(avg_temp-298.0)))
热源受边界平均温度影响，同时热源也影响边界平均温度，是一个耦合的关系。按照一些例子，编了以下程序，但是热源和温度没有耦合在一起，不知哪里出了问题。
#include "udf.h"
real NV_VEC(A);
real avg_temp;
real sum_T_A=0.0;
real sum_A=0.0;
Thread *thread_out;
face_t f;
Domain *domain;

DEFINE_ADJUST(adjust,d)
{
  domain=Get_Domain(2);
  thread_out=Lookup_Thread(domain,6);
  begin_f_loop(f,thread_out)
    {
        F_AREA(A,f,thread_out);
        sum_A+=NV_MAG(A);                 
        sum_T_A+=NV_MAG(A)*F_T(f, thread_out);
     }
      end_f_loop(f,thread_out)
   avg_temp=sum_T_A/sum_A;
}

DEFINE_SOURCE(cell_source,c,t,ds,eqn)
{
real source;
begin_f_loop(f,t)
      {
       source=1000000*(1-0.15*(1-0.0045*(avg_temp-298.0)));
       ds[eqn]=0.0;
       return source;}
    end_f_loop(f,t)

}