Models of electrostatic surfaces in atomic crystals rely on equations involving the Jacobi theta functions. Numerical integration of these is prohibitively time consuming, making it difficult to examine the properties of the fields which give rise to the surfaces. We give simple expressions for the key electrostatic surfaces using Fourier expansions in basis sets of nodal surfaces. Any surface may be computed in seconds in a form ammenable to further analysis. The distribution of the mean and Gaussian curvatures over each surface has been visualised by assigning colours so that the range from minimum to maximum value spans blue to red. We similarly explore the mean and Gaussian scalar fields over a range of triply periodic surfaces of the same morphology.