We present a numerical version of the moment method for the identification of non-radial pulsation modes. The new version requires less computation time than the previous one, allowing users to consider all the information contained in the first three moments of a multiperiodic star and to identify multiple modes simultaneously. This, together with the use of a new discriminant that considers the moments calculated at each time of observation, increases considerably the feasibility and the accuracy of the mode identification for multiperiodic stars. Moreover, the technique is extended to rotating pulsating stars. We apply the new version to three B stars showing multiperiodicity.