运行河包MODFLOW时出现FloPy/Python警告？

Question

我正在尝试向我的模型添加河流，但是在我运行它之后，这条警告恰好出现在那里：

WARNING: Unable to resolve dimension of ('gwf6', 'riv', 'period', 'stress_period_data', 'cellid') based on shape "ncelldim"

我所做的是从栅格中获取数据，以便获得单元格位置(行和列)以及阶段和河床的值。到目前为止，这是我在错误较少的情况下做到这一点的唯一方法(至少我还没有得到异常或类似的东西)。下面是我使用的代码。

df = pd.DataFrame({'Layer': data[0], 'Row': data[1], 'Column': data[2], 'RiverStage': data[3], 'RiverBed': data[4]})

cellid = (dataset['Layer'], dataset['Row'],dataset['Column'])

stress_period_data = [((cellid), dataset['RiverStage'], 10., dataset['RiverBed'])]

river = flopy.mf6.ModflowGwfriv(gwf, stress_period_data = stress_period_data)

另外，我还有两个额外的问题；如果我的模型是稳定状态，我需要添加CHD包吗？另外，我声明我的周期(同样，因为它是稳定的)是1.0，1，1.0。我提到的任何事情都与之前的警告有关吗？

任何帮助都是非常感谢的。

Answer 1

您正在为cellid变量创建一维数组的元组，而不是元组的列表。

您希望每个河流单元都有一个描述层、行和列的元组，如下所示：

cellid = [(0, 0, 0), (0, 1, 1), (0, 1, 2), ...]

不是包含层、行、列值的一维数组的元组，这是您所拥有的：

cellid = (array[0 0 0 0 0 0 0 ...], array[0 1 1 1 1 1 1 1 ...], array[0 1 2 3 4 5 ...])