get_specific_discharge returns discharge at the centers while position='faces'

[OnnoEbbens]

When I call the function flopy.utils.postprocessing.get_specific_discharge() with position='faces' I get the discharge of the cell centers instead of the cell faces. Calling the function with position='centers' I get exactly the same results.

I tried to use this code to get the specific discharge at the cell faces:

cbcobj = flopy.utils.CellBudgetFile(os.path.join(gwf.model_ws, f"{gwf.name}.cbc"))
spdis = cbcobj.get_data(text="DATA-SPDIS")[0]
qx, qy, qz = flopy.utils.postprocessing.get_specific_discharge(spdis, gwf, position="faces")

Looking into the get_specific_discharge function I can see that the position variable is not used when vertices is a recarray. In my case spdis is a recarray.

It would be nice if the function can return the specific discharge on the cell faces when vertices is a recarray. If that is hard or impossible I think an error should be raised when calling this function with a recarray and position='faces'.

My modelgrid is a locally refined grid with type 'vertex', for all that matters.