Question

如何在Flopy MF6的横截面图中绘制排放矢量？

我知道此图用于平面视图：

quiver = mapview.plot_specific_discharge(spdis[0])

我尝试使用以下代码获取特定的放电，但出现错误：

AttributeError: module 'flopy.utils.postprocessing' has no attribute 'get_specific_discharge'

代码：

hds_file = os.path.join(workspace, headfile)
hds = flopy.utils.binaryfile.HeadFile(hds_file)

cbb_file = os.path.join(workspace, budgetfile)
cbb = flopy.utils.CellBudgetFile(fname, precision='double')

q = flopy.utils.postprocessing.get_specific_discharge(gwf,cbb_file, precision='single', kstpkper=(0,1),
                                                  hdsfile=hds_file, position='centers')

Answer 1

对我来说（使用 mf6）在横截面上绘制特定的排放量是这样的：

读取 cbc 输出：

    cbcdobj = flopy.utils.binaryfile.CellBudgetFile(path, precision='double', verbose=True)
    SPDIS = cbcdobj.get_data(kstpkper=kstpkper, text='DATA-SPDIS')[0]

在使用 mf2005 时，您可能需要使用 'verbose=False' 和 'precision=single'。

然后是横截面：

    cros_mp=flopy.plot.PlotCrossSection(model=gwf, line={'row': 200})  
    cros_mp.plot_specific_discharge(SPDIS)

请注意，无法在不规则的横截面（“线”，而不是“行”或“列”）上绘制特定排放量。

软盘MF6在横截面图中绘制比放电

1 个答案: