docs: def get_time_slice(time, z, zdot=None, timeStart=None, timeEnd=None):
      """
      Get slice of time, z and (if provided) zdot from timeStart to timeEnd.
  
      Parameters
      ----------
      time : ndarray
          array of time values 
      z : ndarray
          array of z values
      zdot : ndarray, optional
          array of zdot (velocity) values.
      timeStart : float, optional
          time at which to start the slice.
          Defaults to beginnging of time trace
      timeEnd : float, optional
          time at which to end the slide.
          Defaults to end of time trace
  
      Returns
      -------
      time_sliced : ndarray
          array of time values from timeStart to timeEnd
      z_sliced : ndarray
          array of z values from timeStart to timeEnd
      zdot_sliced : ndarray
          array of zdot values from timeStart to timeEnd.
          None if zdot not provided
  
      """
      if timeStart == None:
          timeStart = time[0]
      if timeEnd == None:
          timeEnd = time[-1]
  
      StartIndex = _np.where(time == take_closest(time, timeStart))[0][0]
      EndIndex = _np.where(time == take_closest(time, timeEnd))[0][0]
  
      time_sliced = time[StartIndex:EndIndex]
      z_sliced = z[StartIndex:EndIndex]
  
      if zdot != None:
          zdot_sliced = zdot[StartIndex:EndIndex]
      else:
          zdot_sliced = None    
      
      return time_sliced, z_sliced, zdot_sliced