outlet pressure
 

Hi all
i'm going to analysis 2d centrifugal pump.in outlet line,i try to define static pressure in a single face of grid. i use this udf. but it doesn't work correctlly.

#include"udf.h"

DEFINE_PROFILE(pressure_outlet, t, i) * this macro is used to define values in BCs*/
{
real pr;
face_t f;
int j;
real p;
int k;
real y;
real x[ND_ND];
if(!Data_Valid_P()) /*check for the availability of data*/
return;
begin_f_loop(f,t)
{
for(k=0;k<12;k++)
p=F_P(f,t); /*F_P read the pressure of each face and assign it to p*/
F_UDMI(f,t,k)= p; /*this macro is used to store/access the values of the user-defined memory on the face.the outlet boundary is devided to 12 subline.so i set the number of memory 12.define>user deine>memory*/

printf("wwwwwwww:%f\n",F_UDMI(f,t,k));
}
end_f_loop(f,t)

begin_f_loop(f, t)
{
for(j=0;j<12;j++)
F_CENTROID(x,f,t);
y=x[0]; /* the outlet line is horizontal and i need the x coordinate*/
if (y<-0.027) /* there is a face centroid in this interior*/
{
F_PROFILE(f,t,i)=100000; /*desired value of the defined face cell*/
}
else
{
pr=F_UDMI(f,t,j); /* access to stored values for each face*/
F_PROFILE(f, t, i) = pr;
}
printf("rrrrrrrrr:%f\n",F_PROFILE(f,t,i));
}
end_f_loop(f, t)
}