اجرای مدل سازی MODFLOW با برنامه نویسی متن باز

۲۶ آبان ۹۶ ، ۲۱:۱۵ نويسنده: بهزاد سرهادی

اجرای مدل سازی

Flopy همچنین می تواند برای اجرای مدل مورد استفاده قرار گیرد. شی مدل (mf در این مثال) یک روش متصل دارد که مدل را اجرا می کند. برای این کار، برنامه MODFLOW باید در جایی در مسیر سیستم یا در داخل دایرکتوری کار قرار بگیرد. در این مثال، ما مشخص کرده ایم که نام برنامه اجرایی 'mf2005' است. برای اجرای این مدل زیر را امتحان کنید:

# Run the MODFLOW model
success, buff = mf.run_model()

در اینجا ما از run_model استفاده کردیم و همچنین می توانیم مقادیر مشخص شده برای کلمات کلیدی اختیاری را ساکت، مکث و گزارش دهیم.

پس پردازش نتایج

حالا که ما مدل MODFLOW را با موفقیت ساختیم و اجرا کردیم، می توانیم نتایج را ببینیم. MODFLOW سرهای شبیه شده را به یک فایل خروجی داده باینری می نویسد. ما نمی توانیم با استفاده از یک ویرایشگر متن به این سر نگاه کنیم، اما فاقد ابزار باینری است که می تواند برای خواندن سر استفاده شود. اظهارات زیر فایل سرفیسم را بخواند و یک طرح از سرهای شبیه سازی شده برای لایه 1 ایجاد کند:

import matplotlib.pyplot as plt
import flopy.utils.binaryfile as bf
plt.subplot(1,1,1,aspect='equal')
hds = bf.HeadFile(modelname+'.hds')
head = hds.get_data(totim=1.0)
levels = np.arange(1,10,1)
extent = (delr/2., Lx - delr/2., Ly - delc/2., delc/2.)
plt.contour(head[0, :, :], levels=levels, extent=extent)
plt.show()