d=load('Movie4level.dat');

t=d(:,1);
z=d(:,2);
t=t*1e6;  % switch to uS
Om1=d(:,3);
Om2=d(:,4);
Om3=d(:,5);
Om_vac=d(:,6);

% figure(1)
% plot3( t, z, Om1);
xlabel('Time (uS)');
ylabel('z (m)');
zlabel('Rabi frequency');


%% End points time dependence

d=load('Movie4level_EndPoints.dat');
t=d(:,1);
t=t*1e6;  % switch to uS
Om1_strt=d(:,2);
Om_vac_strt=d(:,3);
Om1_end=d(:,4);
Om_vac_end=d(:,5);

figure(2)
plot(t, Om1_strt, '-;Om1_{strt};', t, Om_vac_strt, '-;Om_{vac_{strt}};', t,  Om1_end, '-;Om1_{end};', t,Om_vac_end, '-;Om_{vac_{end}};');
xlabel('Time (uS)');
ylabel('Rabi frequency');