hwva for mf6-disu

[hjia1005]

I try to use gridgen to create mf6 with disu grid.

disu6_gridprops = g.get_gridprops_disu6()
disu=flopy.mf6.ModflowGwfdisu(gwf, **disu6_gridprops)

But the hwva data is not correct and the model doesn't run because of that.

This is the first few lines of the hwva in the disu file.

It looks correct after some lines

Any solutions will be appreciated.

[hjia1005]

Found the issue. The default top and bottom for the top layer of the gridgen object is the same. There was a divide by 0 problem when the hwva array was calculated. Just need to reassign the top and bottom of the dis for the gridgen object.