Full zonal mean version with test for volume