How to plot the geomagnetic field intensity at the Earth's surface use pygmt?

Do you have datasets to plot? The steps should be similar to plotting earth relief datasets, which is demonstrated in this tutorial - https://www.pygmt.org/latest/tutorials/advanced/earth_relief.html#sphx-glr-tutorials-advanced-earth-relief-py. To plot multiple datasets, you can also refer to the ‘making subplots’ tutorial - https://www.pygmt.org/latest/tutorials/advanced/subplots.html.

Thanks for your reply :grinning: