Commit 7e281bfb2e8c5f5a4b1f151365e05395d706d498

Authored by Marion Bandet
1 parent 26f5f7e9
Exists in master

update

Showing 116 changed files with 807 additions and 220 deletions   Show diff stats
build/doctrees/20140506_ADCP_V50.doctree
No preview for this file type
build/doctrees/calibPAB2014.doctree 0 → 100644
No preview for this file type
build/doctrees/calibPAO2014.doctree 0 → 100644
No preview for this file type
build/doctrees/calibration.doctree
No preview for this file type
build/doctrees/environment.pickle
No preview for this file type
build/html/20140506_ADCP_V50.html
... ... @@ -121,6 +121,15 @@
121 121 <img alt="_images/ADCP_V50_May2014_19.png" src="_images/ADCP_V50_May2014_19.png" />
122 122 </div>
123 123 </div>
  124 +<div class="section" id="angle-between-eastward-and-northward-velocity-components">
  125 +<h1>angle between eastward and northward velocity components<a class="headerlink" href="#angle-between-eastward-and-northward-velocity-components" title="Permalink to this headline">¶</a></h1>
  126 +<div class="figure align-center">
  127 +<img alt="_images/ADCP_V50_May2014_20.png" src="_images/ADCP_V50_May2014_20.png" />
  128 +</div>
  129 +<div class="figure align-center">
  130 +<img alt="_images/ADCP_V50_May2014_21.png" src="_images/ADCP_V50_May2014_21.png" />
  131 +</div>
  132 +</div>
124 133  
125 134  
126 135 </div>
... ... @@ -140,6 +149,7 @@
140 149 <li><a class="reference internal" href="#averaged-intensity">Averaged Intensity</a></li>
141 150 <li><a class="reference internal" href="#averaged-correlation">Averaged Correlation</a></li>
142 151 <li><a class="reference internal" href="#vertical-profiles-of-averaged-velocity-components">Vertical Profiles of averaged velocity components</a></li>
  152 +<li><a class="reference internal" href="#angle-between-eastward-and-northward-velocity-components">angle between eastward and northward velocity components</a></li>
143 153 </ul>
144 154  
145 155 <h3>This Page</h3>
... ...
build/html/ADCP_V50_May2014.html
... ... @@ -6,7 +6,7 @@
6 6 <!--
7 7 This HTML was auto-generated from MATLAB code.
8 8 To make changes, update the MATLAB code and republish this document.
9   - --><title>ADCP_V50_May2014</title><meta name="generator" content="MATLAB 7.14"><link rel="schema.DC" href="http://purl.org/dc/elements/1.1/"><meta name="DC.date" content="2016-07-12"><meta name="DC.source" content="ADCP_V50_May2014.m"><style type="text/css">
  9 + --><title>ADCP_V50_May2014</title><meta name="generator" content="MATLAB 7.14"><link rel="schema.DC" href="http://purl.org/dc/elements/1.1/"><meta name="DC.date" content="2016-07-18"><meta name="DC.source" content="ADCP_V50_May2014.m"><style type="text/css">
10 10 html,body,div,span,applet,object,iframe,h1,h2,h3,h4,h5,h6,p,blockquote,pre,a,abbr,acronym,address,big,cite,code,del,dfn,em,font,img,ins,kbd,q,s,samp,small,strike,strong,sub,sup,tt,var,b,u,i,center,dl,dt,dd,ol,ul,li,fieldset,form,label,legend,table,caption,tbody,tfoot,thead,tr,th,td{margin:0;padding:0;border:0;outline:0;font-size:100%;vertical-align:baseline;background:transparent}body{line-height:1}ol,ul{list-style:none}blockquote,q{quotes:none}blockquote:before,blockquote:after,q:before,q:after{content:'';content:none}:focus{outine:0}ins{text-decoration:none}del{text-decoration:line-through}table{border-collapse:collapse;border-spacing:0}
11 11  
12 12 html { min-height:100%; margin-bottom:1px; }
... ... @@ -55,7 +55,7 @@ span.syscmd { color:#B28C00 }
55 55 .footer { width:auto; padding:10px 0px; margin:25px 0px 0px; border-top:1px dotted #878787; font-size:0.8em; line-height:140%; font-style:italic; color:#878787; text-align:left; float:none; }
56 56 .footer p { margin:0px; }
57 57  
58   - </style></head><body><div class="content"><h2>Contents</h2><div><ul><li><a href="#2">roll - pitch - heading</a></li><li><a href="#3">orientation - voltage - salinity - pressure sensor depth - temperature -</a></li><li><a href="#4">speed of sound</a></li><li><a href="#5">cell depth</a></li><li><a href="#6">range</a></li><li><a href="#7">PG</a></li><li><a href="#8">CORR</a></li><li><a href="#9">intensity</a></li><li><a href="#10">velocity ENUErr</a></li><li><a href="#11">averaged velocity</a></li><li><a href="#12">intensity</a></li><li><a href="#13">correlation</a></li><li><a href="#14">Vertical profile of averaged velocity</a></li><li><a href="#15">make movie of u and v vs depth</a></li></ul></div><pre class="codeinput"><span class="comment">% ADCP-V50</span>
  58 + </style></head><body><div class="content"><h2>Contents</h2><div><ul><li><a href="#2">roll - pitch - heading</a></li><li><a href="#3">orientation - voltage - salinity - pressure sensor depth - temperature -</a></li><li><a href="#4">speed of sound</a></li><li><a href="#5">cell depth</a></li><li><a href="#6">range</a></li><li><a href="#7">PG</a></li><li><a href="#8">CORR</a></li><li><a href="#9">intensity</a></li><li><a href="#10">velocity ENUErr</a></li><li><a href="#11">averaged velocity</a></li><li><a href="#12">intensity</a></li><li><a href="#13">correlation</a></li><li><a href="#14">Vertical profile of averaged velocity</a></li><li><a href="#15">plot of velocity vector direction vs depth</a></li></ul></div><pre class="codeinput"><span class="comment">% ADCP-V50</span>
59 59 <span class="comment">% first import data (.pd0) using Velocity software dated december 2013 on</span>
60 60 <span class="comment">% windows 8, then export as .mat</span>
61 61  
... ... @@ -390,6 +390,7 @@ grid &lt;span class=&quot;string&quot;&gt;on&lt;/span&gt;
390 390 title(<span class="string">'Averaged velocity (m/s)'</span>)
391 391 ylabel(<span class="string">'depth (m)'</span>)
392 392 ylim([0 20])
  393 +xlim(0.04*[-1 1])
393 394  
394 395 <span class="comment">% compute number of data point</span>
395 396 Nb_pt = sum(~isnan(vb.vel(I,:)));
... ... @@ -416,6 +417,7 @@ grid &lt;span class=&quot;string&quot;&gt;on&lt;/span&gt;
416 417 title(<span class="string">'East velocity (m/s)'</span>)
417 418 ylabel(<span class="string">'depth (m)'</span>)
418 419 ylim([0 20])
  420 +<span class="comment">%xlim(0.3*[-1 1])</span>
419 421  
420 422 subplot(122)
421 423 plot(Nb_pt,depth,<span class="string">'-'</span>)
... ... @@ -471,49 +473,86 @@ set(gca,&lt;span class=&quot;string&quot;&gt;&#39;Ydir&#39;&lt;/span&gt;,&lt;span class=&quot;string&quot;&gt;&#39;reverse&#39;&lt;/span&gt;
471 473 grid <span class="string">on</span>
472 474 title(<span class="string">'Number of data points'</span>)
473 475 ylim([0 20])
  476 +</pre><img vspace="5" hspace="5" src="ADCP_V50_May2014_16.png" alt=""> <img vspace="5" hspace="5" src="ADCP_V50_May2014_17.png" alt=""> <img vspace="5" hspace="5" src="ADCP_V50_May2014_18.png" alt=""> <img vspace="5" hspace="5" src="ADCP_V50_May2014_19.png" alt=""> <h2>plot of velocity vector direction vs depth<a name="15"></a></h2><pre class="codeinput"><span class="comment">%close all</span>
474 477  
475 478  
476   -<span class="keyword">return</span>
477   -</pre><img vspace="5" hspace="5" src="ADCP_V50_May2014_16.png" alt=""> <img vspace="5" hspace="5" src="ADCP_V50_May2014_17.png" alt=""> <img vspace="5" hspace="5" src="ADCP_V50_May2014_18.png" alt=""> <img vspace="5" hspace="5" src="ADCP_V50_May2014_19.png" alt=""> <h2>make movie of u and v vs depth<a name="15"></a></h2><pre class="codeinput"> u = squeeze(wt.vel(I,:,1));
478   - v = squeeze(wt.vel(I,:,2));
479   - w = squeeze(wt.vel(I,:,3));
480   - x = zeros(1,80);
481   - y = zeros(1,80);
482   - z = depth;
  479 +<span class="comment">% time averages</span>
  480 +u_mean = nanmean(squeeze(wt.vel(I,:,1)));
  481 +v_mean = nanmean(squeeze(wt.vel(I,:,2)));
  482 +w_mean = nanmean(squeeze(wt.vel(I,:,3)));
  483 +x = zeros(1,80);
  484 +y = zeros(1,80);
  485 +z = depth;
  486 +
  487 +
  488 +
  489 +clear <span class="string">I</span>
  490 +I = find(depth &lt;= 13); <span class="comment">% limit the depth to 12 m</span>
  491 +
  492 +sc = 1; <span class="comment">% scaling factor to plot vectors</span>
483 493  
484 494 figure
485   -plot3(x,y,z,<span class="string">'-'</span>)
  495 +plot3(x(I),y(I),z(I),<span class="string">'.-'</span>)
486 496 grid <span class="string">on</span>
487   -set(gca,<span class="string">'Zdir'</span>,<span class="string">'reverse'</span>)
488 497 hold <span class="string">on</span>
  498 +plt = quiver3(x(I),y(I),z(I),avg_northvel(I)*100/sc,avg_eastvel(I)*100/sc,avg_upvel(I)/sc,0);
  499 +set(plt,<span class="string">'color'</span>,<span class="string">'b'</span>)
  500 +xlabel(<span class="string">'northward velocity [cm/s]'</span>)
  501 +ylabel(<span class="string">'eastward velocity [cm/s]'</span>)
  502 +zlabel(<span class="string">'depth [m]'</span>)
  503 +zlim([0 14])
  504 +set(gca,<span class="string">'Ydir'</span>,<span class="string">'reverse'</span>)
  505 +set(gca,<span class="string">'Zdir'</span>,<span class="string">'reverse'</span>)
489 506  
490   -<span class="keyword">for</span> i = 1:200<span class="comment">%:length(u)</span>
491   - i
492   - plt = quiver3(x,y,z,u(i,:)/10,v(i,:)/10,w(i,:)/10,0);
493   - set(plt,<span class="string">'color'</span>,<span class="string">'b'</span>)
494   - xlim(0.025*[-1 1])
495   - ylim(0.028*[-1 1])
496   - <span class="comment">% grid on</span>
497   - <span class="comment">% set(gca,'Zdir','reverse')</span>
498   - pause(1)
499   - delete(plt)
500   -<span class="keyword">end</span>
501   -
502   -
503   -
504   -<span class="comment">% save figures</span>
505   -title_fig = {<span class="string">'RollPitchHeading'</span>;<span class="string">'miscel'</span>;<span class="string">'roll'</span>;<span class="string">'pitch'</span>;<span class="string">'heading'</span>;<span class="string">'PressureSensorDepth'</span>;<span class="string">'depthcell'</span>;<span class="string">'range'</span>;<span class="string">'beam_PG'</span>;<span class="string">'beam_correlation'</span>;<span class="string">'beam_intensity'</span>;<span class="string">'beam_velocity'</span>;<span class="string">'Averaged_velocity'</span>;<span class="string">'Averaged_intensity'</span>;<span class="string">'Averaged_correlation'</span>}
506   -fig_w = 15;<span class="comment">%cm</span>
507   -fig_h = 10;<span class="comment">%cm</span>
508   -
509   -<span class="keyword">for</span> i = 1:15
510   - fig_hdl = figure(i);
511   - fname = [title_fig{i} <span class="string">'_'</span> filename(8:22)];
512   - set(fig_hdl,<span class="string">'PaperUnits'</span>,<span class="string">'centimeters'</span>)
513   - set(fig_hdl,<span class="string">'PaperPosition'</span>,[5 5 fig_w fig_h])
514   - print(fig_hdl,<span class="string">'-dpng'</span>,<span class="string">'-r400'</span>,fname)
515   -<span class="keyword">end</span>
516   -</pre><p class="footer"><br>
  507 +view(-116,30)
  508 +
  509 +<span class="comment">% figure</span>
  510 +<span class="comment">% quiver(x(I),z(I),avg_eastvel(I)/sc,avg_northvel(I)/sc,0)</span>
  511 +<span class="comment">% set(gca,'Ydir','reverse')</span>
  512 +<span class="comment">% grid on</span>
  513 +<span class="comment">% xlabel('eastward velocity')</span>
  514 +<span class="comment">% %ylabel('upward velocity')</span>
  515 +<span class="comment">% ylabel('depth [m]')</span>
  516 +
  517 +<span class="comment">% figure</span>
  518 +<span class="comment">% quiver(y(I),z(I),avg_northvel(I)/sc,avg_upvel(I)/sc,0)</span>
  519 +<span class="comment">% set(gca,'Ydir','reverse')</span>
  520 +<span class="comment">% grid on</span>
  521 +<span class="comment">% xlabel('northward velocity')</span>
  522 +<span class="comment">% %ylabel('upward velocity')</span>
  523 +<span class="comment">% ylabel('depth [m]')</span>
  524 +
  525 +<span class="comment">% angle</span>
  526 +theta = atan2(avg_northvel,avg_eastvel);
  527 +
  528 +theta = rad2deg(theta);
  529 +theta = mod(theta,360);
  530 +
  531 +
  532 +figure
  533 +plot(theta(I),depth(I),<span class="string">'.-'</span>)
  534 +grid <span class="string">on</span>
  535 +set(gca,<span class="string">'Ydir'</span>,<span class="string">'reverse'</span>)
  536 +xlabel(<span class="string">'angle btw eastward and northward velocity components wrt True North [deg]'</span>)
  537 +ylabel(<span class="string">'depth [m]'</span>)
  538 +xlim([200 250])
  539 +
  540 +
  541 +<span class="comment">% % save figures</span>
  542 +<span class="comment">% title_fig = {'RollPitchHeading';'miscel';'roll';'pitch';'heading';'PressureSensorDepth';'depthcell';'range';'beam_PG';'beam_correlation';'beam_intensity';'beam_velocity';'Averaged_velocity';'Averaged_intensity';'Averaged_correlation'}</span>
  543 +<span class="comment">% fig_w = 15;%cm</span>
  544 +<span class="comment">% fig_h = 10;%cm</span>
  545 +<span class="comment">%</span>
  546 +<span class="comment">% for i = 1:15</span>
  547 +<span class="comment">% fig_hdl = figure(i);</span>
  548 +<span class="comment">% fname = [title_fig{i} '_' filename(8:22)];</span>
  549 +<span class="comment">% set(fig_hdl,'PaperUnits','centimeters')</span>
  550 +<span class="comment">% set(fig_hdl,'PaperPosition',[5 5 fig_w fig_h])</span>
  551 +<span class="comment">% print(fig_hdl,'-dpng','-r400',fname)</span>
  552 +<span class="comment">% end</span>
  553 +<span class="comment">%</span>
  554 +<span class="comment">%</span>
  555 +</pre><img vspace="5" hspace="5" src="ADCP_V50_May2014_20.png" alt=""> <img vspace="5" hspace="5" src="ADCP_V50_May2014_21.png" alt=""> <p class="footer"><br>
517 556 Published with MATLAB&reg; 7.14<br></p></div><!--
518 557 ##### SOURCE BEGIN #####
519 558 % ADCP-V50
... ... @@ -886,6 +925,7 @@ grid on
886 925 title('Averaged velocity (m/s)')
887 926 ylabel('depth (m)')
888 927 ylim([0 20])
  928 +xlim(0.04*[-1 1])
889 929  
890 930 % compute number of data point
891 931 Nb_pt = sum(~isnan(vb.vel(I,:)));
... ... @@ -912,6 +952,7 @@ grid on
912 952 title('East velocity (m/s)')
913 953 ylabel('depth (m)')
914 954 ylim([0 20])
  955 +%xlim(0.3*[-1 1])
915 956  
916 957 subplot(122)
917 958 plot(Nb_pt,depth,'-')
... ... @@ -969,51 +1010,88 @@ title(&#39;Number of data points&#39;)
969 1010 ylim([0 20])
970 1011  
971 1012  
972   -return
973 1013  
974 1014  
975   -%% make movie of u and v vs depth
976   - u = squeeze(wt.vel(I,:,1));
977   - v = squeeze(wt.vel(I,:,2));
978   - w = squeeze(wt.vel(I,:,3));
979   - x = zeros(1,80);
980   - y = zeros(1,80);
981   - z = depth;
982   -
  1015 +%% plot of velocity vector direction vs depth
  1016 +%close all
  1017 +
  1018 +
  1019 +% time averages
  1020 +u_mean = nanmean(squeeze(wt.vel(I,:,1)));
  1021 +v_mean = nanmean(squeeze(wt.vel(I,:,2)));
  1022 +w_mean = nanmean(squeeze(wt.vel(I,:,3)));
  1023 +x = zeros(1,80);
  1024 +y = zeros(1,80);
  1025 +z = depth;
  1026 +
  1027 +
  1028 +
  1029 +clear I
  1030 +I = find(depth <= 13); % limit the depth to 12 m
  1031 +
  1032 +sc = 1; % scaling factor to plot vectors
  1033 +
983 1034 figure
984   -plot3(x,y,z,'-')
  1035 +plot3(x(I),y(I),z(I),'.-')
985 1036 grid on
986   -set(gca,'Zdir','reverse')
987 1037 hold on
  1038 +plt = quiver3(x(I),y(I),z(I),avg_northvel(I)*100/sc,avg_eastvel(I)*100/sc,avg_upvel(I)/sc,0);
  1039 +set(plt,'color','b')
  1040 +xlabel('northward velocity [cm/s]')
  1041 +ylabel('eastward velocity [cm/s]')
  1042 +zlabel('depth [m]')
  1043 +zlim([0 14])
  1044 +set(gca,'Ydir','reverse')
  1045 +set(gca,'Zdir','reverse')
988 1046  
989   -for i = 1:200%:length(u)
990   - i
991   - plt = quiver3(x,y,z,u(i,:)/10,v(i,:)/10,w(i,:)/10,0);
992   - set(plt,'color','b')
993   - xlim(0.025*[-1 1])
994   - ylim(0.028*[-1 1])
995   - % grid on
996   - % set(gca,'Zdir','reverse')
997   - pause(1)
998   - delete(plt)
999   -end
  1047 +view(-116,30)
1000 1048  
  1049 +% figure
  1050 +% quiver(x(I),z(I),avg_eastvel(I)/sc,avg_northvel(I)/sc,0)
  1051 +% set(gca,'Ydir','reverse')
  1052 +% grid on
  1053 +% xlabel('eastward velocity')
  1054 +% %ylabel('upward velocity')
  1055 +% ylabel('depth [m]')
1001 1056  
  1057 +% figure
  1058 +% quiver(y(I),z(I),avg_northvel(I)/sc,avg_upvel(I)/sc,0)
  1059 +% set(gca,'Ydir','reverse')
  1060 +% grid on
  1061 +% xlabel('northward velocity')
  1062 +% %ylabel('upward velocity')
  1063 +% ylabel('depth [m]')
1002 1064  
1003   -% save figures
1004   -title_fig = {'RollPitchHeading';'miscel';'roll';'pitch';'heading';'PressureSensorDepth';'depthcell';'range';'beam_PG';'beam_correlation';'beam_intensity';'beam_velocity';'Averaged_velocity';'Averaged_intensity';'Averaged_correlation'}
1005   -fig_w = 15;%cm
1006   -fig_h = 10;%cm
  1065 +% angle
  1066 +theta = atan2(avg_northvel,avg_eastvel);
1007 1067  
1008   -for i = 1:15
1009   - fig_hdl = figure(i);
1010   - fname = [title_fig{i} '_' filename(8:22)];
1011   - set(fig_hdl,'PaperUnits','centimeters')
1012   - set(fig_hdl,'PaperPosition',[5 5 fig_w fig_h])
1013   - print(fig_hdl,'-dpng','-r400',fname)
1014   -end
  1068 +theta = rad2deg(theta);
  1069 +theta = mod(theta,360);
1015 1070  
1016 1071  
  1072 +figure
  1073 +plot(theta(I),depth(I),'.-')
  1074 +grid on
  1075 +set(gca,'Ydir','reverse')
  1076 +xlabel('angle btw eastward and northward velocity components wrt True North [deg]')
  1077 +ylabel('depth [m]')
  1078 +xlim([200 250])
  1079 +
  1080 +
  1081 +% % save figures
  1082 +% title_fig = {'RollPitchHeading';'miscel';'roll';'pitch';'heading';'PressureSensorDepth';'depthcell';'range';'beam_PG';'beam_correlation';'beam_intensity';'beam_velocity';'Averaged_velocity';'Averaged_intensity';'Averaged_correlation'}
  1083 +% fig_w = 15;%cm
  1084 +% fig_h = 10;%cm
  1085 +%
  1086 +% for i = 1:15
  1087 +% fig_hdl = figure(i);
  1088 +% fname = [title_fig{i} '_' filename(8:22)];
  1089 +% set(fig_hdl,'PaperUnits','centimeters')
  1090 +% set(fig_hdl,'PaperPosition',[5 5 fig_w fig_h])
  1091 +% print(fig_hdl,'-dpng','-r400',fname)
  1092 +% end
  1093 +%
  1094 +%
1017 1095  
1018 1096 ##### SOURCE END #####
1019 1097 --></body></html>
1020 1098 \ No newline at end of file
... ...
build/html/ADCP_V50_May2014.png

3.28 KB | W: | H:

2.67 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
build/html/ADCP_V50_May2014_01.png

9.09 KB | W: | H:

9.09 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
build/html/ADCP_V50_May2014_02.png

8.75 KB | W: | H:

8.75 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
build/html/ADCP_V50_May2014_03.png

8.59 KB | W: | H:

8.59 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
build/html/ADCP_V50_May2014_04.png

8.37 KB | W: | H:

8.37 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
build/html/ADCP_V50_May2014_05.png

8.84 KB | W: | H:

8.84 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
build/html/ADCP_V50_May2014_06.png

9.27 KB | W: | H:

9.27 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
build/html/ADCP_V50_May2014_07.png

2.98 KB | W: | H:

2.98 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
build/html/ADCP_V50_May2014_08.png

4.19 KB | W: | H:

4.19 KB | W: | H:

  • 2-up
  • Swipe
  • Onion skin
build/html/ADCP_V50_May2014_09.png

18.6 KB | W: | H:

18.6 KB | W: | H: