# --problem size
nlay, nrow, ncol = 1, 41, 40
# --default data if command line argument not defined for variable
head_file = "..\\data\\CoastalAquifer.hds"
# --get available times in the head file
get_cell = 1
headObj = mfb.MODFLOW_Head(nlay, nrow, ncol, head_file)
t = headObj.get_gage(get_cell)
ntimes = t.shape[0]
mf_times = np.zeros((ntimes), np.float)
for i in range(0, ntimes):
mf_times[i] = t[i, 0]
# --create a time-series plot of head
headObj = mfb.MODFLOW_Head(nlay, nrow, ncol, head_file)
get_cell = mfb.icrl_from_kij(1, 18, 31, nlay, nrow, ncol)
t = headObj.get_gage(get_cell)
ztf = figure(figsize=(4, 2), facecolor="w")
ztf.subplots_adjust(wspace=0.2, hspace=0.2, left=0.15, right=0.95, bottom=0.15, top=0.95)
ax = ztf.add_subplot(1, 1, 1)
hp = ax.plot(mf_times / 365.25 - 100.0, t[:, 1])
ax.set_xlim(0, 50)
xlabel("Time since the beginning of groundwater withdrawls (years)")
ylabel("Head (m)")
output_name = "..\\figures.MF\\Cell_{0}.png".format(get_cell)
ztf.savefig(output_name, dpi=300)
# --create a time-series plot of drawdown
h0 = t[:, 1].max()
ztf = figure(figsize=(4, 2), facecolor="w")
ztf.subplots_adjust(wspace=0.2, hspace=0.2, left=0.15, right=0.95, bottom=0.15, top=0.95)
ax = ztf.add_subplot(1, 1, 1)
swr_qaq04[n] = ce1[n,iextbf] #/ 86400. #SWRObj = mfb.SWR_Record(2,SWR_file) #ce1 = SWRObj.get_gage(iobs04) ##print 'qaq04 layer 2\n', ce1 #nt = np.shape(ce1)[0] #for n in range(0,nt): # swr_qaq04[n] = swr_qaq04[n] + ce1[n,iextbf] #/ 86400. #--MODFLOW results #--heads ncol = 6 nrow = 6 nlay = 2 iobs1node = mfb.icrl_from_kij(1,3,4,nlay,nrow,ncol) iobs2node = mfb.icrl_from_kij(1,4,5,nlay,nrow,ncol) iobs3node = mfb.icrl_from_kij(1,5,6,nlay,nrow,ncol) iobs4node = mfb.icrl_from_kij(1,6,5,nlay,nrow,ncol) #print 'MODFLOW observation location', iobs1node, iobs2node, iobs3node, iobs4node HEAD_file = '..\\SWRSample05\\Results\\SWRSample05.hds' HEADObj = mfb.MODFLOW_Head(nlay,ncol,nrow,HEAD_file) mf1 = HEADObj.get_gage(iobs1node) nmft = np.shape(mf1)[0] mf_t = np.zeros(nmft) mf_h01 = np.zeros(nmft) t = np.zeros(nmft) for n in range(0,nmft): mf_t[n] = mf1[n,0]
