bas = fm.ModflowBas(mf, ibound=IBOUND, strt=STRTHD)
lpf = fm.ModflowLpf(mf,
hk=HK,
vka=VKA,
chani=[1.e-20, 1.e-20],
sy=SY,
ss=SS,
laytyp=LAYTYP,
vkcb=VKCB,
ipakcb=53)
ghb = fm.ModflowGhb(mf, stress_period_data=GHB, ipakcb=53)
riv = fm.ModflowRiv(mf, stress_period_data=RIV, ipakcb=53)
drn = fm.ModflowDrn(mf, stress_period_data=DRN, ipakcb=53)
wel = fm.ModflowWel(mf, stress_period_data=SEEP, ipakcb=53)
rch = fm.ModflowRch(mf, nrchop=3, rech=RECH, ipakcb=53)
evt = fm.ModflowEvt(mf, nevtop=3, evtr=EVTR, ipakcb=53)
oc = fm.ModflowOc(mf, stress_period_data=OC, compact=True)
#pcg = fm.ModflowPcg(mf, mxiter=200, iter1=200, hclose=0.001, rclose=0.001)
sms = fm.ModflowSms(mf) #, mxiter=200, iter1=200, hclose=0.001, rclose=0.001)
#%% Write the model input files and run MODFLOW
mf.write_input()
success, mfoutput = mf.run_model(silent=False, pause=False)
print('Running success = {}'.format(success))
if not success:
raise Exception('MODFLOW did not terminate normally.')
#%% SHOWING RESULTS
#%% PLOT AVERAGE HEAD IN THE CROSS SECTION
Nx,
delr=dx,
delc=dy,
top=zTop,
botm=zBot,
laycbd=LAYCBD,
nper=NPER,
perlen=PERLEN,
nstp=NSTP,
steady=STEADY)
bas = fm.ModflowBas(mf, ibound=IBOUND, strt=STRTHD)
lpf = fm.ModflowLpf(mf, hk=HK, vka=VKA, sy=SY, ss=SS, laytyp=LAYTYP, vkcb=VKCB)
ghb = fm.ModflowGhb(mf, stress_period_data=GHB)
wel = fm.ModflowWel(mf, stress_period_data=WEL)
rch = fm.ModflowRch(mf, nrchop=3, rech=RECH)
evt = fm.ModflowEvt(mf, nevtop=3, evtr=EVTR)
oc = fm.ModflowOc(mf, stress_period_data=OC, compact=True)
pcg = fm.ModflowPcg(mf)
#%% Write the model input files and running MODFLOW
mf.write_input()
success, mfoutput = mf.run_model(silent=False, pause=False)
print('Running success = {}'.format(success))
if not success:
raise Exception('MODFLOW did not terminate normally.')
#%% SHOWING RESULTS
def get_package(self, _mf):
content = self.merge()
return mf.ModflowEvt(_mf, **content)