docs: def plot(self, win=None, newfig=True, figsize=None, orientation='hor', topfigfrac=0.8):
          """Plot layout
          
          Parameters
          ----------
      
          win : list or tuple
              [x1, x2, y1, y2]
              
          """
          
          if newfig:
              plt.figure(figsize=figsize)
              ax1 = None
              ax2 = None
              if orientation == 'both':
                  ax1 = plt.axes([0.125, 0.18 + (1 - topfigfrac) * 0.7, (0.9 - 0.125), topfigfrac * 0.7])
                  ax2 = plt.axes([0.125, 0.11, (0.9 - 0.125), (1 - topfigfrac) * 0.7], sharex=ax1)
              elif orientation[:3] == 'hor':
                  ax1 = plt.subplot()
              elif orientation[:3] == 'ver':
                  ax2 = plt.subplot()
          else:
              if orientation == 'both':
                  fig = plt.gcf()
                  ax1 = fig.axes[0]
                  ax2 = fig.axes[1]
              elif orientation[:3] == 'hor':
                  fig = plt.gcf()
                  ax1 = fig.axes[0]
                  ax2 = None
              elif orientation[:3] == 'ver':
                  fig = plt.gcf()
                  ax1 = None
                  ax2 = fig.axes[0]
          if ax1 is not None:
              plt.sca(ax1)
              for e in self.elementlist:
                  e.plot()
              if orientation[:3] == 'hor':
                  plt.axis('scaled')
              elif orientation == 'both':
                  plt.axis('equal')  # cannot be 'scaled' when sharing axes
              if win is not None:
                  plt.axis(win)
          if ax2 is not None:
              plt.sca(ax2)
              for i in range(self.aq.nlayers):
                  if self.aq.ltype[i] == 'l':
                      plt.axhspan(ymin=self.aq.z[i + 1], ymax=self.aq.z[i], color=[0.8, 0.8, 0.8])
              for i in range(1, self.aq.nlayers):
                  if self.aq.ltype[i] == 'a' and self.aq.ltype[i - 1] == 'a':
                      plt.axhspan(ymin=self.aq.z[i], ymax=self.aq.z[i], color=[0.8, 0.8, 0.8])