[jmt]

Hi everyone,

I'm trying to write a small code that creates a mesh grading by manually modifying the location of wall-normal points in the mesh.

The geometry is a backward-facing step with three blocks as shown in the image.



I originally generated the mesh with blockMesh and I can get the points as follows:

pointField newPoints = mesh.points();

The points list is a 1D array which happens to be listed as all block 0 points, all block 1 points, and then all block 2 points. My method is to loop over all points and based on my a priori knowledge of the number of points per block, I can modify the wall-normal coordinate location appropriately.

The code looks like this for just the 0th block. I hard-code the number of points in pointsBlock0 based on what I know, but I wonder if I can retrieve this from the API. Then, I track a global point index to modify the list.

    scalar pointIndex = 0;
    // Block 0  
    vector pointsBlock0(119, 32, 129);

    for (label kk = 0; kk < pointsBlock0.z(); kk++)
    {
        for (label jj = 0; jj < pointsBlock0.y(); jj++)
        {
            Info << "kk: " << kk << " jj: " << jj << " y: " << points[pointIndex][1] << nl;
            for (label ii = 0; ii < pointsBlock0.x(); ii++)
            {
                points[pointIndex][1] = ... // some new calculated y-location to satisfy grading
                pointIndex++;
            }
        }

My question is: Is there some ability in the points API to get the dimensions of each block so I don't have to hard-code these point values?

Or is there some more elegant way to only modify the wall-normal portions of each coordinate? All wall-normal points in a given streamwise-spanwise plane need the same adjustment, but this adjustment is different for each streamwise-spanwise plane.