请教UDF来设定POROSITY随温度的变化，急

vip_lws

各位高手，我想请教一个UDF的问题，
我需要通过UDF来设定一个区域，具体是一个FACE，这个面是POROUS的，我想通过UDF来设定porosity（空隙率）与温度的关系。在边界条件中设定。
我想知道如何通过DEFINE_PROFILE来设定。
这是我写的，想请问大家，问题出在哪里。
#include <udf.h>
#include <math.h>
DEFINE_PROFILE(porosity,thread,position)
{
  face_t f;
real temp=C_T(c,t);
begin_f_loop(f,thread)
{
  F_PROFILE(f,thread,position)=max(0.1,-3.43958+0.024*(temp+273)+5.20e-5*pow((273+temp),2)+4.52e-8*pow((temp+273),3)+1.36e-11*pow((temp+273),4));
/*这是一个空隙率和温度的关系式*/
}
  end_f_loop(f,thread)
}

vip_lws

自己顶一下，希望大家能够帮我解答，我是一个CFD的新手。希望能和大家多交流。
我也想知道，从什么地方学习CFD的UDF比较好。