Unwrap phase tolerance for high phase jumps of moving waves

I want to unwrap my phase but I am not satisfied with my result.

This is the code with the data. The figure is attached.

close all;
clear all;

phase = [134.082437184695;131.369319207920;124.056305092615;116.540978092015;112.526324104033;113.056718979828;116.968087321545;123.995738827168;137.210057313570;146.800175033853;148.916121642133;147.731645164655;145.624232853363;142.910999352223;137.446991994084;125.236287449535;111.465807072913;103.150199921263;98.0321061900557;94.3189471811752;91.8989547166367;90.1309722781345;88.9013493022229;89.3063040271042;90.9468518625895;128.733467921401;93.0322615156086;94.1821339425437;94.6421373969817;92.8233528723675;90.3756373212097;87.6054759585859;85.0790399879507;81.4817459394186;77.2931501364155;72.4943501159715;65.9954392536164;59.2667816057140;53.0772780510376;46.8156719823233;40.4975951438664;33.6488140355637;25.8305297876533;17.7497642454968;10.0447093941267;2.01195686150996;-4.93373373064492;-10.8058476964741;-16.1210410186714;-20.9620668645680;-25.8741662865119;-31.4143243790756;-37.9386960367784;-46.4054314104185;-57.6985446485876;-71.2392737040546;-88.2211341802578;87.0875726325545;-128.700403163520;-147.400646665315;-162.718657833210;-175.985443331330;170.507152091084;157.284618393489;142.141606026656;125.477737739339;106.008289931893;77.7803212724172;42.8308137865843;7.96900785935225;-18.7492759030433;-40.3835610415339;-56.1073993057171;-70.0894379038518;-83.0499549702538;-95.7485694605537;-109.351268931910;-124.175167790244;-141.040343736190;-157.148534839110;-174.455440959497;167.950099474436;151.366245927154;137.377640833089;122.393941546082;106.766670802346;89.6303406191681;68.3853546756871;38.7879856992843;-1.06096691331800;-50.2270290615035;-77.6950685658825;-99.1874910038732;-114.382960942411;-131.368412251259;-153.262330771990;172.677831003320;134.687678481674;104.957799070157;73.6691293830629;39.7386137104618;2.96388294088452;-33.3771786357528;-63.4087720459665;-96.7061015923250;-126.923023240502;-158.008841433759;177.468787434736;149.493546509929;137.143669028367;142.416119400809;148.594469354469;-136.179231721167;-117.902955400767;-128.348268183869;-148.087555054063;-172.697668592545;165.609192745014;147.608220935904;142.729036376178;141.886994177019;176.537944149926;-74.1129490231803;-92.5472376711897;-112.326724867163;-142.913197409068;-170.813432603191;173.682119344360;163.862343736686;163.272323604883;130.692479726823;-87.8287095561921;-56.7016830804346;-50.9671662891455;-74.4042022157955;-118.183961977711;-159.362421028501;-167.883138684716;-162.840111600614;-153.476716454580;-129.584028942079;-74.6467067131749;-35.6386721758194;-97.9246367889239;152.693905380832;-175.914913220363;-171.148447373316;-159.035478834435;-146.256217005272;-140.525979204032;-135.281948003243;123.247915295342;115.034673112608;141.989917436826;163.333486589713;-176.849715088465;-157.489948684471;-136.501231135606;-113.398116414271;-90.0330750499792;46.6871442113775;104.587096461145;129.424194884685;148.001519815181;168.622479089440;-172.833698398609;-153.054969901139;-142.335487445416];
phase = phase*2*pi/360;
unwraped_phase = unwrap(phase);
x= (1:1:168)';

plot(x,phase,'b',x, unwraped_pha

se,'k') xlabel('Distance'); ylabel('Phase in \pi');

Due to the fact that I am measuring moving waves it is not uncommon to meassure this high phases. How can I change the phase jumps in the red ellipses in the other direction. I tried different tolerances like 3*pi or pi/2 but with now actual succes. Is there a way?

Best regards Henrik