Given a partition $\lambda\vdash n$ of $n$, look at its Young diagram $Y_{\lambda}$. Let $a(\lambda)$ be the number of squares (of all sizes) in $Y_{\lambda}$. For example, if $n=4$ then $a(4)=4, a(3,1)=4, a(2,2)=5, a(2,1,1)=4$, and $a(1,1,1,1)=4$. So, we have a total of $4+4+5+4+4=21$ squares.
Together with Richard Stanley (unpublished, Section 9), we proved that $$\sum_{n=0}^{\infty}F_n(q)\, x^n=\sum_{k=0}^{\infty}\frac{x^{k^2}q^{\frac{k(k+1)(2k+1)}6}} {\prod_{j=1}^k\left(1-x^j\,q^{\frac{j(j+1)}2}\right)^2}, \qquad \text{where} \qquad F_n(q)=\sum_{\lambda\vdash n}q^{a(\lambda)}.$$
QUESTION. Now, instead consider $b(\lambda)$ to be the number of rectangles in $Y_{\lambda}$. Can you provide a generating function for $G_n(q)=\sum_{\lambda\vdash n}q^{b(\lambda)}$?