Post-processing for MF6/PESTPP-IES realisations: Would like Cell-based fluxes from binary cbc for Quadtree grid

[eb-oz]

Is there a FloPy command to extract fluxes from a binary cell-by-cell (*.cbc) file produced by MF6 for a DISV (quadtree) grid, and combine it with *.disv information about the cell vertices to produce an array of cellIDs and fluxes? I'm not finding any of the current functions applicable to this situation. When I use the get_data() command, it shows me that the *.cbc contains an array of 2 nodes per flux, which I assume are the two nodes along shared edges between cells, but I would like the fluxes per cellID. I'm trying to make a workflow to do this for 800+realisations. I know an alternative is to use ZBUD6 which produces large .csv files with exactly what I need, but thought I'd ask here, in case there is a way to more directly access the binary cbc with FloPy and maybe combine it with the *.disv file to output, in my case: CHD cellID and its flux. Thank you.

[martclanor]

I think you can derive the "fluxes per cellID" from the "array of 2 nodes per flux" that you already have extracted from the cbc file.

Alternatively, maybe what you need instead is the Observation package with the "chd" observation type. Here's a link to the Observation package tutorial.