Index: /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_a.m
===================================================================
--- /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_a.m	(revision 7937)
+++ /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_a.m	(revision 7938)
@@ -12,6 +12,6 @@
 %--------------------------------------------------------
 
-clear all;
-close all;
+%clear all; %ISSM->THIS LINE HAS BEEN CHANGED
+%close all; %ISSM->THIS LINE HAS BEEN CHANGED
 
 % The names of the folders (also prefix for filenames)
Index: /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_b.m
===================================================================
--- /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_b.m	(revision 7937)
+++ /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_b.m	(revision 7938)
@@ -12,6 +12,6 @@
 %--------------------------------------------------------
 
-clear all;
-close all;
+%clear all; %ISSM->THIS LINE HAS BEEN CHANGED
+%close all; %ISSM->THIS LINE HAS BEEN CHANGED
 
 % The names of the folders (also prefix for filenames)
@@ -167,5 +167,5 @@
 
 % The font to use for the axis labels
-axisLabelFont = 'normal';
+axisLabelFont = 'Arial'; %ISSM->THIS LINE HAS BEEN CHANGED
 % The font size to use for the axis labels
 axisLabelSize = 7;
@@ -246,10 +246,11 @@
 
     % Plot full stokes
-    [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ...
-                       fsFillColor, fsFillColor, 1, 0.5);
-    
-    % Plot non-full stokes
-    [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ...
-                       nfsFillColor, nfsFillColor, 1, 0.35);               
+    %[fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ... %ISSM->THIS LINE HAS BEEN CHANGED
+    %                   fsFillColor, fsFillColor, 1, 0.5);%ISSM->THIS LINE HAS BEEN CHANGED
+    %
+    %% Plot non-full stokes
+    %[fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ...%ISSM->THIS LINE HAS BEEN CHANGED
+    %                   nfsFillColor, nfsFillColor, 1, 0.35);               %ISSM->THIS LINE HAS BEEN CHANGED
+	 [fh1 fh2]=jbfill2(newX,fsMeanVel+fsStdVel,fsMeanVel-fsStdVel,nfsMeanVel+nfsStdVel,nfsMeanVel-nfsStdVel);
            
     % Turn the grid on
@@ -284,10 +285,11 @@
     % These lines plot the mean
     hold on;
-    plot(newX, fsMeanVel, 'k');
-    set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2);
-    plot(newX, nfsMeanVel, 'b');
-    set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2);
-    
-end
+    %plot(newX, fsMeanVel, 'k');%ISSM->THIS LINE HAS BEEN CHANGED 
+    %set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2);
+    %plot(newX, nfsMeanVel, 'b');%ISSM->THIS LINE HAS BEEN CHANGED 
+    %set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2);
+    
+end
+return %ISSM->THIS LINE HAS BEEN CHANGED
 
 % Add the legend to the final subplot (in upper left corner)
@@ -297,5 +299,4 @@
         'portrait','papertype','A4','paperunits', 'centimeters', ...
         'paperposition',[1 4 21 13]);
-return %ISSM->THIS LINE HAS BEEN CHANGED
 saveas(gcf, 'expb.eps', 'psc2');
 
Index: /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_c.m
===================================================================
--- /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_c.m	(revision 7937)
+++ /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_c.m	(revision 7938)
@@ -12,6 +12,6 @@
 %--------------------------------------------------------
 
-clear all;
-close all;
+%clear all; %ISSM->THIS LINE HAS BEEN CHANGED
+%close all; %ISSM->THIS LINE HAS BEEN CHANGED
 
 % The names of the folders (also prefix for filenames)
@@ -175,5 +175,5 @@
 
 % The font to use for the axis labels
-axisLabelFont = 'normal';
+axisLabelFont = 'Arial'; %ISSM->THIS LINE HAS BEEN CHANGED
 % The font size to use for the axis labels
 axisLabelSize = 7;
@@ -254,10 +254,11 @@
 
     % Plot full stokes
-    [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ...
-                       fsFillColor, fsFillColor, 1, 0.5);
-    
-    % Plot non-full stokes
-    [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ...
-                       nfsFillColor, nfsFillColor, 1, 0.35);               
+%    [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ...%ISSM->THIS LINE HAS BEEN CHANGED
+%                       fsFillColor, fsFillColor, 1, 0.5);%ISSM->THIS LINE HAS BEEN CHANGED
+%    
+%    % Plot non-full stokes
+%    [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ...%ISSM->THIS LINE HAS BEEN CHANGED
+%                       nfsFillColor, nfsFillColor, 1, 0.35);               %ISSM->THIS LINE HAS BEEN CHANGED
+		[fh1 fh2]=jbfill2(newX,fsMeanVel+fsStdVel,fsMeanVel-fsStdVel,nfsMeanVel+nfsStdVel,nfsMeanVel-nfsStdVel);
            
     % Turn the grid on
@@ -292,10 +293,11 @@
     % These lines plot the mean
     hold on;
-    plot(newX, fsMeanVel, 'k');
-    set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2);
-    plot(newX, nfsMeanVel, 'b');
-    set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2);
-    
-end
+    %plot(newX, fsMeanVel, 'k');%ISSM->THIS LINE HAS BEEN CHANGED 
+    %set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2);
+    %plot(newX, nfsMeanVel, 'b');%ISSM->THIS LINE HAS BEEN CHANGED 
+    %set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2);
+    
+end
+return %ISSM->THIS LINE HAS BEEN CHANGED
 
 % Add the legend to the final subplot (in upper left corner)
@@ -305,5 +307,4 @@
         'portrait','papertype','A4','paperunits', 'centimeters', ...
         'paperposition',[1 4 21 13]);
-return %ISSM->THIS LINE HAS BEEN CHANGED
 saveas(gcf, 'expc.eps', 'psc2');
 
Index: /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_d.m
===================================================================
--- /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_d.m	(revision 7937)
+++ /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_d.m	(revision 7938)
@@ -12,6 +12,6 @@
 %--------------------------------------------------------
 
-clear all;
-close all;
+%clear all; %ISSM->THIS LINE HAS BEEN CHANGED
+%close all; %ISSM->THIS LINE HAS BEEN CHANGED
 
 % The names of the folders (also prefix for filenames)
@@ -167,5 +167,5 @@
 
 % The font to use for the axis labels
-axisLabelFont = 'normal';
+axisLabelFont = 'Arial'; %ISSM->THIS LINE HAS BEEN CHANGED
 % The font size to use for the axis labels
 axisLabelSize = 7;
@@ -246,10 +246,11 @@
 
     % Plot full stokes
-    [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ...
-                       fsFillColor, fsFillColor, 1, 0.5);
-    
-    % Plot non-full stokes
-    [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ...
-                       nfsFillColor, nfsFillColor, 1, 0.35);               
+%    [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ...%ISSM->THIS LINE HAS BEEN CHANGED
+%                       fsFillColor, fsFillColor, 1, 0.5);%ISSM->THIS LINE HAS BEEN CHANGED
+%    
+%    % Plot non-full stokes
+%    [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ...%ISSM->THIS LINE HAS BEEN CHANGED
+%                       nfsFillColor, nfsFillColor, 1, 0.35);               %ISSM->THIS LINE HAS BEEN CHANGED
+[fh1 fh2]=jbfill2(newX,fsMeanVel+fsStdVel,fsMeanVel-fsStdVel,nfsMeanVel+nfsStdVel,nfsMeanVel-nfsStdVel);
            
     % Turn the grid on
@@ -284,10 +285,11 @@
     % These lines plot the mean
     hold on;
-    plot(newX, fsMeanVel, 'k');
-    set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2);
-    plot(newX, nfsMeanVel, 'b');
-    set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2);
-    
-end
+    %plot(newX, fsMeanVel, 'k');%ISSM->THIS LINE HAS BEEN CHANGED 
+    %set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2);
+    %plot(newX, nfsMeanVel, 'b');%ISSM->THIS LINE HAS BEEN CHANGED 
+    %set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2);
+    
+end
+return %ISSM->THIS LINE HAS BEEN CHANGED
 
 % Add the legend to the final subplot (in upper left corner)
@@ -297,5 +299,4 @@
         'portrait','papertype','A4','paperunits', 'centimeters', ...
         'paperposition',[1 4 21 13]);
-return %ISSM->THIS LINE HAS BEEN CHANGED
 saveas(gcf, 'expd.eps', 'psc2');
 
Index: /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_e.m
===================================================================
--- /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_e.m	(revision 7937)
+++ /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_e.m	(revision 7938)
@@ -12,6 +12,6 @@
 %--------------------------------------------------------
 
-clear all;
-close all;
+%clear all; %ISSM->THIS LINE HAS BEEN CHANGED
+%close all; %ISSM->THIS LINE HAS BEEN CHANGED
 
 % The names of the folders (also prefix for filenames)
@@ -166,5 +166,5 @@
 
 % The font to use for the axis labels
-axisLabelFont = 'normal';
+axisLabelFont = 'Arial'; %ISSM->THIS LINE HAS BEEN CHANGED
 % The font size to use for the axis labels
 axisLabelSize = 7;
@@ -245,10 +245,11 @@
 
     % Plot full stokes
-    [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ...
-                         fsFillColor, fsFillColor, 1, 0.5);
-    
-    % Plot non-full stokes
-    [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ...
-                         nfsFillColor, nfsFillColor, 1, 0.35);               
+%    [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ...%ISSM->THIS LINE HAS BEEN CHANGED
+%                         fsFillColor, fsFillColor, 1, 0.5);%ISSM->THIS LINE HAS BEEN CHANGED
+%    
+%    % Plot non-full stokes
+%    [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ...%ISSM->THIS LINE HAS BEEN CHANGED
+%                         nfsFillColor, nfsFillColor, 1, 0.35);               %ISSM->THIS LINE HAS BEEN CHANGED
+[fh1 fh2]=jbfill2(newX,fsMeanVel+fsStdVel,fsMeanVel-fsStdVel,nfsMeanVel+nfsStdVel,nfsMeanVel-nfsStdVel);
            
     % Turn the grid on
@@ -283,10 +284,11 @@
     % These lines plot the mean
     hold on;
-    plot(newX, fsMeanVel, 'k');
-    set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2);
-    plot(newX, nfsMeanVel, 'b');
-    set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2);
-    
-end
+    %plot(newX, fsMeanVel, 'k');%ISSM->THIS LINE HAS BEEN CHANGED 
+    %set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2);
+    %plot(newX, nfsMeanVel, 'b');%ISSM->THIS LINE HAS BEEN CHANGED 
+    %set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2);
+    
+end
+return %ISSM->THIS LINE HAS BEEN CHANGED
 
 % Add the legend to the final subplot (in upper left corner)
@@ -297,5 +299,4 @@
         'portrait','papertype','A4','paperunits', 'centimeters', ...
         'paperposition',[1 4 11 17]);
-return %ISSM->THIS LINE HAS BEEN CHANGED
 saveas(gcf, 'expe.eps', 'psc2');
 
Index: /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_e_txz.m
===================================================================
--- /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_e_txz.m	(revision 7937)
+++ /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_e_txz.m	(revision 7938)
@@ -12,6 +12,6 @@
 %--------------------------------------------------------
 
-clear all;
-close all;
+%clear all; %ISSM->THIS LINE HAS BEEN CHANGED
+%close all; %ISSM->THIS LINE HAS BEEN CHANGED
 
 % The names of the folders (also prefix for filenames)
@@ -164,5 +164,5 @@
 
 % The font to use for the axis labels
-axisLabelFont = 'normal';
+axisLabelFont = 'Arial';%ISSM->THIS LINE HAS BEEN CHANGED
 % The font size to use for the axis labels
 axisLabelSize = 7;
@@ -243,10 +243,11 @@
 
     % Plot full stokes
-    [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ...
-                         fsFillColor, fsFillColor, 1, 0.5);
-    
-    % Plot non-full stokes
-    [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ...
-                         nfsFillColor, nfsFillColor, 1, 0.35);               
+%    [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ...%ISSM->THIS LINE HAS BEEN CHANGED
+%                         fsFillColor, fsFillColor, 1, 0.5);%ISSM->THIS LINE HAS BEEN CHANGED
+%    
+%    % Plot non-full stokes
+%    [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ...%ISSM->THIS LINE HAS BEEN CHANGED
+%                         nfsFillColor, nfsFillColor, 1, 0.35);               %ISSM->THIS LINE HAS BEEN CHANGED
+[fh1 fh2]=jbfill2(newX,fsMeanVel+fsStdVel,fsMeanVel-fsStdVel,nfsMeanVel+nfsStdVel,nfsMeanVel-nfsStdVel);
            
     % Turn the grid on
@@ -281,10 +282,11 @@
     % These lines plot the mean
     hold on;
-    plot(newX, fsMeanVel, 'k');
-    set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2);
-    plot(newX, nfsMeanVel, 'b');
-    set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2);
+    %plot(newX, fsMeanVel, 'k');%ISSM->THIS LINE HAS BEEN CHANGED 
+    %set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2);
+    %plot(newX, nfsMeanVel, 'b');%ISSM->THIS LINE HAS BEEN CHANGED 
+    %set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2);
     
 end
+return %ISSM->THIS LINE HAS BEEN CHANGED
 
 % Add the legend to the final subplot (in upper left corner)
@@ -295,4 +297,3 @@
         'portrait','papertype','A4','paperunits', 'centimeters', ...
         'paperposition',[1 4 11 17]);
-return %ISSM->THIS LINE HAS BEEN CHANGED
 saveas(gcf, 'expetxz.eps', 'psc2');
Index: /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_f.m
===================================================================
--- /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_f.m	(revision 7937)
+++ /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_f.m	(revision 7938)
@@ -12,6 +12,6 @@
 %--------------------------------------------------------
 
-clear all;
-close all;
+%clear all; %ISSM->THIS LINE HAS BEEN CHANGED
+%close all; %ISSM->THIS LINE HAS BEEN CHANGED
 
 % The names of the folders (also prefix for filenames)
@@ -201,5 +201,5 @@
 
 % The font to use for the axis labels
-axisLabelFont = 'normal';
+axisLabelFont = 'Arial'; %ISSM->THIS LINE HAS BEEN CHANGED
 % The font size to use for the axis labels
 axisLabelSize = 7;
@@ -272,5 +272,5 @@
     
     % Tell MATLAB which subplot we are currently on
-    subplot(2, 1, i);
+    subplot(2, 2, 2*i-1); %ISSM->THIS LINE HAS BEEN CHANGED
     
     % The fill color (rgb) for full stokes
@@ -280,10 +280,11 @@
 
     % Plot full stokes
-    [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ...
-                         fsFillColor, fsFillColor, 1, 0.5);
-    
-    % Plot non-full stokes
-    [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ...
-                         nfsFillColor, nfsFillColor, 1, 0.35);               
+%    [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ...%ISSM->THIS LINE HAS BEEN CHANGED
+%                         fsFillColor, fsFillColor, 1, 0.5);%ISSM->THIS LINE HAS BEEN CHANGED
+%    
+%    % Plot non-full stokes
+%    [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ...%ISSM->THIS LINE HAS BEEN CHANGED
+%                         nfsFillColor, nfsFillColor, 1, 0.35);               %ISSM->THIS LINE HAS BEEN CHANGED
+[fh1 fh2]=jbfill2(newX,fsMeanVel+fsStdVel,fsMeanVel-fsStdVel,nfsMeanVel+nfsStdVel,nfsMeanVel-nfsStdVel);
            
     % Turn the grid on
@@ -318,8 +319,8 @@
     % These lines plot the mean
     hold on;
-    plot(newX, fsMeanVel, 'k');
-    set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2);
-    plot(newX, nfsMeanVel, 'b');
-    set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2);
+    %plot(newX, fsMeanVel, 'k');%ISSM->THIS LINE HAS BEEN CHANGED 
+    %set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2);
+    %plot(newX, nfsMeanVel, 'b');%ISSM->THIS LINE HAS BEEN CHANGED 
+    %set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2);
     
     if (i==1)
@@ -327,4 +328,5 @@
     end
 end
+return %ISSM->THIS LINE HAS BEEN CHANGED
 
 % Add the legend to the final subplot (in upper left corner)
@@ -335,5 +337,4 @@
         'portrait','papertype','A4','paperunits', 'centimeters', ...
         'paperposition',[1 4 11 17]);
-return %ISSM->THIS LINE HAS BEEN CHANGED
 saveas(gcf, 'expf.eps', 'psc2');
 
Index: /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_fs.m
===================================================================
--- /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_fs.m	(revision 7937)
+++ /issm/trunk/test/IsmipHomResults/plotting_scripts/plot_exp_fs.m	(revision 7938)
@@ -12,6 +12,6 @@
 %--------------------------------------------------------
 
-clear all;
-close all;
+%clear all; %ISSM->THIS LINE HAS BEEN CHANGED
+%close all; %ISSM->THIS LINE HAS BEEN CHANGED
 
 % The names of the folders (also prefix for filenames)
@@ -200,5 +200,5 @@
 
 % The font to use for the axis labels
-axisLabelFont = 'normal';
+axisLabelFont = 'Arial'; %ISSM->THIS LINE HAS BEEN CHANGED
 % The font size to use for the axis labels
 axisLabelSize = 7;
@@ -271,5 +271,5 @@
     
     % Tell MATLAB which subplot we are currently on
-    subplot(2, 1, i);
+    subplot(2, 2,2*i);  %ISSM->THIS LINE HAS BEEN CHANGED
     
     % The fill color (rgb) for full stokes
@@ -279,10 +279,11 @@
 
     % Plot full stokes
-    [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ...
-                         fsFillColor, fsFillColor, 1, 0.5);
-    
-    % Plot non-full stokes
-    [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ...
-                         nfsFillColor, nfsFillColor, 1, 0.35);               
+%    [fh1, msg1] = jbfill(newX, fsMeanVel+fsStdVel, fsMeanVel-fsStdVel, ...%ISSM->THIS LINE HAS BEEN CHANGED
+%                         fsFillColor, fsFillColor, 1, 0.5);%ISSM->THIS LINE HAS BEEN CHANGED
+%    
+%    % Plot non-full stokes
+%    [fh2, msg2] = jbfill(newX, nfsMeanVel+nfsStdVel, nfsMeanVel-nfsStdVel, ...%ISSM->THIS LINE HAS BEEN CHANGED
+%                         nfsFillColor, nfsFillColor, 1, 0.35);               %ISSM->THIS LINE HAS BEEN CHANGED
+[fh1 fh2]=jbfill2(newX,fsMeanVel+fsStdVel,fsMeanVel-fsStdVel,nfsMeanVel+nfsStdVel,nfsMeanVel-nfsStdVel);
            
     % Turn the grid on
@@ -317,8 +318,8 @@
     % These lines plot the mean
     hold on;
-    plot(newX, fsMeanVel, 'k');
-    set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2);
-    plot(newX, nfsMeanVel, 'b');
-    set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2);
+    %plot(newX, fsMeanVel, 'k');%ISSM->THIS LINE HAS BEEN CHANGED 
+    %set(findobj(gca,'Type','line','Color',[0 0 0]),'Color',fsFillColor,'LineWidth',2);
+    %plot(newX, nfsMeanVel, 'b');%ISSM->THIS LINE HAS BEEN CHANGED 
+    %set(findobj(gca,'Type','line','Color',[0 0 1]),'Color',nfsFillColor,'LineWidth',2);
     
     if (i==1)
@@ -326,4 +327,5 @@
     end
 end
+return %ISSM->THIS LINE HAS BEEN CHANGED
 
 % Add the legend to the final subplot (in upper left corner)
@@ -334,5 +336,4 @@
         'portrait','papertype','A4','paperunits', 'centimeters', ...
         'paperposition',[1 4 11 17]);
-return %ISSM->THIS LINE HAS BEEN CHANGED
 saveas(gcf, 'expfs.eps', 'psc2');
 
