[求助]请教关于UDF中源项的问题

gauss-55

下面是一个简化了的模块，首先想利用lookup来找到靠近边界上的一层流场区域中，本来想用F_c0 来获取，结果出错，再简化成下面的，FT=F_T(f,t);FP=F_P(f, t)+Patm; 编译时提示未初始化！
DEFINE_SOURCE(mass,cell,thread,dS,eqn)
#define Tign=500;
{
real source;
        cell_t c;
real r, mflux;
real c0,FT,FP,FR;
        face_t f;
Thread *t;
          Domain* domain_root=Get_Domain(ROOT_DOMAIN_ID);
          Thread* thread_f4=Lookup_Thread(domain_root,4);
                  FT=F_T(f,t);
                  FP=F_P(f, t)+Patm;           
  if(FT<Tign)
     {
                       source=0;         
                       dS[eqn]=0;
                     }
  
                  else
                     {
                        source=FP;
                        dS[eqn]=0;
                     }
        
                        return source;
}
帮我看看是哪个问题，我主要是想在靠近边界上的一层流体单元添加源项，该如何做？

gauss-55

那位帮忙看一下啊，
对了，还有一个小问题，就是如何获得某个cell的ID阿，是那条命令？欧没查到

liuguangyizju

在可以加源项的区域自己建一个domain。