Simulating horizontal crustal motions of glacial isostatic adjustment using compressible cartesian models