diff --git a/moored_adcp_proc/adcp_check_surface.m b/moored_adcp_proc/adcp_check_surface.m
index 129b514d3e40da70b20f5d7da09cd3bd94295e0f..4ef402b7aa7cddd7f59ba767a038c619d510eb6f 100644
--- a/moored_adcp_proc/adcp_check_surface.m
+++ b/moored_adcp_proc/adcp_check_surface.m
@@ -7,6 +7,7 @@ orient tall;
 subplot 321;
 imagesc(time,bin(bins),u(bins,:),[-1 1]);
 set(gca,'YTick',bins);
+ylabel('Bins');
 hold on;
 gregtick;
 title('Meridional velocity before correction');
@@ -15,6 +16,7 @@ subplot 323;
 imagesc(time,bin(bins),u1(bins,:),[-1 1]);
 set(gca,'YTick',bins);
 hold on;
+ylabel('Bins');
 gregtick;
 title('Meridional velocity after correction');
 
@@ -22,13 +24,15 @@ subplot 322;
 imagesc(time,bin(bins),v(bins,:),[-1 1]);
 set(gca,'YTick',bins);
 hold on;
+ylabel('Bins');
 gregtick;
-title('Zonal velocity after correction');
+title('Zonal velocity before correction');
 
 subplot 324;
 imagesc(time,bin(bins),v1(bins,:),[-1 1]);
 set(gca,'YTick',bins);
 hold on;
+ylabel('Bins');
 gregtick;
 title('Zonal velocity after correction');
 
@@ -39,6 +43,7 @@ for i=1:length(bins);
 end
 set(gca,'XLim',[time(1) time(end)]);
 gregtick;
+ylabel('Bins');
 title('Bin depth after correction');
 
 subplot 326;
@@ -48,4 +53,5 @@ for i=1:length(bins);
 end
 set(gca,'XLim',[time(1) time(end)]);
 gregtick;
+ylabel('Bins');
 title('Bin depth after correction');
diff --git a/moored_adcp_proc/adcp_filt_sub.m b/moored_adcp_proc/adcp_filt_sub.m
index a1764eb2c00686ff8cdaf2fa4f7f83cfe03ff61d..65c211d469531c2cf5d7fce4f34118fb0cf8808e 100644
--- a/moored_adcp_proc/adcp_filt_sub.m
+++ b/moored_adcp_proc/adcp_filt_sub.m
@@ -63,30 +63,30 @@ for iidep = 1:length(intdepvec)
       
 end
 
-
-figure;
-subplot 122;
-imagesc(data.time,intdepvec,uintfilt,[-1 1]);
-set(gca,'YLim',[min(intdepvec)-10 max(intdepvec)+10]);
-gregtick;
-title('data filtered and subsampled');
-
+hf=figure('position', [0, 0, 1400, 1000]);
 subplot 121;
 imagesc(data.time,intdepvec,ui,[-1 1]);
-set(gca,'YLim',[min(intdepvec)-10 max(intdepvec)+10]);
+%set(gca,'YLim',[min(intdepvec)-10 max(intdepvec)+10]);
 gregtick;
-title('data raw');
-
-
-figure;
+ylabel('Bins');
+title('U field - data raw');
 subplot 122;
-imagesc(data.time,intdepvec,vintfilt,[-1 1]);
-set(gca,'YLim',[min(intdepvec)-10 max(intdepvec)+10]);
+imagesc(data.time,intdepvec,uintfilt,[-1 1]);
+%set(gca,'YLim',[min(intdepvec)-10 max(intdepvec)+10]);
 gregtick;
-title('data filtered and subsampled');
+ylabel('Bins');
+title('U field - data interpolated, filtered and subsampled');
 
+hf=figure('position', [0, 0, 1400, 1000]);
 subplot 121;
 imagesc(data.time,intdepvec,vi,[-1 1]);
-set(gca,'YLim',[min(intdepvec)-10 max(intdepvec)+10])
+%set(gca,'YLim',[min(intdepvec)-10 max(intdepvec)+10])
+gregtick;
+ylabel('Bins');
+title('V field - data raw');
+subplot 122;
+imagesc(data.time,intdepvec,vintfilt,[-1 1]);
+%set(gca,'YLim',[min(intdepvec)-10 max(intdepvec)+10]);
 gregtick;
-title('data raw');
+ylabel('Bins');
+title('V field - data interpolated, filtered and subsampled');
diff --git a/moored_adcp_proc/adcp_surface_fit.m b/moored_adcp_proc/adcp_surface_fit.m
index ea4d9774d15f8d19967e563f5672db939a3b5912..9461dfd84100310c44a9386b97a269ca8364205b 100644
--- a/moored_adcp_proc/adcp_surface_fit.m
+++ b/moored_adcp_proc/adcp_surface_fit.m
@@ -46,6 +46,8 @@ function [zbins,zadcp1,offset,x_null]=adcp_surface_fit(zadcp,ea,surface_bins,ble
     figure(1);
     bar(count,ncount);
     title('Histogram of differences between original and reconstructed depth record');
+    ylabel('Frequency');
+    xlabel('Difference of Depth (m)');
     
     figure(2);
     plot(zadcp,'b');
@@ -75,12 +77,18 @@ function [zbins,zadcp1,offset,x_null]=adcp_surface_fit(zadcp,ea,surface_bins,ble
     text(300, max(zadcp),['Offset applied: ' num2str(offset) ' m']);
     %,'fonts',12,'fontw','bold','backgroundc','w');
     legend('Original','Reconstructed from surface reflection','Offset applied');
+    ylabel('Depth (m)');
+    xlabel('Time index');
     
     %print -dpng surface_fit;
     
     figure(3);
     pcolor([1:length(x_null)],zbins,ea); shading flat;
-    
+    set(gca,'ydir', 'reverse');
+    ylabel('Corrected Depth (m)');
+    xlabel('Time index');
+    title('Amplitude of the bins with corrected depth');
+
     %print -dpng surface_ea;
 
 end
\ No newline at end of file