The study of plasma and conductive fluid motions in presence of electro-magnetic fields is based on the magnetohydrodynamic balance equations. Depending on several parameters such as the kinetic and magnetic Reynolds numbers, as well as the magnetic Prandtl number, a very wide class of phenomena can be described by these equations. For instance, almost the same formalism may be used to investigate fairly different systems such as astrophysical media, planetary cores, fusion plasmas and liquid metal flows. This has motivated the development of a unified numerical tool capable of describing these phenomena, at least in the limit of the available computational power. We briefly describe the strategy adopted in developing this code as well as various examples in which it has been exploited.