[Aj_45]

I want to calculate the liquid phase volume of a droplet in a 2D axisymmetric domain. I have written following UDF for it. But I get segmentation fault error after the "I am here" message at the end of a time step.

DEFINE_EXECUTE_AT_END(execute_at_end)
{
Domain *domain;
Thread *t;
Thread **pt;

real volume = 0.;
cell_t c;
domain = Get_Domain(1);

mp_thread_loop_c(t, domain, pt)
Message("I am here");
if (FLUID_THREAD_P(t))
Message("Hey there!");
{
begin_c_loop(c, t)
{
volume += 2*3.1415*C_VOLUME(c, t)*C_VOF(c, pt[1]);
}
end_c_loop(c, t)
}
Message("Volume = %f", volume);
}

How do I make sure a which thread is being passed?
Thanks!

[AlexanderZ]

was

Message("I am here");
if (FLUID_THREAD_P(t))
Message("Hey there!");
{

to be

{
Message("I am here");
if (FLUID_THREAD_P(t))
{
Message("Hey there!");

[Aj_45]

Thanks for the help. Although, I used the volume integral option in Monitors tab. That solved the issue.