[PierceH]

Hello,


I want to specify a certain flux for my user defined scalar. My idea is to create a profile that I can then implement as the flux condition on the given wall. In this example I am trying to set a boundary that would recover pressure. and So the pressure BC is nabla(P) dot n. My code as such would ideally read as follows, where F_UDSI(f,t,12) is dp/dx and 13 is dp/dy.


DEFINE_PROFILE(pressure_profile,t,i)
{
real x[ND_ND];
real y;
real A[ND_ND];
face_t f;
begin_f_loop(f,t)
{
F_AREA(A, f, t);
F_PROFILE(f,t,i) = A[0]*F_UDSI(f,t,12)+A[1]*F_UDSI(f,t,13);
end_f_loop(f,t)
}
}

It appears to be working in the sense that it will hook but is not doing as desired, I'm unsure If this is incorrect from a computational standpoint or I am missing something here.

Thanks for your help.