[Bibill]

Dear all,

I ran a first simulation on a mesh and i have exported a quantity of interest (QoI) at each cell center of my CFD domain in a ASCII file type. It is thus a very large file (millions of lines):

cellNumber, x-coordinate, y-coordinate, z-coordinate, QoI
.
.
.


I would like to re-use this result for a second simulation on the same mesh. How can I read the file and store the QoI inside a C_UDMI(c,t) ?

I guess the following option is naive

DEFINE_ON_DEMAND(read_QoI)
{

#if !RP_NODE /* serial or host process */

char line[100];
int index;
float f1, f2, f3, f4;
int number_cell = 13786084;

float *QoI;
QoI = (float *)malloc((number_cell -1)*sizeof(float));

FILE* file_r;
	file_r = fopen("QoI_ascii", "r");
	if(file_r==NULL)
	{
		Message("\n");
		Message("Error opening file \n");
		Message(" \n");
	}


    /* Skip one line by reading one line and ignoring */
    fgets(line, 100, file_r);
	
	
    while (fscanf(file_r, "%d%f%f%f%f", &index, &f1, &f2, &f3, &f4) > 0)
    {
		QoI[index-1] = f4;
    }

fclose(file_r);


Domain *d;
d = Get_Domain(1);
Thread *t;
cell_t c;

thread_loop_c(t,d)
{
    begin_c_loop(c,t)
    {
	int ID = C_ID(c,t);
        C_UDMI(c,t,0) = QoI[ID-1];
    }
    end_c_loop(c,t)
}

free(signal);

#endif


}

I would be highly interested if anyone could help!

Thanks in advance !

[AlexanderZ]

1. you may write profile and later read it without any UDF

2. on the other hand for single core your code should work well in case you are storing the same cell index as C_ID(c,t) macro returns.

[Bibill]

Originally Posted by AlexanderZ

1. you may write profile and later read it without any UDF

2. on the other hand for single core your code should work well in case you are storing the same cell index as C_ID(c,t) macro returns.

Thanks AlexanderZ!

1. Can the profile tool also be used with field data? Or should it be used for surface ?

2. Ok thanks a lot, I'll give a try with a smaller file then!