def METHOD_NAME(self, parameters):
    sys.path.insert(
        1,
        parameters.parthenon_path
        + "/scripts/python/packages/parthenon_tools/parthenon_tools",
    )
    try:
        import phdf
        from phdf_diff import compare
    except ModuleNotFoundError:
        print("Couldn't find modules to load Parthenon hdf5.")
        return False
    all_pass = True
    # Reflection: First comparing initial condition to final state
    res = compare(
        ["advection.reflecting.00000.phdf", "advection.reflecting.final.phdf"],
        check_metadata=False,
        quiet=True,
    )
    if res != 0:
        print("Double reflection test failed: output start != end ")
        all_pass = False
    # Compare that data in between is actually different.
    res = compare(
        ["advection.reflecting.00000.phdf", "advection.reflecting.00002.phdf"],
        check_metadata=False,
        quiet=True,
    )
    if res == 0:
        print("Double reflection test failed: data identical during sim.")
        all_pass = False
    # Periodic: First comparing initial condition to final state
    res = compare(
        ["advection.periodic.00000.phdf", "advection.periodic.final.phdf"],
        check_metadata=False,
        quiet=True,
    )
    if res != 0:
        print("Fully periodic test failed: output start != end ")
        all_pass = False
    # Compare that data in between is actually different.
    res = compare(
        ["advection.periodic.00000.phdf", "advection.periodic.00001.phdf"],
        check_metadata=False,
        quiet=True,
    )
    if res == 0:
        print("Fully periodic test failed: data identical during sim.")
        all_pass = False
    # Outflow: Check there's sth in the beginning...
    outflow_data_initial = phdf.phdf("advection.outflow.00000.phdf")
    advected_initial = outflow_data_initial.Get("advected")
    if np.sum(advected_initial) != 60.0:
        print("Unexpected initial data for outflow test.")
        all_pass = False
    # ... and nothing at the end
    outflow_data_final = phdf.phdf("advection.outflow.final.phdf")
    advected_final = outflow_data_final.Get("advected")
    if np.sum(advected_final) != 0.0:
        print("Some 'advected' did not leave the box in outflow test.")
        all_pass = False
    return all_pass