[luzikato]

Greetings all,

I'm trying to compute average volume temperature of a region. This is my code.

/* Calculate Volume-weighted Average of Module Temperature. */
real volume;
real T_AVG_TEMP;

#if !RP_HOST
	thread = Lookup_Thread(domain, TEM_ID);
	begin_c_loop(c, thread)
	{
		volume += C_VOLUME(c, thread);
		T_AVG_TEM += C_T(c, thread)*C_VOLUME(c, thread);
	}
	end_c_loop(c, thread)

#if RP_NODE
	volume = PRF_GRSUM1(volume);
	T_AVG_TEM = PRF_GRSUM1(T_AVG_TEM);
#endif /* RP_NODE */

#endif /* RP_HOST */

I insert this into a DEFINE_ADJUST macro so the volume this will run at the start of every iterations. I am also using DEFINE_AT_END macro to print total volume calculated every iteration in console to check the value.
After running the simulation, I compared the total volume I calculated with UDF and using Volume Integral with Reports feature.

I found that by using UDF, I got the amount of total volume.

Total volume: 6.18199e-05

But then, using Reports I got.

                    Total Volume                 (m3)
-------------------------------- --------------------
                       tem_solid        5.9999975e-05
                ---------------- --------------------
                             Net        5.9999975e-05

I know I could just ignore it. However, in this case, unfortunately, I really can't tolerate any value error.
Is there something I need to fix in my code that could results difference in total volume? Any insight, would be highly appreciated.

Thanks in advance.
Regards,

A.

[AlexanderZ]

define volume = 0 explicitly before summation as well as T_AVG_TEMP = 0