diff --git a/Utilities/Matlab/scripts/atmospheric_boundary_layer.m b/Utilities/Matlab/scripts/atmospheric_boundary_layer.m
index 5469e4f4ed0..79e43a7c5d1 100644
--- a/Utilities/Matlab/scripts/atmospheric_boundary_layer.m
+++ b/Utilities/Matlab/scripts/atmospheric_boundary_layer.m
@@ -48,7 +48,7 @@
 
 z = [z_0{i} 10*z_0{i} 1 2 3 4 5 6 7 8 9 10 15 20 25 30 50 100];
 
-figure(1)
+f1=figure;
 set(gca,'Units',Plot_Units)
 set(gca,'Position',[Plot_X Plot_Y Plot_Width Plot_Height])
 
@@ -68,7 +68,7 @@
 T = T + (theta_0-T(12));
 
 plot(u,z,'ko'); hold on
-plot(M2.data(:,2),M2.data(:,1),'k-'); hold on
+plot(M2.data(:,2),M2.data(:,1),'k-'); hold off
 set(gca,'FontName',Font_Name)
 set(gca,'FontSize',Title_Font_Size)
 axis([0 u_high{i} 0 100])
@@ -88,21 +88,19 @@
 git_file = [datafile, '_git.txt'];
 addverstr(gca,git_file,'linear')
 
-set(gcf,'Visible',Figure_Visibility);
-set(gcf,'Units',Paper_Units);
-set(gcf,'PaperUnits',Paper_Units);
-set(gcf,'PaperSize',[Paper_Width Paper_Height]);
-set(gcf,'Position',[0 0 Paper_Width Paper_Height]);
-print(gcf,'-dpdf',[outfile '_vel'])
+set(f1,'Visible',Figure_Visibility);
+set(f1,'Units',Paper_Units);
+set(f1,'PaperUnits',Paper_Units);
+set(f1,'PaperSize',[Paper_Width Paper_Height]);
+set(f1,'Position',[0 0 Paper_Width Paper_Height]);
+print(f1,'-dpdf',[outfile '_vel'])
 
-clf('reset')
-
-figure(2)
+f2=figure;
 set(gca,'Units',Plot_Units)
 set(gca,'Position',[Plot_X Plot_Y Plot_Width Plot_Height])
 
 plot(T-273.15,z,'ko'); hold on
-plot(M2.data(:,3),M2.data(:,1),'k-'); hold on
+plot(M2.data(:,3),M2.data(:,1),'k-'); hold off
 set(gca,'FontName',Font_Name)
 set(gca,'FontSize',Title_Font_Size)
 axis([T_low{i} T_high{i} 0 100])
@@ -114,13 +112,11 @@
 
 addverstr(gca,git_file,'linear')
 
-set(gcf,'Visible',Figure_Visibility);
-set(gcf,'Units',Paper_Units);
-set(gcf,'PaperUnits',Paper_Units);
-set(gcf,'PaperSize',[Paper_Width Paper_Height]);
-set(gcf,'Position',[0 0 Paper_Width Paper_Height]);
-print(gcf,'-dpdf',[outfile '_tmp'])
-
-clf('reset')
+set(f2,'Visible',Figure_Visibility);
+set(f2,'Units',Paper_Units);
+set(f2,'PaperUnits',Paper_Units);
+set(f2,'PaperSize',[Paper_Width Paper_Height]);
+set(f2,'Position',[0 0 Paper_Width Paper_Height]);
+print(f2,'-dpdf',[outfile '_tmp'])
 
 end