Raise error when basic half cell model is used with simulations

pybamm/simulation.py

@@ -97,6 +97,10 @@ def __init__(
     ):
         self.parameter_values = parameter_values or model.default_parameter_values
 
+        if isinstance(model, pybamm.lithium_ion.BasicDFNHalfCell):
+            if experiment is not None:
+                raise NotImplementedError("BasicDFNHalfCell is not compatible with experiment simulations yet.")
+
         if experiment is None:
             # Check to see if the current is provided as data (i.e. drive cycle)
             current = self._parameter_values.get("Current function [A]")

tests/unit/test_models/test_full_battery_models/test_lithium_ion/test_basic_models.py

@@ -36,6 +36,19 @@ def test_dfn_half_cell_well_posed(self):
         copy = model.new_copy()
         copy.check_well_posedness()
 
+    def test_dfn_half_cell_simulation_with_experiment_error(self):
+        options = {"working electrode": "negative"}
+        model = pybamm.lithium_ion.BasicDFNHalfCell(options=options)
+        experiment = pybamm.Experiment([("Discharge at C/10 for 10 hours or until 3.5 V")])
+        with self.assertRaisesRegex(
+            NotImplementedError, "BasicDFNHalfCell is not compatible with experiment simulations yet."):
+            pybamm.Simulation(model, experiment=experiment)
+
+    def basic_dfn_half_cell_simulation(self):
+        model = pybamm.lithium_ion.BasicDFNHalfCell(options={"working electrode": "positive"})
+        sim = pybamm.Simulation(model=model)
+        sim.solve([0, 3600])
+
     def test_dfn_half_cell_defaults(self):
         # test default geometry
         var = half_cell_spatial_vars