مدل سازی آبخوان با کد اسکریپت پایتون
بسته پایه
در گام بعد ما می توانیم جسم FloPy را ایجاد کنیم که MODFLOW Basic Pack را نشان می دهد. جزئیات در کلاس BAS فلاپی در: flopy.modflow.mfbas موجود است. برای این مدل ساده، مقادیر سطح پایه 10 و 0 به ترتیب به ستون اول و آخر مدل (در همه لایه ها) اختصاص داده می شود. کد پایتون برای انجام این کار:
# Variables for the BAS package ibound = np.ones((nlay, nrow, ncol), dtype=np.int32) ibound[:, :, 0] = -1 ibound[:, :, -1] = -1 strt = np.ones((nlay, nrow, ncol), dtype=np.float32) strt[:, :, 0] = 10. strt[:, :, -1] = 0. bas = flopy.modflow.ModflowBas(mf, ibound=ibound, strt=strt)
Packet Flow Layer Property
اطلاعات در مورد کلاس LPF مدل FloPy در: flopy.modflow.mflpf قرار دارد. مقادیر 10 برای هدایت هیدرولیکی افقی و عمودی اختصاص داده شده است:
# Add LPF package to the MODFLOW model lpf = flopy.modflow.ModflowLpf(mf, hk=10., vka=10.)
از آنجا که ما یک مقدار برای laytyp مشخص نمی کنیم، Flopy از مقدار پیش فرض 0 استفاده می کند، یعنی این مدل محدود خواهد شد.
کنترل خروجی
جزئیات کلاس OC مدل FloPy در: flopy.modflow.mfoc موجود است. در اینجا می توانیم از تنظیمات پیش فرض OC با استفاده از موارد زیر استفاده کنیم:
# Add OC package to the MODFLOW model oc = flopy.modflow.ModflowOc(mf)
Packaged Gradient Conjugate Preconditioned
اطلاعات در مورد کلاس PCG مئل FloPy در: flopy.modflow.mfpcg قرار دارد. تنظیمات پیش فرض استفاده شده توسط FloPy با مشخص کردن دستورات زیر استفاده می شود:
نوشتن فایل های داده MODFLOW
فایل های داده ورودی MODFLOW به سادگی با انتشار موارد زیر نوشته می شوند:
# Add PCG package to the MODFLOW model pcg = flopy.modflow.ModflowPcg(mf)
ادامه دارد...
بخش بعدی را اینجا بیابید.
شناسه تلگرام مدیر سایت: SubBasin@
نشانی ایمیل: behzadsarhadi@gmail.com
(سوالات تخصصی را در گروه تلگرام ارسال کنید)
_______________________________________________________
نظرات (۰)