def METHOD_NAME(m):
    m.fs.feed = Feed(property_package=m.fs.stoich_softening_thermo_params)
    comp_list = ["Na", "Ca", "Mg", "SO4", "Cl"]
    feed_flow_mass = 1  # kg/s
    feed_mass_frac_comp = {
        "Na": 11122e-6,
        "Ca": 382e-6,
        "Mg": 1394e-6,
        "SO4": 2136e-6,
        "Cl": 20316.88e-6,
    }
    feed_mass_frac_comp["H2O"] = 1 - sum(x for x in feed_mass_frac_comp.values())
    mw_comp = {
        "H2O": 18.015e-3,
        "Na": 22.990e-3,
        "Ca": 40.078e-3,
        "Mg": 24.305e-3,
        "SO4": 96.06e-3,
        "Cl": 35.453e-3,
    }
    feed_flow_mol_comp = {}
    for j in feed_mass_frac_comp:
        feed_flow_mol_comp[j] = feed_flow_mass * feed_mass_frac_comp[j] / mw_comp[j]
    feed_mole_frac_comp = {}
    for j in feed_flow_mol_comp:
        feed_mole_frac_comp[j] = feed_flow_mol_comp[j] / sum(
            feed_flow_mol_comp[jj] for jj in feed_flow_mol_comp
        )
    m.fs.feed.properties[0].mole_frac_comp["Ca(HCO3)2"].fix(feed_mole_frac_comp["Ca"])
    m.fs.feed.properties[0].mole_frac_comp["Ca(OH)2"].fix(0)
    m.fs.feed.properties[0].mole_frac_comp["CaCO3"].fix(0)
    m.fs.feed.properties[0].mole_frac_comp["H2O"].fix(feed_mole_frac_comp["H2O"])
    m.fs.feed.properties[0].mole_frac_comp["Mg(HCO3)2"].fix(feed_mole_frac_comp["Mg"])
    m.fs.feed.properties[0].mole_frac_comp["Mg(OH)2"].fix(0)
    m.fs.feed.properties[0].mole_frac_comp["NaCl"].fix(feed_mole_frac_comp["Na"])
    m.fs.feed.properties[0].mole_frac_comp["SO4_2-"].fix(feed_mole_frac_comp["SO4"])
    m.fs.feed.properties[0].mole_frac_comp["Cl_-"].fix(
        feed_mole_frac_comp["Cl"] - feed_mole_frac_comp["Na"]
    )
    m.fs.feed.properties[0].flow_mol.fix(55.56)
    m.fs.feed.properties[0].pressure.fix(101325)
    m.fs.feed.properties[0].temperature.fix(298)