Implement Rusanov MHD solver