3 files changed, 3 insertions, 83 deletions
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp
index a0cd31f..e11c569 100644..120000
--- a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I2.gp
@@ -1,15 +1 @@
-set terminal postscript portrait enhanced color solid size 5,3.5
-set output 'fields_propagation_I2.eps'
-set dgrid3d 100,100 qnorm 4
-set pm3d map
-#set contour
-set hidden3d
-set palette rgb 10,13,31 negative
-
-
-set xlabel "z (cm)"
-set ylabel "t ({/Symbol m}S)"
-set zlabel "I_2 (1/S)"
-set nokey
-#set view map
-splot [0:][-0.2:0.4] 'I2.dat' 
+../plot_fields_propagation_I2.gp
+\ No newline at end of file
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I4.gp b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I4.gp
new file mode 120000
index 0000000..134a380
--- /dev/null
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/plot_fields_propagation_I4.gp
@@ -0,0 +1 @@
+../plot_fields_propagation_I4.gp
+\ No newline at end of file
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m
index ecfdbbb..4a1ed9b 100644..120000
--- a/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm/pulse_split/pp_I2.m
@@ -1,68 +1 @@
-Nlevels_no_dopler_with_z_4wm
-
-%% field propagation
-z_1=z_1*100; % z in cm
-t_1=t_1*1e6; % time now measured in uS
-figure(1)
-set(gca,'fontsize',20); 
-imagesc(z_1, t_1, I2_out_1); colorbar
-xlabel('z (cm)')
-ylabel('t (uS)')
-zlabel('I_2')
-title('I_2')
-
-xskip=1;
-yskip=10;
-map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip);
-
-
-
-print('-color','fields_propagation_I2.eps')
-
-
-
-%% fields before and after the cell
-figure(2)
-%set(gca,'fontsize',30); 
-plot( ...
-	t_1,I2_out_1(:,1),'.-;before;', "linewidth", 4,  ...
-	t_1,I2_out_1(:,end), '-;after;', "linewidth", 4 ...
-	)
-xlabel('t (uS)')
-ylabel('I_2 (1/s)^2')
-title('I_2 before and after cell')
-legend('location', 'northeast');
-
-[b, a]=butter(3, 0.05);
-I2_out_after=I2_out_1(:,end);
-I2_out_after_filtered=filtfilt(b,a,I2_out_after);
-settling_time=0.01; %uS
-t_good_indx=t_1> min(t_1 + settling_time);
-[m,max_pos_before]=max(I2_out_1(t_good_indx,1) ); [m,max_pos_after]=max(I2_out_after_filtered(t_good_indx));
-delay_time=t_1(max_pos_after)-t_1(max_pos_before);
-printf('Second field delay time  = %f uS\n',delay_time);
-
-%set(gca,'fontsize',40); 
-set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts
-print('-color','fields_before_after_cell_I2.eps')
-
-figure(4)
-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;
-tmax=0.05;
-indx=(t_1>=tmin & t_1<=tmax);  % soom in in time to this region
-plot( ...
-	t_1(indx),I2_in_norm(indx),'.-;before;', "linewidth", 4,  ...
-	t_1(indx),I2_out_norm(indx), '-;after;', "linewidth", 4 ...
-	)
-xlim([tmin,tmax],'manual');	
-xlabel('t (uS)')
-ylabel('I_2')
-title('I_2 before and after cell normalized')
-legend('location', 'northeast');
-set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts
-print('-color','probe_before_after_cell_I2_normalized.eps')
-
+../pp_I2.m
+\ No newline at end of file