ایجاد مدل MODFLOW با کد پایتون

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

ایجاد مدل MODFLOW

یکی از موارد مفید در مورد ایجاد مدل در پایتون این است که فرآیندی بسیار آسان داشته که با تغییر یک یا دو جزء، به طور کلی شبکه برای مدل شما را تغییر می دهد. بنابراین در این مثال، ما اسکریپت پایتون را طراحی می کنیم تا تعداد لایه ها، ستون ها و ردیف ها به راحتی تغییر کند.

ما می توانیم مدل MODFLOW بسیار ساده ای را تهیه کنیم که یک بسته پایه (BAS)، دیسک ورودی دیسکاوری (DIS)، بسته LPF یا package property، کنترل خروجی (OC) و حل کننده موتور گرادیان (PCG) را ایجاد کند. هر یک از این موارد فایل ورودی خود را دارد، که به صورت خودکار بوسیله FloPy ایجاد می شود.

گسسته سازی

با ایجاد شیء مدل FloPy ما به صورت زیر شروع می کنیم:

# Assign name and create modflow model object
modelname = 'tutorial1'
mf = flopy.modflow.Modflow(modelname, exe_name='mf2005')

در گام بعد، اجازه دهید با تعیین دامنه مدل ما با ایجاد یک شبکه MODFLOW به توسعه دامنه ادامه دهیم:

# Model domain and grid definition
Lx = 1000.
Ly = 1000.
ztop = 0.
zbot = -50.
nlay = 1
nrow = 10
ncol = 10
delr = Lx/ncol
delc = Ly/nrow
delv = (ztop - zbot) / nlay
botm = np.linspace(ztop, zbot, nlay + 1)

با استفاده از این اطلاعات، اکنون می توانیم شیء گسسته را با وارد کردن موارد زیر ایجاد کنیم:

# Create the discretization object
dis = flopy.modflow.ModflowDis(mf, nlay, nrow, ncol, delr=delr, delc=delc,
                               top=ztop, botm=botm[1:])

سوال واضح در این مرحله این است که چگونه می توانیم از این مفهوم ناشناخته به نام flopy.modflow.ModflowDis تشخیص دهیم که چه استدلال هایی لازم است؟ خوشبختانه یک صفحه راهنما برای هر یک از اشیاء مدل وجود دارد. صفحه برای فایل ورودی DIS در flopy.modflow.mfdis قرار دارد.

ادامه دارد...

بخش بعدی را اینجا بیابید.