The Fermilab Muon g-2 Experiment (E989) contains flat-plate electrostatic quadrupoles, curved with the reference trajectory as defined by the constant, uniform magnetic dipole field. To understand the beam behaviour at a sufficient level, we require fast, high-accuracy particle tracking methods for this layout. Standard multipole fits to numerically calculated 2D transverse electric field maps have provided a first approximation to the electric field within the main part of the quadrupole, but cannot model the longitudinal curvature or extended fringe fields of the electrostatic plates. Expressions for curvilinear multipoles can be fit to a 2D transverse slice taken from the central point of a numerically calculated 3D electric field map of the quadrupole, providing a curved-multipole description. Generalised gradients can be used to model the fringe field regions. We present the results of curvilinear multipole and generalised gradient fits to the curved quadrupole fields, and the differences in tracking using these fields over 200 turns of a model of the storage ring in BMAD.