求助：固体的温度沿着x轴方向改变的udf

zhysh2009

我模型是一个扁的长方体，我想让长方体（该长方体是多孔介质）的温度沿着x轴线性变化，但是编译没问题，但是初始化出现错误是怎么回事？各位大神帮忙看一下~
源程序：
#include "udf.h"
DEFINE_INIT(my_init,d)
{
        cell_t c;
        Thread *t;
        real xc[ND_ND];
        thread_loop_c(t,d)
        {
        if(FLUID_THREAD_P(t)&&THREAD_VAR(t).fluid.porous)
                {
                begin_c_loop(c,t)
                {
                        C_CENTROID(xc,c,t);
                        y=xc[0];

                                C_T(c,t)=1386.-y*0.02;
                }
                                end_c_loop(c,t)
               
                }       
                else C_T(c,t)=300.;
        }

        }

[ 本帖最后由 zhysh2009 于 2011-9-30 10:40 编辑 ]