udf编写湍流入口边界条件时如何调用雷诺应力？

tianfluid

采用如下代码，compile没问题，但是hook时提示"violate access",请各位大侠执教！感谢。
DEFINE_PROFILE(k_profile, t, i)
{
  real y, del, h, x[ND_ND], ustar; /* variable declarations */
  /*Thread *tc;*/
  cell_t c;
  face_t f;
  begin_f_loop(f, t)
    {
   c=F_C0(f,t);
   /*tc=THREAD_T0(t);*/
                  F_CENTROID(x,f,t);
                  y=x[1];
                  F_PROFILE(f,t,i)=1.5*sqrt(C_RUU(c,t));
    }
  end_f_loop(f, t)
}

gallongallon

include "sg_mem.h"看看