走过路过的大侠快来看看啊

xiangxiang

DEFINE_PROFILE(time_temprature,thread,position)
{
   face_t f;
   begin_f_loop(f,thread)
{
  real t,tem;
  real t=RP_Get_Real("flow-time");
  if (t<=0.4)
  tem=512.12654+7716.04938*t;
  else
  tem=0;
}
  end_f_loop(f,thread)
}
  上面的程序是某个面的温度随时间的增加程序，请各位看看有什么问题吗？我的计算过程和现实情况不符，不知是不是程序的问题？请各位指点