diff --git a/user-doc/tutorials/belfast-4.txt b/user-doc/tutorials/belfast-4.txt
index 3a5841b87775163338fa842cf43be8f6d158cbe7..3531956e6a18c6cc342554fcef3f1f07ae31f907 100644
--- a/user-doc/tutorials/belfast-4.txt
+++ b/user-doc/tutorials/belfast-4.txt
@@ -425,7 +425,7 @@ can be mixed and it is not necessary to keep track of which restraint was used t
 every single frame. Let's get the concatenated trajectory
 
 \verbatim
-trjcat -cat -f traj*.xtc -o alltraj.xtc
+gmx_mpi trjcat -cat -f traj*.xtc -o alltraj.xtc
 \endverbatim
 
 Now we should compute the value of each of the bias potentials on the entire (concatenated) trajectory
diff --git a/user-doc/tutorials/belfast-7.txt b/user-doc/tutorials/belfast-7.txt
index 1cd873611feecfca77b41158ef96582e1929a64f..b81abd651fd9751472c1e852a54e96533b9542c3 100644
--- a/user-doc/tutorials/belfast-7.txt
+++ b/user-doc/tutorials/belfast-7.txt
@@ -177,7 +177,7 @@ so there seems not to be any problem on this side.
 We can use the second file to reconstruct the continuous trajectories of each replica in temperature:
 
 \verbatim
-trjcat_mpi -f traj0.trr traj1.trr -demux replica_index.xvg
+gmx_mpi trjcat -f traj0.trr traj1.trr -demux replica_index.xvg
 \endverbatim
 
 and the following PLUMED input file (plumed_demux.dat) to recalculate the value of the CVs on the demuxed trajectories,