code: def truncate_array(array1, array2, position):
      
  
      slices = []
      for i in range(array1.ndim):
          xmin = 0
          xmax = array1.shape[i]
          dxlo = array1.shape[i] // 2
          dxhi = array1.shape[i] - dxlo
          if position[i] - dxlo < 0:
              xmin = max(dxlo - position[i], 0)
  
          if position[i] + dxhi > array2.shape[i]:
              xmax = array1.shape[i] - (position[i] + dxhi - array2.shape[i])
              xmax = max(xmax, 0)
          slices += [slice(xmin, xmax)]
  
      return array1[slices]
docs: Truncate array1 by finding the overlap with array2 when the
      array1 center is located at the given position in array2.