Skip to content
Snippets Groups Projects
Commit 846b4e16 authored by Pierre Rousselot's avatar Pierre Rousselot
Browse files

Version reviewed

parent 6c112ee1
No related branches found
No related tags found
No related merge requests found
...@@ -7,6 +7,7 @@ orient tall; ...@@ -7,6 +7,7 @@ orient tall;
subplot 321; subplot 321;
imagesc(time,bin(bins),u(bins,:),[-1 1]); imagesc(time,bin(bins),u(bins,:),[-1 1]);
set(gca,'YTick',bins); set(gca,'YTick',bins);
ylabel('Bins');
hold on; hold on;
gregtick; gregtick;
title('Meridional velocity before correction'); title('Meridional velocity before correction');
...@@ -15,6 +16,7 @@ subplot 323; ...@@ -15,6 +16,7 @@ subplot 323;
imagesc(time,bin(bins),u1(bins,:),[-1 1]); imagesc(time,bin(bins),u1(bins,:),[-1 1]);
set(gca,'YTick',bins); set(gca,'YTick',bins);
hold on; hold on;
ylabel('Bins');
gregtick; gregtick;
title('Meridional velocity after correction'); title('Meridional velocity after correction');
...@@ -22,13 +24,15 @@ subplot 322; ...@@ -22,13 +24,15 @@ subplot 322;
imagesc(time,bin(bins),v(bins,:),[-1 1]); imagesc(time,bin(bins),v(bins,:),[-1 1]);
set(gca,'YTick',bins); set(gca,'YTick',bins);
hold on; hold on;
ylabel('Bins');
gregtick; gregtick;
title('Zonal velocity after correction'); title('Zonal velocity before correction');
subplot 324; subplot 324;
imagesc(time,bin(bins),v1(bins,:),[-1 1]); imagesc(time,bin(bins),v1(bins,:),[-1 1]);
set(gca,'YTick',bins); set(gca,'YTick',bins);
hold on; hold on;
ylabel('Bins');
gregtick; gregtick;
title('Zonal velocity after correction'); title('Zonal velocity after correction');
...@@ -39,6 +43,7 @@ for i=1:length(bins); ...@@ -39,6 +43,7 @@ for i=1:length(bins);
end end
set(gca,'XLim',[time(1) time(end)]); set(gca,'XLim',[time(1) time(end)]);
gregtick; gregtick;
ylabel('Bins');
title('Bin depth after correction'); title('Bin depth after correction');
subplot 326; subplot 326;
...@@ -48,4 +53,5 @@ for i=1:length(bins); ...@@ -48,4 +53,5 @@ for i=1:length(bins);
end end
set(gca,'XLim',[time(1) time(end)]); set(gca,'XLim',[time(1) time(end)]);
gregtick; gregtick;
ylabel('Bins');
title('Bin depth after correction'); title('Bin depth after correction');
...@@ -63,30 +63,30 @@ for iidep = 1:length(intdepvec) ...@@ -63,30 +63,30 @@ for iidep = 1:length(intdepvec)
end end
hf=figure('position', [0, 0, 1400, 1000]);
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');
subplot 121; subplot 121;
imagesc(data.time,intdepvec,ui,[-1 1]); 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; gregtick;
title('data raw'); ylabel('Bins');
title('U field - data raw');
figure;
subplot 122; subplot 122;
imagesc(data.time,intdepvec,vintfilt,[-1 1]); imagesc(data.time,intdepvec,uintfilt,[-1 1]);
set(gca,'YLim',[min(intdepvec)-10 max(intdepvec)+10]); %set(gca,'YLim',[min(intdepvec)-10 max(intdepvec)+10]);
gregtick; 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; subplot 121;
imagesc(data.time,intdepvec,vi,[-1 1]); 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; gregtick;
title('data raw'); ylabel('Bins');
title('V field - data interpolated, filtered and subsampled');
...@@ -46,6 +46,8 @@ function [zbins,zadcp1,offset,x_null]=adcp_surface_fit(zadcp,ea,surface_bins,ble ...@@ -46,6 +46,8 @@ function [zbins,zadcp1,offset,x_null]=adcp_surface_fit(zadcp,ea,surface_bins,ble
figure(1); figure(1);
bar(count,ncount); bar(count,ncount);
title('Histogram of differences between original and reconstructed depth record'); title('Histogram of differences between original and reconstructed depth record');
ylabel('Frequency');
xlabel('Difference of Depth (m)');
figure(2); figure(2);
plot(zadcp,'b'); plot(zadcp,'b');
...@@ -75,12 +77,18 @@ function [zbins,zadcp1,offset,x_null]=adcp_surface_fit(zadcp,ea,surface_bins,ble ...@@ -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']); text(300, max(zadcp),['Offset applied: ' num2str(offset) ' m']);
%,'fonts',12,'fontw','bold','backgroundc','w'); %,'fonts',12,'fontw','bold','backgroundc','w');
legend('Original','Reconstructed from surface reflection','Offset applied'); legend('Original','Reconstructed from surface reflection','Offset applied');
ylabel('Depth (m)');
xlabel('Time index');
%print -dpng surface_fit; %print -dpng surface_fit;
figure(3); figure(3);
pcolor([1:length(x_null)],zbins,ea); shading flat; 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; %print -dpng surface_ea;
end end
\ No newline at end of file
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment