def METHOD_NAME(filename, cmap, axes1, Obs_xml, Obs_data):
    # Instantiate the analytic solution
    mymodel = model_linear_head_head_1d.createFromXML(filename)
    table_values = []
    # Create a set of points to plot the solution
    x = numpy.linspace(mymodel.x_0,mymodel.x_1,11)
    coords = numpy.zeros((11,2))
    coords[:,0] = x
    # Plot a line for each z-coordinate in the observations
    for (z_val, color) in cmap.items():
        coords[:,1] = z_val
        pressure = mymodel.pressure(coords)
        axes1.plot(x,pressure,color,label='$z = %0.2f $'%z_val)
        axes1.legend(loc="upper right" , fancybox = True , shadow = True)