Commit d001201d5613043170ab35d516b11ca0d6abfe03

Authored by Marion Bandet
1 parent f1d33131
Exists in master

update

build/doctrees/Comparison_WERA-ADCP_Sept2013.doctree 0 → 100644
No preview for this file type
build/doctrees/CurrentsCombination.doctree
No preview for this file type
build/doctrees/WERAprocessing.doctree
No preview for this file type
build/doctrees/environment.pickle
No preview for this file type
build/doctrees/index.doctree
No preview for this file type
build/html/Comparison_WERA-ADCP_Sept2013.html 0 → 100644
... ... @@ -0,0 +1,163 @@
  1 +<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN"
  2 + "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
  3 +
  4 +
  5 +<html xmlns="http://www.w3.org/1999/xhtml">
  6 + <head>
  7 + <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
  8 +
  9 + <title>Comparison WERA and ADCP radial velocities September 2013 &mdash; documentation_HFR_analysis 0.0.1 documentation</title>
  10 +
  11 + <link rel="stylesheet" href="_static/default.css" type="text/css" />
  12 + <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
  13 +
  14 + <script type="text/javascript">
  15 + var DOCUMENTATION_OPTIONS = {
  16 + URL_ROOT: './',
  17 + VERSION: '0.0.1',
  18 + COLLAPSE_INDEX: false,
  19 + FILE_SUFFIX: '.html',
  20 + HAS_SOURCE: true
  21 + };
  22 + </script>
  23 + <script type="text/javascript" src="_static/jquery.js"></script>
  24 + <script type="text/javascript" src="_static/underscore.js"></script>
  25 + <script type="text/javascript" src="_static/doctools.js"></script>
  26 + <link rel="top" title="documentation_HFR_analysis 0.0.1 documentation" href="index.html" />
  27 + <link rel="next" title="Current combination and OGSL data feed" href="CurrentsCombination.html" />
  28 + <link rel="prev" title="WERA radial currents processing" href="WERAprocessing.html" />
  29 + </head>
  30 + <body>
  31 + <div class="related">
  32 + <h3>Navigation</h3>
  33 + <ul>
  34 + <li class="right" style="margin-right: 10px">
  35 + <a href="genindex.html" title="General Index"
  36 + accesskey="I">index</a></li>
  37 + <li class="right" >
  38 + <a href="CurrentsCombination.html" title="Current combination and OGSL data feed"
  39 + accesskey="N">next</a> |</li>
  40 + <li class="right" >
  41 + <a href="WERAprocessing.html" title="WERA radial currents processing"
  42 + accesskey="P">previous</a> |</li>
  43 + <li><a href="index.html">documentation_HFR_analysis 0.0.1 documentation</a> &raquo;</li>
  44 + </ul>
  45 + </div>
  46 +
  47 + <div class="document">
  48 + <div class="documentwrapper">
  49 + <div class="bodywrapper">
  50 + <div class="body">
  51 +
  52 + <div class="section" id="comparison-wera-and-adcp-radial-velocities-september-2013">
  53 +<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>
  54 +<p>Matlab script used to do the analysis is <strong>RadialVeloCompa_Sept2013_HFR_ADCP_WERA_ADCPcorrected_v2.m</strong></p>
  55 +<p>Steps of the algorithm:</p>
  56 +<p>1- we choose a site to analyze, either PAO (site 1) or PAB (site 2)</p>
  57 +<p>2- we load the bathymetry and the HFR sites&#8217;coordinates</p>
  58 +<p>3- we specify if we use the &#8220;old&#8221; or &#8220;new&#8221; calibrations. Usually &#8220;New&#8221;.</p>
  59 +<p>4- we load IML-4 ADCP data (2013) and gps data</p>
  60 +<p>5- we make a list of all the HFR files available for the site of interest.</p>
  61 +<p>6- we analyze the ADCP data (bin 2 only to get near-surface data):</p>
  62 +<blockquote>
  63 +<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>
  64 +<p>For each index, we look for the ADCP position (xadcp,yadcp)(via the gps coordinates).</p>
  65 +<p>We then look for the positions 15min before and 15min after.</p>
  66 +<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>
  67 +<p>We find the radial orientation at the ADCP position given by xadcp,yadcp: azimuth in degrees from HFR site to adcp position</p>
  68 +<p>We project U_mean_adcp onto the radial direction to get the radial current components (ur_adcp_proj, vr_adcp_proj).</p>
  69 +<p>We store all the meaningful parameters into a structure named adcp.</p>
  70 +</div></blockquote>
  71 +<p>7- We analyze the HFR data:</p>
  72 +<blockquote>
  73 +<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>
  74 +<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>
  75 +</div></blockquote>
  76 +<div class="figure align-center">
  77 +<img alt="_images/Radial_ADCPbin2_HFR_PAO_version5_Sept2013.png" src="_images/Radial_ADCPbin2_HFR_PAO_version5_Sept2013.png" />
  78 +<p class="caption">Hourly radial currents from ADCP bin 2 and HFR WERA at PAO - September 2013</p>
  79 +<div class="legend">
  80 +Superimposed are the 48-h averages</div>
  81 +</div>
  82 +<p>[[Fichier:gridwithnan_adcp.png]]</p>
  83 +<blockquote>
  84 +<div><p>We find the HFR Radial velocity at this position and project it onto zonal and meridional components.</p>
  85 +<p>We store all the meaningful parameters into a structure named radar.</p>
  86 +</div></blockquote>
  87 +<p>This analysis allows us to compare radial velocity components from each HFR site with the ADCP components.</p>
  88 +</div>
  89 +<div class="section" id="hfr-direction">
  90 +<h1>HFR direction<a class="headerlink" href="#hfr-direction" title="Permalink to this headline">¶</a></h1>
  91 +<p>To verify that the HFR looks in the right direction:</p>
  92 +<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>
  93 +<p>[[Fichier:rangeADCP.png]]</p>
  94 +<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>
  95 +<p>We store all the meaningful parameters into a structure named range.</p>
  96 +<p>Results from the whole analysis are stored, for example, under:
  97 +/run/media/bandma01/SeagateWireless/MEOPAR/Matlab/Analysis/Radial_Comparison/201309/ADCP-PAO/ADCP_PAO_version5adcpCorrected_newcalibrations_v2.mat</p>
  98 +<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>
  99 +<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>
  100 +<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>
  101 +</div>
  102 +
  103 +
  104 + </div>
  105 + </div>
  106 + </div>
  107 + <div class="sphinxsidebar">
  108 + <div class="sphinxsidebarwrapper">
  109 + <h3><a href="index.html">Table Of Contents</a></h3>
  110 + <ul>
  111 +<li><a class="reference internal" href="#">Comparison WERA and ADCP radial velocities September 2013</a></li>
  112 +<li><a class="reference internal" href="#hfr-direction">HFR direction</a></li>
  113 +</ul>
  114 +
  115 + <h4>Previous topic</h4>
  116 + <p class="topless"><a href="WERAprocessing.html"
  117 + title="previous chapter">WERA radial currents processing</a></p>
  118 + <h4>Next topic</h4>
  119 + <p class="topless"><a href="CurrentsCombination.html"
  120 + title="next chapter">Current combination and OGSL data feed</a></p>
  121 + <h3>This Page</h3>
  122 + <ul class="this-page-menu">
  123 + <li><a href="_sources/Comparison_WERA-ADCP_Sept2013.txt"
  124 + rel="nofollow">Show Source</a></li>
  125 + </ul>
  126 +<div id="searchbox" style="display: none">
  127 + <h3>Quick search</h3>
  128 + <form class="search" action="search.html" method="get">
  129 + <input type="text" name="q" />
  130 + <input type="submit" value="Go" />
  131 + <input type="hidden" name="check_keywords" value="yes" />
  132 + <input type="hidden" name="area" value="default" />
  133 + </form>
  134 + <p class="searchtip" style="font-size: 90%">
  135 + Enter search terms or a module, class or function name.
  136 + </p>
  137 +</div>
  138 +<script type="text/javascript">$('#searchbox').show(0);</script>
  139 + </div>
  140 + </div>
  141 + <div class="clearer"></div>
  142 + </div>
  143 + <div class="related">
  144 + <h3>Navigation</h3>
  145 + <ul>
  146 + <li class="right" style="margin-right: 10px">
  147 + <a href="genindex.html" title="General Index"
  148 + >index</a></li>
  149 + <li class="right" >
  150 + <a href="CurrentsCombination.html" title="Current combination and OGSL data feed"
  151 + >next</a> |</li>
  152 + <li class="right" >
  153 + <a href="WERAprocessing.html" title="WERA radial currents processing"
  154 + >previous</a> |</li>
  155 + <li><a href="index.html">documentation_HFR_analysis 0.0.1 documentation</a> &raquo;</li>
  156 + </ul>
  157 + </div>
  158 + <div class="footer">
  159 + &copy; Copyright 2015, Marion Bandet.
  160 + Created using <a href="http://sphinx-doc.org/">Sphinx</a> 1.2.3.
  161 + </div>
  162 + </body>
  163 +</html>
0 164 \ No newline at end of file
... ...
build/html/CurrentsCombination.html
... ... @@ -6,7 +6,7 @@
6 6 <head>
7 7 <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
8 8  
9   - <title>Current combination &mdash; documentation_HFR_analysis 0.0.1 documentation</title>
  9 + <title>Current combination and OGSL data feed &mdash; documentation_HFR_analysis 0.0.1 documentation</title>
10 10  
11 11 <link rel="stylesheet" href="_static/default.css" type="text/css" />
12 12 <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
... ... @@ -24,8 +24,8 @@
24 24 <script type="text/javascript" src="_static/underscore.js"></script>
25 25 <script type="text/javascript" src="_static/doctools.js"></script>
26 26 <link rel="top" title="documentation_HFR_analysis 0.0.1 documentation" href="index.html" />
27   - <link rel="next" title="OGSL quasi-realtime data feed" href="OGSL.html" />
28   - <link rel="prev" title="WERA POST=PROCESSING" href="WERA_processing.html" />
  27 + <link rel="next" title="Experiment on May 15, 2013" href="20130515.html" />
  28 + <link rel="prev" title="Comparison WERA and ADCP radial velocities September 2013" href="Comparison_WERA-ADCP_Sept2013.html" />
29 29 </head>
30 30 <body>
31 31 <div class="related">
... ... @@ -35,10 +35,10 @@
35 35 <a href="genindex.html" title="General Index"
36 36 accesskey="I">index</a></li>
37 37 <li class="right" >
38   - <a href="OGSL.html" title="OGSL quasi-realtime data feed"
  38 + <a href="20130515.html" title="Experiment on May 15, 2013"
39 39 accesskey="N">next</a> |</li>
40 40 <li class="right" >
41   - <a href="WERA_processing.html" title="WERA POST=PROCESSING"
  41 + <a href="Comparison_WERA-ADCP_Sept2013.html" title="Comparison WERA and ADCP radial velocities September 2013"
42 42 accesskey="P">previous</a> |</li>
43 43 <li><a href="index.html">documentation_HFR_analysis 0.0.1 documentation</a> &raquo;</li>
44 44 </ul>
... ... @@ -49,35 +49,49 @@
49 49 <div class="bodywrapper">
50 50 <div class="body">
51 51  
52   - <div class="section" id="current-combination">
53   -<h1>Current combination<a class="headerlink" href="#current-combination" title="Permalink to this headline">¶</a></h1>
  52 + <div class="section" id="current-combination-and-ogsl-data-feed">
  53 +<h1>Current combination and OGSL data feed<a class="headerlink" href="#current-combination-and-ogsl-data-feed" title="Permalink to this headline">¶</a></h1>
54 54 <p>Radial currents from WERA and CODAR systems are combined together on a grid to provide total currents approximately every hour.
55 55 The processing is conducted on the WERA server (ssh <a class="reference external" href="mailto:wera&#37;&#52;&#48;132&#46;215&#46;11&#46;14">wera<span>&#64;</span>132<span>&#46;</span>215<span>&#46;</span>11<span>&#46;</span>14</a>) and the resulting hourly currents are exported to .csv files to be displayed on he OGSL website (<a class="reference external" href="http://ogsl.ca/observations/">http://ogsl.ca/observations/</a>).</p>
56 56 <p>All the matlab scripts are located in ~/projets/hfr_data_processing/hfr_processing.</p>
57 57 <p>A crontab runs the analysis (crontab).</p>
58 58 <p>The processing runs as follows:</p>
59   -<p>&#8216;&#8217;&#8216;1- download the last available data at WERA sites&#8217;&#8216;&#8217; (radial currents - ascii format)</p>
60   -<p>For PAO, data is available at HH+00 and HH+30</p>
61   -<p>For PAB, data is available at HH+15 et HH+45</p>
  59 +<p>1- Download the last available data at each WERA sites (radial currents - ascii format)</p>
  60 +<blockquote>
  61 +<div><p>For PAO, data is available at HH+00 and HH+30.</p>
  62 +<p>For PAB, data is available at HH+15 et HH+45.</p>
62 63 <p>Every 10 min the crontab tries to download data via a shell script:</p>
  64 +</div></blockquote>
63 65 <p>5,15,25,35,45,50 * * * * /home/wera/projets/hfr-data-processing/hfr_process/download_wera_rt.sh</p>
64   -<p>The data is located in: /home/wera/codar_wera_ogsl/archives/data</p>
65   -<p>&#8216;&#8217;&#8216;2- download the last available data at CODAR sites&#8217;&#8216;&#8217; (.ruv)</p>
66   -<p>For PAP and STF, data is available at HH+00</p>
  66 +<blockquote>
  67 +<div>The data is located in: /home/wera/codar_wera_ogsl/archives/data</div></blockquote>
  68 +<p>2- Download the last available data at each CODAR sites (.ruv)</p>
  69 +<blockquote>
  70 +<div><p>For PAP and STF, data is available at HH+00</p>
67 71 <p>Every 5 minutes, the crontab tries to download data via a shell script:</p>
  72 +</div></blockquote>
68 73 <p>0,5,10,15,20,25,30,35,40,45,50,55 * * * * /home/wera/projets/hfr-data-processing/hfr_process/download_codar_rt.sh</p>
69   -<p>The data is located in: /home/wera/codar_wera_ogsl/archives/data</p>
70   -<p>&#8216;&#8217;&#8216;3-&#8216;&#8217;&#8217; at H+55min the crontab runs to &#8216;&#8217;&#8216;combine currents for the last hour&#8217;&#8216;&#8217; via a shell script:</p>
  74 +<blockquote>
  75 +<div>The data is located in: /home/wera/codar_wera_ogsl/archives/data</div></blockquote>
  76 +<p>3- at H+55min the crontab runs to combine currents for the last hour via a shell script:</p>
71 77 <p>55 * * * * /home/wera/projets/hfr-data-processing/hfr_process/ogsl_hfr_treatment.sh</p>
72 78 <p>What does the shell do?</p>
73   -<p>1- copy the .ruv and ascii files located in /home/wera/codar_wera_ogsl/archives/data in the working directory /home/wera/codar_wera_ogsl/travail</p>
74   -<p>2- unzip .gz files</p>
75   -<p>3- convert .ruv to .mat format (ruv2mat)</p>
76   -<p>4- convert ascii to .mat format (ascii2mat)</p>
77   -<p>5- generate (combine) total currents (rad2tot)</p>
  79 +<blockquote>
  80 +<div><ol class="loweralpha simple">
  81 +<li>copy the .ruv and ascii files located in /home/wera/codar_wera_ogsl/archives/data in the working directory /home/wera/codar_wera_ogsl/travail</li>
  82 +<li>unzip .gz files</li>
  83 +<li>convert .ruv to .mat format (<strong>ruv2mat.m</strong>)</li>
  84 +<li>convert ascii to .mat format (<strong>ascii2mat.m</strong>)</li>
  85 +<li>generate (combine) total currents (<strong>rad2tot.m</strong>)</li>
  86 +</ol>
  87 +</div></blockquote>
78 88 <p>At this stage a .mat file is generated (i.e. HFR_LSLE_20150817T170000Z.mat) for each hour.</p>
79   -<p>6- the routine currentsHfrWithAverageExport.m computes the 25-hour current average and exports the hourly currents into a .csv file to be displayed on the OGSL website.</p>
80   -<p>7- the .csv files are moved to the webserver (/srv/www/htdocs/wera/hfrdata) and the .mat are moved in archives (/home/wera/codar_wera_ogsl/archives/matfiles and /srv/www/htdocs/wera/hfrmat) and logs are kept in /home/wera/codar_wera_ogsl/logs (download_codar_rt.log, download_wera_rt.log, ogsl_hfr_treatment.log, reprocess_ogsl_hfr_treatment.log)</p>
  89 +<blockquote>
  90 +<div><ol class="loweralpha simple" start="6">
  91 +<li>the routine <strong>currentsHfrWithAverageExport.m</strong> computes the 25-hour current average and exports the hourly currents into a .csv file to be displayed on the OGSL website.</li>
  92 +<li>the .csv files are moved to the webserver (/srv/www/htdocs/wera/hfrdata) and the .mat are moved in archives (/home/wera/codar_wera_ogsl/archives/matfiles and /srv/www/htdocs/wera/hfrmat) and logs are kept in /home/wera/codar_wera_ogsl/logs (download_codar_rt.log, download_wera_rt.log, ogsl_hfr_treatment.log, reprocess_ogsl_hfr_treatment.log)</li>
  93 +</ol>
  94 +</div></blockquote>
81 95 </div>
82 96  
83 97  
... ... @@ -87,11 +101,11 @@ The processing is conducted on the WERA server (ssh &lt;a class=&quot;reference external
87 101 <div class="sphinxsidebar">
88 102 <div class="sphinxsidebarwrapper">
89 103 <h4>Previous topic</h4>
90   - <p class="topless"><a href="WERA_processing.html"
91   - title="previous chapter">WERA POST=PROCESSING</a></p>
  104 + <p class="topless"><a href="Comparison_WERA-ADCP_Sept2013.html"
  105 + title="previous chapter">Comparison WERA and ADCP radial velocities September 2013</a></p>
92 106 <h4>Next topic</h4>
93   - <p class="topless"><a href="OGSL.html"
94   - title="next chapter">OGSL quasi-realtime data feed</a></p>
  107 + <p class="topless"><a href="20130515.html"
  108 + title="next chapter">Experiment on May 15, 2013</a></p>
95 109 <h3>This Page</h3>
96 110 <ul class="this-page-menu">
97 111 <li><a href="_sources/CurrentsCombination.txt"
... ... @@ -121,10 +135,10 @@ The processing is conducted on the WERA server (ssh &lt;a class=&quot;reference external
121 135 <a href="genindex.html" title="General Index"
122 136 >index</a></li>
123 137 <li class="right" >
124   - <a href="OGSL.html" title="OGSL quasi-realtime data feed"
  138 + <a href="20130515.html" title="Experiment on May 15, 2013"
125 139 >next</a> |</li>
126 140 <li class="right" >
127   - <a href="WERA_processing.html" title="WERA POST=PROCESSING"
  141 + <a href="Comparison_WERA-ADCP_Sept2013.html" title="Comparison WERA and ADCP radial velocities September 2013"
128 142 >previous</a> |</li>
129 143 <li><a href="index.html">documentation_HFR_analysis 0.0.1 documentation</a> &raquo;</li>
130 144 </ul>
... ...
build/html/WERAprocessing.html
... ... @@ -6,7 +6,7 @@
6 6 <head>
7 7 <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
8 8  
9   - <title>WERA POST=PROCESSING &mdash; documentation_HFR_analysis 0.0.1 documentation</title>
  9 + <title>WERA radial currents processing &mdash; documentation_HFR_analysis 0.0.1 documentation</title>
10 10  
11 11 <link rel="stylesheet" href="_static/default.css" type="text/css" />
12 12 <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
... ... @@ -24,7 +24,7 @@
24 24 <script type="text/javascript" src="_static/underscore.js"></script>
25 25 <script type="text/javascript" src="_static/doctools.js"></script>
26 26 <link rel="top" title="documentation_HFR_analysis 0.0.1 documentation" href="index.html" />
27   - <link rel="next" title="Current combination" href="CurrentsCombination.html" />
  27 + <link rel="next" title="Comparison WERA and ADCP radial velocities September 2013" href="Comparison_WERA-ADCP_Sept2013.html" />
28 28 <link rel="prev" title="CODAR radial currents processing" href="CODARprocessing.html" />
29 29 </head>
30 30 <body>
... ... @@ -35,7 +35,7 @@
35 35 <a href="genindex.html" title="General Index"
36 36 accesskey="I">index</a></li>
37 37 <li class="right" >
38   - <a href="CurrentsCombination.html" title="Current combination"
  38 + <a href="Comparison_WERA-ADCP_Sept2013.html" title="Comparison WERA and ADCP radial velocities September 2013"
39 39 accesskey="N">next</a> |</li>
40 40 <li class="right" >
41 41 <a href="CODARprocessing.html" title="CODAR radial currents processing"
... ... @@ -49,8 +49,8 @@
49 49 <div class="bodywrapper">
50 50 <div class="body">
51 51  
52   - <div class="section" id="wera-post-processing">
53   -<h1>WERA POST=PROCESSING<a class="headerlink" href="#wera-post-processing" title="Permalink to this headline">¶</a></h1>
  52 + <div class="section" id="wera-radial-currents-processing">
  53 +<h1>WERA radial currents processing<a class="headerlink" href="#wera-radial-currents-processing" title="Permalink to this headline">¶</a></h1>
54 54 <p>Post-processing of WERA radial currents is much simpler than for the CODAR system.</p>
55 55 <p>The WERA system outputs data in folders identified by the year and Julian day (for example 2014101 for April 11, 2014).</p>
56 56 <p>In each folder we find .RFI and .SORT files for every hour (should find 24 files of each extention per day).</p>
... ... @@ -71,8 +71,8 @@
71 71 <p class="topless"><a href="CODARprocessing.html"
72 72 title="previous chapter">CODAR radial currents processing</a></p>
73 73 <h4>Next topic</h4>
74   - <p class="topless"><a href="CurrentsCombination.html"
75   - title="next chapter">Current combination</a></p>
  74 + <p class="topless"><a href="Comparison_WERA-ADCP_Sept2013.html"
  75 + title="next chapter">Comparison WERA and ADCP radial velocities September 2013</a></p>
76 76 <h3>This Page</h3>
77 77 <ul class="this-page-menu">
78 78 <li><a href="_sources/WERAprocessing.txt"
... ... @@ -102,7 +102,7 @@
102 102 <a href="genindex.html" title="General Index"
103 103 >index</a></li>
104 104 <li class="right" >
105   - <a href="CurrentsCombination.html" title="Current combination"
  105 + <a href="Comparison_WERA-ADCP_Sept2013.html" title="Comparison WERA and ADCP radial velocities September 2013"
106 106 >next</a> |</li>
107 107 <li class="right" >
108 108 <a href="CODARprocessing.html" title="CODAR radial currents processing"
... ...
build/html/_images/Radial_ADCPbin2_HFR_PAO_version5_Sept2013.png 0 → 100644

22.3 KB

build/html/_sources/Comparison_WERA-ADCP_Sept2013.txt 0 → 100644
... ... @@ -0,0 +1,87 @@
  1 +Comparison WERA and ADCP radial velocities September 2013
  2 +===============================================================
  3 +
  4 +Matlab script used to do the analysis is **RadialVeloCompa_Sept2013_HFR_ADCP_WERA_ADCPcorrected_v2.m**
  5 +
  6 +Steps of the algorithm:
  7 +
  8 +1- we choose a site to analyze, either PAO (site 1) or PAB (site 2)
  9 +
  10 +2- we load the bathymetry and the HFR sites'coordinates
  11 +
  12 +3- we specify if we use the "old" or "new" calibrations. Usually "New".
  13 +
  14 +4- we load IML-4 ADCP data (2013) and gps data
  15 +
  16 +5- we make a list of all the HFR files available for the site of interest.
  17 +
  18 +6- we analyze the ADCP data (bin 2 only to get near-surface data):
  19 +
  20 + We look for the indices of the ADCP time vector corresponding to the time period of interest (here 01 to 30 September 2013).
  21 +
  22 + For each index, we look for the ADCP position (xadcp,yadcp)(via the gps coordinates).
  23 +
  24 + We then look for the positions 15min before and 15min after.
  25 +
  26 + We estimate the mean drift of the ADCP within these 30 minutes (Umean_adcp_u, Umean_adcp_v and U_mean_adcp (total)).
  27 +
  28 + We find the radial orientation at the ADCP position given by xadcp,yadcp: azimuth in degrees from HFR site to adcp position
  29 +
  30 + We project U_mean_adcp onto the radial direction to get the radial current components (ur_adcp_proj, vr_adcp_proj).
  31 +
  32 + We store all the meaningful parameters into a structure named adcp.
  33 +
  34 +7- We analyze the HFR data:
  35 +
  36 + We look for the indices of the HFR time vector corresponding to the time period of interest (here 01 to 30 September 2013)
  37 +
  38 + 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).
  39 +
  40 +.. figure:: ../../figures/RadialVelocities/ADCP/Radial_ADCPbin2_HFR_PAO_version5_Sept2013.png
  41 + :align: center
  42 +
  43 + Hourly radial currents from ADCP bin 2 and HFR WERA at PAO - September 2013
  44 +
  45 + Superimposed are the 48-h averages
  46 +
  47 +
  48 +.. [Fichier:Exple_allgridpt.png]] [[Fichier:Gridwithnan.png]]
  49 +
  50 +
  51 + We then look for the HFR grid point nearest the ADCP position at the time of interest.
  52 +
  53 +[[Fichier:gridwithnan_adcp.png]]
  54 +
  55 + We find the HFR Radial velocity at this position and project it onto zonal and meridional components.
  56 +
  57 + We store all the meaningful parameters into a structure named radar.
  58 +
  59 +
  60 +This analysis allows us to compare radial velocity components from each HFR site with the ADCP components.
  61 +
  62 +
  63 +HFR direction
  64 +===============
  65 +
  66 +To verify that the HFR looks in the right direction:
  67 +
  68 +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.
  69 +
  70 +[[Fichier:rangeADCP.png]]
  71 +
  72 +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).
  73 +
  74 +We store all the meaningful parameters into a structure named range.
  75 +
  76 +Results from the whole analysis are stored, for example, under:
  77 +/run/media/bandma01/SeagateWireless/MEOPAR/Matlab/Analysis/Radial_Comparison/201309/ADCP-PAO/ADCP_PAO_version5adcpCorrected_newcalibrations_v2.mat
  78 +
  79 +
  80 +[[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]]
  81 +
  82 +
  83 +[[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]]
  84 +
  85 +
  86 +[[Fichier:ADCPbin2vsPAOSept2013.png|200px|thumb|left|radial velocities from ADCP bin 2 vs HFR WERA at PAO site September 2013 - in red: 48h averages]]]]
  87 +
... ...
build/html/_sources/CurrentsCombination.txt
1   -Current combination
2   -===================
  1 +Current combination and OGSL data feed
  2 +==========================================
3 3  
4 4 Radial currents from WERA and CODAR systems are combined together on a grid to provide total currents approximately every hour.
5 5 The processing is conducted on the WERA server (ssh wera@132.215.11.14) and the resulting hourly currents are exported to .csv files to be displayed on he OGSL website (http://ogsl.ca/observations/).
... ... @@ -10,49 +10,50 @@ A crontab runs the analysis (crontab).
10 10  
11 11 The processing runs as follows:
12 12  
13   -'''1- download the last available data at WERA sites''' (radial currents - ascii format)
  13 +1- Download the last available data at each WERA sites (radial currents - ascii format)
14 14  
15   -For PAO, data is available at HH+00 and HH+30
  15 + For PAO, data is available at HH+00 and HH+30.
16 16  
17   -For PAB, data is available at HH+15 et HH+45
  17 + For PAB, data is available at HH+15 et HH+45.
18 18  
19   -Every 10 min the crontab tries to download data via a shell script:
  19 + Every 10 min the crontab tries to download data via a shell script:
20 20  
21 21 5,15,25,35,45,50 * * * * /home/wera/projets/hfr-data-processing/hfr_process/download_wera_rt.sh
22 22  
23   -The data is located in: /home/wera/codar_wera_ogsl/archives/data
  23 + The data is located in: /home/wera/codar_wera_ogsl/archives/data
24 24  
  25 +..
25 26  
26   -'''2- download the last available data at CODAR sites''' (.ruv)
  27 +2- Download the last available data at each CODAR sites (.ruv)
27 28  
28   -For PAP and STF, data is available at HH+00
  29 + For PAP and STF, data is available at HH+00
29 30  
30   -Every 5 minutes, the crontab tries to download data via a shell script:
  31 + Every 5 minutes, the crontab tries to download data via a shell script:
31 32  
32 33 0,5,10,15,20,25,30,35,40,45,50,55 * * * * /home/wera/projets/hfr-data-processing/hfr_process/download_codar_rt.sh
33 34  
34   -The data is located in: /home/wera/codar_wera_ogsl/archives/data
  35 + The data is located in: /home/wera/codar_wera_ogsl/archives/data
35 36  
  37 +..
36 38  
37   -
38   -'''3-''' at H+55min the crontab runs to '''combine currents for the last hour''' via a shell script:
  39 +3- at H+55min the crontab runs to combine currents for the last hour via a shell script:
39 40  
40 41 55 * * * * /home/wera/projets/hfr-data-processing/hfr_process/ogsl_hfr_treatment.sh
41 42  
42 43 What does the shell do?
43 44  
44   -1- copy the .ruv and ascii files located in /home/wera/codar_wera_ogsl/archives/data in the working directory /home/wera/codar_wera_ogsl/travail
  45 + a) copy the .ruv and ascii files located in /home/wera/codar_wera_ogsl/archives/data in the working directory /home/wera/codar_wera_ogsl/travail
45 46  
46   -2- unzip .gz files
  47 + b) unzip .gz files
47 48  
48   -3- convert .ruv to .mat format (ruv2mat)
  49 + c) convert .ruv to .mat format (**ruv2mat.m**)
49 50  
50   -4- convert ascii to .mat format (ascii2mat)
  51 + d) convert ascii to .mat format (**ascii2mat.m**)
51 52  
52   -5- generate (combine) total currents (rad2tot)
  53 + e) generate (combine) total currents (**rad2tot.m**)
53 54  
54 55 At this stage a .mat file is generated (i.e. HFR_LSLE_20150817T170000Z.mat) for each hour.
55 56  
56   -6- the routine currentsHfrWithAverageExport.m computes the 25-hour current average and exports the hourly currents into a .csv file to be displayed on the OGSL website.
  57 + f) the routine **currentsHfrWithAverageExport.m** computes the 25-hour current average and exports the hourly currents into a .csv file to be displayed on the OGSL website.
57 58  
58   -7- the .csv files are moved to the webserver (/srv/www/htdocs/wera/hfrdata) and the .mat are moved in archives (/home/wera/codar_wera_ogsl/archives/matfiles and /srv/www/htdocs/wera/hfrmat) and logs are kept in /home/wera/codar_wera_ogsl/logs (download_codar_rt.log, download_wera_rt.log, ogsl_hfr_treatment.log, reprocess_ogsl_hfr_treatment.log)
  59 + g) the .csv files are moved to the webserver (/srv/www/htdocs/wera/hfrdata) and the .mat are moved in archives (/home/wera/codar_wera_ogsl/archives/matfiles and /srv/www/htdocs/wera/hfrmat) and logs are kept in /home/wera/codar_wera_ogsl/logs (download_codar_rt.log, download_wera_rt.log, ogsl_hfr_treatment.log, reprocess_ogsl_hfr_treatment.log)
... ...
build/html/_sources/WERAprocessing.txt
1   -WERA POST=PROCESSING
2   -====================
  1 +WERA radial currents processing
  2 +==================================
3 3  
4 4 Post-processing of WERA radial currents is much simpler than for the CODAR system.
5 5  
... ...
build/html/_sources/index.txt
... ... @@ -20,8 +20,8 @@ Documentation about the steps of the analysis
20 20  
21 21 CODARprocessing
22 22 WERAprocessing
  23 + Comparison_WERA-ADCP_Sept2013
23 24 CurrentsCombination
24   - OGSL
25 25  
26 26  
27 27  
... ...
build/html/index.html
... ... @@ -59,9 +59,10 @@
59 59 <li class="toctree-l2"><a class="reference internal" href="CODARprocessing.html#check-up">CHECK-UP</a></li>
60 60 </ul>
61 61 </li>
62   -<li class="toctree-l1"><a class="reference internal" href="WERAprocessing.html">WERA POST=PROCESSING</a></li>
63   -<li class="toctree-l1"><a class="reference internal" href="CurrentsCombination.html">Current combination</a></li>
64   -<li class="toctree-l1"><a class="reference internal" href="OGSL.html">OGSL quasi-realtime data feed</a></li>
  62 +<li class="toctree-l1"><a class="reference internal" href="WERAprocessing.html">WERA radial currents processing</a></li>
  63 +<li class="toctree-l1"><a class="reference internal" href="Comparison_WERA-ADCP_Sept2013.html">Comparison WERA and ADCP radial velocities September 2013</a></li>
  64 +<li class="toctree-l1"><a class="reference internal" href="Comparison_WERA-ADCP_Sept2013.html#hfr-direction">HFR direction</a></li>
  65 +<li class="toctree-l1"><a class="reference internal" href="CurrentsCombination.html">Current combination and OGSL data feed</a></li>
65 66 </ul>
66 67 </div>
67 68 </div>
... ...
build/html/searchindex.js
1   -Search.setIndex({envversion:42,terms:{all:[10,2],identifi:5,spectratoprocess:10,move:[10,2],lst:5,do_process_sort_pol:5,marion:[],estuari:6,follow:[10,2],disk:10,find:[0,5],paramet:10,cell:[],onli:10,postprocess:[],locat:2,hfr_lsle_20150817t170000z:2,radial2uv:10,css_:10,pab:2,dispers:0,instrument:6,theme:[],matlab:2,radialconfig:10,everi:[5,2],lsle:[8,5,6],main:0,local:10,approxim:2,lawrenc:6,patt_pap1_2014_01_14_1154_nosmooth:10,unzip:2,march:10,variou:6,get:10,than:5,conduct:[0,8,2],"export":2,dai:[10,5],projet:2,volum:10,ssh:[5,2],adcp:0,radiallst2uv:5,tri:2,measur:10,csv:2,ascii2mat:2,hour:[5,2],organ:0,hfrdata:2,hfrmat:2,should:5,spectraprocess:10,server:[10,2],measpattern:10,diagram:10,yyyyjjj:5,each:[5,2,6],output:[0,5],reprocess_ogsl_hfr_treat:2,page:0,xlp:10,www:2,transform:10,codarprocess:[],ftp_transfer:5,direct:10,crontab:2,past:10,sure:10,sampl:6,result:[5,2],year:5,download:[5,2],home:[5,2],"55min":2,webserv:2,index:[],what:[10,2],total:2,miss:10,section:0,per:5,research:[],categori:[],hourli:2,between:0,content:10,element:[],run:[10,2],codar2:10,codar1:10,advanc:10,spot:6,gener:[5,2],radial:2,htdoc:2,water:6,step:0,extent:5,found:5,org:5,come:10,rfi:5,care:10,comparison:0,convert:[5,2],last:2,ogsl_hfr_treat:2,row:[],codar_wera_ogsl:2,forecast:0,grid:2,idealpattern:10,filter:10,april:5,srv:2,span:[],improv:0,chang:10,origin:10,routin:2,via:2,pap1_css_:10,onc:10,modul:[],automat:10,system:[5,2],header:10,weather:[],drogu:6,txt:10,open:10,stf:[10,2],simpler:5,avail:2,from:[5,2,10],log:2,coupl:0,script:[5,2],transfer:10,loggin:10,currentshfrwithaverageexport:2,add:10,been:10,websit:2,start:10,much:5,music:10,interest:10,julian:5,spectraseri:10,stage:2,rad2tot:2,store:10,more:0,sort:5,download_wera_rt:2,setup:10,shell:[5,2],travail:2,saint:6,delet:10,copi:[10,2],search:[],back:10,flag:10,part:0,modifi:[],togeth:2,line:10,folder:[10,5],had:6,info:0,codardata:10,save:10,retriev:10,sometim:10,currentscombin:[],provid:2,remov:10,work:[0,2],smooth:10,displai:2,project:0,column:6,below:0,can:10,were:6,crad_pol_lst:5,about:0,minut:2,radialconfigs_pap1_20130201_bapm_music402002_nosmooth:10,observ:2,loop:10,deg:10,hfr:[0,2],comput:2,file:[5,2,10],creat:10,desktop:10,hfr_data_process:2,site:[0,5,2,10],subfold:5,heur:5,ruv2mat:2,archiv:[10,2],"20d":10,have:10,field:0,need:10,ismer:0,doe:[10,2],pap1:10,check:[],download_codar_rt:2,fill:10,uncheck:10,"_ascii":5,date:5,thei:5,antenna:10,min:2,check_up:[],ascii:2,make:10,format:2,dyndn:5,also:10,ideal:10,other:0,option:[],uqar:[0,10],verifi:5,you:[0,10],trash:10,config:10,pao:[5,2],pap:[10,2],ruv:[10,2],relat:0,http:2,mat:[5,2,10],after:10,drift:0,rsync:5,seasond:10,marin:0,befor:10,meopar:0,attent:10,deploi:6,phase:10,radialconfigs_pap1_20130201_bapm_music402002_20dsmooth:10,contain:[],data:[],model:0,averag:2,lower:6,appropri:10,kept:2,hfr_process:2,directori:2,spectra:10,season:[],buoi:[0,6],depth:6,exampl:[10,5],matfil:2,thi:[0,2,6],cs4:10,everyth:10,bodi:[],bandet:[],order:6,weraprocess:[]},objtypes:{},objnames:{},filenames:["index","20130826","CurrentsCombination","20130830","20130515","WERAprocessing","20140506","OGSL","experiments","postprocessing","CODARprocessing"],titles:["Welcome to the documentation of High Frequency Radars (CODAR and WERA)","Experiment on August 26, 2013","Current combination","Experiment on August 20, 2013","Experiment on May 15, 2013","WERA POST=PROCESSING","Experiment on May 06, 2014","OGSL quasi-realtime data feed","&lt;no title&gt;","CODAR processing","CODAR radial currents processing"],objects:{},titleterms:{feed:[9,7],process:[10,5,9],indic:[],high:0,radar:0,tabl:[],check:10,mai:[4,6],welcom:0,current:[10,2],quasi:7,postprocess:0,document:0,codar:[0,9,10],analysi:[],august:[1,3],radial:10,ogsl:[9,7],documentation_hfr_analysi:[],post:5,data:[9,7],experi:[0,1,4,6,3],realtim:7,list:[],wera:[0,5,9],combin:2,frequenc:0}})
2 1 \ No newline at end of file
  2 +Search.setIndex({envversion:42,terms:{all:[0,10,3],correl:0,marion:[],follow:[10,3],disk:10,row:[],content:10,matlab:[0,3],under:0,unzip:3,everi:[6,3],projet:3,volum:10,spectratoprocess:10,cyan:0,tri:3,hour:[6,3],hfrmat:3,measpattern:10,direct:10,past:10,download:[6,3],index:0,what:[10,3],compar:0,section:1,advanc:10,cell:[],hourli:[0,3],conduct:[1,8,3],"new":0,whatev:[],codar2:10,codar1:10,gener:[6,3],here:0,water:7,modifi:[],"200px":0,search:[],idealpattern:10,hfr_data_process:3,chang:10,via:[0,3],pap1_css_:10,lawrenc:7,modul:[],"_ascii":6,total:[0,3],select:0,from:[0,6,3,10],websit:3,music:10,postprocess:[],stage:3,rad2tot:3,more:1,sort:6,relat:1,flag:10,retriev:10,currentscombin:[],setup:10,work:[1,3],can:10,deg:10,heur:6,download_codar_rt:3,uncheck:10,surfac:0,radial_comparison:0,check_up:[],radialvelocompa_sept2013_hfr_adcp_wera_adcpcorrected_v2:0,uqar:[1,10],verifi:[0,6],csv:3,config:10,mat:[0,6,3,10],after:[0,10],spot:7,diagram:10,befor:[0,10],meopar:[0,1],attent:10,date:6,data:[],averag:[0,3],youwish:[],spectra:10,correspond:0,element:[],allow:0,cs4:10,order:7,move:[10,3],lst:6,same:0,estuari:7,paramet:[0,10],thumb:0,hfr_lsle_20150817t170000z:3,iml:0,onc:10,lsle:[8,6,7],main:1,thei:6,dai:[10,6],name:0,each:[0,6,3,7],found:6,reprocess_ogsl_hfr_treat:3,xlp:10,crontab:3,year:6,extract:0,orient:0,open:10,research:[],categori:[],vector:0,red:0,bathymetri:0,gridwithnan_adcp:0,org:6,care:10,angl:0,filter:10,caption:[],onto:0,origin:10,rang:0,do_process_sort_pol:6,currentshfrwithaverageexport:3,sometim:10,compsnewcalib_radialvelo_adcp_pao_sept2013:0,miss:10,given:0,script:[0,6,3],system:[6,3],spectraseri:10,saint:7,store:[0,10],shell:[6,3],travail:3,option:[],season:[],copi:[10,3],specifi:0,part:1,than:[0,6],png:0,stf:[10,3],provid:3,dispers:1,remov:10,structur:0,project:[0,1],were:7,posit:0,minut:[0,3],hfr:[],comput:[0,3],adcp_pao_version5adcpcorrected_newcalibrations_v2:0,julian:6,"20d":10,have:10,need:10,ismer:1,min:3,superimpos:0,also:10,ideal:10,sure:10,seasond:10,marin:1,ogsl_hfr_treat:3,deploi:7,phase:10,everyth:10,doe:[10,3],azimuth:0,left:0,radial2uv:10,find:[0,1,6],fichier:0,onli:[0,10],locat:3,delet:10,should:6,theme:[],analyz:0,folder:[10,6],adcpbin2vspaosept2013:0,local:10,info:1,variou:7,get:[0,10],yadcp:0,ssh:[6,3],adcp:[],ascii2mat:3,organ:1,contain:0,rfi:6,ftp_transfer:6,result:[0,6,3],"55min":3,forecast:1,figur:0,between:1,theta:0,betwen:0,extent:6,subfold:6,bodi:[],last:3,umean_adcp_v:0,umean_adcp_u:0,april:6,"15min":0,improv:1,whole:0,load:0,point:0,ur_adcp_proj:0,period:0,header:10,drogu:7,diff:0,coupl:1,come:10,loggin:10,been:10,compon:0,much:6,interest:[0,10],download_wera_rt:3,imag:[],convert:[6,3],coordin:0,togeth:3,look:0,observ:3,loop:10,site:[0,1,6,3,10],archiv:[10,3],zonal:0,ascii:3,pap:[10,3],media:0,make:[0,10],dyndn:6,pab:[0,3],meridion:0,pao:[0,6,3],http:3,webserv:3,closest:0,radialconfigs_pap1_20130201_bapm_music402002_20dsmooth:10,lower:7,appropri:10,kept:3,hfr_process:3,exampl:[0,10,6],thi:[0,1,3,7],choos:0,model:1,bandet:[],usual:0,distanc:0,identifi:6,calibr:0,css_:10,had:7,desktop:10,instrument:7,add:10,radialconfig:10,els:[],save:10,approxim:3,bin:0,march:10,format:3,grid:[0,3],weraprocess:[],radiallst2uv:6,measur:10,like:[],hfrdata:3,rangeadcp:0,spectraprocess:10,server:[10,3],either:0,output:[1,6],page:1,www:3,right:0,old:0,codarprocess:[],back:10,sampl:7,home:[6,3],rsync:6,txt:10,per:6,seagatewireless:0,estim:0,run:[0,10,3],radial:[],step:[0,1],ruv:[10,3],"20132440000_pao":0,simpler:6,comparison:[],about:1,bandma01:0,column:7,codar_wera_ogsl:3,degre:0,srv:3,mean:0,routin:3,within:0,automat:10,vr_adcp_proj:0,weather:[],span:[],log:3,transfer:10,transform:10,avail:[0,3],start:10,yyyyjjj:6,line:10,u_mean_adcp:0,patt_pap1_2014_01_14_1154_nosmooth:10,codardata:10,smooth:10,displai:3,ruv2mat:3,below:1,crad_pol_lst:6,"export":3,radialconfigs_pap1_20130201_bapm_music402002_nosmooth:10,creat:10,"48h":0,file:[0,6,3,10],pap1:10,check:[],fill:10,antenna:10,field:1,other:1,you:[1,10],trash:10,drift:[0,1],meaning:0,xadcp:0,algorithm:0,directori:3,htdoc:3,orientationpaosept2013:0,buoi:[1,7],depth:7,matfil:3,time:0},objtypes:{},objnames:{},filenames:["Comparison_WERA-ADCP_Sept2013","index","20130826","CurrentsCombination","20130830","20130515","WERAprocessing","20140506","experiments","postprocessing","CODARprocessing"],titles:["Comparison WERA and ADCP radial velocities September 2013","Welcome to the documentation of High Frequency Radars (CODAR and WERA)","Experiment on August 26, 2013","Current combination and OGSL data feed","Experiment on August 20, 2013","Experiment on May 15, 2013","WERA radial currents processing","Experiment on May 06, 2014","&lt;no title&gt;","CODAR processing","CODAR radial currents processing"],objects:{},titleterms:{feed:[3,9],process:[10,6,9],direct:0,indic:[],high:1,radar:1,tabl:[],check:10,mai:[5,7],welcom:1,current:[6,3,10],quasi:[],postprocess:1,document:1,codar:[1,9,10],analysi:[],august:[2,4],realtim:[],ogsl:[3,9],adcp:0,documentation_hfr_analysi:[],post:[],data:[3,9],comparison:0,experi:[1,2,5,7,4],frequenc:1,veloc:0,hfr:0,list:[],wera:[0,1,6,9],combin:3,radial:[0,10,6],septemb:0}})
3 3 \ No newline at end of file
... ...
source/Comparison_WERA-ADCP_Sept2013.rst
... ... @@ -37,7 +37,15 @@ Steps of the algorithm:
37 37  
38 38 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).
39 39  
40   -[[Fichier:Exple_allgridpt.png]] [[Fichier:Gridwithnan.png]]
  40 +.. figure:: ../../figures/RadialVelocities/ADCP/Radial_ADCPbin2_HFR_PAO_version5_Sept2013.png
  41 + :align: center
  42 +
  43 + Hourly radial currents from ADCP bin 2 and HFR WERA at PAO - September 2013
  44 +
  45 + Superimposed are the 48-h averages
  46 +
  47 +
  48 +.. [Fichier:Exple_allgridpt.png]] [[Fichier:Gridwithnan.png]]
41 49  
42 50  
43 51 We then look for the HFR grid point nearest the ADCP position at the time of interest.
... ... @@ -52,8 +60,10 @@ Steps of the algorithm:
52 60 This analysis allows us to compare radial velocity components from each HFR site with the ADCP components.
53 61  
54 62  
  63 +HFR direction
  64 +===============
55 65  
56   -To verify that the HFR look in the right direction:
  66 +To verify that the HFR looks in the right direction:
57 67  
58 68 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.
59 69  
... ...
source/figures/Radial_ADCPbin2_HFR_PAO_version5_Sept2013.png 0 → 100644

22.3 KB