We suggest a simple and an efficient way of selecting a suitable set of interpolation points for the well-known rational Krylov based model order reduction techniques. To do this, some sampling points from the frequency response of the transfer function are taken. These points correspond to the places where the sign of the numerical derivation of transfer function changes. The suggested method requires a set of linear system's solutions several times. But, they can be computed concurrently by different processors in a parallel computing environment. Serial performance of the method is compared to the well-known H-2 optimal method for several benchmark examples. The method achieves acceptable accuracies (the same order of magnitude) compared to that of H-2 optimal methods and has a better performance than the common selection procedures such as linearly distributed points.