Zonal harmonic model of Saturn's magnetic field from Voyager 1 and 2 observations