This paper presents nonlinear gyrokinetic simulations of electron temperature gradient (ETG) turbulence based on spherical tokamak (ST) parameters. Most significantly the simulations include the strong toroidal flow and flow shear present in STs that suppress ion-scale turbulence while using kinetic ions at full mass ratio (mi/me=3600). The flow shear provides a physical long-wavelength cutoff mechanism that aids saturation of the simulations, which has previously been demonstrated to be problematic depending on magnetic shear. As magnetic shear varies widely in STs we systematically demonstrate saturation and convergence of the ETG simulations with respect to grid resolution, physical domain size, and boundary conditions. While using reduced ion mass or adiabatic ions can lessen computational expense they do not always provide reliable results. The resulting spectra from converged simulations are anisotropic everywhere in contrast to previous ETG simulations without flow shear. These results have implications for interpreting turbulence measurements, and represent an important step in determining when and where ETG turbulence is expected to be relevant in ST plasmas. They are also important in the context of validating simulations with both experimental transport analysis and turbulence measurements.