Ab initio calculation of magnetic susceptibilities and screening constants using gauge invariant Gaussian orbitals