Parabolic approximation method for fast magnetosonic wave propagation in tokamaks