[chenxizh]

Hi, I have a problem about the format of data. The problem is as follows: I have a bed as 0~Xmax, 0~Ymax, the size of grid is dx, dy. as in the SIMPLE scheme, the pressure is stored in cell center, so in order to visual the data in paraview, I do the following output
1. first do the loop to output the cell center coordinate, i*dx , i*dy,
2. then I want to output the cell center data of pressure P[i][j],
So the PROBLEM is that should I specific the type of pressure data to be CELL_DATA or POINT_DATA in Paraview?
can any one help me?
Thanks a lot!

The code is as follows,

// the version
    fprintf(output, "# vtk DataFile Version 2.0\n");

    // the title
    fprintf(output, "Fluid Flied\n");

    // the file format
    fprintf(output, "ASCII\n");

    // the datast structure
    fprintf(output, "DATASET STRUCTURED_GRID\n");

    // the datast dimensions 
    fprintf(output, "DIMENSIONS %d %d %d\n", NX, NY, NZ);  
  // points info
    fprintf(output, "POINTS %d int\n", NX*NY*NZ);  // this is the position
    for (z=0; z<NZ; z++)
    {
        for (y=0; y<NY; y++)
        {
            for (x=0; x<NX; x++e)
            {
                fprintf(output, "%5d%5d%5d\n", x, y, z);
            }
        }
    }
    
    // the dataset attributes
    // point id 
    fprintf(output, "POINT_DATA %d\n", NX*NY*NZ);  // Or This should be CELL_DATA:confused:, the number should be (nx-1)*(ny-1)*(nz-1) and start from 1
    fprintf(output, "SCALARS presuss double %d\n", 1);
    fprintf(output, "LOOKUP_TABLE default\n");

    for (z=0; z<NZ; z++)
    {
        for (y=0; y<NY; y++)
        {
            for (x=0; x<NX; x++)
            {
                fprintf(output, "%lf\n", press[z][y][x]);
            }
        }
    }