|
|

楼主 |
发表于 2010-12-19 10:09:55
|
显示全部楼层
#include "udf.h"
DEFINE_SOURCE(cell_T_source, cell, thread, dS, eqn)
{
real source;
real x[ND_ND]; /* this will hold the position vector */
real y;
face_t f;
{
F_CENTROID(x,f,thread);
if (x[0]<.1&&x[0]>-.1)
{
source = -3E6;
dS[eqn] = 0;
}
else
{source=0;dS[eqn] = 0;}
}
return source;
} |
|