在flopy中创建MF-USG模型的示例

import flopy
import numpy as np

modelname = "example_model"
model = flopy.modflow.Modflow(modelname, exe_name="mfusg")

nlay = 1  # 模型的层数
nrow = 10  # 模型的行数
ncol = 10  # 模型的列数
delr = 100.0  # 列宽
delc = 100.0  # 行宽
top = 0.0  # 顶部高程
botm = -10.0  # 底部高程
dis = flopy.modflow.ModflowDis(model, nlay=nlay, nrow=nrow, ncol=ncol, delr=delr, delc=delc, top=top, botm=botm)

hk = 10.0  # 渗透系数
laytyp = 1  # 层类型
upw = flopy.modflow.ModflowUpw(model, laytyp=laytyp, hk=hk)

ibound = np.ones((nlay, nrow, ncol), dtype=np.int32)  # 边界条件
strt = np.ones((nlay, nrow, ncol), dtype=np.float32)  # 初始水位
bas = flopy.modflow.ModflowBas(model, ibound=ibound, strt=strt)

nper = 1  # 模拟的时间步数
perlen = [1.0]  # 时间步长
nstp = [1]  # 每个时间步的计算步数
stress_period_data = {}  # 应力周期数据
for kper in range(nper):
    stress_period_data[kper] = [(0, 'ET', 0.1)]  # 定义蒸散发边界条件
bas = flopy.modflow.ModflowBas(model, ibound=ibound, strt=strt)

model.write_input()
model.run_model()