Nlevels_no_dopler_no_z;

figure(1)
imagesc(t_1, d2_1, r23_imOut_1)
colorbar;
xlabel('t')
ylabel('d_2')
zlabel('r23_imOut_1')
title('Imag(rho_{23}) vs time and d_2 detuning')


figure(2)
plot(d2_1,r23_imOut_1(:,end),'-')
xlabel('d_2')
title('Imag(rho_{23}) vs  d_2 detuning')
figure(3)
plot(d2_1,r23_rlOut_1(:,end),'-')
xlabel('d_2')
title('Real(rho_{23}) vs  d_2 detuning')