diff --git a/examples/symba_gr_test/init_cond.py b/examples/symba_gr_test/init_cond.py
index e59294190..a88ed26fb 100755
--- a/examples/symba_gr_test/init_cond.py
+++ b/examples/symba_gr_test/init_cond.py
@@ -61,7 +61,7 @@
 GMcb = sim.ds.isel(id=0)['Gmass'].values
 GU = swiftest.GC / (sim.param['DU2M']**3 / (sim.param['MU2KG'] * sim.param['TU2S']**2))
 dens = 3000.0 / (sim.param['MU2KG'] / sim.param['DU2M']**3) # Assume a bulk density of 3 g/cm^3
-GM_pl = 1e-11
+GM_pl = sim.ds.sel(id=1)['Gmass'].values * 0.01
 M_pl = GM_pl / GU
 R_pl = (3 * M_pl / (4 * np.pi * dens))**(1.0 / 3.0)
 Rh_pl = Me_a * (GM_pl / (3 * GMcb))**(1.0/3.0)
diff --git a/examples/symba_gr_test/pl.swiftest.in b/examples/symba_gr_test/pl.swiftest.in
index cd36bda61..1e4f398e9 100644
--- a/examples/symba_gr_test/pl.swiftest.in
+++ b/examples/symba_gr_test/pl.swiftest.in
@@ -31,7 +31,7 @@ Neptune 0.0020336100526728302319 0.7749718408665498732
 0.000164587904124493665
 30.038959912561129073 0.008955570159296157365 1.7711193542961420899
 131.82211597488270627 284.47484279411258967 308.45137222693909962
-Planetesimal 9.999999999999999395e-12 1.698243864025464137e-05
-2.2876635862715332271e-07
+Planetesimal 6.553709809565314101e-08 0.00031780626665496970825
+4.281092009908029071e-06
 0.38709858430953941744 0.20562340108970009189 7.0033025080013837638
 48.296118373786072198 29.20442403952453958 339.33948746828792764