From 846b4e16e2b99b27c7dfc8e90f83383de75ca27c Mon Sep 17 00:00:00 2001
From: Pierre Rousselot <pier.rousselot@gmail.com>
Date: Fri, 12 May 2017 16:16:00 +0200
Subject: [PATCH] Version reviewed

---
 moored_adcp_proc/adcp_check_surface.m |  8 +++++-
 moored_adcp_proc/adcp_filt_sub.m      | 36 +++++++++++++--------------
 moored_adcp_proc/adcp_surface_fit.m   | 10 +++++++-
 3 files changed, 34 insertions(+), 20 deletions(-)

diff --git a/moored_adcp_proc/adcp_check_surface.m b/moored_adcp_proc/adcp_check_surface.m
index 129b514..4ef402b 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 a1764eb..65c211d 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 ea4d977..9461dfd 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
-- 
GitLab