【求助】版主请进，表面反应速率的UDF编译通过，但是没有计算过程没有发生反应

340440776

#include "udf.h"
#define C1 28.773
#define C2 2233.9
#define C3 22.1117
#define C4 -2
#define C5 0.021347058
#define C6 -3.6571
#define C7 1
#define C8 0.94
#define C9 -0.06
real Tsat(real press)
{
  return C1-C2/(log(press)-C3);
}
real Tweifen(real press)
{
  return C2*pow(log(press)-C3,C4)/press;
}
real Rate1(real press,real temper)
{
return exp(C6*pow(temper/Tsat(press)-C7,C8));
}
real Rate2(real press,real temper)
{
return C6*C8*pow(temper/Tsat(press)-C7,C9);
}
DEFINE_SR_RATE(my_rate1,f,t,r,mw,yi,rr)
{
Thread *t0=t->t0;
cell_t c0=F_C0(f,t);
real temper   = C_T(c0,t0);
real press    = C_P(c0,t0);
real press_d  = (press-C_P_M1(c0,t0))/CURRENT_TIMESTEP;
real temper_d = (temper-C_T_M1(c0,t0))/CURRENT_TIMESTEP;
*rr = C5*Rate1(press,temper)*Rate2(press,temper)*(temper_d/Tsat(press)+temper*Tweifen(press)*press_d*pow(Tsat(press),C4));
}
编译过程通过了，但是在非稳态过程，明明有温度变化，却没有发生反应的迹象。如果需要，我可以把case和data文件传一下，求版主帮忙看看，实在是不知道错在哪里？

[ 本帖最后由 340440776 于 2010-5-9 18:21 编辑 ]

030301029

程序怎么看着那么别扭。