Blame view

build/html/Comparison_WERA-ADCP_Sept2013.html 11.6 KB
d001201d   Marion Bandet   update
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN"
  "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">


<html xmlns="http://www.w3.org/1999/xhtml">
  <head>
    <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
    
    <title>Comparison WERA and ADCP radial velocities September 2013 &mdash; documentation_HFR_analysis 0.0.1 documentation</title>
    
    <link rel="stylesheet" href="_static/default.css" type="text/css" />
    <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
    
    <script type="text/javascript">
      var DOCUMENTATION_OPTIONS = {
        URL_ROOT:    './',
        VERSION:     '0.0.1',
        COLLAPSE_INDEX: false,
        FILE_SUFFIX: '.html',
        HAS_SOURCE:  true
      };
    </script>
    <script type="text/javascript" src="_static/jquery.js"></script>
    <script type="text/javascript" src="_static/underscore.js"></script>
    <script type="text/javascript" src="_static/doctools.js"></script>
    <link rel="top" title="documentation_HFR_analysis 0.0.1 documentation" href="index.html" />
ee74eb7f   Marion Bandet   update fig ADCP V...
27
    <link rel="next" title="Comparison WERA and ADCP radial velocities May-June 2013" href="Comparison_WERA-ADCP_MayJune2013.html" />
64977cd9   Marion Bandet   check Rx antennas...
28
    <link rel="prev" title="Outputs of WERAcalibration.m for PAB site" href="calibPAB.html" /> 
d001201d   Marion Bandet   update
29
30
31
32
33
34
35
36
37
  </head>
  <body>
    <div class="related">
      <h3>Navigation</h3>
      <ul>
        <li class="right" style="margin-right: 10px">
          <a href="genindex.html" title="General Index"
             accesskey="I">index</a></li>
        <li class="right" >
ee74eb7f   Marion Bandet   update fig ADCP V...
38
          <a href="Comparison_WERA-ADCP_MayJune2013.html" title="Comparison WERA and ADCP radial velocities May-June 2013"
d001201d   Marion Bandet   update
39
40
             accesskey="N">next</a> |</li>
        <li class="right" >
64977cd9   Marion Bandet   check Rx antennas...
41
          <a href="calibPAB.html" title="Outputs of WERAcalibration.m for PAB site"
d001201d   Marion Bandet   update
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
             accesskey="P">previous</a> |</li>
        <li><a href="index.html">documentation_HFR_analysis 0.0.1 documentation</a> &raquo;</li> 
      </ul>
    </div>  

    <div class="document">
      <div class="documentwrapper">
        <div class="bodywrapper">
          <div class="body">
            
  <div class="section" id="comparison-wera-and-adcp-radial-velocities-september-2013">
<h1>Comparison WERA and ADCP radial velocities September 2013<a class="headerlink" href="#comparison-wera-and-adcp-radial-velocities-september-2013" title="Permalink to this headline">¶</a></h1>
<p>Matlab script used to do the analysis is <strong>RadialVeloCompa_Sept2013_HFR_ADCP_WERA_ADCPcorrected_v2.m</strong></p>
<p>Steps of the algorithm:</p>
<p>1- we choose a site to analyze, either PAO (site 1) or PAB (site 2)</p>
2b149458   Marion Bandet   mise a jour
57
58
<div class="figure align-center">
<img alt="_images/sites.png" src="_images/sites.png" />
54645489   Marion Bandet   update HFR direction
59
<p class="caption"><em>location of HFRs domain in the St Lawrence estuary</em></p>
2b149458   Marion Bandet   mise a jour
60
</div>
d001201d   Marion Bandet   update
61
62
63
64
65
66
<p>2- we load the bathymetry and the HFR sites&#8217;coordinates</p>
<p>3- we specify if we use the &#8220;old&#8221; or &#8220;new&#8221; calibrations. Usually &#8220;New&#8221;.</p>
<p>4- we load IML-4 ADCP data (2013) and gps data</p>
<p>5- we make a list of all the HFR files available for the site of interest.</p>
<p>6- we analyze the ADCP data (bin 2 only to get near-surface data):</p>
<blockquote>
2b149458   Marion Bandet   mise a jour
67
68
69
70
71
72
73
74
75
<div><ul class="simple">
<li>We look for the indices of the ADCP time vector corresponding to the time period of interest (here 01 to 30 September 2013).</li>
<li>For each index, we look for the ADCP position (xadcp,yadcp)(via the gps coordinates).</li>
<li>We then look for the positions 15min before and 15min after.</li>
<li>We estimate the mean drift of the ADCP within these 30 minutes (Umean_adcp_u, Umean_adcp_v and U_mean_adcp (total)).</li>
<li>We find the radial orientation at the ADCP position given by xadcp,yadcp: azimuth in degrees from HFR site to adcp position</li>
<li>We project U_mean_adcp onto the radial direction to get the radial current components (ur_adcp_proj, vr_adcp_proj).</li>
<li>We store all the meaningful parameters into a structure named adcp.</li>
</ul>
d001201d   Marion Bandet   update
76
77
78
</div></blockquote>
<p>7- We analyze the HFR data:</p>
<blockquote>
2b149458   Marion Bandet   mise a jour
79
80
81
82
<div><ul class="simple">
<li>We look for the indices of the HFR time vector corresponding to the time period of interest (here 01 to 30 September 2013)</li>
<li>For each index, we load the corresponding HFR file (i.e. 20132440000_pao.mat). Each file contains a grid (X,Y) with radial currents at each grid point (Ur), a corresponding time vector (t) and the angle and distance from the site (theta,r).</li>
</ul>
d001201d   Marion Bandet   update
83
84
</div></blockquote>
<div class="figure align-center">
2b149458   Marion Bandet   mise a jour
85
86
87
<img alt="_images/explecoveragegridallpoint.png" src="_images/explecoveragegridallpoint.png" />
</div>
<div class="figure align-center">
d001201d   Marion Bandet   update
88
<img alt="_images/Radial_ADCPbin2_HFR_PAO_version5_Sept2013.png" src="_images/Radial_ADCPbin2_HFR_PAO_version5_Sept2013.png" />
54645489   Marion Bandet   update HFR direction
89
<p class="caption"><em>Hourly radial currents from ADCP bin 2 and HFR WERA at PAO - September 2013 - superimposed are the 48-h averages</em></p>
d001201d   Marion Bandet   update
90
<div class="legend">
54645489   Marion Bandet   update HFR direction
91
92
93
94
<ul class="simple">
<li>We then look for the HFR grid point nearest the ADCP position at the time of interest.</li>
</ul>
</div>
d001201d   Marion Bandet   update
95
</div>
2b149458   Marion Bandet   mise a jour
96
97
<div class="figure align-center">
<img alt="_images/explecoveragegrid2.png" src="_images/explecoveragegrid2.png" />
54645489   Marion Bandet   update HFR direction
98
<p class="caption"><em>example of grid coverage</em></p>
2b149458   Marion Bandet   mise a jour
99
<div class="legend">
54645489   Marion Bandet   update HFR direction
100
101
102
103
104
<ul class="simple">
<li>We find the HFR Radial velocity at this position and project it onto zonal and meridional components.</li>
<li>We store all the meaningful parameters into a structure named radar.</li>
</ul>
</div>
2b149458   Marion Bandet   mise a jour
105
</div>
d001201d   Marion Bandet   update
106
107
108
109
110
<p>This analysis allows us to compare radial velocity components from each HFR site with the ADCP components.</p>
</div>
<div class="section" id="hfr-direction">
<h1>HFR direction<a class="headerlink" href="#hfr-direction" title="Permalink to this headline">¶</a></h1>
<p>To verify that the HFR looks in the right direction:</p>
54645489   Marion Bandet   update HFR direction
111
<p># For each time step of the time period of interest, we compute the distance from the HFR site to each grid point. We select the grid points that are at the same range from the HFR site than the grid point closest to the ADCP.</p>
2b149458   Marion Bandet   mise a jour
112
113
114
115
116
117
<div class="figure align-center">
<img alt="_images/explecoveragegrid_rangeADCP.png" src="_images/explecoveragegrid_rangeADCP.png" />
</div>
<div class="figure align-center">
<img alt="_images/explecoveragegrid_rangeADCP2.png" src="_images/explecoveragegrid_rangeADCP2.png" />
</div>
54645489   Marion Bandet   update HFR direction
118
119
<p># For each of these grid point (in cyan in the figure), we compute the azimuth the the site coordinates and we extract the radial velocity (total) and its components (zonal and meridional).</p>
<p>#We store all the meaningful parameters into a structure named range.</p>
d001201d   Marion Bandet   update
120
121
<p>Results from the whole analysis are stored, for example, under:
/run/media/bandma01/SeagateWireless/MEOPAR/Matlab/Analysis/Radial_Comparison/201309/ADCP-PAO/ADCP_PAO_version5adcpCorrected_newcalibrations_v2.mat</p>
54645489   Marion Bandet   update HFR direction
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
<p><strong>RESULTS WITH OLD CALIBRATIONS</strong></p>
<p>generated with PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2.m</p>
<div class="figure align-center">
<img alt="_images/PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2_oldcalib_01.png" src="_images/PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2_oldcalib_01.png" />
<p class="caption"><em>Comparison of radial velocities from ADCP bin 2 and HFR WERA at PAO site September 2013</em></p>
</div>
<div class="figure align-center">
<img alt="_images/PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2_oldcalib_02.png" src="_images/PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2_oldcalib_02.png" />
</div>
<div class="figure align-center">
<img alt="_images/PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2_oldcalib_03.png" src="_images/PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2_oldcalib_03.png" />
</div>
<div class="figure align-center">
<img alt="_images/PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2_oldcalib_04.png" src="_images/PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2_oldcalib_04.png" />
</div>
<div class="figure align-center">
<img alt="_images/Plot_WERA_direction_Sept2013_ADCP_oldcalib_03.png" src="_images/Plot_WERA_direction_Sept2013_ADCP_oldcalib_03.png" />
</div>
<p><strong>RESULTS WITH NEW CALIBRATIONS</strong></p>
<p>generated with PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2.m</p>
<div class="figure align-center">
<img alt="_images/PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2_newcalib_01.png" src="_images/PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2_newcalib_01.png" />
</div>
<div class="figure align-center">
<img alt="_images/PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2_newcalib_02.png" src="_images/PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2_newcalib_02.png" />
</div>
<div class="figure align-center">
<img alt="_images/PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2_newcalib_03.png" src="_images/PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2_newcalib_03.png" />
</div>
<div class="figure align-center">
<img alt="_images/PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2_newcalib_04.png" src="_images/PlotRadialVeloCompa_Sept2013_HFR_ADCP_WERA_v2_newcalib_04.png" />
</div>
<div class="figure align-center">
<img alt="_images/Plot_WERA_direction_Sept2013_ADCP_newcalib_03.png" src="_images/Plot_WERA_direction_Sept2013_ADCP_newcalib_03.png" />
<p class="caption"><em>RMS diff and correlation from radial velocities compared betwen 1/ ADCP bin 2 and HFR WERA at PAO site September 2013 (HFR grid point closest to ADCP location) and 2/ using model data at ADCP location and at HFR grid point closest to ADCP</em></p>
</div>
d001201d   Marion Bandet   update
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
</div>


          </div>
        </div>
      </div>
      <div class="sphinxsidebar">
        <div class="sphinxsidebarwrapper">
  <h3><a href="index.html">Table Of Contents</a></h3>
  <ul>
<li><a class="reference internal" href="#">Comparison WERA and ADCP radial velocities September 2013</a></li>
<li><a class="reference internal" href="#hfr-direction">HFR direction</a></li>
</ul>

  <h4>Previous topic</h4>
64977cd9   Marion Bandet   check Rx antennas...
173
174
  <p class="topless"><a href="calibPAB.html"
                        title="previous chapter">Outputs of WERAcalibration.m for PAB site</a></p>
d001201d   Marion Bandet   update
175
  <h4>Next topic</h4>
ee74eb7f   Marion Bandet   update fig ADCP V...
176
177
  <p class="topless"><a href="Comparison_WERA-ADCP_MayJune2013.html"
                        title="next chapter">Comparison WERA and ADCP radial velocities May-June 2013</a></p>
d001201d   Marion Bandet   update
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
  <h3>This Page</h3>
  <ul class="this-page-menu">
    <li><a href="_sources/Comparison_WERA-ADCP_Sept2013.txt"
           rel="nofollow">Show Source</a></li>
  </ul>
<div id="searchbox" style="display: none">
  <h3>Quick search</h3>
    <form class="search" action="search.html" method="get">
      <input type="text" name="q" />
      <input type="submit" value="Go" />
      <input type="hidden" name="check_keywords" value="yes" />
      <input type="hidden" name="area" value="default" />
    </form>
    <p class="searchtip" style="font-size: 90%">
    Enter search terms or a module, class or function name.
    </p>
</div>
<script type="text/javascript">$('#searchbox').show(0);</script>
        </div>
      </div>
      <div class="clearer"></div>
    </div>
    <div class="related">
      <h3>Navigation</h3>
      <ul>
        <li class="right" style="margin-right: 10px">
          <a href="genindex.html" title="General Index"
             >index</a></li>
        <li class="right" >
ee74eb7f   Marion Bandet   update fig ADCP V...
207
          <a href="Comparison_WERA-ADCP_MayJune2013.html" title="Comparison WERA and ADCP radial velocities May-June 2013"
d001201d   Marion Bandet   update
208
209
             >next</a> |</li>
        <li class="right" >
64977cd9   Marion Bandet   check Rx antennas...
210
          <a href="calibPAB.html" title="Outputs of WERAcalibration.m for PAB site"
d001201d   Marion Bandet   update
211
212
213
214
215
216
217
218
219
220
             >previous</a> |</li>
        <li><a href="index.html">documentation_HFR_analysis 0.0.1 documentation</a> &raquo;</li> 
      </ul>
    </div>
    <div class="footer">
        &copy; Copyright 2015, Marion Bandet.
      Created using <a href="http://sphinx-doc.org/">Sphinx</a> 1.2.3.
    </div>
  </body>
</html>