def METHOD_NAME():
    #-------- EmbeddeFrac File path
    hdf5File1Path = "displacementJump_embeddedFrac.hdf5"
    # Read HDF5
    hf = h5py.File(hdf5File1Path, 'r')
    jump = hf.get('displacementJump')
    jump = np.array(jump)
    aperture_EmbeddeFrac = jump[0, :, 0]
    x = hf.get('displacementJump elementCenter')
    loc_EmbeddeFrac = x[0, :, 1]
    #-------- Lagrange Contact File path
    hdf5File2Path = "displacementJump_contactMechanics.hdf5"
    # Read HDF5
    hf = h5py.File(hdf5File2Path, 'r')
    jump = hf.get('displacementJump')
    jump = np.array(jump)
    aperture_Contact = jump[0, :, 0]
    x = hf.get('displacementJump elementCenter')
    loc_Contact = x[0, :, 1]
    #-------- HydroFrac File path
    hdf5File3Path = "displacementJump_hydroFrac.hdf5"
    # Read HDF5
    hf = h5py.File(hdf5File3Path, 'r')
    jump = hf.get('elementAperture')
    jump = np.array(jump)
    aperture_HydroFrac = jump[0, :]
    x = hf.get('elementAperture elementCenter')
    loc_HydroFrac = x[0, :, 1]
    #-------- Extract info from XML
    xmlFilePath = "../../../../../../../inputFiles/efemFractureMechanics/Sneddon_embeddedFrac_base.xml"
    mechanicalParameters = getMechanicalParametersFromXML(xmlFilePath)
    appliedPressure = getFracturePressureFromXML(xmlFilePath)
    # Get length of the fracture
    length, origin = getFractureLengthFromXML(xmlFilePath)
    # Initialize Sneddon's analytical solution
    sneddonAnalyticalSolution = Sneddon(mechanicalParameters, length, appliedPressure)
    # Plot analytical (continuous line) and numerical (markers) aperture solution
    x_analytical = np.linspace(-length, length, 501, endpoint=True)
    aperture_analytical = np.empty(len(x_analytical))
    i = 0
    for xCell in x_analytical:
        aperture_analytical[i] = sneddonAnalyticalSolution.computeAperture(xCell)
        i += 1
    fsize = 30
    msize = 15
    lw = 6
    fig, ax = plt.subplots(1, figsize=(16, 12))
    cmap = plt.get_cmap("tab10")
    N1 = 1
    ax.plot(x_analytical, aperture_analytical * 1.0e3, color='k', label='Analytical Solution', lw=lw)
    ax.plot(loc_EmbeddeFrac[0::N1],
            aperture_EmbeddeFrac[0::N1] * 1.0e3,
            'o',
            color=cmap(0),
            label='Embedded Fracture',
            markersize=msize * 0.6,
            alpha=0.8)
    ax.plot(loc_Contact[0::N1],
            aperture_Contact[0::N1] * 1.0e3,
            's',
            color=cmap(4),
            label='Lagrange Contact',
            markersize=msize * 0.6,
            alpha=0.6)
    ax.plot(loc_HydroFrac[0::N1],
            aperture_HydroFrac[0::N1] * 1.0e3,
            'X',
            color=cmap(1),
            label='HydroFrac Solver',
            markersize=msize * 0.6,
            alpha=0.4)
    ax.grid()
    ax.set_xlabel('Fracture Length [m]', size=fsize, weight="bold")
    ax.set_ylabel('Fracture Aperture [mm]', size=fsize, weight="bold")
    ax.legend(bbox_to_anchor=(0.5, 0.2), loc='center', borderaxespad=0., fontsize=fsize)
    ax.xaxis.set_tick_params(labelsize=fsize)
    ax.yaxis.set_tick_params(labelsize=fsize)
    plt.show()