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

[LiHai]

Thanks for your reply