Homoclinic tangency of codimension 2 and historic wandering domains

We construct a four-dimensional diffeomorphism exhibiting a homoclinic tangency of the largest codimension, which admits a historic wandering domain of positive Lebesgue measure. Every orbit in this wandering domain exhibits historic behavior, in the sense that time averages do not converge. This example shows that homoclinic tangencies of the largest codimension can still give rise to positive Lebesgue measure sets with non-convergent statistical behavior.