Consider a finite field $\mathbb{F}_p$ such that $p \equiv 1 \ (\mathrm{mod} \ 3)$ and its element $\zeta \neq 1$, $\zeta^3 = 1$.
Also, let $E\!: y^2 = x^3 + b$ be an elliptic curve of $j$-invariant $0$, where $b \in \mathbb{F}_p^* \setminus (\mathbb{F}_p^*)^3.$ This curve has the order $3$ automorphism $$[\zeta]\!: (x,y) \mapsto (\zeta x, y).$$
Consider the diagonal matrices $$A := \mathrm{diag}(1, \zeta, \zeta^2),\qquad B := \mathrm{diag}(\zeta^2, \zeta, 1).$$ They generate the subgroup $$G := \langle A, B\rangle \subset \mathrm{SL}(3, \mathbb{Z}[\zeta]),$$ which is isomorphic to $(\mathbb{Z}/3)^2$. Any element $\mathrm{diag}(\alpha,\beta,\gamma) \in G$ naturally acts on the abelian threefold $E^3$: $$ (P, Q, R) \mapsto ([\alpha]P, [\beta]Q, [\gamma]R). $$
What is the geometric quotient $E^3\!/G$? Could you explicitly write a (affine) model for this variety?