Hi,
You may use FDE, which is specfically developed to calculate the modes of a waveguide. Also, the mode source in FDTD has an inbuilt eignemode solver, which can also calculate the modes. But I would recommend FDE if you have the option.
Regards,
Amrita