docs: def plotcapzone(self, nt=10, zstart=None, hstepmax=20, vstepfrac=0.2,
                      tmax=365, nstepmax=100, silent='.', color=None, orientation='hor',
                     win=[-1e30, 1e30, -1e30, 1e30], newfig=False, figsize=None):
          """Plot a capture zone
          
          Parameters
          ----------
          nt : int
              number of path lines
          zstart : scalar
              starting elevation of the path lines
          hstepmax : scalar
              maximum step in horizontal space
          vstepfrac : float
              maximum fraction of aquifer layer thickness during one step
          tmax : scalar
              maximum time
          nstepmax : scalar(int)
              maximum number of steps
          silent : boolean or string
              True (no messages), False (all messages), or '.'
              (print dot for each path line)
          color : color
          orientation : string
              'hor' for horizontal, 'ver' for vertical, or 'both' for both
          win : array_like (length 4)
              [xmin, xmax, ymin, ymax]
          newfig : boolean (default False)
              boolean indicating if new figure should be created
          figsize : tuple of integers, optional, default: None
              width, height in inches.
              
          """
          
          xstart, ystart, zstart = self.capzonestart(nt, zstart)
          self.model.tracelines(xstart, ystart, zstart, hstepmax=-abs(hstepmax), vstepfrac=vstepfrac,
                     tmax=tmax, nstepmax=nstepmax, silent=silent, color=color, orientation=orientation,
                     win=win, newfig=newfig, figsize=figsize)