Question about error using interpolate_with_sva_pilotpoints_2d() provided grid specification file from modflow-usg (unstructured grid)

[DataGeoranger]

Hello!

I am trying to use interpolate_with_sva_pilotpoints_2d() function with unstructured grid providing the grid specification file associated with modflow-usg. Please see below the screenshot of the code and error message:

I also tried to use the function provided by pypestutils: 'get_grid_info_from_gridspec. Then I realized that it only works for structured grid type.

Could you please provide any suggestions on reading in the grid spec file from modflow-usg for this function or constructing SpatialReference object in general? Any tips are much appreciated!

[hjia1005]

According to the above, you don't have to provide the gsf file for an unstructured grid, a pandas dataframe with x,y, and layer will be accepted. I normally use this way for unstructured grid. Hope it helps.

[DataGeoranger]

@hjia1005 Thanks much for the tips! It works when passing x, y, layer as df.

However, the interpolate_with_sva_pilotpoints_2d() broke again, when I tried to do interpolation just for a particular K zone.
For illustration, it works when I have pilot points covering the entire model area, and it stops working when I interpolate a sub zone.

My guess is the correlation length for max distance for kriging probably breaks when having sub-zone pp. Another issue is related to grid_info file, there is not column for zone number for each node, I don't see get_2d_grid_info_from_file() is trying to read in zone info as a input parameter. I am not sure where to provide that, which is sort like the node_property that plproc takes, so it allows me to perform kriging just within a defined area within the model. Maybe this should be a new topic for discussion. Thanks.

[hjia1005]

@DataGeoranger In that case, I would only use the grid information (x and y) in that specific zone other than the entire model. You can do the same for multiple zones, just like what you do when using plproc in PEST utilities.

[DataGeoranger]

@hjia1005. Gotcha. So it's up to the user to provide the selected cells compared to the selection in plproc. Thanks for the tips!