[sakura006]

I have problem with parallel mode. I defined a source in UDF. If I run in serial mode, everything is fine. However, if I run in parallel mode with the same code, the region of the source is incomplete. That means part of the source is missing in the region where it should be. I am not sure if it was my parallelization problem.
Here is the code.

DEFINE_SOURCE(Wsource,c,tc,dsource,Vequation)
{
#if !RP_HOST
real x[ND_ND];
real x_c;
real y_c;
real z_c;
real wsource=0.0;
real volume2;
C_CENTROID(x,c,tc);

x_c = x[0];
y_c = x[1];
z_c = x[2];

volume2 = 0.029890004736;
dsource[Vequation] = 0.0;

if ((x_c >= 4.572 && x_c <= 5.1816) && (y_c >= 0.0254 && y_c <= 0.508) && (z_c >= -1.016 && z_c <= -0.9144))
{
wsource = -1.0 * (rho_air * server_area * mv_hots1 * ((mv_hots1 / porosity) - mv_hots1)) / volume2;
}
else if ((x_c >= 4.572 && x_c <= 5.1816) && (y_c >= 0.508 && y_c <= 0.9906) && (z_c >= -1.016 && z_c <= -0.9144))
{
wsource = -1.0 * (rho_air * server_area * mv_hots2 * ((mv_hots2 / porosity) - mv_hots2)) / volume2;
}
else if ((x_c >= 4.572 && x_c <= 5.1816) && (y_c >= 0.9906 && y_c <= 1.4732) && (z_c >= -1.016 && z_c <= -0.9144))
{
wsource = -1.0 * (rho_air * server_area * mv_hots3* ((mv_hots3 / porosity) - mv_hots3)) / volume2;
}
else if ((x_c >= 4.572 && x_c <= 5.1816) && (y_c >= 1.4732 && y_c <= 1.9558) && (z_c >= -1.016 && z_c <= -0.9144))
{
wsource = -1.0 * (rho_air * server_area * mv_hots4 * ((mv_hots4 / porosity) - mv_hots4)) / volume2;
}
else
{
wsource = 0.0;
}

C_UDMI(c,tc,1) = wsource;
return wsource;
#endif

}

Here is the source in serial mode
serial source.PNG
And source in parallel mode
parallel source.PNG

[sakura006]

Can someone help me?