[om1234]

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, but when i try to use it with if else ,answers are correct. i myself think that the problem is that the define-profile macro is a array of data and i must define the pressure in all of the cell faces, but i'm not sure.I wish any one can help me.

DEFINE_PROFILE(outlet_p, thread, position)
{
real x;
real y[ND_ND];
face_t f;

begin_f_loop(f, thread)
{
F_CENTROID(y,f,thread);
if (0.004<x<0.006) /* in this interior there is a face centroid or if (0.004<x & x<0.006)*/
F_PROFILE(f, thread, position) = 1e5;
}
end_f_loop(f, thread)
}