hfb6.export function issue?

[Matt-Webb1]

Hi,

I wondered if anyone could tell me if the hfb6.export function is working correctly? I can use .export on a number of the packages without issue (e.g. DIS/BAS/RIV/UPW/GHB). It is listed as an available function of hfb6 but I can't get it to export anything. E.g.:

model_ws = os.path.join("model_data")
output_folder = os.path.join('output')
ml = flopy.modflow.Modflow.load(model.nam, model_ws=model_ws, check=False)

# This creates the RIV folder and exports the River BCs to a .vtu file
river_elev_output = os.path.join(output_folder, 'RIV')
ml.riv.export(river_elev_output, fmt='vtk')

# This creates the folder HFB but does not export anything
hfb_output = os.path.join(output_folder, 'HFB')
ml.hfb6.export(hfb_output, fmt='vtk')

Thanks!

Matt

[christianlangevin]

Sorry, but it does not look like HFB information can be written to vtk yet. Ideally, the vtk representation should be represented with planes that correspond to the shared cell face where the hfb is located. Will keep this issue open as a reminder to implement.

[Matt-Webb1]

Thanks! I will try to see if I can use the HFB data to generate the required vertices etc to create the VTK file.

[jdhughes-dev]

The VTK export functionality is currently being revised. We may be able to implement this as part of that effort.

Also see #873

[jlarsen-usgs]

@jdhughes-usgs, @Matt-Webb1

I am working on a hfb export method for VTK as part of the Vtk class overhaul.