def METHOD_NAME(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)