The analysis of electromagnetic (EM) wave propagation in an inhomogeneous, magnetically confined plasma is essential for nuclear fusion, particularly for ion and electron heating. This paper presents EMWAVE, a finite element code that simulates EM wave propagation in an anisotropic plasma by solving the Helmholtz equation in the frequency domain. The solver supports linear and quadratic elements and efficiently models complex domains. Various dielectric media and tokamak relevant configurations have been benchmarked, demonstrating the code’s reliability. EMWAVE incorporates the cold plasma permittivity tensor, showing excellent agreement with the well established ERMES code and theoretical predictions for tokamak cutoffs. Developed as a standalone tool and integrated into the HPC framework Alya, it aims to support reactor-scale simulations via a digital twin approach. Future work includes implementing a hot plasma permittivity tensor to capture wave particle interactions and extending the model to 3D reactor geometries.