Cannot export river package as shapefile properly

[sebmeng]

Hello community:
I have some troubles to export the river package of a model I created. If I use the function riv.export(os.path.join(model_ws, "riv.shp")) I get a shapefile of the whole model, but the only attributes that it has is the node number.

If I use flopy.plot.PlotMapView I get a nice view of the cells that has the river boundary condition:

#plot river BC using flopy
fig, ax = plt.subplots(1, 1, figsize=(10, 8), constrained_layout=True)
mv = flopy.plot.PlotMapView(model=gwf, ax=ax, layer=6) #in layer 7 the bottom is present
mv.plot_bc("RIV",plotAll=True)

and if I export the model to a VTK file and add the river package, I can also see them quite easily:

But I need them in shp format to make some checking with GIS. Maybe I am doing something wrong in the export process?

Thanks in advance for your help!
Cheers,
Sebastian

[martclanor]

Hi Sebastian,

Trying to look around here. Do you maybe happen to have multiple stress periods in your model while your river package's stress period data is defined only on a subset of your stress periods, (perhaps the initial one)?

If it's anything similar to that, then my guess is that it's because the MFPandasTransientList does not re-use the data from the preceding stress periods when it is empty (which is probably intended as mf6 would handle that logic already). In that case, an error is caught on the following code block, and thus exporting of the stress period data is skipped.

flopy/flopy/export/shapefile_utils.py

Lines 380 to 383 in 181e101

	try:
	list(a.masked_4D_arrays_itr())
	except:
	continue

I can try to submit a PR on this, but maybe the faster solution for you is to just assign the river stress period data for each stress period.

Regards,
Mart

[sebmeng]

Hi Mart,
That is indeed how I managed the river data, since MF6 copies the information of the previous defined stress-period I just have the information of the first SP defined in my RIV package. I will add the info for the rest of the SP and see how it goes,
Thanks!

Sebastian