0

我正在尝试将河流添加到我的模型中，但是在我运行它之后会出现这个警告：

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。我提到的是否与之前的警告有关？

很感谢任何形式的帮助。

4

1 回答 1

1

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

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

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 ...])
于 2020-07-13T13:52:28.263 回答