You need ax_gs = Axis(f[1, 1], spinewidth=2.0pt, aspect=DataAspect())
ax_gs = Axis(f[1, 1], spinewidth=2.0pt, aspect=DataAspect())