In Abaqus the default output points through the thickness of a shell section are the points that are on the bottom and top surfaces of the shell section (for integration with Simpson's rule) or the points that are closest to the bottom and top surfaces (for Gauss quadrature). For example, if five integration points are used through a single layer shell, output will be provided for section points 1 (bottom) and 5 (top). So, how to get the through thickness distribution?