Adapt `plotYieldObservedVsModel()` and `matchYield()` to work with multiple gears

[gustavdelius]

These functions currently assume without checking that a species is selected by a single gear only and that the observed yield is stored in the species_params data frame rather than the gear_params data frame.