def METHOD_NAME(self):
    # Test isocontour properties for AMRGridData
    fields = ["density", "cell_mass"]
    units = ["g/cm**3", "g"]
    ds = fake_amr_ds(fields=fields, units=units, particles=16**3)
    dd = ds.all_data()
    q = dd.quantities["WeightedAverageQuantity"]
    rho = q(("gas", "density"), weight=("gas", "cell_mass"))
    dd.extract_isocontours(("gas", "density"), rho, "triangles.obj", True)
    dd.calculate_isocontour_flux(
        ("gas", "density"),
        rho,
        ("index", "x"),
        ("index", "y"),
        ("index", "z"),
        ("index", "dx"),
    )
    # Test error in case of ParticleData
    ds = fake_particle_ds()
    dd = ds.all_data()
    q = dd.quantities["WeightedAverageQuantity"]
    rho = q(("all", "particle_velocity_x"), weight=("all", "particle_mass"))
    with assert_raises(NotImplementedError):
        dd.extract_isocontours("density", rho, sample_values="x")