[aditya.pandare]

Hello
I am trying to simulate cavitation using a UDF for the source term of vapor phase in the mixture model.
I am having trouble accessing the volume fraction of the secondary phase. this is the UDF I wrote. I am new to UDFs, so please help me.

#include "udf.h"

DEFINE_SOURCE(re,c,t,dS,eqn)
{
real source, prsr, alpha, pv, zoneid;
zoneid = 11;
Thread *mixture_thread = Lookup_Thread(domain,zid);
Thread **pt;
pt = THREAD_SUB_THREADS(mixture_thread);
prsr = C_P(c,t);
alpha = C_VOF(c,pt[3]);
pv = -97000;
source = 0.02 * (sqrt((2/3)*(pv - prsr)/998.2)) * (1-alpha) * 998.2 * 0.5542 / 0.0717;
dS[eqn] = - 0.02 * (sqrt((2/3)*(pv - prsr)/998.2)) * 998.2 * 0.5542 / 0.0717;
return source;
}

It gives a parse error for the 2 lines beginning with "Thread"
Any help in understanding the working of the THREAD syntaxes!