Reproducible B0 shim matrix calibration at 7T and dynamic range analysis of spherical harmonics for residual compensation