[hussam_alziadi86]

I have two matrices ones is the u-component velocity and the other is v-component velocity,I am trying to read the values of these two matrices and put them as an initial conditions in the domain of my model. I use the code below but I do not know how to read the values from the matrices and then write them into the model.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
#include "udf.h"
#define pi 4.*atan(1.)

DEFINE_INIT(uv_intial,d)
{

cell_t c;
Thread *t;
real xc[ND_ND];
FILE *fp;

real r1;
real xo1 = 0.00;
real yo1 = 0.00;

/* loop over all cell threads in the domain */
thread_loop_c(t,d)
{
/* loop over all cells */
begin_c_loop_all(c,t)
{
C_CENTROID(xc,c,t);

r1=sqrt(ND_SUM(pow(xc[0] - xo1,2.),
pow(xc[1] - yo1,2.),
pow(xc[2] - 0,2.)));
/* Note xc[0] is x coordinate, xc[1] is y and xc[2] is z */

if ( r1 != 0.0)
{
fp=fopen("u.txt","r");

fscanf(fp,"lf% %lf",&r1,&u);

DEFINE_RW_FILE(uv_intial,fp)


C_U(c,t) =
C_V(c,t) =
}
else
{
C_U(c,t) = 0.0;
C_V(c,t) = 0.0;
}

}

end_c_loop_all(c,t)
fclose(fp);
}
}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


thanks in advance for your help