Convert from structured grid to disv - how to get cell IDs?

[jkennedy-usgs]

I'm trying to convert an existing model from a structured grid to a DISV grid. I can import the model fine, and create the DISV grid following the example here but all of the packages (MAW, for example) retain the original (layer, row, column) coordinates. What (I think) I need is a routine that converts (layer, row, column) to cellid. Surely something like this has been done before? Thanks

[jlarsen-usgs]

@jkennedy-usgs

If the node ordering of the new grid is consistent with the node ordering of the original grid, then you could use the original modelgrid object to create a map of layer, row, col to node number. The get_node() method on StructuredGrid objects converts layer, row, column to node number. Alternatively, the get_lrc() method will convert node numbers into layer, row, column.

If the node ordering of the new grid is not consistent with the original grid, you could use intersect the cell centers of the original grid with the new VertexGrid object to create a mapping from the original layer, row, column to the node numbers. The method call for that approach is: modelgrid.intersect(x, y, [z])