code: def test_multi_panel_single_band(one_band_3dims):
      """Test that multi panel works with single band arr."""
  
      title1 = "Title axis one"
      title2 = "Title axis two"
      f, (ax1, ax2) = plt.subplots(2, 1)
      ep.plot_bands(one_band_3dims, title=title1, ax=ax1)
      ep.plot_bands(one_band_3dims, title=title2, ax=ax2)
  
      # get all axis subplot elements - note this returns subplots and axes
      all_axes = f.axes
  
      assert len(all_axes) == 4
      assert all_axes[0].get_title() == title1
      assert all_axes[1].get_title() == title2
  
imports: import pytest
  import matplotlib.pyplot
  import matplotlib.colors
  import earthpy.plot
  import numpy
label_desc: Use image node_id 16 for calling a global function with example usage: test_multi_panel_single_band(one_band_3dims) without return types