|
|
还没有人注意啊？
不知道这样行不行？求指教啊
DEFINE_ADJUST(coor_adjust,domain)
{
Thread *t;
cell_t c;
thread_loop_c(t,domain)
{
real x[ND_ND];
real xi,yi,zi;
face_t f;
real cs;
cs=3/pow(R0,2);
begin_c_loop(c,t)
{
real temp=C_T(c,t);
if(temp=3375)
{
C_CENTROID(x,c,t);
xi=x[0];
yi=x[1];
zi=x[2];
C_UDMI(c,t,1)=cs/(3.14159*H*(1-1/exp(3)))*exp(-(pow(xi,2)+pow(yi,2))*cs/log10(H/zi+0.0001));
}
else C_UDMI(c,t,1)=0;
}
end_c_loop(c,t)
}
}
DEFINE_SOURCE(energy_source,c,t,dS,eqn)
{
real source;
real xi,yi,zi,x[ND_ND];
C_CENTROID(x,c,t);
xi=x[0];
yi=x[1];
zi=x[2];
source=C_UDMI(c,t,1);
dS[eqn]=0;
return source;
} |
|
回复此楼
seazsh
