Plotting Impact Map Returns Array Instead of Axes

[DanielHerzbergMeteoSwiss]

When calling plot_imp_map function of climada.engine.forecast module, it does not return a GeoAxes object as documented (and expected) but an array containing that object.

This happens in _plot_imp_map function with the following code:

    if not isinstance(axis_sub, np.ndarray):
        axis_sub = np.array([[axis_sub]])

The return value of the function should a GeoAxes object analog to the other plotting functions.

[chahank]

Thanks for noting! Would it be possible for you to make a pull request with a bugfix?

[peanutfun]

@DanielHerzbergMeteoSwiss Thanks for the report! I think the idea of the code was to return what plt.subplots returns. But the implementation does not really hit the mark. plt.subplots returns a squeezed numpy array of the plots, e.g. an array [[ax1, ax2], [ax3, ax4]] for two rows and columns, and array [ax1, ax2] for one row, two columns, and ax (no numpy array) for a single row and column. A suggestion for improvement would be very welcome!