def METHOD_NAME(
    before: MeasurementData, after: MeasurementData
) -> float:
    """Calculate volume of water."""
    # TODO: actually calculate volume
    avg_env = get_average_reading([before.environment, after.environment])
    water_density_at_this_temperature_kg_m3 = sum(
        [
            CONSTANT_WATER_DENSITY["1"],
            CONSTANT_WATER_DENSITY["2"] * pow(10, -2) * pow(avg_env.celsius_liquid, 1),
            CONSTANT_WATER_DENSITY["3"] * pow(10, -3) * pow(avg_env.celsius_liquid, 2),
            CONSTANT_WATER_DENSITY["4"] * pow(10, -5) * pow(avg_env.celsius_liquid, 3),
            CONSTANT_WATER_DENSITY["5"] * pow(10, -7) * pow(avg_env.celsius_liquid, 4),
        ]
    )
    liquid_density_kg_m3 = (
        RELATIVE_DENSITY_WATER * water_density_at_this_temperature_kg_m3
    )
    # equations in ISO use hPa, so sticking with that
    air_pressure_hpa = avg_env.pascals_air / 100
    air_density_kg_m3 = (
        (CONSTANT_AIR_DENSITY_K["1"] * air_pressure_hpa)
        + (
            avg_env.humidity_air
            * (
                (CONSTANT_AIR_DENSITY_K["2"] * avg_env.celsius_air)
                + CONSTANT_AIR_DENSITY_K["3"]
            )
        )
    ) / (avg_env.celsius_air + CONSTANT_TEMPERATURE_TA0)
    z_factor = (
        (1.0 / CONSTANT_SCALE_CALIBRATED_DENSITY_KG_M3)
        * (
            (CONSTANT_SCALE_CALIBRATED_DENSITY_KG_M3 - air_density_kg_m3)
            / (liquid_density_kg_m3 - air_density_kg_m3)
        )
        * 1000.0
    )
    liquid_grams = abs(before.grams_average - after.grams_average)
    liquid_micro_liters = liquid_grams * z_factor * 1000.0
    return liquid_micro_liters