diff --git a/Trajectories/trajectories.py b/Trajectories/trajectories.py
index afd8efc102f43a087210236da114fe6deba883b4..8cb7b8a9498b1002165cd90ea8f808f712a23de4 100755
--- a/Trajectories/trajectories.py
+++ b/Trajectories/trajectories.py
@@ -29,6 +29,7 @@ def new_traj(ind_traj, traj_prop, n, traj_vert_ind):
 
 if len(sys.argv) != 2: sys.exit("Required argument: input-graph")
 g = graph_tool.load_graph(sys.argv[1])
+print("Loading done")
 
 if "cost_function" in g.edge_properties:
     closest_succ = g.new_vertex_property('int')
@@ -61,6 +62,8 @@ traj_vert_ind = []
 # An element of traj_vert_ind is a trajectory, represented by a list
 # of vertex indices.
 
+print("Starting loop to define traj_prop...")
+
 for n in topology.topological_sort(g):
     # First, we check whether we have a phantom pattern. Phantom
     # eddies may appear in altimetric data because of
@@ -166,6 +169,8 @@ for t in traj_vert_ind:
     traj_segm.append(list_segm)
     expanded_traj.append(list_eddies_traj)
 
+print("Saving...")
+
 with open("traj_segm.json", "w") as outfile:
     json.dump(traj_segm, outfile, indent = 4)