plots to compare fits

We should develop some plotting tools to facilitate model comparisons. A good starting point would be allowing each panel of the annual plots to include multiple fits. `plot_biomass` and `plot_rec` can do this using facet_wrap so just need to extend this capability to single panels. Would be handy for SIs too