[sinaited]

Hi everyone, hope you be healthy!
i wrote this udf for calculating the average pressure on a wall, but i think it only can be interpret in serial and i should rewrite this code for parallel or linux node, anyone can help me or rerwite it?




#include "udf.h"
real pavg;
DEFINE_EXECUTE_AT_END(average_pressure)
{
Domain *d;
face_t f;
real pressure = 0.0;
real A[ND_ND];
real area = 0.0;
real area_tot = 0.0;
int ID = 7; /*this is the ID of the boundary wall that I want to get the temperature from*/
Thread *t;
int zone_ID;
d = Get_Domain(1);
t = Lookup_Thread(d,ID);
pavg = 0.0;
begin_f_loop(f,t)
{
F_AREA(A,f,t);
area = NV_MAG(A)*2.0*M_PI; /*AREA*/
area_tot += area;
pressure = F_P(f,t); pavg += pressure*area;
Message ("Text %g \n", pavg);
}
end_f_loop(f,t)
pavg /= area_tot;
printf("Pavg = %g area_tot = %g\n",pavg,area_tot);
Message ("Text %g \n", pavg);
}

DEFINE_CG_MOTION(move_body, dt, vel,omega, time, dtime)
{
real ome,Ra,A,B;
/* reset velocities */
NV_S(vel,=,0.0);
NV_S(omega,=,0.0);
/* motion */
ome=pavg / 1000;
vel[0] = 0.0;
vel[1] = ome;
vel[2] = 0.0;
omega[0] = 0.0;
omega[1] = 0.0;
omega[2] = 0.0;
}

[AlexanderZ]

compile the code

#include "udf.h"
real pavg;
DEFINE_EXECUTE_AT_END(average_pressure)
{
Domain *d;
face_t f;
real pressure = 0.0;
real A[ND_ND];
real area = 0.0;
real area_tot = 0.0;
int ID = 7; /*this is the ID of the boundary wall that I want to get the temperature from*/
Thread *t;
int zone_ID;
d = Get_Domain(1);
t = Lookup_Thread(d,ID);
pavg = 0.0;
begin_f_loop(f,t)
{
F_AREA(A,f,t);
area = NV_MAG(A)*2.0*M_PI; /*AREA*/
area_tot += area;
pressure = F_P(f,t); pavg += pressure*area;
Message0("Text %g \n", pavg);
}
end_f_loop(f,t)
#if RP_NODE
	area_tot = PRF_GRSUM1(area_tot);
#endif
pavg /= area_tot;
Message0("Pavg = %g area_tot = %g\n",pavg,area_tot);
Message0("Text %g \n", pavg);
}

DEFINE_CG_MOTION(move_body, dt, vel,omega, time, dtime)
{
real ome,Ra,A,B;
/* reset velocities */
NV_S(vel,=,0.0);
NV_S(omega,=,0.0);
/* motion */
ome=pavg / 1000;
vel[0] = 0.0;
vel[1] = ome;
vel[2] = 0.0;
omega[0] = 0.0;
omega[1] = 0.0;
omega[2] = 0.0;
}

[sinaited]

Thanks a lot alexander, it worked for area and based on that, i corrected it for pressure