Nlevels_no_dopler_no_z;

figure(1)
imagesc(t_1, d1_1, r13_imOut_1)
colorbar;
xlabel('t (uS)')
ylabel('d_1 (MHz)')
zlabel('r13_imOut_1')
title('Imag(rho_{13}) vs time and d_2 detuning')


figure(2)
plot(d1_1, r13_imOut_1(:,end),'-')
xlabel('d_1 (MHz)')
ylim ([-6e-4,0])
title('Imag(rho_{13}) vs  d_1 detuning')
print('imag_r13.pdf');

figure(3)
plot(d1_1,r13_rlOut_1(:,end),'-')
xlabel('d_1 (MHz)')
title('Real(rho_{13}) vs  d_1 detuning')
print('real_r13.pdf');