Add Tolerance to S-DTW

[Steenbuck]

Hey,

for my masterthesis on I want to implement a tolerance proposed in this paper:
Time-series Averaging Using Constrained Dynamic Time Warping with Tolerance

The aim is: "To better represent a set of signals, we added to its modelling the tolerance that depicts
admissible values around the average signal. This tolerance is computed for each time
step of the average signal. It corresponds to the standard deviation of the sets of signals
aligned on each point of the average time-series. The tolerance of the CDBA at time
step k is thus computed as the standard deviation of the points aligned on this time step
k of the average time-series."
The pseudo-code looks like this:

In order to produce a figure like this we need to add the standard deviation for the list of points which are regarded for computing each respective point of the centroid:

Can someone point me in the right direction where to start? Like in which function are the centroids computed and which variables I need to mess with?

Thank you very much