files renamed to better reflect the processes

10 files changed, 61 insertions, 61 deletions

diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/map2dat.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/map2dat.m
deleted file mode 100644
index 969b6dc..0000000
--- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/map2dat.m
+++ /dev/null
@@ -1,33 +0,0 @@
-function map2dat(outfile, x,y,z, xskip, yskip)
-% saves 3D data in suitable way to be drawn by gnuplot
-% x,y - vectors of x,y values
-% z map of z values as used by Octave/Matlab
-% xskip, yskip - skip paprameters
-%  only every, xskip, yskip point will be written
-
-
-
-Nx=length(x);   
-Ny=length(y);   
-Nxs=Nx/xskip;
-Nys=Ny/yskip;
-points=zeros(1,3*Nxs*Nys);
-%points=[];
-tic; 
-for i=1:Nxs
-	for k=1:Nys
-		%points=[points x(i*xskip) y(k*yskip) z(k*yskip,i*xskip)];
-		points((i-1)*(Nys-1)*3+3*(k-1)+1) =   x(i*xskip);
-		points((i-1)*(Nys-1)*3+3*(k-1)+2) = y(k*yskip);
-		points((i-1)*(Nys-1)*3+3*(k-1)+3) = z(k*yskip,i*xskip);
-	end
-end
-disp('=== points formation complete ===');
-toc
-tic;
-%points
-fd = fopen(outfile, "wt");
-fprintf (fd, "%g %g %g\n", points);
-fclose(fd);
-disp('=== points saving complete ===');
-toc;
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile
index 95ebb0a..76d3f5c 100644
--- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/Makefile
@@ -9,11 +9,9 @@ GNUPLOT_FILES = $(wildcard *.gp)
 XSIL2GRAPHICS = xsil2graphics
 
 # fast light
-# PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4
+PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4
 # slow light EIT
 #PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=0 --E4o=0
-#Fast light to Slow light switch
-PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=2e7 --E2o=3e3 --E3o=6e6 --E4o=3e2
 
 all: $(XSIL_FILES) Nlevels_no_dopler_with_z_4wm.xsil $(M_FILES) plot png
 
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/README b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/README
index 1de063d..1de063d 100644
--- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/README
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/README
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/map2dat.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/map2dat.m
index 14fae30..969b6dc 120000..100644
--- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/map2dat.m
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/map2dat.m
@@ -1 +1,33 @@
-../map2dat.m
\ No newline at end of file
+function map2dat(outfile, x,y,z, xskip, yskip)
+% saves 3D data in suitable way to be drawn by gnuplot
+% x,y - vectors of x,y values
+% z map of z values as used by Octave/Matlab
+% xskip, yskip - skip paprameters
+%  only every, xskip, yskip point will be written
+
+
+
+Nx=length(x);   
+Ny=length(y);   
+Nxs=Nx/xskip;
+Nys=Ny/yskip;
+points=zeros(1,3*Nxs*Nys);
+%points=[];
+tic; 
+for i=1:Nxs
+	for k=1:Nys
+		%points=[points x(i*xskip) y(k*yskip) z(k*yskip,i*xskip)];
+		points((i-1)*(Nys-1)*3+3*(k-1)+1) =   x(i*xskip);
+		points((i-1)*(Nys-1)*3+3*(k-1)+2) = y(k*yskip);
+		points((i-1)*(Nys-1)*3+3*(k-1)+3) = z(k*yskip,i*xskip);
+	end
+end
+disp('=== points formation complete ===');
+toc
+tic;
+%points
+fd = fopen(outfile, "wt");
+fprintf (fd, "%g %g %g\n", points);
+fclose(fd);
+disp('=== points saving complete ===');
+toc;
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp
index a0cd31f..8591d87 100644
--- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/plot_fields_propagation_I2.gp
@@ -12,4 +12,4 @@ 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' 
+splot [0:1.5][-0.2:0.2] 'I2.dat' 
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m
index ecfdbbb..0d46ec7 100644
--- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_to_slow_switch/pp_I2.m
@@ -4,8 +4,11 @@ Nlevels_no_dopler_with_z_4wm
 z_1=z_1*100; % z in cm
 t_1=t_1*1e6; % time now measured in uS
 figure(1)
-set(gca,'fontsize',20); 
+%set(gca,'fontsize',20); 
 imagesc(z_1, t_1, I2_out_1); colorbar
+tmin=-0.4;
+tmax= 0.4;
+ylim([tmin,tmax],'manual');	
 xlabel('z (cm)')
 ylabel('t (uS)')
 zlabel('I_2')
@@ -17,7 +20,7 @@ map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip);
 
 
 
-print('-color','fields_propagation_I2.eps')
+print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I2.eps')
 
 
 
@@ -43,26 +46,26 @@ 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')
+%set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts
+print('-color','-depsc2', '-tight','-S200,120', '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;
+tmin=-.05;
+tmax=.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 ...
 	)
+legend('location', 'southeast');
 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')
+%set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts
+print('-color','-depsc2', '-tight','-S200,120', 'probe_before_after_cell_I2_normalized.eps')
 
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/Makefile b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/Makefile
index 76d3f5c..95ebb0a 100644
--- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/Makefile
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/Makefile
@@ -9,9 +9,11 @@ GNUPLOT_FILES = $(wildcard *.gp)
 XSIL2GRAPHICS = xsil2graphics
 
 # fast light
-PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4
+# PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=3.8e7 --E4o=6.3e4
 # slow light EIT
 #PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=1.9e7 --E2o=3.1e5 --E3o=0 --E4o=0
+#Fast light to Slow light switch
+PARAMS =  --delta1=0 --delta2=0 --delta3=0 --E1o=2e7 --E2o=3e3 --E3o=6e6 --E4o=3e2
 
 all: $(XSIL_FILES) Nlevels_no_dopler_with_z_4wm.xsil $(M_FILES) plot png
 
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/map2dat.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/map2dat.m
new file mode 120000
index 0000000..14fae30
--- /dev/null
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/map2dat.m
@@ -0,0 +1 @@
+../map2dat.m
\ No newline at end of file
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/plot_fields_propagation_I2.gp b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/plot_fields_propagation_I2.gp
index 8591d87..a0cd31f 100644
--- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/plot_fields_propagation_I2.gp
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/plot_fields_propagation_I2.gp
@@ -12,4 +12,4 @@ set ylabel "t ({/Symbol m}S)"
 set zlabel "I_2 (1/S)"
 set nokey
 #set view map
-splot [0:1.5][-0.2:0.2] 'I2.dat' 
+splot [0:][-0.2:0.4] 'I2.dat' 
diff --git a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/pp_I2.m b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/pp_I2.m
index 0d46ec7..ecfdbbb 100644
--- a/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/fast_light/pp_I2.m
+++ b/xmds2/Nlevels_no_dopler_with_z_4wm_for_irina_pqe/pulse_split/pp_I2.m
@@ -4,11 +4,8 @@ Nlevels_no_dopler_with_z_4wm
 z_1=z_1*100; % z in cm
 t_1=t_1*1e6; % time now measured in uS
 figure(1)
-%set(gca,'fontsize',20); 
+set(gca,'fontsize',20); 
 imagesc(z_1, t_1, I2_out_1); colorbar
-tmin=-0.4;
-tmax= 0.4;
-ylim([tmin,tmax],'manual');	
 xlabel('z (cm)')
 ylabel('t (uS)')
 zlabel('I_2')
@@ -20,7 +17,7 @@ map2dat('I2.dat',z_1,t_1, I2_out_1, xskip, yskip);
 
 
 
-print('-color','-depsc2', '-tight', '-S200,120', 'fields_propagation_I2.eps')
+print('-color','fields_propagation_I2.eps')
 
 
 
@@ -46,26 +43,26 @@ 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','-depsc2', '-tight','-S200,120', 'fields_before_after_cell_I2.eps')
+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=-.05;
-tmax=.05;
+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 ...
 	)
-legend('location', 'southeast');
 xlim([tmin,tmax],'manual');	
 xlabel('t (uS)')
 ylabel('I_2')
 title('I_2 before and after cell normalized')
-%set (gcf,'paperposition',[0.5 0 2.5,1.5]); % IMPORTANT to shrink eps size for readable fonts
-print('-color','-depsc2', '-tight','-S200,120', 'probe_before_after_cell_I2_normalized.eps')
+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')