We present an analysis of the anisotropic flow harmonics in Pb+Pb collisions at beam momenta of 30A GeV/c collected by the NA61/SHINE experiment in the year 2016. Directed and elliptic flow coefficients are measured relative to the spectator plane estimated with the Projectile Spectators Detector (PSD). The flow coefficients are reported as a function of transverse momentum in different classes of collision centrality. The results are compared with a new analysis of the NA49 data for Pb+Pb collisions at 40A GeV using forward calorimeters (VCal and RCal) for event plane estimation.