各位高手帮忙看看错误吧！

sj51603

#include "udf.h"
DEFINE_PROFILE(unsteady_velocity, thread, position)
{
  face_t f;
begin_f_loop(f, thread)
    {
      real t = RP_Get_Real("flow-time");
  if(t<=1.57)
            F_PROFILE(f, thread, position) =0.2+0.1*sin(4*t);
  else
             F_PROFILE(f, thread, position)=0.2;
    }
  end_f_loop(f, thread)
}
一个小程序，想在入口处定义边界条件，在1，57秒前加一个小的扰动，以后就是常数，但是运行是出错误，各位高手帮忙看一下，小弟谢谢了

午夜流星

没有问题吧

wangwenquanw

real t = RP_Get_Real("flow-time");最好放在前面，且应包括&#35;include "unsteady.h"

invader

只是调整了下位置，提高下代码执行效率，其他貌似没问题，

#include "udf.h"
DEFINE_PROFILE(unsteady_velocity, thread, position)
{
        face_t f;
        real t = RP_Get_Real("flow-time");
        begin_f_loop(f, thread)
        {                  
                if(t<=1.57)
                        F_PROFILE(f, thread, position) =0.2+0.1*sin(4*t);
                else
                        F_PROFILE(f, thread, position)=0.2;
        }
        end_f_loop(f, thread)
}