4th field plotted

1 files changed, 26 insertions, 9 deletions

diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pp.m b/xmds2/Nlevels_no_dopler_with_z_4wm/pp.m
index ab536be..186d822 100644
--- a/xmds2/Nlevels_no_dopler_with_z_4wm/pp.m
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pp.m
@@ -4,21 +4,27 @@ Nlevels_no_dopler_with_z
 z_1=z_1*100; % z in cm
 t_1=t_1*1e6; % time now measured in uS
 figure(1)
-subplot(1,3,1); imagesc(z_1, t_1, I1_out_1); colorbar
+subplot(2,2,1); imagesc(z_1, t_1, I1_out_1); colorbar
 xlabel('z (cm)')
 ylabel('t (uS)')
 zlabel('I_1')
 title('I_1')
-subplot(1,3,2); imagesc(z_1, t_1, I2_out_1); colorbar
+subplot(2,2,2); imagesc(z_1, t_1, I2_out_1); colorbar
 xlabel('z (cm)')
 ylabel('t (uS)')
 zlabel('I_2')
 title('I_2')
-subplot(1,3,3); imagesc(z_1, t_1, I3_out_1); colorbar
+subplot(2,2,3); imagesc(z_1, t_1, I3_out_1); colorbar
 xlabel('z (cm)')
 ylabel('t (uS)')
 zlabel('I_3')
 title('I_3')
+subplot(2,2,4); imagesc(z_1, t_1, I4_out_1); colorbar
+xlabel('z (cm)')
+ylabel('t (uS)')
+zlabel('I_4')
+title('I_4')
+
 
 print('-color','fields_propagation.eps')
 
@@ -26,7 +32,7 @@ print('-color','fields_propagation.eps')
 
 %% fields before and after the cell
 figure(2)
-subplot(1,3,1);
+subplot(2,2,1);
 plot( ...
 	t_1,I1_out_1(:,1),'-;before;', "linewidth", 4,  ...
 	t_1,I1_out_1(:,end), '-;after;', "linewidth", 4 ...
@@ -34,7 +40,7 @@ plot( ...
 xlabel('t (uS)')
 ylabel('I_1 (1/s)^2')
 title('I_1 before and after cell')
-subplot(1,3,2);
+subplot(2,2,2);
 plot( ...
 	t_1,I2_out_1(:,1),'-;before;', "linewidth", 4,  ...
 	t_1,I2_out_1(:,end), '-;after;', "linewidth", 4 ...
@@ -42,7 +48,7 @@ plot( ...
 xlabel('t (uS)')
 ylabel('I_2 (1/s)^2')
 title('I_2 before and after cell')
-subplot(1,3,3);
+subplot(2,2,3);
 plot( ...
 	t_1,I3_out_1(:,1),'-;before;', "linewidth", 4,  ...
 	t_1,I3_out_1(:,end), '-;after;', "linewidth", 4 ...
@@ -60,9 +66,20 @@ printf('Second field delay time  = %f uS\n',delay_time);
 
 print('-color','fields_before_after_cell.eps')
 
+subplot(2,2,4);
+plot( ...
+	t_1,I4_out_1(:,1),'-;before;', "linewidth", 4,  ...
+	t_1,I4_out_1(:,end), '-;after;', "linewidth", 4 ...
+	)
+xlabel('t (uS)')
+ylabel('I_3 (1/s)^2')
+title('I_3 before and after cell')
+
 figure(4)
-I2_max_in=max(I2_out_1(:,1));
-I2_max_out=max(I2_out_1(:,end));
+settling_time=0.5; %uS
+t_good_indx=t_1> min(t_1 + 0.5);
+I2_max_in=max(I2_out_1(t_good_indx,1));
+I2_max_out=max(I2_out_1(t_good_indx,end));
 I2_in_norm=(I2_out_1(:,1))/I2_max_in;
 I2_out_norm=(I2_out_1(:,end))/I2_max_out;
 tmin=-0.05;
@@ -78,7 +95,7 @@ ylabel('I_2')
 title('I_2 before and after cell normalized')
 print('-color','probe_before_after_cell.eps')
 
-%return;
+return;
 
 %% all density matrix elements in one plot
 % diagonal populations,