We continue our research on Fourier restriction for hyperbolic surfaces, by studying local perturbations of the hyperbolic paraboloid z = x y which are of the form z = x y + h(y), where h(y) is a smooth function which is flat at the origin. The case of perturbations of finite type had already been handled before, but the flat case imposes several new obstacles. By means of a decomposition into intervals on which |h | is of a fixed size λ, we can apply methods devised in preceding papers, but since we lose control on higher order derivatives of h we are forced to rework the bilinear method for wave packets that are only slowly decaying. Another problem lies in the passage from bilinear estimates to linear estimates, for which we need to require some monotonicity of h .