From 7ef443b1de470f830c3610ec4cee6c30277593a1 Mon Sep 17 00:00:00 2001 From: Federico Esteban Date: Thu, 11 May 2023 08:56:23 -0300 Subject: [PATCH] Update gravity_disturbance_topofree.py At the beginnig of the scripts says: "The gravity disturbance is the observed gravity minus the normal gravity (boule.Ellipsoid.normal_gravity). The signal that is left is assumed to be due to the differences between the actual Earth and the reference ellipsoid." So, the value 0 has a special meaning. But, in the figure the hinge (white color) of the CPT is around 100. I think that it should be in 0. To fix it, I propose to force the hinge at 0 (+h0). --- examples/gravity_disturbance_topofree.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/examples/gravity_disturbance_topofree.py b/examples/gravity_disturbance_topofree.py index 737fbf6f7..888136bea 100644 --- a/examples/gravity_disturbance_topofree.py +++ b/examples/gravity_disturbance_topofree.py @@ -61,7 +61,7 @@ # Make a plot of data using PyGMT fig = pygmt.Figure() -pygmt.grd2cpt(grid=disturbance_topofree, cmap="vik", continuous=True) +pygmt.grd2cpt(grid=disturbance_topofree, cmap="vik+h0", continuous=True) title = "Topography-free (Bouguer) gravity disturbance of the Earth"