flopy - 如何在软盘中获得带有BC输出的cbc
问题描述
我正在尝试在“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)
解决方案
Ok, given your reply to my comment let me expand a bit here. Given the following name file:
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 will storge the GHB
package in memory on unit 23 and output to the cell-budget file on unit 53. The IGHBCB
variable of the GHB
package is used to push output to the cell-budget file; from the online guide to MODFLOW:
Therefore, inside of your GHB
package file, the second entry on the first line should also be "53" to send the output to the cell-budget file:
# GHB package for MODFLOW-NWT, generated by Flopy.
43183 53 NOPRINT
43183 0 # stress period 1
推荐阅读
- node.js - 将密钥文件与在 Kubernetes 集群上运行的应用程序一起使用
- testing - clojure.spec 检查是否会产生错误的输入?
- reactjs - 无法加载 PNG 图片 | NextJS
- firebase - 如何在 Firebase 中构造数据?
- python - 用相同的时间过滤一个带有itens的列表,并在Python中只放一个
- python-3.x - 结合 pandas df 列和行以获取字符串中的对 - python
- php - 网站联系表和 Google Workspace
- python - 在 Panda Python 中删除行中的非数字
- c# - 如何在列表中使用 IndexOf
- asp.net-core - 如何在 .Net 5 OAuth API 上授权不同的角色