def METHOD_NAME(self, field='velocity', index=0, times=[0., 100]):
    '''
    Provide field time average
    field - the field to time-average
    index - the component index (for example: velocity has 3 components)
    times - the times to average
    '''
    # Filter the field based on the times
    filt = ((self.time[:] >= times[0]) & (self.time[:] <= times[1]))
    # Filtered time
    t = self.time[filt]
    # The total time
    dt = np.amax(t) - np.amin(t)
    # Filtered field
    f = self.abl_stats[field][filt,:,index]
    # Compute the time average as an integral
    integral = np.trapz(f, x=t, axis=0) / dt
    return integral