docs: def merge_limits(axes, xlim=True, ylim=True):
      '''Set maximum and minimum limits from list of axis objects to each axis
  
      Args
      ----
      axes: iterable
          list of `matplotlib.pyplot` axis objects whose limits should be modified
      xlim: bool
          Flag to set modification of x axis limits
      ylim: bool
          Flag to set modification of y axis limits
      '''
  
      # Compile lists of all x/y limits
      xlims = list()
      ylims = list()
      for ax in axes:
          [xlims.append(lim) for lim in ax.get_xlim()]
          [ylims.append(lim) for lim in ax.get_ylim()]
  
      # Iterate over axes objects and set limits
      for ax in axes:
          if xlim:
              ax.set_xlim(min(xlims), max(xlims))
          if ylim:
              ax.set_ylim(min(ylims), max(ylims))
  
      return None