udf求修改

#include "udf.h"
DEFINE_PROPERTY(partical_diameter, cell, thread)
{
cell_t c;
Thread *t;
Thread *cell_thread;
real xc[ND_ND];
real y, u, v, s, a, diam;
{
begin_c_loop(cell,cell_thread)
{
C_CENTROID(xc, cell,cell_thread)
y = xc[1];
}
end_c_loop(cell,cell_thread)
}
u = C_U(cell, thread);
v = C_V(cell, thread);
s = sqrt(u*u + v*v);
a = C_VOF(cell, thread);
diam = 0.0000206*pow(s, 2.0)*pow(a, -4.8) + 0.06701*pow(s - 0.14, 0.294)*pow(y, 0.288)*pow(a, 0.6) - sqrt(0.17*pow(s, 4.0)*pow(a, -9.6) + 790548.61*pow(s - 0.14, 0.784)*pow(y, 0.768)*pow(a, 1.6) + 706.742*pow(s - 0.14, 0.392)*pow(y, 0.384)*pow(s, 2.0)*pow(a, -4.0) - 2.251e-16) / 19360.918;
return diam;
}