Comparison_WERA-ADCP_Sept2013.html 8.52 KB
Newer Older
Marion Bandet's avatar
update  
Marion Bandet committed
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 27
<!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" />
    <link rel="next" title="Current combination and OGSL data feed" href="CurrentsCombination.html" />
Marion Bandet's avatar
Marion Bandet committed
28
    <link rel="prev" title="Outputs of WERAcalibration.m for PAB site" href="calibPAB.html" /> 
Marion Bandet's avatar
update  
Marion Bandet committed
29 30 31 32 33 34 35 36 37 38 39 40
  </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" >
          <a href="CurrentsCombination.html" title="Current combination and OGSL data feed"
             accesskey="N">next</a> |</li>
        <li class="right" >
Marion Bandet's avatar
Marion Bandet committed
41
          <a href="calibPAB.html" title="Outputs of WERAcalibration.m for PAB site"
Marion Bandet's avatar
update  
Marion Bandet committed
42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115
             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>
<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>
<div><p>We look for the indices of the ADCP time vector corresponding to the time period of interest (here 01 to 30 September 2013).</p>
<p>For each index, we look for the ADCP position (xadcp,yadcp)(via the gps coordinates).</p>
<p>We then look for the positions 15min before and 15min after.</p>
<p>We estimate the mean drift of the ADCP within these 30 minutes (Umean_adcp_u, Umean_adcp_v and U_mean_adcp (total)).</p>
<p>We find the radial orientation at the ADCP position given by xadcp,yadcp: azimuth in degrees from HFR site to adcp position</p>
<p>We project U_mean_adcp onto the radial direction to get the radial current components (ur_adcp_proj, vr_adcp_proj).</p>
<p>We store all the meaningful parameters into a structure named adcp.</p>
</div></blockquote>
<p>7- We analyze the HFR data:</p>
<blockquote>
<div><p>We look for the indices of the HFR time vector corresponding to the time period of interest (here 01 to 30 September 2013)</p>
<p>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).</p>
</div></blockquote>
<div class="figure align-center">
<img alt="_images/Radial_ADCPbin2_HFR_PAO_version5_Sept2013.png" src="_images/Radial_ADCPbin2_HFR_PAO_version5_Sept2013.png" />
<p class="caption">Hourly radial currents from ADCP bin 2 and HFR WERA at PAO - September 2013</p>
<div class="legend">
Superimposed are the 48-h averages</div>
</div>
<p>[[Fichier:gridwithnan_adcp.png]]</p>
<blockquote>
<div><p>We find the HFR Radial velocity at this position and project it onto zonal and meridional components.</p>
<p>We store all the meaningful parameters into a structure named radar.</p>
</div></blockquote>
<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>
<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>
<p>[[Fichier:rangeADCP.png]]</p>
<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>
<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>
<p>[[Fichier:CompsNewCalib_RadialVelo_ADCP_PAO_Sept2013.png|200px|thumb|left|Comparison of radial velocities from ADCP bin 2 and HFR WERA at PAO site September 2013]]</p>
<p>[[Fichier:orientationPAOSept2013.png|200px|thumb|left|RMS diff and correlation from radial velocities compared betwen ADCP bin 2 and HFR WERA at PAO site September 2013]]</p>
<p>[[Fichier:ADCPbin2vsPAOSept2013.png|200px|thumb|left|radial velocities from ADCP bin 2 vs HFR WERA at PAO site September 2013 - in red: 48h averages]]]]</p>
</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>
Marion Bandet's avatar
Marion Bandet committed
116 117
  <p class="topless"><a href="calibPAB.html"
                        title="previous chapter">Outputs of WERAcalibration.m for PAB site</a></p>
Marion Bandet's avatar
update  
Marion Bandet committed
118 119 120 121 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
  <h4>Next topic</h4>
  <p class="topless"><a href="CurrentsCombination.html"
                        title="next chapter">Current combination and OGSL data feed</a></p>
  <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" >
          <a href="CurrentsCombination.html" title="Current combination and OGSL data feed"
             >next</a> |</li>
        <li class="right" >
Marion Bandet's avatar
Marion Bandet committed
153
          <a href="calibPAB.html" title="Outputs of WERAcalibration.m for PAB site"
Marion Bandet's avatar
update  
Marion Bandet committed
154 155 156 157 158 159 160 161 162 163
             >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>