def plotTheisAnalytic(filename, cmap, axes1, Obs_xml, Obs_data):
mymodel = model_theis_isotropic_1d.createFromXML(filename)
tindex = numpy.arange(125)
times = []
table_values = []
press_amanzi = []
press_analytic = []
axes1.set_ylabel('Drawdown [m]')
axes1.set_xlabel('Time after pumping [s]')
axes1.set_title('Drawdown vs Time after Pumping')
for i in tindex:
times.append(1 + math.exp(float(i) * (i + 1) / (10.3 * len(tindex))))
for rad in mymodel.r:
r = abs(rad)
color = cmap[rad]
drawdown = mymodel.runForFixedRadius(times, r)
axes1.plot(times, drawdown, label='$r=%0.1f m$' % r, c=color)
press_analytic.append([r, drawdown])
for obs in Obs_data.observations.values():
press_amanzi.append(obs.data)
for press1, press2 in zip(press_analytic, press_amanzi):
table_values.append([press1[0], press1[1], press2])
axes1.legend(title='Theis Solution',
loc='lower right',
fancybox=True,
shadow=True)
def MakeTable(Obs_data,Obs_xml,filename):
#### ==== Initial height of water table is set to 20.0 meters ==== ####
#### ==== The Table is set to generate at a coordinate of (55,y,z) *only* ==== ####
#### ==== Modify if statement of another coordinate is desired ==== ####
drawdown_amanzi = []
coordinates = []
mymodel = model_theis_isotropic_1d.createFromXML(filename)
for obs in Obs_data.observations.itervalues():
if obs.coordinate[0] == -55.0:
coordinates.append([abs(obs.coordinate[0]), obs.coordinate[2]])
for press in obs.data:
pres0 = 101325 -9806.65 * (obs.coordinate[2] - 20)
pres_drop = (pres0 - press)
drawdown = pres_drop / 9806.65
drawdown_amanzi.append(drawdown)
drawdown_analytic = mymodel.runForFixedRadius(mymodel.times,coordinates[0][0])
x = prettytable.PrettyTable(["time [s]","r [m]", "z [m]", "Analytic [m]","Amanzi [m]"])
x.padding_width = 2
x.hrules = 1
for time,d_analytic, d_amanzi in zip(mymodel.times,drawdown_analytic, drawdown_amanzi):
x.add_row([time,coordinates[0][0],coordinates[0][1],"%.10f" % d_analytic,"%.10f" % d_amanzi])
if os.path.exists("table_values_theis.txt"):
os.remove("table_values_theis.txt")
table_file = open("table_values_theis.txt", "w+")
table_file.write(x.get_string())
table_file.close()