Commit e48dcce599f612b64e7930ad54ffc15baadd6af3

Authored by Marion Bandet
1 parent 037a0eaa
Exists in master

update YSI data analysis May 2014

Showing 51 changed files with 1005 additions and 60 deletions   Show diff stats
build/doctrees/20140506.doctree
No preview for this file type
build/doctrees/20140506_ADCP_V50.doctree
No preview for this file type
build/doctrees/20140506_YSI_CastAway_CTD.doctree 0 → 100644
No preview for this file type
build/doctrees/YSI_CastAway_CTD.doctree 0 → 100644
No preview for this file type
build/doctrees/environment.pickle
No preview for this file type
build/html/20140506.html
... ... @@ -55,12 +55,13 @@
55 55 GPS were programmed to sample every 10 seconds and spots every 5 min.</p>
56 56 <p>The buoys were deployed by pairs of matching drogue&#8217;s depth, however due to satellite communication issues, not all GPSs provided data right after deployment. A delay of up to ?? min was observed.</p>
57 57 <p>A RDInstrument 300 kHz Acoustic Doppler Current Profiler (ADCP) was deployed for four (4) hours on May 7, 2014 while the buoys were drifting. The ADCP was mounted at the center of a donut-shaped float on a swivel ring to minimize the effects of pitch and roll. The ADCP was programmed to sample the water column using twenty-six (26) 0.5m bins. Vertical velocity profiles were exploited to compare with other estimates of velocity components.</p>
58   -<p>A portable YSI CTD was used to sample vertical profiles of temperature and conductivity at 10 locations in the study area. The data were used to identify different water masses.</p>
  58 +<p>To view data from the ADCP V50, go here: <a class="reference external" href="20140506_ADCP_V50.html">20140506_ADCP_V50</a></p>
  59 +<p>A portable YSI CTD was used to sample vertical profiles of temperature and conductivity at 10 locations in the study area. The data was used to identify different water masses.</p>
  60 +<p>To view data from the YSI CTD, go here: <a class="reference external" href="20140506_YSI_CastAway_CTD.html">YSI_CastAway_CTD</a></p>
59 61 <p>In addition, velocity field from MOGSL model output was obtained for comparison purposes.</p>
60 62 <p>Tides were obtained from <a class="reference external" href="http://tides.gc.ca/eng/data/table/2014/">http://tides.gc.ca/eng/data/table/2014/</a> at Rimouski.
61 63 Hourly wind data were downloaded from <a class="reference external" href="http://climat.meteo.gc.ca/">http://climat.meteo.gc.ca/</a> at station Pointe-au-Père, QC.
62 64 Simulated wave data were obtained from the wave model WAVEWATCH III (<a class="reference external" href="http://polar.ncep.noaa.gov/waves/index2.shtml">http://polar.ncep.noaa.gov/waves/index2.shtml</a>) outputs run at UQAR for the area of interest.</p>
63   -<p>To view data from the ADCP V50, go here: &#8216;20140506_ADCP_V50.rst`_</p>
64 65 </div>
65 66  
66 67  
... ...
build/html/20140506_ADCP_V50.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>Data from ADCP V50 (from routine &mdash; documentation_HFR_analysis 0.0.1 documentation</title>
  9 + <title>roll - pitch - heading &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" />
... ... @@ -41,9 +41,10 @@
41 41 <div class="bodywrapper">
42 42 <div class="body">
43 43  
44   - <div class="section" id="data-from-adcp-v50-from-routine">
45   -<h1>Data from ADCP V50 (from routine<a class="headerlink" href="#data-from-adcp-v50-from-routine" title="Permalink to this headline">¶</a></h1>
46   -<p>roll - pitch - heading</p>
  44 + <p>Data from ADCP V50 (from routine <a class="reference external" href="ADCP_V50_May2014.html">ADCP_V50_may2014.m</a>)</p>
  45 +<hr class="docutils" />
  46 +<div class="section" id="roll-pitch-heading">
  47 +<h1>roll - pitch - heading<a class="headerlink" href="#roll-pitch-heading" title="Permalink to this headline">¶</a></h1>
47 48 <div class="figure align-center">
48 49 <img alt="_images/ADCP_V50_May2014_01.png" src="_images/ADCP_V50_May2014_01.png" />
49 50 </div>
... ... @@ -56,39 +57,57 @@
56 57 <div class="figure align-center">
57 58 <img alt="_images/ADCP_V50_May2014_05.png" src="_images/ADCP_V50_May2014_05.png" />
58 59 </div>
59   -<p>range</p>
  60 +</div>
  61 +<div class="section" id="range">
  62 +<h1>range<a class="headerlink" href="#range" title="Permalink to this headline">¶</a></h1>
60 63 <div class="figure align-center">
61 64 <img alt="_images/ADCP_V50_May2014_08.png" src="_images/ADCP_V50_May2014_08.png" />
62 65 </div>
63   -<p>Percent Good</p>
  66 +</div>
  67 +<div class="section" id="percent-good">
  68 +<h1>Percent Good<a class="headerlink" href="#percent-good" title="Permalink to this headline">¶</a></h1>
64 69 <div class="figure align-center">
65 70 <img alt="_images/ADCP_V50_May2014_09.png" src="_images/ADCP_V50_May2014_09.png" />
66 71 </div>
67   -<p>Correlation</p>
  72 +</div>
  73 +<div class="section" id="correlation">
  74 +<h1>Correlation<a class="headerlink" href="#correlation" title="Permalink to this headline">¶</a></h1>
68 75 <div class="figure align-center">
69 76 <img alt="_images/ADCP_V50_May2014_10.png" src="_images/ADCP_V50_May2014_10.png" />
70 77 </div>
71   -<p>Intensity</p>
  78 +</div>
  79 +<div class="section" id="intensity">
  80 +<h1>Intensity<a class="headerlink" href="#intensity" title="Permalink to this headline">¶</a></h1>
72 81 <div class="figure align-center">
73 82 <img alt="_images/ADCP_V50_May2014_11.png" src="_images/ADCP_V50_May2014_11.png" />
74 83 </div>
75   -<p>Velocity (ENUErr)</p>
  84 +</div>
  85 +<div class="section" id="velocity-enuerr">
  86 +<h1>Velocity (ENUErr)<a class="headerlink" href="#velocity-enuerr" title="Permalink to this headline">¶</a></h1>
76 87 <div class="figure align-center">
77 88 <img alt="_images/ADCP_V50_May2014_12.png" src="_images/ADCP_V50_May2014_12.png" />
78 89 </div>
79   -<p>Averaged velocity</p>
  90 +</div>
  91 +<div class="section" id="averaged-velocity">
  92 +<h1>Averaged velocity<a class="headerlink" href="#averaged-velocity" title="Permalink to this headline">¶</a></h1>
80 93 <div class="figure align-center">
81 94 <img alt="_images/ADCP_V50_May2014_13.png" src="_images/ADCP_V50_May2014_13.png" />
82 95 </div>
83   -<p>Averaged Intensity</p>
  96 +</div>
  97 +<div class="section" id="averaged-intensity">
  98 +<h1>Averaged Intensity<a class="headerlink" href="#averaged-intensity" title="Permalink to this headline">¶</a></h1>
84 99 <div class="figure align-center">
85 100 <img alt="_images/ADCP_V50_May2014_14.png" src="_images/ADCP_V50_May2014_14.png" />
86 101 </div>
87   -<p>Averaged Correlation</p>
  102 +</div>
  103 +<div class="section" id="averaged-correlation">
  104 +<h1>Averaged Correlation<a class="headerlink" href="#averaged-correlation" title="Permalink to this headline">¶</a></h1>
88 105 <div class="figure align-center">
89 106 <img alt="_images/ADCP_V50_May2014_15.png" src="_images/ADCP_V50_May2014_15.png" />
90 107 </div>
91   -<p>Vertical Profiles of averaged velocity components</p>
  108 +</div>
  109 +<div class="section" id="vertical-profiles-of-averaged-velocity-components">
  110 +<h1>Vertical Profiles of averaged velocity components<a class="headerlink" href="#vertical-profiles-of-averaged-velocity-components" title="Permalink to this headline">¶</a></h1>
92 111 <div class="figure align-center">
93 112 <img alt="_images/ADCP_V50_May2014_16.png" src="_images/ADCP_V50_May2014_16.png" />
94 113 </div>
... ... @@ -109,6 +128,20 @@
109 128 </div>
110 129 <div class="sphinxsidebar">
111 130 <div class="sphinxsidebarwrapper">
  131 + <h3><a href="index.html">Table Of Contents</a></h3>
  132 + <ul>
  133 +<li><a class="reference internal" href="#">roll - pitch - heading</a></li>
  134 +<li><a class="reference internal" href="#range">range</a></li>
  135 +<li><a class="reference internal" href="#percent-good">Percent Good</a></li>
  136 +<li><a class="reference internal" href="#correlation">Correlation</a></li>
  137 +<li><a class="reference internal" href="#intensity">Intensity</a></li>
  138 +<li><a class="reference internal" href="#velocity-enuerr">Velocity (ENUErr)</a></li>
  139 +<li><a class="reference internal" href="#averaged-velocity">Averaged velocity</a></li>
  140 +<li><a class="reference internal" href="#averaged-intensity">Averaged Intensity</a></li>
  141 +<li><a class="reference internal" href="#averaged-correlation">Averaged Correlation</a></li>
  142 +<li><a class="reference internal" href="#vertical-profiles-of-averaged-velocity-components">Vertical Profiles of averaged velocity components</a></li>
  143 +</ul>
  144 +
112 145 <h3>This Page</h3>
113 146 <ul class="this-page-menu">
114 147 <li><a href="_sources/20140506_ADCP_V50.txt"
... ...
build/html/20140506_YSI_CastAway_CTD.html 0 → 100644
... ... @@ -0,0 +1,128 @@
  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>Casts location &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 + </head>
  28 + <body>
  29 + <div class="related">
  30 + <h3>Navigation</h3>
  31 + <ul>
  32 + <li class="right" style="margin-right: 10px">
  33 + <a href="genindex.html" title="General Index"
  34 + accesskey="I">index</a></li>
  35 + <li><a href="index.html">documentation_HFR_analysis 0.0.1 documentation</a> &raquo;</li>
  36 + </ul>
  37 + </div>
  38 +
  39 + <div class="document">
  40 + <div class="documentwrapper">
  41 + <div class="bodywrapper">
  42 + <div class="body">
  43 +
  44 + <p>Data from YSI CastAway (from routine <a class="reference external" href="YSI_CastAway_CTD_May2014.html">YSI_CastAway_CTD_May2014.m</a>)</p>
  45 +<hr class="docutils" />
  46 +<div class="section" id="casts-location">
  47 +<h1>Casts location<a class="headerlink" href="#casts-location" title="Permalink to this headline">¶</a></h1>
  48 +<div class="figure align-center">
  49 +<img alt="_images/YSI_CastAway_CTD_May2014_01.png" src="_images/YSI_CastAway_CTD_May2014_01.png" />
  50 +</div>
  51 +</div>
  52 +<div class="section" id="temperature">
  53 +<h1>Temperature<a class="headerlink" href="#temperature" title="Permalink to this headline">¶</a></h1>
  54 +<div class="figure align-center">
  55 +<img alt="_images/YSI_CastAway_CTD_May2014_02.png" src="_images/YSI_CastAway_CTD_May2014_02.png" />
  56 +</div>
  57 +</div>
  58 +<div class="section" id="conductivity">
  59 +<h1>Conductivity<a class="headerlink" href="#conductivity" title="Permalink to this headline">¶</a></h1>
  60 +<div class="figure align-center">
  61 +<img alt="_images/YSI_CastAway_CTD_May2014_03.png" src="_images/YSI_CastAway_CTD_May2014_03.png" />
  62 +</div>
  63 +</div>
  64 +<div class="section" id="salinity">
  65 +<h1>Salinity<a class="headerlink" href="#salinity" title="Permalink to this headline">¶</a></h1>
  66 +<div class="figure align-center">
  67 +<img alt="_images/YSI_CastAway_CTD_May2014_04.png" src="_images/YSI_CastAway_CTD_May2014_04.png" />
  68 +</div>
  69 +</div>
  70 +<div class="section" id="density">
  71 +<h1>Density<a class="headerlink" href="#density" title="Permalink to this headline">¶</a></h1>
  72 +<div class="figure align-center">
  73 +<img alt="_images/YSI_CastAway_CTD_May2014_05.png" src="_images/YSI_CastAway_CTD_May2014_05.png" />
  74 +</div>
  75 +</div>
  76 +
  77 +
  78 + </div>
  79 + </div>
  80 + </div>
  81 + <div class="sphinxsidebar">
  82 + <div class="sphinxsidebarwrapper">
  83 + <h3><a href="index.html">Table Of Contents</a></h3>
  84 + <ul>
  85 +<li><a class="reference internal" href="#">Casts location</a></li>
  86 +<li><a class="reference internal" href="#temperature">Temperature</a></li>
  87 +<li><a class="reference internal" href="#conductivity">Conductivity</a></li>
  88 +<li><a class="reference internal" href="#salinity">Salinity</a></li>
  89 +<li><a class="reference internal" href="#density">Density</a></li>
  90 +</ul>
  91 +
  92 + <h3>This Page</h3>
  93 + <ul class="this-page-menu">
  94 + <li><a href="_sources/20140506_YSI_CastAway_CTD.txt"
  95 + rel="nofollow">Show Source</a></li>
  96 + </ul>
  97 +<div id="searchbox" style="display: none">
  98 + <h3>Quick search</h3>
  99 + <form class="search" action="search.html" method="get">
  100 + <input type="text" name="q" />
  101 + <input type="submit" value="Go" />
  102 + <input type="hidden" name="check_keywords" value="yes" />
  103 + <input type="hidden" name="area" value="default" />
  104 + </form>
  105 + <p class="searchtip" style="font-size: 90%">
  106 + Enter search terms or a module, class or function name.
  107 + </p>
  108 +</div>
  109 +<script type="text/javascript">$('#searchbox').show(0);</script>
  110 + </div>
  111 + </div>
  112 + <div class="clearer"></div>
  113 + </div>
  114 + <div class="related">
  115 + <h3>Navigation</h3>
  116 + <ul>
  117 + <li class="right" style="margin-right: 10px">
  118 + <a href="genindex.html" title="General Index"
  119 + >index</a></li>
  120 + <li><a href="index.html">documentation_HFR_analysis 0.0.1 documentation</a> &raquo;</li>
  121 + </ul>
  122 + </div>
  123 + <div class="footer">
  124 + &copy; Copyright 2015, Marion Bandet.
  125 + Created using <a href="http://sphinx-doc.org/">Sphinx</a> 1.2.3.
  126 + </div>
  127 + </body>
  128 +</html>
0 129 \ No newline at end of file
... ...
source/figures/20140506/html/ADCP_V50_May2014.html renamed to build/html/ADCP_V50_May2014.html
source/figures/20140506/html/ADCP_V50_May2014.png renamed to build/html/ADCP_V50_May2014.png

3.28 KB

source/figures/20140506/html/ADCP_V50_May2014_01.png renamed to build/html/ADCP_V50_May2014_01.png

9.09 KB

source/figures/20140506/html/ADCP_V50_May2014_02.png renamed to build/html/ADCP_V50_May2014_02.png

8.75 KB

source/figures/20140506/html/ADCP_V50_May2014_03.png renamed to build/html/ADCP_V50_May2014_03.png

8.59 KB

source/figures/20140506/html/ADCP_V50_May2014_04.png renamed to build/html/ADCP_V50_May2014_04.png

8.37 KB

source/figures/20140506/html/ADCP_V50_May2014_05.png renamed to build/html/ADCP_V50_May2014_05.png

8.84 KB

source/figures/20140506/html/ADCP_V50_May2014_06.png renamed to build/html/ADCP_V50_May2014_06.png

9.27 KB

source/figures/20140506/html/ADCP_V50_May2014_07.png renamed to build/html/ADCP_V50_May2014_07.png

2.98 KB

source/figures/20140506/html/ADCP_V50_May2014_08.png renamed to build/html/ADCP_V50_May2014_08.png

4.19 KB

source/figures/20140506/html/ADCP_V50_May2014_09.png renamed to build/html/ADCP_V50_May2014_09.png

18.6 KB

source/figures/20140506/html/ADCP_V50_May2014_10.png renamed to build/html/ADCP_V50_May2014_10.png

7.72 KB

source/figures/20140506/html/ADCP_V50_May2014_11.png renamed to build/html/ADCP_V50_May2014_11.png

19.9 KB

source/figures/20140506/html/ADCP_V50_May2014_12.png renamed to build/html/ADCP_V50_May2014_12.png

10.4 KB

source/figures/20140506/html/ADCP_V50_May2014_13.png renamed to build/html/ADCP_V50_May2014_13.png

9.9 KB

source/figures/20140506/html/ADCP_V50_May2014_14.png renamed to build/html/ADCP_V50_May2014_14.png

13.1 KB

source/figures/20140506/html/ADCP_V50_May2014_15.png renamed to build/html/ADCP_V50_May2014_15.png

33.9 KB

source/figures/20140506/html/ADCP_V50_May2014_16.png renamed to build/html/ADCP_V50_May2014_16.png

7.47 KB

source/figures/20140506/html/ADCP_V50_May2014_17.png renamed to build/html/ADCP_V50_May2014_17.png

7.62 KB

source/figures/20140506/html/ADCP_V50_May2014_18.png renamed to build/html/ADCP_V50_May2014_18.png

7.89 KB

source/figures/20140506/html/ADCP_V50_May2014_19.png renamed to build/html/ADCP_V50_May2014_19.png

7.79 KB

build/html/YSI_CastAway_CTD.html 0 → 100644
... ... @@ -0,0 +1,107 @@
  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>Temperature &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 + </head>
  28 + <body>
  29 + <div class="related">
  30 + <h3>Navigation</h3>
  31 + <ul>
  32 + <li class="right" style="margin-right: 10px">
  33 + <a href="genindex.html" title="General Index"
  34 + accesskey="I">index</a></li>
  35 + <li><a href="index.html">documentation_HFR_analysis 0.0.1 documentation</a> &raquo;</li>
  36 + </ul>
  37 + </div>
  38 +
  39 + <div class="document">
  40 + <div class="documentwrapper">
  41 + <div class="bodywrapper">
  42 + <div class="body">
  43 +
  44 + <p>Data from YSI CTD (from routine <a class="reference external" href="YSI_CastAway_CTD_May2014.html">YSI_CastAway_CTD_May2014.m</a>)</p>
  45 +<hr class="docutils" />
  46 +<div class="section" id="temperature">
  47 +<h1>Temperature<a class="headerlink" href="#temperature" title="Permalink to this headline">¶</a></h1>
  48 +<div class="figure align-center">
  49 +<img alt="_images/YSI_CastAway_CTD_May2014_01.png" src="_images/YSI_CastAway_CTD_May2014_01.png" />
  50 +</div>
  51 +</div>
  52 +<div class="section" id="conductivity">
  53 +<h1>Conductivity<a class="headerlink" href="#conductivity" title="Permalink to this headline">¶</a></h1>
  54 +<div class="figure align-center">
  55 +<img alt="_images/YSI_CastAway_CTD_May2014_02.png" src="_images/YSI_CastAway_CTD_May2014_02.png" />
  56 +</div>
  57 +</div>
  58 +
  59 +
  60 + </div>
  61 + </div>
  62 + </div>
  63 + <div class="sphinxsidebar">
  64 + <div class="sphinxsidebarwrapper">
  65 + <h3><a href="index.html">Table Of Contents</a></h3>
  66 + <ul>
  67 +<li><a class="reference internal" href="#">Temperature</a></li>
  68 +<li><a class="reference internal" href="#conductivity">Conductivity</a></li>
  69 +</ul>
  70 +
  71 + <h3>This Page</h3>
  72 + <ul class="this-page-menu">
  73 + <li><a href="_sources/YSI_CastAway_CTD.txt"
  74 + rel="nofollow">Show Source</a></li>
  75 + </ul>
  76 +<div id="searchbox" style="display: none">
  77 + <h3>Quick search</h3>
  78 + <form class="search" action="search.html" method="get">
  79 + <input type="text" name="q" />
  80 + <input type="submit" value="Go" />
  81 + <input type="hidden" name="check_keywords" value="yes" />
  82 + <input type="hidden" name="area" value="default" />
  83 + </form>
  84 + <p class="searchtip" style="font-size: 90%">
  85 + Enter search terms or a module, class or function name.
  86 + </p>
  87 +</div>
  88 +<script type="text/javascript">$('#searchbox').show(0);</script>
  89 + </div>
  90 + </div>
  91 + <div class="clearer"></div>
  92 + </div>
  93 + <div class="related">
  94 + <h3>Navigation</h3>
  95 + <ul>
  96 + <li class="right" style="margin-right: 10px">
  97 + <a href="genindex.html" title="General Index"
  98 + >index</a></li>
  99 + <li><a href="index.html">documentation_HFR_analysis 0.0.1 documentation</a> &raquo;</li>
  100 + </ul>
  101 + </div>
  102 + <div class="footer">
  103 + &copy; Copyright 2015, Marion Bandet.
  104 + Created using <a href="http://sphinx-doc.org/">Sphinx</a> 1.2.3.
  105 + </div>
  106 + </body>
  107 +</html>
0 108 \ No newline at end of file
... ...
build/html/YSI_CastAway_CTD_May2014.html 0 → 100755
... ... @@ -0,0 +1,532 @@
  1 +
  2 +<!DOCTYPE html
  3 + PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
  4 +<html><head>
  5 + <meta http-equiv="Content-Type" content="text/html; charset=utf-8">
  6 + <!--
  7 +This HTML was auto-generated from MATLAB code.
  8 +To make changes, update the MATLAB code and republish this document.
  9 + --><title>YSI_CastAway_CTD_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-14"><meta name="DC.source" content="YSI_CastAway_CTD_May2014.m"><style type="text/css">
  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 +
  12 +html { min-height:100%; margin-bottom:1px; }
  13 +html body { height:100%; margin:0px; font-family:Arial, Helvetica, sans-serif; font-size:10px; color:#000; line-height:140%; background:#fff none; overflow-y:scroll; }
  14 +html body td { vertical-align:top; text-align:left; }
  15 +
  16 +h1 { padding:0px; margin:0px 0px 25px; font-family:Arial, Helvetica, sans-serif; font-size:1.5em; color:#d55000; line-height:100%; font-weight:normal; }
  17 +h2 { padding:0px; margin:0px 0px 8px; font-family:Arial, Helvetica, sans-serif; font-size:1.2em; color:#000; font-weight:bold; line-height:140%; border-bottom:1px solid #d6d4d4; display:block; }
  18 +h3 { padding:0px; margin:0px 0px 5px; font-family:Arial, Helvetica, sans-serif; font-size:1.1em; color:#000; font-weight:bold; line-height:140%; }
  19 +
  20 +a { color:#005fce; text-decoration:none; }
  21 +a:hover { color:#005fce; text-decoration:underline; }
  22 +a:visited { color:#004aa0; text-decoration:none; }
  23 +
  24 +p { padding:0px; margin:0px 0px 20px; }
  25 +img { padding:0px; margin:0px 0px 20px; border:none; }
  26 +p img, pre img, tt img, li img { margin-bottom:0px; }
  27 +
  28 +ul { padding:0px; margin:0px 0px 20px 23px; list-style:square; }
  29 +ul li { padding:0px; margin:0px 0px 7px 0px; }
  30 +ul li ul { padding:5px 0px 0px; margin:0px 0px 7px 23px; }
  31 +ul li ol li { list-style:decimal; }
  32 +ol { padding:0px; margin:0px 0px 20px 0px; list-style:decimal; }
  33 +ol li { padding:0px; margin:0px 0px 7px 23px; list-style-type:decimal; }
  34 +ol li ol { padding:5px 0px 0px; margin:0px 0px 7px 0px; }
  35 +ol li ol li { list-style-type:lower-alpha; }
  36 +ol li ul { padding-top:7px; }
  37 +ol li ul li { list-style:square; }
  38 +
  39 +.content { font-size:1.2em; line-height:140%; padding: 20px; }
  40 +
  41 +pre, tt, code { font-size:12px; }
  42 +pre { margin:0px 0px 20px; }
  43 +pre.error { color:red; }
  44 +pre.codeinput { padding:10px; border:1px solid #d3d3d3; background:#f7f7f7; }
  45 +pre.codeoutput { padding:10px 11px; margin:0px 0px 20px; color:#4c4c4c; }
  46 +
  47 +@media print { pre.codeinput, pre.codeoutput { word-wrap:break-word; width:100%; } }
  48 +
  49 +span.keyword { color:#0000FF }
  50 +span.comment { color:#228B22 }
  51 +span.string { color:#A020F0 }
  52 +span.untermstring { color:#B20000 }
  53 +span.syscmd { color:#B28C00 }
  54 +
  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 +.footer p { margin:0px; }
  57 +
  58 + </style></head><body><div class="content"><h2>Contents</h2><div><ul><li><a href="#2">place casts on a map of SL</a></li><li><a href="#3">compute mean profiles</a></li><li><a href="#4">plots</a></li></ul></div><pre class="codeinput"><span class="comment">% YSI_CastAway_CTD</span>
  59 +
  60 +clear <span class="string">all</span>
  61 +close <span class="string">all</span>
  62 +
  63 +addpath (genpath(<span class="string">'/home/bandma01/Dropbox/MEOPAR_Marion/Matlab/mfiles'</span>));
  64 +addpath (genpath(<span class="string">'/home/bandma01/Dropbox/UQAR/MEOPAR_Marion/Matlab_Analysis/'</span>));
  65 +addpath (genpath(<span class="string">'/home/bandma01/Dropbox/Matlab_Tools/'</span>));
  66 +
  67 +
  68 +<span class="comment">%cd('H:\MEOPAR\DATA\ADCP-V50')</span>
  69 +cd(<span class="string">'/run/media/bandma01/SeagateWireless/MEOPAR/DATA/field_data/2014_05_07/YSI_CastAway_CTD'</span>)
  70 +
  71 +d = dir(<span class="string">'*.mat'</span>);
  72 +
  73 +
  74 +<span class="comment">% load bathy</span>
  75 +cd(<span class="string">'/home/bandma01/hfr-data-processing'</span>)
  76 +load <span class="string">CHS_500m_gulf_grid.mat</span>
  77 +
  78 +<span class="comment">% load radars sites coordinates</span>
  79 +load <span class="string">sites_coord</span>
  80 +
  81 +
  82 +<span class="comment">% initialize matrices</span>
  83 +Press_down = NaN(800,length(d));
  84 +Temp_down = Press_down;
  85 +Cond_down = Press_down;
  86 +Sal_down = Press_down;
  87 +Dens_down = Press_down;
  88 +Press_up = Press_down;
  89 +Temp_up = Press_down;
  90 +Cond_up = Press_down;
  91 +Sal_up = Press_down;
  92 +Dens_up = Press_down;
  93 +
  94 +
  95 +<span class="keyword">for</span> i = 1:length(d)
  96 +
  97 + filename = d(i).name;
  98 +
  99 + load(filename)
  100 +
  101 + <span class="comment">% compute Conductivity ratio</span>
  102 + Cndr = Conductivity/1000/sw_c3515; <span class="comment">%[no units] conductivity in mS/cm</span>
  103 +
  104 + <span class="comment">% compute Salinity</span>
  105 + Salinity = sw_salt(Cndr,Temperature,Pressure);
  106 +
  107 + <span class="comment">% compute density</span>
  108 + Density = sw_dens(Salinity,Temperature,Pressure);
  109 +
  110 +
  111 +
  112 + <span class="comment">% get rid of data below 1m depth</span>
  113 + clear <span class="string">I</span>
  114 + Presssure_threshold = 1; <span class="comment">%[m]</span>
  115 + I = find(Pressure &lt; 1);
  116 + Pressure(I) = [];
  117 + Temperature(I) = [];
  118 + Conductivity(I) = [];
  119 + Salinity(I) = [];
  120 + Density(I) = [];
  121 +
  122 + <span class="comment">% compute means downcast and upcast</span>
  123 + clear <span class="string">J</span>
  124 + J = find(Pressure == max(abs(Pressure)));
  125 +
  126 + Press_down(1:J,i) = Pressure(1:J);
  127 + Temp_down(1:J,i) = Temperature(1:J);
  128 + Cond_down(1:J,i) = Conductivity(1:J);
  129 + Sal_down(1:J,i) = Salinity(1:J);
  130 + Dens_down(1:J,i) = Density(1:J);
  131 +
  132 + Press_up(1:length(Pressure)-J,i) = Pressure(J+1:end);
  133 + Temp_up(1:length(Pressure)-J,i) = Temperature(J+1:end);
  134 + Cond_up(1:length(Pressure)-J,i) = Conductivity(J+1:end);
  135 + Sal_up(1:length(Pressure)-J,i) = Salinity(J+1:end);
  136 + Dens_up(1:length(Pressure)-J,i) = Density(J+1:end);
  137 +
  138 +
  139 + <span class="comment">% compute drift</span>
  140 + clear <span class="string">E</span> <span class="string">W</span>
  141 + [E,W] = lonlat2km(LongitudeStart,LatitudeStart,LongitudeEnd,LatitudeEnd);
  142 + drift(i) = sqrt(E^2+W^2);
  143 +
  144 + Long_start(i) = LongitudeStart;
  145 + Lat_start(i) = LatitudeStart;
  146 +
  147 +<span class="keyword">end</span>
  148 +</pre><h2>place casts on a map of SL<a name="2"></a></h2><pre class="codeinput">mapSLE
  149 +
  150 +hold <span class="string">on</span>
  151 +m_plot(Long_start,Lat_start,<span class="string">'or'</span>,<span class="string">'MarkerFaceColor'</span>,<span class="string">'r'</span>)
  152 +</pre><img vspace="5" hspace="5" src="YSI_CastAway_CTD_May2014_01.png" alt=""> <h2>compute mean profiles<a name="3"></a></h2><pre class="codeinput">press_min = 1; <span class="comment">% [m]</span>
  153 +press_max = 100; <span class="comment">% [m]</span>
  154 +
  155 +press = press_min:.1:press_max;
  156 +
  157 +<span class="comment">% interpolate temperature and conductivity data</span>
  158 +<span class="keyword">for</span> i = 1:length(d)
  159 +
  160 + k = ~isnan(Press_down(:,i));
  161 + Tempdwn_interp(:,i) = interp1(Press_down(k,i),Temp_down(k,i),press');
  162 + Conddwn_interp(:,i) = interp1(Press_down(k,i),Cond_down(k,i),press');
  163 + Saldwn_interp(:,i) = interp1(Press_down(k,i),Sal_down(k,i),press');
  164 + Densdwn_interp(:,i) = interp1(Press_down(k,i),Dens_down(k,i),press');
  165 +
  166 + m = ~isnan(Press_up(:,i));
  167 + Tempup_interp(:,i) = interp1((Press_up(m,i)),(Temp_up(m,i)),press');
  168 + Condup_interp(:,i) = interp1(Press_up(m,i),Cond_up(m,i),press');
  169 + Salup_interp(:,i) = interp1(Press_up(m,i),Sal_up(m,i),press');
  170 + Densup_interp(:,i) = interp1(Press_up(m,i),Dens_up(m,i),press');
  171 +
  172 + <span class="comment">% figure(i+10)</span>
  173 + <span class="comment">% plot(Press_down,Temp_downp,'.-')</span>
  174 + <span class="comment">% hold on</span>
  175 + <span class="comment">% plot(press,Tempdwn_interp,'.r')</span>
  176 +
  177 + <span class="comment">% figure(i+20)</span>
  178 + <span class="comment">% plot(Press_up,Temp_up,'.-')</span>
  179 + <span class="comment">% hold on</span>
  180 + <span class="comment">% plot(press,Tempup_interp,'.r')</span>
  181 +
  182 +<span class="keyword">end</span>
  183 +</pre><h2>plots<a name="4"></a></h2><pre class="codeinput"><span class="comment">% temperature</span>
  184 +figure
  185 +subplot(121)
  186 +plot(Temp_down,Press_down,<span class="string">'-'</span>)
  187 +grid <span class="string">on</span>
  188 +set(gca,<span class="string">'Ydir'</span>,<span class="string">'reverse'</span>)
  189 +xlabel(<span class="string">'Temperature (deg C)'</span>)
  190 +ylabel(<span class="string">'Pressure (dbar)'</span>)
  191 +<span class="comment">%title(['CAST # ' num2str(i) ' - Starts at (' num2str(LongitudeStart) ',' num2str(LatitudeStart) ')'])</span>
  192 +title(<span class="string">'down casts'</span>)
  193 +
  194 +hold <span class="string">on</span>
  195 +plot(nanmean(Tempdwn_interp,2),press,<span class="string">'-k'</span>,<span class="string">'linewidth'</span>,2)
  196 +
  197 +subplot(122)
  198 +plot(Temp_up,Press_up,<span class="string">'-'</span>)
  199 +grid <span class="string">on</span>
  200 +set(gca,<span class="string">'Ydir'</span>,<span class="string">'reverse'</span>)
  201 +xlabel(<span class="string">'Temperature (deg C)'</span>)
  202 +ylabel(<span class="string">'Pressure (dbar)'</span>)
  203 +title(<span class="string">' up casts'</span>)
  204 +
  205 +hold <span class="string">on</span>
  206 +plot(nanmean(Tempup_interp,2),press,<span class="string">'-k'</span>,<span class="string">'linewidth'</span>,2)
  207 +
  208 +
  209 +
  210 +<span class="comment">% Conductivity</span>
  211 +figure
  212 +subplot(121)
  213 +plot(Cond_down,Press_down,<span class="string">'-'</span>)
  214 +grid <span class="string">on</span>
  215 +set(gca,<span class="string">'Ydir'</span>,<span class="string">'reverse'</span>)
  216 +xlabel(<span class="string">'Conductivity (S/m)'</span>)
  217 +ylabel(<span class="string">'Pressure (dbar)'</span>)
  218 +<span class="comment">%title(['CAST # ' num2str(i) ' - Starts at (' num2str(LongitudeStart) ',' num2str(LatitudeStart) ')'])</span>
  219 +title(<span class="string">'down casts'</span>)
  220 +
  221 +hold <span class="string">on</span>
  222 +plot(nanmean(Conddwn_interp,2),press,<span class="string">'-k'</span>,<span class="string">'linewidth'</span>,2)
  223 +
  224 +subplot(122)
  225 +plot(Cond_up,Press_up,<span class="string">'-'</span>)
  226 +grid <span class="string">on</span>
  227 +set(gca,<span class="string">'Ydir'</span>,<span class="string">'reverse'</span>)
  228 +xlabel(<span class="string">'Conductivity (S/m)'</span>)
  229 +ylabel(<span class="string">'Pressure (dbar)'</span>)
  230 +title(<span class="string">' up casts'</span>)
  231 +
  232 +hold <span class="string">on</span>
  233 +plot(nanmean(Condup_interp,2),press,<span class="string">'-k'</span>,<span class="string">'linewidth'</span>,2)
  234 +
  235 +<span class="comment">% Salinity</span>
  236 +figure
  237 +subplot(121)
  238 +plot(Sal_down,Press_down,<span class="string">'-'</span>)
  239 +grid <span class="string">on</span>
  240 +set(gca,<span class="string">'Ydir'</span>,<span class="string">'reverse'</span>)
  241 +xlabel(<span class="string">'Salinity (psu)'</span>)
  242 +ylabel(<span class="string">'Pressure (dbar)'</span>)
  243 +<span class="comment">%title(['CAST # ' num2str(i) ' - Starts at (' num2str(LongitudeStart) ',' num2str(LatitudeStart) ')'])</span>
  244 +title(<span class="string">'down casts'</span>)
  245 +
  246 +hold <span class="string">on</span>
  247 +plot(nanmean(Saldwn_interp,2),press,<span class="string">'-k'</span>,<span class="string">'linewidth'</span>,2)
  248 +
  249 +subplot(122)
  250 +plot(Sal_up,Press_up,<span class="string">'-'</span>)
  251 +grid <span class="string">on</span>
  252 +set(gca,<span class="string">'Ydir'</span>,<span class="string">'reverse'</span>)
  253 +xlabel(<span class="string">'Salinity (psu)'</span>)
  254 +ylabel(<span class="string">'Pressure (dbar)'</span>)
  255 +title(<span class="string">' up casts'</span>)
  256 +
  257 +hold <span class="string">on</span>
  258 +plot(nanmean(Salup_interp,2),press,<span class="string">'-k'</span>,<span class="string">'linewidth'</span>,2)
  259 +
  260 +
  261 +
  262 +<span class="comment">% Density</span>
  263 +figure
  264 +subplot(121)
  265 +plot(Dens_down,Press_down,<span class="string">'-'</span>)
  266 +grid <span class="string">on</span>
  267 +set(gca,<span class="string">'Ydir'</span>,<span class="string">'reverse'</span>)
  268 +xlabel(<span class="string">'Density (g/m^3)'</span>)
  269 +ylabel(<span class="string">'Pressure (dbar)'</span>)
  270 +<span class="comment">%title(['CAST # ' num2str(i) ' - Starts at (' num2str(LongitudeStart) ',' num2str(LatitudeStart) ')'])</span>
  271 +title(<span class="string">'down casts'</span>)
  272 +
  273 +hold <span class="string">on</span>
  274 +plot(nanmean(Densdwn_interp,2),press,<span class="string">'-k'</span>,<span class="string">'linewidth'</span>,2)
  275 +
  276 +subplot(122)