0

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

1 回答 1

0

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:

enter image description here

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
于 2019-07-31T22:23:36.147 回答