"""
supports multi-well extraction & with constant head boundaries from your input
"""
from mymf_v3 import mymf

# single wells
init_head = 1 # initial head distribution = 1 m
model = mymf(init_head)
well_rcs = [[4,4]] # nrow = 10, ncol = 10
Qs = [-903.166236] # we go ahead and use the results from optimization
model.run(well_rcs,Qs)
model.plot()

print('min head is %f' % (model.minhead()))