Finite element derivative superconvergent points for harmonic functions under local rectangular mesh are investigated. All superconvergent points for the finite element space of any order that is contained in the tensorproduct space and contains the intermediate family can be predicted. In the case of the serendipity family, results are given for finite element spaces of order below 6. The results justify the computer findings of Babuska, et al.