[p.norouzi.gv]

I am trying to solve a UDS in a particular part of a pipe (attached file). I prepared a UDF for the source term. Now I am trying to define a face-loop for boundary condition in this section. I was wondering how I can make a loop for face boundary conditions in internal faces. I wrote the UDF, but it belongs to the external boundary. I have to define C0 and C1 as far as I know, but I could not define them.

I appreciate any help you can provide.


DEFINE_PROFILE(Boundray_FI,t,nv)
{

face_t f;
Thread *t0;
cell_t c, c0;
real x[ND_ND];
real time;
time=CURRENT_TIME;

begin_f_loop(f,t)
{
F_CENTROID(x,f,t);
c0=F_C0(f,t);
t0=THREAD_T0(t);
F_PROFILE(f, t, nv) = C_UDSI(c0, t0, 0)+30*time;
}
end_f_loop(f,t)
}