[blowback]

When I compile below function for steady-state simulation, everything is fine. In case of transient simulation I got an error message when journal file set up wall BC with heat flux on it.

After this command:
/def/boundary/wall selection 0 n 0 n n y y "udf" "p1::libudf" n n n n 0 n 0
I got an error (only if transient is set up) - errorlog in enclosed file

http://sendfile.pl/194535/fluenterror.log

#include "udf.h"
#include <math.h>
#define Vp 0.815
#define Vk 0.20375
#define dT 3.0
int Npts = 99 ; 
float *mean_power;

float get_q( int i)
{
  float ydata;
  ydata = mean_power[i];
  return ydata;
}
DEFINE_ON_DEMAND(read_mean_power)
{
int i=0;
  FILE* plik;
  plik = fopen("moce99.dat","r");
  mean_power = (float *) malloc(Npts*sizeof(float));
   for ( i=0;i<Npts;i++ )
    {
      fscanf(plik,"%f\n",&mean_power[i]);
    }
  fclose(plik);
}
DEFINE_PROFILE(velocity, ft, var)
{
    real time;
    face_t f;
    time = CURRENT_TIME;
    begin_f_loop(f,ft)
    {
        if (dT>time)
        {
            F_PROFILE(f,ft,var) = Vp - (Vp-Vk)*time/dT;
        }
        else
        {
            F_PROFILE(f,ft,var) = Vk;
        }
    }
    end_f_loop(f,ft)
}

DEFINE_PROFILE(p1, ft, var)
{
    float x[3];
    float z;
    face_t f;
    int i = 0;
    float q0;
q0 = get_q(i);
    
    begin_f_loop(f,ft)
    {
        F_CENTROID(x,f,ft);
        z = x[2];
        F_PROFILE(f,ft,var)= q0*sin(((z-0.231)/0.65)*M_PI);
    }
    end_f_loop(f,ft)
}
...
...
...
DEFINE_PROFILE(p99, ft, var)
{
    float x[3];
    float z;
    face_t f;
    int i = 98;
    float q0;
q0 = get_q(i);
    
    begin_f_loop(f,ft)
    {
        F_CENTROID(x,f,ft);
        z = x[2];
        F_PROFILE(f,ft,var)= q0*sin(((z-0.231)/0.65)*M_PI);
    }
    end_f_loop(f,ft)
}