def METHOD_NAME(self, variables):
    domain, Domain = self.domain_Domain
    eps_solid_init = self.domain_param.prim.epsilon_s
    if self.x_average is True:
        eps_solid_xav = variables[
            f"X-averaged {domain} electrode active material volume fraction"
        ]
        self.initial_conditions = {eps_solid_xav: pybamm.x_average(eps_solid_init)}
    else:
        eps_solid = variables[f"{Domain} electrode active material volume fraction"]
        self.initial_conditions = {eps_solid: eps_solid_init}
    lli_due_to_lam = variables[
        "Loss of lithium due to loss of active material "
        f"in {domain} electrode [mol]"
    ]
    self.initial_conditions[lli_due_to_lam] = pybamm.Scalar(0)