“…Thus, equations (1) and (2) turn into scalar differential equations with complex exponential functions as solutions, representing the upward and downward traveling waves. In the spectrum domain, to ensure numerical stability, we adopt the idea of Das (1995) and write the solutions as A are the spectral counterparts of A z,l and F z,l of the space domain (two-dimensional Fourier transformation). A l + and A l -are the amplitudes of the downgoing and upgoing TM mode wave potentials, while F l + and F l -are the amplitudes for the TE mode.…”