如何在 FloPy Modflow6 中输出所有时间步长的 MAW 水头值？

如何解决如何在 FloPy Modflow6 中输出所有时间步长的 MAW 水头值？

我正在创建一个 MAW 井，并希望将其用作观察井，以便稍后将其与现场数据进行比较，应该对其进行多层筛选。但是，我只在输出文件的最后一个时间步的井中获得水头值。关于如何在输出中获取所有时间步长的任何想法？

FloPy 手册说它需要在输出控制中，但我不知道如何做到这一点：

print_head (boolean) –   print_head (boolean) keyword to indicate that the list of multi-aquifer well heads will be printed to the listing file for every stress period in which “HEAD PRINT” is specified in Output Control. If there is no Output Control option and PRINT_HEAD is specified,then heads are printed for the last time step of each stress period.

在 MODFLOW6 手册中，我看到可以进行连续输出： modflow6

我的 MAW 定义如下：

maw = flopy.mf6.ModflowGwfmaw(gwf,nmawwells=1,packagedata=[0,Rwell,minbot,wellhead,'MEAN',OBS1welllayers],connectiondata=OBS1connectiondata,perioddata=[(0,'STATUS','ACTIVE')],flowing_wells=False,save_flows=True,mover=True,flow_correction=True,budget_filerecord='OBS1wellbudget',print_flows=True,print_head=True,head_filerecord='OBS1wellhead',)

我的输出控件如下所示：

oc = flopy.mf6.ModflowGwfoc(gwf,budget_filerecord=budget_file,head_filerecord=head_file,saverecord=[('HEAD','ALL'),('BUDGET',],)

希望这一切都清楚了，有人可以帮助我，谢谢！

解决方法

您需要初始化 MAW 观测文件……这不是在 OC 包中完成的。

您可以在此处的 MF6 文档中找到三个 MAW 示例的脚本： https://github.com/MODFLOW-USGS/modflow6-examples/tree/master/notebooks

看起来像这样：

obs_file = "{}.maw.obs".format(name)
csv_file = obs_file + ".csv"
obs_dict = {csv_file: [
                ("head","head",(0,)),("Q1","maw",),("Q2",(1,("Q3",(2,]}
maw.obs.initialize(filename=obs_file,digits=10,print_input=True,continuous=obs_dict)