accessing gradient in boundary for define_proile

Asghari

Hi all;

I want accessing to pressure reconstruction gradient or usual gradient in bondary using my define_profile. My define_profile for pressure outlet is in following: DEFINE_PROFILE(pressure_outlet, thread, index)

{

real A[ND_ND],dr0[ND_ND],es[ND_ND],d1,dr,A_by_es ;

face_t f;

begin_f_loop(f, thread) /* loops over all faces in the thread passed in the DEFINE macro argument */

{ /** identify the cell thread adjacent to the face thread **/ Thread *t0=thread->t0;

cell_t C0=F_C0(f,thread);

BOUNDARY_FACE_GEOMETRY(f,thread,A,dr,es,A_by_es,dr 0) ; if (NNULLP(T_STORAGE_R_NV(t0,SV_P_RG)) )

{

d1=dr0[0]*C_P_RG(C0,t0)[0]+dr0[1]*C_P_RG(C0,t0)[1] ;

F_PROFILE(f, thread, index) = C_P(C0,t0)+d1;

}

else

{

F_PROFILE(f, thread, index) = -500.0;

}

}

end_f_loop(f, thread)

} But pressure in outlet is -500 at all and is'nt entered to first loop in if command. How can i solve this problem ?

Thanks alot for your help.