Problem with UDF in VOF model

vijay

Hi,

I am doing simulation using VOF model. I am using a UDF to add a source term in y-momentum equation. But when I click "iterate" it is giving some error. Can anybody help me in identifying the mistake in this UDF.

DEFINE_SOURCE(cell_y_source, cell, thread, dS, eqn) { real source;

Thread *t_phase2 = THREAD_SUB_THREAD(THREAD_T0(thread),1);

source=0.0; if (C_VOF(cell,t_phase2)<1)

{

source = 0.0;

dS[eqn] = 0.0;

} else

/* source term */

source = C_R(cell,thread)*29.43;

/* derivative of source term w.r.t. x-velocity. */

dS[eqn]=0.0;

return source; }