Fast wave (FW) studies of mode conversion (MC) processes at the ion-ion hybrid layer in toroidal plasmas must capture the disparate scales of the FW and mode converted ion Bernstein and ion cyclotron waves. Correct modeling of the MC layer requires resolving wavelengths on the order of k(perpendicular to)rho(i)similar to1 which leads to a scaling of the maximum poloidal mode number, M-max, proportional to 1/rho(*) (rho(*)equivalent torho(i)/L). The computational resources needed scale with the number of radial (N-r), poloidal (N-theta), and toroidal (N-phi) elements as N-r * N-phi * N-theta(3). Two full wave codes, a massively-parallel-processor (MPP) version of the TORIC-2D finite Larmor radius code [M. Brambilla, Plasma Phys. Controlled Fusion 41, 1 (1999)] and also an all orders spectral code AORSA2D [E. F. Jaeger , Phys. Plasmas 9, 1873 (2002)], have been developed which for the first time are capable of achieving the resolution and speed necessary to address mode conversion phenomena in full two-dimensional (2-D) toroidal geometry. These codes have been used in conjunction with theory and experimental data from the Alcator C-Mod [I. H. Hutchinson , Phys. Plasmas 1, 1511 (1994)] to gain new understanding into the nature of FWMC in tokamaks. The massively-parallel-processor version of TORIC is also now capable of running with sufficient resolution to model planned lower hybrid range of frequencies experiments in the Alcator C-Mod. (C) 2004 American Institute of Physics.