level 2
#include "udf.h"
#define CELL_ID 64
DEFINE_SOURCE(energy_sourced, c, t, dS, eqn)
{
real source;
real x = CURRENT_TIME;
Domain *d;
real temperature_sum, count, y;
d = Get_Domain(1);
t = Lookup_Thread(d, CELL_ID);
thread_loop_c(t, d)
{
begin_c_loop(c, t)
{
real temperature = C_T(c, t);
temperature_sum += temperature;
count++;
}
end_c_loop(c, t)
}
temperature_sum /= count;
y = temperature_sum;
if (x >= 0 && x < 1220)
{
source = 1000*y;
}
else if (x >= 1220 && x <= 1540)
{
source = 100000;
}
return source;
}
2023年11月08日 03点11分
