Heat Flux UDF not working
 

#include "udf.h"
#include "mem.h"


DEFINE_PROFILE(HeatFlux,t,i)
{

real temp, press, temp_grad, a, b, cz, Psat, denominator, mass_flux, heat_flux;
real vel,tc,tp,dt,den,dtdxw;

cell_t c,c0;
face_t f;
Thread *t0 = t->t0;

if (!Data_Valid_P()) return;

begin_f_loop (f,t)
{

c0 = F_C0(f,t);
t0 = F_C0_THREAD(f,t);



press = C_P(c0,t0);
temp = F_T(f,t);


a=0.1;
b=10.0896-2926.6/temp;
cz=2*3.14*65*temp;
Psat=133.3*pow(2.71,b);

denominator=pow(cz,0.5);

mass_flux=a*(Psat-press)/denominator;
heat_flux=mass_flux*(574000.0);


F_PROFILE(f,t,i) = 574000.0*mass_flux;

}
end_f_loop(f,t)

}


i am trying to use heat flux on the wall but im getting error ''compute process interrupted". (density based solver)

When i am using a lower value(574*mass_flux) then it is working. I have tried replacing 'real' with 'float' but no luck.

Thread *t0;