GRMHD simulations of accreting neutron stars with non-dipole fields