|
|
DEFINE_ADJUST(my_adjust,domain)
{
Thread *t;
cell_t c;
d=Get_Domain(1);
thread_loop_c(t,domain)
{
if(C_VOF(c,t)>0&&C_VOF(c,t)<1)
{
begin_c_loop(c,t)
{
real X[ND_ND];
C_CENTROID(X,c,t);
real z;
real Zmin=0;
z=X[2];
if(z<Zmin||Zmin==0.) Zmin=z;
C_UDMI(c,t,0)=Zmin;
}
end_c_loop(c,t)
}
}
} |
|