c_face_loop的用法

jyy162

我得udf如下：
其中在不加
c_face_loop(cell, mix_th, n)      
           {
          f = C_FACE(cell,mix_th,n);
          tf = C_FACE_THREAD(cell,mix_th,n);
          if(F_T(f,tf)>t0)
          t0=F_T(f,tf);
          }
这一段在后面程序中把t0设为常数时计算很好。
但是加上这段后歘现
Error: C:\water\m-tran2.c: line 28: non-integer subscript expression: unsigned char
不解应该改正那里？求高手指点迷津！！！
#include "udf.h"
#include "metric.h"
#include "mem.h"
DEFINE_MASS_TRANSFER(liq_gas_source,cell,mix_th,from_index,from_species_index,to_index,t
o_species_index)
{
   real m_lg;
   real Na;
   real freq;
   real xc[ND_ND];
   real t0;
   face_t f;
   Thread *tf;
   int n;
   
   real T_SAT = 273.15;
   real Db=1e-3;
   real area=3.9275e-6;
      
   Thread *liq = THREAD_SUB_THREAD(mix_th, from_index);
   Thread *gas = THREAD_SUB_THREAD(mix_th, to_index);
              
   m_lg=0.;
   C_CENTROID(xc,cell,mix_th);
   if (sqrt(ND_SUM(pow(xc[0]-xc[0],2.),pow(xc[1],2.),pow(xc[2],2.))) >0.0048)
     {
          t0=273.15;
          c_face_loop(cell, mix_th, n)      
           {
          f = C_FACE(cell,mix_th,n);
          tf = C_FACE_THREAD(cell,mix_th,n);
          if(F_T(f,tf)>t0)
          t0=F_T(f,tf);
          }
          Na=pow(210*(t0-T_SAT),1.805);
          freq=113.7;
          if(C_T(cell,mix_th) >= 275.15)
              {      
              m_lg=3.412*Na*freq*pow(Db,3)*area/6;   
              }
          else if(C_VOF(cell,gas)>0.0)
              {
              m_lg=3.412*Na*freq*pow(Db,3)*area/6;
              }
        
      }  
    return (m_lg);
}[USECHGFONTE]

cocowei

c_face_loop循环后面好像要加一个结束循环语句end_face_loop

cocowei

不好意思，我错了，当我没说

sdlff

我也要在cell中寻找各个面，能不能用c_face_loop？请高手指点！

jyy162

错误的提示的意思是下标错误，那个量表示这个？请大家给点意见！呵呵！

sdlff

请问在cell的宏中能不能用c_face_loop？例如：DEFINE_MASS_TRANSFER( name, c, mixture_thread, from_phase_index, from_species_index,
to_phase_index, to_species_index) 其中c是cell_t c。
在类似这种宏中能不能用c_face_loop？请高手指点！

a2883524

F_T(f,tf)>t0，这个tf是不是应该指向相的指针而不是面指针？

yfcao

楼主您好 请问你的这个问题解决了吗？因为我也遇到了类似的问题