Kinematic data of NGC1023. Top: fitted kinematic parameters in elliptical bins as a function of major axis distance a for NGC1023 (see Eq. (1)). Results on the Atlas3D velocity field are shown in red, and those on the ePN.S smoothed velocity field are in green. Blue symbols show rotational velocities and kinematic position angles from kinemetry on SLUGGS data (Foster et al. 2016). These agree within errors with the V_rot and PA_kin values measured from Atlas3D. ePN.S values also agree within errors except for the first PN velocity point at 70 arcsec, where the PN systematically underestimate the mean velocity (Pulsoni et al. 2018, Sect. 7.2). In the radial range marked with vertical lines, we therefore use the V_rot(a) and PA_kin(a) profiles from SLUGSS, and V_s3(a) and V_c3(a) from the linear interpolation of the values fitted on the Atlas3D and ePN.S velocity fields. Bottom: the two-dimensional rotation velocity field of NGC1023 reconstructed from Atlas3D, SLUGGS, and ePN.S., using Eq. (1). The velocity field is rotated so that the photometric major axis (PA_phot = 83.3 deg) is horizontal. The final mean velocity field is given by the Atlas3D data in the center (within the black square), the reconstructed two-dimensional field from SLUGGS at intermediate radii (within the black ellipse), and the ePN.S data at large radii. The PN positions are highlighted by open black circles.