Comments from code review in #20:
I wouldn't have simulate() return results at all, since that encourages users to separate the results from the originating simulation. I would encourage the use of...
>>> sim.simulate()
>>> sim.results.plot_something()
Or, if the result was a xarray.Dataset then one could do...
>>> sim.simulate()
>>> sim.results.T_e.plot()