[Shoh]

Hello. I would really appreciate if you could give me an answer.

I'd like to apply a boundary condition so that the UDS of the present surface is 10% of the sum of UDS of the other surface of which id is 7.

Here is my udf code.
-------------------------------------------
#include "udf.h"
#define idoutlet 7

DEFINE_PROFILE(udf_filter,t,i)
{
real sumudf;
face_t f;

sumudf=0;
begin_f_loop(f,idoutlet)
{
sumudf+=F_UDSI(f,idoutlet,0);
}
end_f_loop(f,idoutlet)

begin_f_loop(f,t)
{
F_PROFILE(f,t,i)=sumudf*0.1;
}
end_f_loop(f,t)

}
-------------------------------------------

But, when the calculation starts, the solver stops and the program is shut down.

Could you tell me where is the problem in my code?

[AlexanderZ]

#include "udf.h"
#define idoutlet 7
real sumudf;

DEFINE_PROFILE(udf_filter,t,i)
{
	face_t f;
	begin_f_loop(f,t)
	{
		F_PROFILE(f,t,i)=sumudf*0.1;
	}
	end_f_loop(f,t)
}

DEFINE_ADJUST(surf_avg,d)
{
	Thread *t;
	face_t f;
	t = Lookup_Thread(d,idoutlet);
	sumudf=0.0;
	thread_loop_f(t,d)
	{
		begin_f_loop(f,t)
		{
			sumudf+=F_UDSI(f,t,0);
		}
		end_f_loop(f,t)
	}
	# if RP_NODE
		sumudf = PRF_GRSUM1(sumudf);
	# endif
	Message0("**** sumudf = %f ****\n",sumudf);
}