Get xyz from vertices in LGRs and 8-corners-of-refined-cells #792
Conversation
| // | ||
| std::vector<std::vector<std::array<double,3>>> xyz_for_zcorn(grid.currentData()[level]->size(0)); // size total amount of cells in the level refined grid. | ||
| for (const auto& element : elements(grid.levelGridView(level))) { | ||
| xyz_for_zcorn[element.index()].reserve(8); |
There was a problem hiding this comment.
Maybe also get the local global index. Similar to above. That together with the order you documented above will be all the information needed.
There was a problem hiding this comment.
I will add that here then. Thanks for your feedback!
| // Expected values in each LGR: 8*refined_cells_in_that_lgr | ||
| // Each refined cell has its 8 corners stored in <refined_level_grid>.cell_to_point_[ refined element index ] | ||
| // in the following order: | ||
| // BOTTOM FACE TOP FACE | ||
| // 2 --- 3 6 --- 7 | ||
| // / / / / | ||
| // 0 --- 1 4 --- 5 |
There was a problem hiding this comment.
It typically does not matter, but for what it's worth, the default coordinate system is right-handed with the origin in the upper, left, back corner and Z axis pointing down. In other words, the typical node ordering is as outlined in the picture in MRST's processGRDECL function.
There was a problem hiding this comment.
I see, in the grid, we follow the described ordering. Therefore, for output, we need to re-order them according to your description. Thanks for your feedback!
There was a problem hiding this comment.
Therefore, for output, we need to re-order them according to your description.
Sorry for being too brief: It's not that we need to do anything in particular to the order–it's actually the same in both views–it's just that in terms of labelling the top face is normally comprised of vertices 0..3 while the bottom face normally comprises vertices 4..7.
There was a problem hiding this comment.
Thanks for clarifying! I added your comment into the documentation.
|
Printed values: coord_zcorn.log |
|
jenkins build this please |
This PR aims to contribute to the generation of output files for simulations involving Corner Point Grids (CpGrid) with Local Grid Refinement (LGR). It demonstrates how to extract the necessary values for defining COORD and ZCORN in CpGrids with LGRs. Note that the ordering presented here does not follow the "pillars" construction. To correctly output COORD and ZCORN, a reordering step must be implemented.
Not relevant for the Reference Manual.
This commit translates one of the test cases from OPM/opm-common#4255 using the grid interface.