#include "udf.h"
DEFINE_PROFILE(topcell_temp, thread, index)
{
face_t f;
int curr_ts;
curr_ts= RP_Get_Integer("time-step");
begin_f_loop(f, thread)
{
if (curr_ts>= 0 && curr_ts<= 4 )
{
F_PROFILE(f, thread, index) = 305.7;
}
else if (curr_ts>= 5 && curr_ts<= 8 )
{
F_PROFILE(f, thread, index) = 301;
}
else if (curr_ts>= 9 && curr_ts<=12 )
{
F_PROFILE(f, thread, index) = 306;
}
else if (curr_ts>=13 && curr_ts<=16)
{
F_PROFILE(f, thread, index) = 305;
}
else if (curr_ts>=17 && curr_ts<=22)
{
F_PROFILE(f, thread, index) = 311;
}
else if (curr_ts>=23 && curr_ts<=26)
{
F_PROFILE(f, thread, index) =306;
}
else if (curr_ts>=27 && curr_ts<=32)
{
F_PROFILE(f, thread, index) = 320;
}
}
end_f_loop(f, thread)
}
DEFINE_PROFILE(sidewin_temp, thread, index)
{
face_t f;
int curr_ts;
curr_ts= RP_Get_Integer("time-step");
begin_f_loop(f, thread)
{
if (curr_ts>= 0 && curr_ts<= 4 )
{
F_PROFILE(f, thread, index) = 320.18;
}
else if (curr_ts>= 5 && curr_ts<= 8 )
{
F_PROFILE(f, thread, index) = 306;
}
else if (curr_ts>= 9 && curr_ts<=12 )
{
F_PROFILE(f, thread, index) = 320.75;
}
else if (curr_ts>=13 && curr_ts<=16)
{
F_PROFILE(f, thread, index) = 310.4;
}
else if (curr_ts>=17 && curr_ts<=22)
{
F_PROFILE(f, thread, index) = 337;
}
else if (curr_ts>=23 && curr_ts<=26)
{
F_PROFILE(f, thread, index) =302;
}
else if (curr_ts>=27 && curr_ts<=32)
{
F_PROFILE(f, thread, index) = 344;
}
}
end_f_loop(f, thread)
}
DEFINE_PROFILE(backwin_temp, thread, index)
{
face_t f;
int curr_ts;
curr_ts= RP_Get_Integer("time-step");
begin_f_loop(f, thread)
{
if (curr_ts>= 0 && curr_ts<= 4 )
{
F_PROFILE(f, thread, index) =323.7;
}
else if (curr_ts>= 5 && curr_ts<= 8 )
{
F_PROFILE(f, thread, index) = 309;
}
else if (curr_ts>= 9 && curr_ts<=12 )
{
F_PROFILE(f, thread, index) = 324.28;
}
else if (curr_ts>=13 && curr_ts<=16)
{
F_PROFILE(f, thread, index) = 312.57;
}
else if (curr_ts>=17 && curr_ts<=22)
{
F_PROFILE(f, thread, index) = 343;
}
else if (curr_ts>=23 && curr_ts<=26)
{
F_PROFILE(f, thread, index) =303;
}
else if (curr_ts>=27 && curr_ts<=32)
{
F_PROFILE(f, thread, index) = 349;
}
}
end_f_loop(f, thread)
}
DEFINE_PROFILE(frontwin_temp, thread, index)
{
face_t f;
int curr_ts;
curr_ts= RP_Get_Integer("time-step");
begin_f_loop(f, thread)
{
if (curr_ts>= 0 && curr_ts<= 4 )
{
F_PROFILE(f, thread, index) =366.38;
}
else if (curr_ts>= 5 && curr_ts<= 8 )
{
F_PROFILE(f, thread, index) = 355;
}
else if (curr_ts>= 9 && curr_ts<=12 )
{
F_PROFILE(f, thread, index) = 366.76;
}
else if (curr_ts>=13 && curr_ts<=16)
{
F_PROFILE(f, thread, index) = 357.5;
}
else if (curr_ts>=17 && curr_ts<=22)
{
F_PROFILE(f, thread, index) = 381;
}
else if (curr_ts>=23 && curr_ts<=26)
{
F_PROFILE(f, thread, index) =349;
}
else if (curr_ts>=27 && curr_ts<=32)
{
F_PROFILE(f, thread, index) = 386;
}
}
end_f_loop(f, thread)
}
DEFINE_PROFILE(midleftlet_temp, thread, index)
{
face_t f;
int t;
begin_f_loop(f, thread)
{
t= RP_Get_Integer("time-step");
if (t >= 0 && t<=8)
{
F_PROFILE(f, thread, index) =273+21.5-0.7 * (t*5 - 0) + 0.0216 * (t*5-0) *(t*5- 9)-0.000486* (t*5 - 0) *(t*5- 9)*(t*5- 19);
}
if (t >= 9 && t<=16)
{
F_PROFILE(f, thread, index) =273+ 13.8+ 0.18*((t-8)*5 - 0)-0.012 * ((t-8)*5 - 0) *((t-8)*5- 10)
+0.000364 * ((t-8)*5 - 0) *((t-8)*5- 10)*((t-8)*5- 22);
}
if (t >= 17 && t<=22)
{
F_PROFILE(f, thread, index) = 273+ 15.7+0.325 * ((t-16)*5-0)-0.000357* ((t-16)*5 - 0) *((t-16)*5- 4)+ 0.000355* ((t-16)*5 - 0) *((t-16)*5- 4)*((t-16)*5- 14);
}
if (t >= 23 && t<=26)
{
F_PROFILE(f, thread, index) =273+ 23.6-1.47 * ((t-22)*5 - 0) +0.106*((t-22)*5- 0) *((t-22)*5- 3)-0.00449 * ((t-22)*5 - 0) *((t-22)*5- 3)*((t-22)*5- 9);
}
if (t >= 27 && t<=32)
{
F_PROFILE(f, thread, index) =273+ 20.8+ 1.075 * ((t-26)*5 - 0) -0.068 * ((t-26)*5 - 0) *((t-26)*5- 4)+ 0.00203 * ((t-26)*5 - 0) *((t-26)*5- 4)*((t-26)*5- 12);
}
}
end_f_loop(f, thread)
}
DEFINE_PROFILE(frmidleftlet_temp, thread, index)
{
face_t f;
int t;
begin_f_loop(f, thread)
{
t= RP_Get_Integer("time-step");
if (t >= 0 && t<=8)
{
F_PROFILE(f, thread, index) =273+21.3-0.73 * (t*5-0)+0.022* (t*5-0)*(t*5-9)-0.00045 * (t*5-0) *(t*5-9)*(t*5-19);
}
if (t >= 9 && t<=16)
{
F_PROFILE(f, thread, index) = 273+16.9+0.16 * ((t-8)*5 - 0)-0.00917*((t-8)*5-0)*((t-8)*5-10)+ 0.000245 * ((t-8)*5 - 0) *((t-8)*5- 10)*((t-8)*5- 22);
}
if (t >= 17 && t<=22)
{
F_PROFILE(f, thread, index) = 273+15.5+0.225 * ((t-16)*5 - 0)+0.00321*((t-16)*5-0) *((t-16)*5- 4)-0.0000489 * ((t-16)*5 - 0) *((t-16)*5- 4)*((t-16)*5- 14);
}
if (t >= 23 && t<=26)
{
F_PROFILE(f, thread, index) =273+ 23-0.57 * ((t-22)*5 - 0) + 0.0204 * ((t-22)*5 - 0) *((t-22)*5- 3)-0.0000823 * ((t-22)*5 - 0) *((t-22)*5- 3)*((t-22)*5- 9);
}
if (t >= 27 && t<=32)
{
F_PROFILE(f, thread, index) =273+ 21.6+ 0.725 * ((t-26)*5 - 0)-0.0395 * ((t-26)*5 - 0)
*((t-26)*5- 4)+ 0.00121 * ((t-26)*5 - 0) *((t-26)*5- 4)*((t-26)*5- 12);
}
}
end_f_loop(f, thread)
}
DEFINE_PROFILE(frrightlet_temp, thread, index)
{
face_t f;
int t;
begin_f_loop(f, thread)
{
t= RP_Get_Integer("time-step");
if (t >= 0 && t<=8)
{
F_PROFILE(f, thread, index) =273+25.7-0.87 * (t*5 - 0)+0.0277 * (t*5 - 0) *(t*5- 9)-0.000614 * (t*5 - 0) *(t*5- 9)*(t*5- 19);
}
if (t >= 9 && t<=16)
{
F_PROFILE(f, thread, index) =273+17.7+0.24 * ((t-8)*5 - 0) -0.0124*((t-8)*5-0) *((t-8)*5-10)+ 0.000329 * ((t-8)*5 - 0) *((t-8)*5- 10)*((t-8)*5- 22);
}
if (t >= 17 && t<=22)
{
F_PROFILE(f, thread, index) = 273+ 17.6+0.15 * ((t-16)*5-0)+0.00785 * ((t-16)*5-0) *((t-16)*5-4)-0.00021 * ((t-16)*5 - 0) *((t-16)*5- 4)*((t-16)*5- 14);
}
if (t >= 23 && t<=26)
{
F_PROFILE(f, thread, index) =273+ 24.9-0.6 * ((t-22)*5-0) + 0.0204* ((t-22)*5 - 0) *((t-22)*5- 3);
}
if (t >= 27 && t<=32)
{
F_PROFILE(f, thread, index) =273+ 22.8+0.8*((t-26)*5-0)-0.0458 * ((t-26)*5-0) *((t-26)*5-4)+ 0.00151* ((t-26)*5- 0) *((t-26)*5- 4)*((t-26)*5- 12);
}
}
end_f_loop(f, thread)
}
DEFINE_PROFILE(frleftlet_temp, thread, index)
{
face_t f;
int t;
begin_f_loop(f, thread)
{
t= RP_Get_Integer("time-step");
if (t >= 0 && t<=8)
{
F_PROFILE(f, thread, index) =273+23.2-0.7 * (t*5 - 0) + 0.021 * (t*5 - 0) *(t*5- 9)-0.000472* (t*5 - 0) *(t*5- 9)*(t*5- 19);
}
if (t >= 9 && t<=16)
{
F_PROFILE(f, thread, index) =273+18+ 0.2 * ((t-8)*5 - 0)-0.0106*((t-8)*5-0) *((t-8)*5-10)
+ 0.000269 * ((t-8)*5 - 0) *((t-8)*5- 10)*((t-8)*5- 22);
}
if (t >= 17 && t<=22)
{
F_PROFILE(f, thread, index) = 273+ 16.6+0.175*((t-16)*5-0)+0.00678* ((t-16)*5 - 0) *((t-16)*5- 4)
-0.000165 * ((t-16)*5 - 0) *((t-16)*5- 4)*((t-16)*5- 14);
}
if (t >= 23 && t<=26)
{
F_PROFILE(f, thread, index) =273+ 24.9-0.63 * ((t-22)*5 - 0) + 0.0241 * ((t-22)*5 - 0) *((t-22)*5- 3)-0.000205 * ((t-22)*5- 0) *((t-22)*5- 3)*((t-22)*5- 9);
}
if (t >= 27 && t<=32)
{
F_PROFILE(f, thread, index) =273+22.6+0.775*((t-26)*5-0)-0.040625*((t-26)*5-0) *((t-26)*5- 4)+ 0.00117 * ((t-26)*5- 0) *((t-26)*5- 4)*((t