def METHOD_NAME(Obs_data,Obs_xml,filename):
    pressure_amanzi = []
    coordinates = []
    mymodel = model_linear_head_head_1d.createFromXML(filename)
    for obs in Obs_data.observations.values():
        coordinates.append([obs.coordinate[0], obs.coordinate[2]])
        pressure_amanzi.append(str(obs.data).rstrip(']').lstrip('['))
    pressure_analytic = list(mymodel.pressure(numpy.array(coordinates)))
    x = prettytable.PrettyTable(["x [m]", "z [m]", "Analytic [Pa]","Amanzi [Pa]"])
    x.padding_width = 2
    x.hrules = 1
    x.horizontal_header_char="="
    for coords, p_analytic, p_amanzi in zip(coordinates,pressure_analytic,pressure_amanzi):
        x.add_row([coords[0],coords[1],"%.4f" % float(p_analytic),"%.4f" % float(p_amanzi)])
    if os.path.exists("table_values.txt"):
        os.remove("table_values.txt")
    table_file = open("table_values.txt", "w+")
    table_file.write('.. tabularcolumns:: ' + '|R|C|C|C|' + '\n\n')
    table_file.write(x.get_string(sortby="x [m]"))
    table_file.write('\n')
    table_file.close()