pybamm-team#4123 fix

pybamm/models/submodels/interface/sei/sei_growth.py

@@ -71,7 +71,8 @@ def get_fundamental_variables(self):
 
         L_inner, L_outer = Ls
 
-        SEI_option = getattr(self.options, domain)["SEI"]
+        domain_options = getattr(self.options, domain)
+        SEI_option = getattr(domain_options, self.phase)["SEI"]
         if SEI_option.startswith("ec reaction limited"):
             L_inner = 0 * L_inner  # Set L_inner to zero, copying domains
 
@@ -83,7 +84,8 @@ def get_coupled_variables(self, variables):
         param = self.param
         phase_param = self.phase_param
         domain, Domain = self.domain_Domain
-        SEI_option = getattr(self.options, domain)["SEI"]
+        domain_options = getattr(self.options, domain)
+        SEI_option = getattr(domain_options, self.phase)["SEI"]
         T = variables[f"{Domain} electrode temperature [K]"]
         # delta_phi = phi_s - phi_e
         if self.reaction_loc == "interface":
@@ -256,7 +258,8 @@ def set_rhs(self, variables):
         )
 
         # we have to add the spreading rate to account for cracking
-        SEI_option = getattr(self.options, domain)["SEI"]
+        domain_options = getattr(self.options, domain)
+        SEI_option = getattr(domain_options, self.phase)["SEI"]
         if SEI_option.startswith("ec reaction limited"):
             self.rhs = {L_outer: -dLdt_SEI_outer + spreading_outer}
         else:
@@ -284,7 +287,8 @@ def set_initial_conditions(self, variables):
         else:
             L_inner_0 = self.phase_param.L_inner_0
             L_outer_0 = self.phase_param.L_outer_0
-        SEI_option = getattr(self.options, domain)["SEI"]
+        domain_options = getattr(self.options, domain)
+        SEI_option = getattr(domain_options, self.phase)["SEI"]
         if SEI_option.startswith("ec reaction limited"):
             self.initial_conditions = {L_outer: L_inner_0 + L_outer_0}
         else: