A multipole reaction-field model for gauge-origin independent magnetic properties of solvated molecules