GAUZZ (pronounced "Gauss") : GPU-Accelerated nested UZawa solver for Stokes flow in geodynamics. The name combines the solver's two algorithmic pillars (GPU and Uzawa iteration) Source code and benchmark input scripts for the GPU-accelerated finite-element Stokes solver presented in Lee et al. (2026), submitted to Geoscientific Model Development. The code couples FEniCS with PETSc's CUDA backend for solving incompressible Stokes equations in geodynamics. Benchmarks include a manufactured solution, the SolCx viscosity-discontinuity benchmark, sticky-air free-surface tests, 2D/3D Rayleigh-Taylor instabilities, 3D thermo-mechanical convection, and a 2D subduction model with nonlinear viscosity.