各位udf高手请进来看看，救救我吧！

linken

我的热源是一个移动的高斯热源 加载在一个表面！儿热源的加载方式是脉冲式的 及在某时间作用 在另一段时间 不作用！如此循环，按照这种思路我编了一个udf程序但是 编译连接的时候能成功 初始化就出问题了！希望各位给予指点！附udf：（我的udf的思路是在小于T0时间加载热源在T0到T这段时间不加载，如此往复）
#include "udf.h"
DEFINE_PROFILE(temperature_profile, thread, position)
{
real m[ND_ND];
real x;
real y;
face_t f;
real x0,y0,R;
real v=0.00167;
real T=0.03;
real T0=0.006;
real r=0.0003;
real J=1.5;
real t=RP_Get_Real("flow-time");
real i=0;
do
{
if(t<=i*T+T0)
  {
    begin_f_loop(f,thread)
   {
    F_CENTROID(m,f,thread);
    x=m[0];
    y=m[1];
    x0=v*t;
    y0=0;
    R=sqrt(x-x0)*(x-x0)+(y-y0)*(y-y0);
    if(R<=r)
    F_PROFILE(f,thread,position)=3*1.5*0.8/(3.14159*0.006*r*r)*exp(-3*R*R/r*r);
    else
    F_PROFILE(f, thread, position)=0;
   }
   end_f_loop(f, thread)
   printf("a");
   printf("%f",t);
   }
else if(T+i*T>t>i*T+T0)
  {
  begin_f_loop(f,thread)
   {
    F_CENTROID(m,f,thread);
    x=m[0];
    y=m[1];
    F_PROFILE(f, thread, position)=0;
   }
   end_f_loop(f, thread)
   printf("2");
   printf("%f",t);
   }
else
{
i++;
printf("3");
printf("%f",i);
}
}while(i<=10);
}

[ 本帖最后由 linken 于 2010-12-13 17:46 编辑 ]