如何在flopy中使用BC输出获取cbc

Question

我正在尝试在“ FloPy”（v'3.2.12'）中建立模型，但是我在modflow输出文件上遇到了一些麻烦。遵循开发者网站（https://modflowpy.github.io/flopydoc/tutorial2.html）上提供的“ turorial2”时，一切似乎正常。但是，cbc文件不包含有关常规头部边界条件的信息。我在排水边界条件上也遇到了同样的问题。运行模型之前，我是否需要更改或设置一些其他设置？还是有其他方法可以从GHB边界条件中读取费率信息？

仅供参考- 到目前为止，我已经尝试过更改ModflowOc设置，并且还使用其他求解器来运行模型。

感谢您的帮助，

有关设置的一些其他信息（应该与教程相同）：

mf = flopy.modflow.Modflow(modelname, exe_name='mf2005')


ghb = flopy.modflow.ModflowGhb(mf, stress_period_data=stress_period_data)


oc = flopy.modflow.ModflowOc(mf, stress_period_data=stress_period_data,
                             compact=True)

Answer 1

好吧，鉴于您对我的评论的答复，让我在这里扩大一下。给定以下名称文件：

LIST               2  ss-kper60.list
DIS               11  ss-kper60.dis 
OC                14  ss-kper60.oc 
BAS6              13  ss-kper60.bas6 
NWT               32  ss-kper60.nwt 
GHB               23  ss-kper60.ghb 
RCH               19  ss-kper60.rch 
UPW               31  ss-kper60.upw 
DATA(BINARY)      51  ss-kper60.hds REPLACE
DATA(BINARY)      53  ss-kper60.cbc REPLACE

MODFLOW将在单元23的存储器中存储GHB包，并在单元53上输出到单元预算文件。IGHBCB包的GHB变量用于推送输出到单元预算文件；来自online guide to MODFLOW：

因此，在您的GHB软件包文件中，第一行的第二个条目也应为“ 53”，以将输出发送到单元预算文件：

# GHB package for  MODFLOW-NWT, generated by Flopy.
     43183        53  NOPRINT
     43183         0 # stress period 1