如何在软盘中获得带有BC输出的cbc答案

【问题标题】：how to get cbc with BC output in flopy如何在软盘中获得带有BC输出的cbc
【发布时间】：2019-07-29 19:10:09
【问题描述】：

我正在尝试在“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)

【问题讨论】：

GHB 和 DRN 包文件中的输出单元号是否与单元预算文件的输出单元号匹配？如果您可以发布您的 OC 包、名称文件和相关边界包文件的示例，它将帮助其他人解决您的问题。
您好，感谢您的快速回复！不，他们不匹配。因此，对于教程 2 中的示例，.cbc 文件的单元为 53，.ghb 的单元为 23（基于列表文件）。就像在线示例一样，脚本中没有指定单元编号。我接受它，它们应该被指定？

标签： flopy

【解决方案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 会将GHB 包存储在单元23 的内存中，并输出到单元53 的单元预算文件。GHB 包的IGHBCB 变量用于将输出推送到单元预算文件;来自online guide to MODFLOW：

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

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

【讨论】：