Как загрузить существующую модель MODFLOW-USG с помощью flopy - PullRequest
Новая
       10

Как загрузить существующую модель MODFLOW-USG с помощью flopy

0 голосов
/ 07 февраля 2019

Я пытаюсь выучить флопи, и я попытался загрузить один из тестовых случаев MODFLOW-USG.Тестовый пример - это однослойная модель с вложенной сеткой и пакетом GNC (пример 01B).Я не изменял файлы тестовых примеров, распространяемые с помощью MODFLOW-USG, и тестовый пример хорошо работал с файлом run.bat.

Сначала я попытался следовать примеру для загрузки файла имени с помощью модуль flopy.modflow.mf : 

m = flopy.modflow.Modflow.load('flow.nam', version='mfusg')

Но я получил следующее ValueError: 

File "C:\Users\Kent\Anaconda3\lib\site-packages\flopy\modflow\mf.py", line 758, in load
  ml.check(f='{}.chk'.format(ml.name), verbose=ml.verbose, level=0)

File "C:\Users\Kent\Anaconda3\lib\site-packages\flopy\mbase.py", line 1121, in check
  level=level - 1)

File "C:\Users\Kent\Anaconda3\lib\site-packages\flopy\modflow\mfbas.py", line 185, in check
  neighbors = get_neighbors(self.ibound.array)

File "C:\Users\Kent\Anaconda3\lib\site-packages\flopy\utils\check.py", line 590, in get_neighbors
  nk, ni, nj = a.shape

ValueError: not enough values to unpack (expected 3, got 1)

Затем я попытался загрузить пакет DISU, следуя примеру для модуль flopy.modflow.mfdisu : 

m = flopy.modflow.Modflow(version='mfusg')
disu = flopy.modflow.ModflowDisU.load('flow.disu', m)

Но это дало мне следующее TypeError: 

File "C:\Users\Kent\Anaconda3\lib\site-packages\flopy\modflow\mfdisu.py", line 622, in load
  fahl = Util2d.load(f, model, (n,), np.float32, 'fahl', ext_unit_dict)

File "C:\Users\Kent\Anaconda3\lib\site-packages\flopy\utils\util_array.py", line 2678, in load
  ext_unit = ext_unit_dict[cr_dict['nunit']]

TypeError: 'NoneType' object is not subscriptable

Как правильно загрузить существующийМодель MODFLOW-USG с дискетой?

...