docs: def plot_dives(dv0, dv1, p, dp, t_on, t_off):
      '''Plots depths and delta depths with dive start stop markers
  
      Args
      ----
      dv0: int
          Index position of dive start in cue array
      dv1: int
          Index position of dive stop in cue array
      p: ndarray
          Depth values
      dp: ndarray
          Delta depths
      t_on: ndarray
          Cue array with start index position of dives
      t_off: ndarray
          Cue array with stop index postition of dives
      '''
  
      fig, (ax1, ax2) = plt.subplots(2, 1, sharex=True)
  
      x0   = t_on[dv0:dv1] - t_on[dv0]
      x1   = t_off[dv0:dv1] - t_on[dv0]
  
      # Extract start end depths
      y0_p = p[t_on[dv0:dv1]]
      y1_p = p[t_off[dv0:dv1]]
  
      # Extract start end delta depths
      y0_dp = dp[t_on[dv0:dv1]]
      y1_dp = dp[t_off[dv0:dv1]]
  
      start = t_on[dv0]
      stop  = t_off[dv1]
  
      ax1.title.set_text('Dives depths')
      ax1.plot(range(len(p[start:stop])), p[start:stop])
      ax1.scatter(x0, y0_p, label='start')
      ax1.scatter(x1, y1_p, label='stop')
      ax1.set_ylabel('depth (m)')
  
      ax1.title.set_text('Depth rate of change')
      ax2.plot(range(len(dp[start:stop])), dp[start:stop])
      ax2.scatter(x0, y0_dp, label='start')
      ax2.scatter(x1, y1_dp, label='stop')
      ax2.set_ylabel('depth (dm/t)')
      ax2.set_xlabel('sample')
  
      for ax in [ax1, ax2]:
          ax.legend(loc='upper right')
          ax.set_xlim([-50, len(dp[start:stop])+50])
  
      plt.show()
  
      return None