Download the seismological Moho depth grid from the supplementary material of the paper and convert it to netCDF for easier manipulation and loading with GMT and xarray.
Szwillus, W., Afonso, J. C. C., Ebbing, J., & Mooney, W. D. (2019). Global crustal thickness and velocity structure from geostatistical analysis of seismic data. Journal of Geophysical Research: Solid Earth, 124, 1626– 1652. https://doi.org/10.1029/2018JB016593
All source code is made available under a MIT license. You can freely use and modify the code, without warranty, so long as you provide attribution to the original authors. See LICENSE for the full license text.