def METHOD_NAME(self, variables):
    domain = self.domain
    delta_phi = variables[
        f"X-averaged {domain} electrode surface potential difference [V]"
    ]
    delta_phi_init = self.domain_param.prim.U_init
    self.initial_conditions = {delta_phi: delta_phi_init}