[lav1]

Hello! I am building a DPM model with a UDF. In that UDF, at some point, I need to use the property of a cell adjacent to the wall. How can I do that? I was thinking of using macros THREAD_T(t) and F_C0(f, t), but I am not sure on how to define the face. Setting it as “face_t f;” seems not enough.

Does anyone have a hint on how to approach it? Any reply will be much appreciated!

[Siba11]

Originally Posted by lav1

Hello! I am building a DPM model with a UDF. In that UDF, at some point, I need to use the property of a cell adjacent to the wall. How can I do that? I was thinking of using macros THREAD_T(t) and F_C0(f, t), but I am not sure on how to define the face. Setting it as “face_t f;” seems not enough.

Does anyone have a hint on how to approach it? Any reply will be much appreciated!

Hi. You can use something like this.

c_face_loop(cell,t,n)
	{
		f = C_FACE(cell,t,n);
		tf = C_FACE_THREAD(cell,t,n);
	if(!BOUNDARY_FACE_THREAD_P(tf))
		{
		c0 = F_C0(f,tf);
		c1 = F_C1(f,tf);
                }
        }