diff --git a/examples/telemac3d/particles/doc/particles.tex b/examples/telemac3d/particles/doc/particles.tex
index c31430d38b05d216dfe96a5a45b411e882683965..886f55c05fafb14165bd6353a5e2c5af875a01cd 100644
--- a/examples/telemac3d/particles/doc/particles.tex
+++ b/examples/telemac3d/particles/doc/particles.tex
@@ -73,9 +73,18 @@ The time step is 5~s for a simulated period of  2~h (= 7,200~s).
 
 The non-hydrostatic version is used.
 
-To solve the advection, the method of characteristics
-is used for the velocities (scheme 1).
-The GMRES is used for solving the propagation and PPE steps (option 7).
+To solve the advection of velocities, the default choice (LIPS) is used rather
+than the method of characteristics (scheme 1) so that the differences
+between sequential and parallel runs are decreased.
+NERD scheme (\# 13) can also be used for that purpose.
+The GMRES is used for solving the propagation (option 7 = default choice)
+whereas the conjugate gradient (option 1) combined with preconditioning = 34
+= 2 $times$ 17 is used for solving PPE as more efficient than GMRES for this
+example.
+Moreover, preconditioning for velocities is also set to 34 (a little bit
+more efficient).
+With this setting, accuracies for propagation, PPE and diffusion of velocities
+can be let to default values (10$^{-8}$).
 
 A maximum of 100 drogues are released.
 %and control sections are calculated between two couples of nodes.
diff --git a/examples/telemac3d/particles/f3d_particles.slf b/examples/telemac3d/particles/f3d_particles.slf
index e15df02db307a54452cfbdf798a5d4c43aa6e86e..303ce58f9a181fa749b122422cbdbfd0257c46e2 100644
--- a/examples/telemac3d/particles/f3d_particles.slf
+++ b/examples/telemac3d/particles/f3d_particles.slf
@@ -1,3 +1,3 @@
 version https://git-lfs.github.com/spec/v1
-oid sha256:23a2e3103c5d0c45429464dd426dfb7327ceaa815f9cb82c165cdc62c85ed839
-size 1714096
+oid sha256:92a93b00054273dfe02fef8280dff1d45733cc4c42ebd27dd0c3916f4b3371d2
+size 1387844
diff --git a/examples/telemac3d/particles/t3d_particles.cas b/examples/telemac3d/particles/t3d_particles.cas
index 3c9b2383341e8e5c6d894a4fefddfb26a61551da..82bfe761b8045dfc2dde89c48643aa00dd6591dc 100644
--- a/examples/telemac3d/particles/t3d_particles.cas
+++ b/examples/telemac3d/particles/t3d_particles.cas
@@ -36,7 +36,6 @@ HORIZONTAL TURBULENCE MODEL = 1
 VERTICAL TURBULENCE MODEL = 2
 MIXING LENGTH MODEL       = 1   / OLD DEFAULT VALUE UNTIL V8P5
 /
-IMPLICITATION FOR DEPTH = 1.
 PRESCRIBED ELEVATIONS   :  265.;0.
 PRESCRIBED FLOWRATES    :  0.;700.
 /
@@ -48,12 +47,10 @@ MASS-LUMPING FOR DIFFUSION : 1.
 /
 /  SYSTEMES LINEAIRES DE L'OPTION NON-HYDROSTATIQUE
 /
-ACCURACY FOR PPE = 1.E-5
+SOLVER FOR PPE = 1
+PRECONDITIONING FOR PPE = 34
+/
+PRECONDITIONING FOR DIFFUSION OF VELOCITIES = 34
 /
 MAXIMUM NUMBER OF DROGUES : 100
 PRINTOUT PERIOD FOR DROGUES :  36
-/
-/ DEFAULT VALUES UNTIL V8P0 KEPT FOR NON REGRESSION
-SCHEME FOR ADVECTION OF VELOCITIES = 1
-IMPLICITATION FOR VELOCITIES = 1.
-ACCURACY FOR PROPAGATION = 1.E-6
diff --git a/examples/telemac3d/particles/vnv_particles.py b/examples/telemac3d/particles/vnv_particles.py
index fc8f75006277190bdb1bc3665cdfed425354f1a8..6b0df4402fb3a57100b46cbad8da3e7de3c0e640 100644
--- a/examples/telemac3d/particles/vnv_particles.py
+++ b/examples/telemac3d/particles/vnv_particles.py
@@ -49,17 +49,17 @@ class VnvStudy(AbstractVnvStudy):
         # Comparison with the last time frame of the reference file.
         self.check_epsilons('vnv_1:T3DRES',
                             'f3d_particles.slf',
-                            eps=[7e-5, 0.013, 7e-3, 3e-3])
+                            eps=[1e-8])
 
         # Comparison with the last time frame of the reference file.
         self.check_epsilons('vnv_2:T3DRES',
                             'f3d_particles.slf',
-                            eps=[4e-5, 0.01, 7e-3, 3e-3])
+                            eps=[1e-8])
 
         # Comparison between sequential and parallel run.
         self.check_epsilons('vnv_1:T3DRES',
                             'vnv_2:T3DRES',
-                            eps=[7e-5, 0.014, 7e-3, 3e-3])
+                            eps=[1e-8])
 
 
     def _post(self):