good job!
But you know that ‘w’ is only the wrapped phase? If you want to have metric velocities you have to use ‘v’. I am referring to it because you asked this recently.
You find a list of all commands and variables here: https://github.com/dbekaert/StaMPS/blob/master/matlab/ps_plot.m