udf求助！！急急！！！！！！

wade19890130

#include "udf.h"
/*profile for x_velocity*/
DEFINE_PROFILE(x_velocity,thread,index)
{
  real x[3];
  real z;
  face_t f;
  begin_f_loop(f,thread)
  {
  F_CENTROID(x,f,thread);
  z=x[2];
  F_PROFILE(f,thread,index)=10*pow(z/10,0.22);
  }
  end_f_loop(f,thread)
}

/*profile for intensity*/
DEFINE_PROFILE(intensity_profile,thread,index)
{
  real x[3];
  real z;
  face_t f;
  begin_f_loop(f,thread)
  {
  F_CENTROID(x,f,thread);
  z=x[2];
  F_PROFILE(f,thread,index)=0.23*pow(10/z,0.22);
  }
  end_f_loop(f,thread)
}
已经编出x-velocity即U1=10*(z/10)^0.22的程序和Iz=（0.23*(z/10)^(-0.22)湍流强度的程序，接下来想编写湍流动能K的程序，我该怎么编呢？？因为U1和Iz都是变量。
其中K=0.5*(U1^2+U2^2+U3^2)*Iz^2，U1=x-velocity,U2=0.86*U1,U3=0.45*U1
怎么编K的程序啊，求大侠指导啊。因为以前使用CFX，刚学FLUENT，得用到大论文啊！！重谢！！！QQ:352516839