def METHOD_NAME(b, t):
    return b.MgCl2_flowrate[t] == (
        b.magnesium_chloride_dosage
        * pyunits.convert(
            b.properties_byproduct[t].flow_mass_comp["phosphorus"],
            to_units=pyunits.kg / pyunits.hour,
        )
    )