In this paper, we develop numerical methods based on the weighted Birkhoff average for studying two-dimensional invariant tori for volume-preserving maps. The methods do not rely on symmetries, such as time-reversal symmetry, nor on approximating tori by periodic orbits. The rate of convergence of the average gives a sharp distinction between chaotic and regular dynamics and allows accurate computation of rotation vectors for regular orbits. Resonant and rotational tori are distinguished by computing the resonance order of the rotation vector to a given precision. Critical parameter values, where tori are destroyed, are computed by a sharp decrease in convergence rate of the Birkhoff average. We apply these methods for a three-dimensional generalization of Chirikov’s standard map: an angle–action map with two angle variables. Computations on grids in frequency and perturbation amplitude allow estimates of the critical set. We also use continuation to follow tori with fixed rotation vectors. We test three conjectures for cubic fields that have been proposed to give locally robust invariant tori, but are not able to provide compelling evidence that one of these three fields is more robust than the other two.