Runge Kutta methods are well known high order methods for ODEs. In scalar advection problems with a single family of characteristics, any high order Runge-Kutta method can be used to compute high order solutions that don't diffuse by following characteristics. This is known as the method of lines. When the advection equation is no longer scalar, two (or more) families of characteristics may appear. By putting the unknowns at the intersection between the characteristics of these two different families, designing a Two-Dimensional Runge-Kutta method of order $2$ is hardly more difficult. But going beyond order $2$, designing Two-Dimensional Runge-Kutta methods of order $3$ or above is far more difficult.