Numerical calculations using the full MHD equations in toroidal geometry