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

Question

我正在创建一个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 = 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', 'ALL'), ],
                        )

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

Answer 1

你需要初始化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", (0,), (0,)),
                ("Q2", "maw", (0,), (1,)),
                ("Q3", "maw", (0,), (2,)),
                ]}
maw.obs.initialize(filename=obs_file, digits=10, print_input=True, continuous=obs_dict)