diff --git a/doc/examples/6_inversion/plot_3_petro_joint_inv.py b/doc/examples/6_inversion/plot_3_petro_joint_inv.py
index f147cf136..0bd1bf470 100644
--- a/doc/examples/6_inversion/plot_3_petro_joint_inv.py
+++ b/doc/examples/6_inversion/plot_3_petro_joint_inv.py
@@ -62,7 +62,7 @@
 saturation = pg.solver.parseArgToArray([[1, 0.6], [2, 0.9], [3, 0.3]],
                                        mMesh.cellCount(), mMesh)
 
-satKW = dict(cMin=0.3, cMax=0.9, logScale=False, cMap="plasma")
+satKW = dict(cMin=0.3, cMax=0.9, logScale=False, cMap="plasma_r")
 ax, _ = pg.show(mMesh, saturation, **satKW,
              showMesh=True, label=r'Saturation (${\tt petro}$)')