Commit f45a21332ece50b8c0d39e0dee01602ebb29dc1f

Authored by Marion Bandet
1 parent 516e2d17
Exists in master

update calibration outputs PAO and PAB

Showing 109 changed files with 798 additions and 617 deletions   Show diff stats

Too many changes.

To preserve performance only 100 of 109 files displayed.

build/doctrees/CurrentsCombination.doctree
No preview for this file type
build/doctrees/calibPAB.doctree 0 → 100644
No preview for this file type
build/doctrees/calibPAO.doctree 0 → 100644
No preview for this file type
build/doctrees/calibration.doctree
No preview for this file type
build/doctrees/environment.pickle
No preview for this file type
build/doctrees/index.doctree
No preview for this file type
build/html/WERA_processing.html deleted
... ... @@ -1,118 +0,0 @@
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>WERA POST=PROCESSING &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" href="CurrentsCombination.html" />
28   - <link rel="prev" title="CODAR post-processing" href="CODAR_processing.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"
39   - accesskey="N">next</a> |</li>
40   - <li class="right" >
41   - <a href="CODAR_processing.html" title="CODAR post-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="wera-post-processing">
53   -<h1>WERA POST=PROCESSING<a class="headerlink" href="#wera-post-processing" title="Permalink to this headline">¶</a></h1>
54   -<p>Post-processing of WERA radial currents is much simpler than for the CODAR system.</p>
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   -<p>In each folder we find .RFI and .SORT files for every hour (should find 24 files of each extention per day).</p>
57   -<p>Use the shell script do_process_SORT_pol to generate .crad_pol_lst files. They will be found in the subfolder YYYYJJJ/results/.</p>
58   -<p><a href="#id1"><span class="problematic" id="id2">**</span></a>** Use radiallst2uv.m to convert lst 2 mat files (To be verified).****</p>
59   -<p>To download radial currents from the sites:</p>
60   -<p>rsync -v -e ssh <a class="reference external" href="mailto:wera&#37;&#52;&#48;wera-lsle-pao&#46;dyndns&#46;org">wera<span>&#64;</span>wera-lsle-pao<span>&#46;</span>dyndns<span>&#46;</span>org</a>:/home/wera/data/ftp_transfer/${DATE}${HEURE}*_ascii* .</p>
61   -<p>with i.e DATE=2015121 (yyyyJJJ) and HEURE =0100</p>
62   -</div>
63   -
64   -
65   - </div>
66   - </div>
67   - </div>
68   - <div class="sphinxsidebar">
69   - <div class="sphinxsidebarwrapper">
70   - <h4>Previous topic</h4>
71   - <p class="topless"><a href="CODAR_processing.html"
72   - title="previous chapter">CODAR post-processing</a></p>
73   - <h4>Next topic</h4>
74   - <p class="topless"><a href="CurrentsCombination.html"
75   - title="next chapter">Current combination</a></p>
76   - <h3>This Page</h3>
77   - <ul class="this-page-menu">
78   - <li><a href="_sources/WERA_processing.txt"
79   - rel="nofollow">Show Source</a></li>
80   - </ul>
81   -<div id="searchbox" style="display: none">
82   - <h3>Quick search</h3>
83   - <form class="search" action="search.html" method="get">
84   - <input type="text" name="q" />
85   - <input type="submit" value="Go" />
86   - <input type="hidden" name="check_keywords" value="yes" />
87   - <input type="hidden" name="area" value="default" />
88   - </form>
89   - <p class="searchtip" style="font-size: 90%">
90   - Enter search terms or a module, class or function name.
91   - </p>
92   -</div>
93   -<script type="text/javascript">$('#searchbox').show(0);</script>
94   - </div>
95   - </div>
96   - <div class="clearer"></div>
97   - </div>
98   - <div class="related">
99   - <h3>Navigation</h3>
100   - <ul>
101   - <li class="right" style="margin-right: 10px">
102   - <a href="genindex.html" title="General Index"
103   - >index</a></li>
104   - <li class="right" >
105   - <a href="CurrentsCombination.html" title="Current combination"
106   - >next</a> |</li>
107   - <li class="right" >
108   - <a href="CODAR_processing.html" title="CODAR post-processing"
109   - >previous</a> |</li>
110   - <li><a href="index.html">documentation_HFR_analysis 0.0.1 documentation</a> &raquo;</li>
111   - </ul>
112   - </div>
113   - <div class="footer">
114   - &copy; Copyright 2015, Marion Bandet.
115   - Created using <a href="http://sphinx-doc.org/">Sphinx</a> 1.2.3.
116   - </div>
117   - </body>
118   -</html>
119 0 \ No newline at end of file
build/html/WERAcalibration.html
... ... @@ -6,7 +6,7 @@
6 6 <!--
7 7 This HTML was auto-generated from MATLAB code.
8 8 To make changes, update the MATLAB code and republish this document.
9   - --><title>WERAcalibration</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="2015-10-06"><meta name="DC.source" content="WERAcalibration.m"><style type="text/css">
  9 + --><title>WERAcalibration</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="2015-10-07"><meta name="DC.source" content="WERAcalibration.m"><style type="text/css">
10 10 html,body,div,span,applet,object,iframe,h1,h2,h3,h4,h5,h6,p,blockquote,pre,a,abbr,acronym,address,big,cite,code,del,dfn,em,font,img,ins,kbd,q,s,samp,small,strike,strong,sub,sup,tt,var,b,u,i,center,dl,dt,dd,ol,ul,li,fieldset,form,label,legend,table,caption,tbody,tfoot,thead,tr,th,td{margin:0;padding:0;border:0;outline:0;font-size:100%;vertical-align:baseline;background:transparent}body{line-height:1}ol,ul{list-style:none}blockquote,q{quotes:none}blockquote:before,blockquote:after,q:before,q:after{content:'';content:none}:focus{outine:0}ins{text-decoration:none}del{text-decoration:line-through}table{border-collapse:collapse;border-spacing:0}
11 11  
12 12 html { min-height:100%; margin-bottom:1px; }
... ... @@ -55,19 +55,21 @@ span.syscmd { color:#B28C00 }
55 55 .footer { width:auto; padding:10px 0px; margin:25px 0px 0px; border-top:1px dotted #878787; font-size:0.8em; line-height:140%; font-style:italic; color:#878787; text-align:left; float:none; }
56 56 .footer p { margin:0px; }
57 57  
58   - </style></head><body><div class="content"><h2>Contents</h2><div><ul><li><a href="#3">phases and time of calibration</a></li><li><a href="#4">separate between calibrations made with Helzel's system and Pierre's system</a></li><li><a href="#5">theoretical values</a></li></ul></div><pre class="codeinput">clear <span class="string">all</span>
  58 + </style></head><body><div class="content"><pre class="codeinput">clear <span class="string">all</span>
59 59 close <span class="string">all</span>
60 60  
61 61 <span class="comment">%cd('/home/bandma01/Documents/WERAs/')</span>
62 62 cd(<span class="string">'/home/bandma01/hfr-data-processing/calibrations/'</span>)
63 63  
64   -<span class="comment">% site = 'PAO';</span>
65   -site = <span class="string">'PAB'</span>;
  64 +site = <span class="string">'PAO'</span>;
  65 +<span class="comment">% site = 'PAB';</span>
  66 +
66 67 cd(site)
67 68  
68 69 <span class="comment">% system</span>
69 70 <span class="keyword">for</span> sys = 1:2
70   -</pre><pre class="codeinput"> <span class="keyword">if</span> sys == 1
  71 +
  72 + <span class="keyword">if</span> sys == 1
71 73 systm = <span class="string">'Helzel'</span>;
72 74 <span class="keyword">elseif</span> sys == 2
73 75 systm = <span class="string">'Flament'</span>;
... ... @@ -90,7 +92,10 @@ cd(site)
90 92 <span class="keyword">end</span>
91 93  
92 94 <span class="keyword">end</span>
93   -</pre><h2>phases and time of calibration<a name="3"></a></h2><pre class="codeinput"> <span class="keyword">if</span> strcmp(site,<span class="string">'PAO'</span>) == 1
  95 +
  96 +
  97 + <span class="comment">% phases and time of calibration</span>
  98 + <span class="keyword">if</span> strcmp(site,<span class="string">'PAO'</span>) == 1
94 99 cd(<span class="string">'2014302'</span>)
95 100 <span class="keyword">elseif</span> strcmp(site,<span class="string">'PAB'</span>) == 1
96 101 cd(<span class="string">'2014303'</span>)
... ... @@ -115,6 +120,8 @@ cd(site)
115 120 <span class="keyword">end</span>
116 121 clear <span class="string">phase</span>
117 122  
  123 + <span class="comment">% for each time stamp (column) compute phase difference between 2</span>
  124 + <span class="comment">% consecutive antennas (row)</span>
118 125 delta_phase_obs = modsym(diff(Phase,1),180);
119 126  
120 127 To = datenum(2014,01,01,0,0,0);
... ... @@ -123,30 +130,37 @@ cd(site)
123 130 calib_time = To + calib_jd -1 + calib_hr./24 + calib_min/24/60 + 20/24/3600;
124 131 <span class="comment">% datevec(calib_time)</span>
125 132  
  133 + clear <span class="string">calib_year</span> <span class="string">calib_jd</span> <span class="string">calib_hr</span> <span class="string">calib_min</span>
  134 +
126 135 <span class="keyword">if</span> strcmp(site,<span class="string">'PAO'</span>) == 1
127 136 <span class="comment">% calib time 2min34sec late wrt GPS</span>
128 137 calib_time = calib_time + 2/24/60 + 34/24/3600;
129 138 <span class="keyword">end</span>
130   -</pre><h2>separate between calibrations made with Helzel's system and Pierre's system<a name="4"></a></h2><p>date of new calibration files Helzel 2014-10-30 16:52 to 17:54 Flament 2014-10-30 17:55 to 18:20</p><pre class="codeinput"> T1 = datenum(2014,10,30,17,59,00);
  139 + disp([<span class="string">'calibrations started at '</span> datestr(calib_time(1)) <span class="string">' and ended at '</span> datestr(calib_time(end))])
131 140  
  141 + <span class="comment">% separate between calibrations made with Helzel's system and Flament's system</span>
  142 + <span class="comment">% date of new calibration files</span>
  143 + <span class="comment">% Helzel 2014-10-30 16:52 to 17:54</span>
  144 + <span class="comment">% Flament 2014-10-30 17:55 to 18:20</span>
  145 + T1 = datenum(2014,10,30,17,55,00);
132 146  
133   - clear <span class="string">J</span>
134   - <span class="keyword">if</span> strcmp(systm,<span class="string">'Helzel'</span>) == 1
135   - J = find(calib_time &gt; T1);
136   - <span class="keyword">elseif</span> strcmp(systm,<span class="string">'Flament'</span>) == 1
137   - J = find(calib_time &gt; T1);
138   - <span class="keyword">end</span>
  147 + clear <span class="string">J</span>
  148 + <span class="keyword">if</span> strcmp(systm,<span class="string">'Helzel'</span>) == 1
  149 + J = find(calib_time &gt; T1);
  150 + <span class="keyword">elseif</span> strcmp(systm,<span class="string">'Flament'</span>) == 1
  151 + J = find(calib_time &lt; T1);
  152 + <span class="keyword">end</span>
139 153  
140   - <span class="comment">%remove unwanted data to avoid confusion!</span>
141   - calib_year(J)=[];
142   - calib_jd(J)=[];
143   - calib_hr(J)=[];
144   - calib_min(J)=[];
145   - Phase(:,J)=[];
146   - delta_phase_obs(:,J)=[];
147   - calib_time(J)=[];
148   - d(J)=[];
149   - clear <span class="string">c</span> <span class="string">Nfiles</span>
  154 + <span class="comment">% remove unwanted data to avoid confusion!</span>
  155 +<span class="comment">% calib_year(J)=[];</span>
  156 +<span class="comment">% calib_jd(J)=[];</span>
  157 +<span class="comment">% calib_hr(J)=[];</span>
  158 +<span class="comment">% calib_min(J)=[];</span>
  159 + Phase(:,J)=[];
  160 + delta_phase_obs(:,J)=[];
  161 + calib_time(J)=[];
  162 + d(J)=[];
  163 + clear <span class="string">c</span> <span class="string">Nfiles</span>
150 164  
151 165 <span class="comment">%+++++++++++++++++++++++++++++++++++++++++</span>
152 166 <span class="comment">% for each calibration file (series# 1)</span>
... ... @@ -190,8 +204,16 @@ cd(site)
190 204 hold <span class="string">on</span>
191 205 plot(lon,lat,<span class="string">'x'</span>)
192 206 grid <span class="string">on</span>
  207 + legend(<span class="string">'gps'</span>,<span class="string">'antennas'</span>,<span class="string">'location'</span>,<span class="string">'southeast'</span>)
  208 + xlabel(<span class="string">'longitude'</span>)
  209 + ylabel(<span class="string">'latitude'</span>)
193 210 title([<span class="string">'site = '</span> site ])
194   -</pre><img vspace="5" hspace="5" src="WERAcalibration_01.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_16.png" alt=""> <h2>theoretical values<a name="5"></a></h2><p>wavelength</p><pre class="codeinput"> lambda = 3E8/16.15E6;
  211 +
  212 +
  213 +
  214 + <span class="comment">% compute theoretical phase values</span>
  215 + <span class="comment">% wavelength</span>
  216 + lambda = 3E8/16.15E6;
195 217  
196 218  
197 219 <span class="comment">% initialize variables</span>
... ... @@ -209,8 +231,7 @@ cd(site)
209 231 <span class="keyword">end</span>
210 232  
211 233  
212   -
213   - <span class="comment">%for each pair of antennas</span>
  234 + <span class="comment">% for each pair of antennas</span>
214 235 <span class="keyword">for</span> p = 1:11
215 236  
216 237 <span class="comment">%delta_r</span>
... ... @@ -248,14 +269,8 @@ cd(site)
248 269 theta_pairs(p,pt) = modsym(zeta(p)-phi(p)+90,180);
249 270 <span class="comment">%rad2deg(theta_pairs)</span>
250 271 <span class="keyword">end</span>
251   -
252   -
253 272 <span class="keyword">end</span>
254 273  
255   - figure(1)
256   - legend(<span class="string">'gps'</span>,<span class="string">'antennas'</span>,<span class="string">'location'</span>,<span class="string">'southeast'</span>)
257   - xlabel(<span class="string">'longitude'</span>)
258   - ylabel(<span class="string">'latitude'</span>)
259 274  
260 275  
261 276 <span class="comment">% for each pair of antenna</span>
... ... @@ -268,8 +283,6 @@ cd(site)
268 283 xlabel(<span class="string">'incidence (degrees)'</span>)
269 284 ylabel(<span class="string">'phase difference (degrees)'</span>)
270 285 legend(<span class="string">'theoretical'</span>,<span class="string">'observed'</span>,<span class="string">'location'</span>,<span class="string">'best'</span>)
271   - title([site <span class="string">' - antennas '</span> num2str(p) <span class="string">'-'</span> num2str(p+1)])
272   -
273 286 <span class="keyword">if</span> strcmp(systm,<span class="string">'Helzel'</span>) == 1
274 287 set(h2,<span class="string">'color'</span>,<span class="string">'b'</span>)
275 288 set(h3,<span class="string">'color'</span>,<span class="string">'r'</span>,<span class="string">'marker'</span>,<span class="string">'o'</span>)
... ... @@ -277,6 +290,7 @@ cd(site)
277 290 set(h2,<span class="string">'color'</span>,<span class="string">'k'</span>,<span class="string">'marker'</span>,<span class="string">'*'</span>)
278 291 set(h3,<span class="string">'color'</span>,<span class="string">'r'</span>,<span class="string">'marker'</span>,<span class="string">'sq'</span>,<span class="string">'markerfacecolor'</span>,<span class="string">'r'</span>)
279 292 <span class="keyword">end</span>
  293 + title([site <span class="string">' - antennas '</span> num2str(p) <span class="string">'-'</span> num2str(p+1)])
280 294  
281 295  
282 296 <span class="comment">% compute calibration value for antennas pairs</span>
... ... @@ -288,9 +302,9 @@ cd(site)
288 302  
289 303 <span class="keyword">end</span>
290 304  
291   - P=ones(16,3);
292   - P(13:16,:)=0;
293   - P(1:12,3)=phi_ant(:,1);
  305 + P = ones(16,3);
  306 + P(13:16,:) = 0;
  307 + P(1:12,3) = phi_ant(:,1);
294 308  
295 309  
296 310 <span class="comment">% save file (we only keep results from Helzel</span>
... ... @@ -303,30 +317,30 @@ cd(site)
303 317  
304 318 figure(13)
305 319 hold <span class="string">on</span>
306   - h4=plot(calib_value,<span class="string">'.-'</span>);
  320 + h4 = plot(calib_value,<span class="string">'.-'</span>);
307 321 grid <span class="string">on</span>
308 322 xlabel(<span class="string">'pairs of antennas'</span>)
309   - ylabel({<span class="string">'phase difference between'</span>,<span class="string">'theoretical and observed values (degrees)'</span>})
310   - title(site)
311   -
  323 + ylabel([<span class="string">'phase differences between theoretical'</span>,10 ,<span class="string">'and observed values (degrees)'</span>])
312 324 <span class="keyword">if</span> strcmp(systm,<span class="string">'Helzel'</span>) == 1
313 325 set(h4,<span class="string">'color'</span>,<span class="string">'b'</span>)
  326 + legend(<span class="string">'Helzel'</span>)
314 327 <span class="keyword">elseif</span> strcmp(systm,<span class="string">'Flament'</span>) == 1
315 328 set(h4,<span class="string">'color'</span>,<span class="string">'k'</span>)
316 329 legend(<span class="string">'Helzel'</span>,<span class="string">'Flament'</span>,<span class="string">'location'</span>,<span class="string">'southeast'</span>)
317 330 <span class="keyword">end</span>
  331 + title(site)
318 332  
319   -
320   - figure(14)
321   - hold <span class="string">on</span>
322   - h5=plot([1:12],phi_ant(:,sys),<span class="string">'r'</span>);
323   - grid <span class="string">on</span>
324   - <span class="keyword">if</span> strcmp(systm,<span class="string">'Helzel'</span>) == 1
325   - set(h5,<span class="string">'marker'</span>,<span class="string">'o'</span>)
326   - <span class="keyword">elseif</span> strcmp(systm,<span class="string">'Flament'</span>) == 1
327   - set(h5,<span class="string">'marker'</span>,<span class="string">'sq'</span>,<span class="string">'markerfacecolor'</span>,<span class="string">'r'</span>)
328   - legend(<span class="string">'Helzel'</span>,<span class="string">'Flament'</span>,<span class="string">'location'</span>,<span class="string">'southeast'</span>)
329   - <span class="keyword">end</span>
  333 + figure(14)
  334 + hold <span class="string">on</span>
  335 + h5 = plot([1:12],phi_ant(:,sys),<span class="string">'r'</span>);
  336 + grid <span class="string">on</span>
  337 + <span class="keyword">if</span> strcmp(systm,<span class="string">'Helzel'</span>) == 1
  338 + set(h5,<span class="string">'marker'</span>,<span class="string">'o'</span>)
  339 + legend(<span class="string">'Helzel'</span>)
  340 + <span class="keyword">elseif</span> strcmp(systm,<span class="string">'Flament'</span>) == 1
  341 + set(h5,<span class="string">'marker'</span>,<span class="string">'sq'</span>,<span class="string">'markerfacecolor'</span>,<span class="string">'r'</span>)
  342 + legend(<span class="string">'Helzel'</span>,<span class="string">'Flament'</span>,<span class="string">'location'</span>,<span class="string">'southeast'</span>)
  343 + <span class="keyword">end</span>
330 344 xlabel(<span class="string">'antenna number'</span>)
331 345 ylabel(<span class="string">'calibration phase to be applied (degrees)'</span>)
332 346 title(site)
... ... @@ -334,9 +348,10 @@ cd(site)
334 348 <span class="keyword">if</span> strcmp(site,<span class="string">'PAO'</span>) == 1
335 349 <span class="keyword">break</span>
336 350 <span class="keyword">end</span>
337   -</pre><img vspace="5" hspace="5" src="WERAcalibration_02.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_03.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_04.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_05.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_06.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_07.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_08.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_09.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_10.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_11.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_12.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_13.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_14.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_15.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_17.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_18.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_19.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_20.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_21.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_22.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_23.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_24.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_25.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_26.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_27.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_28.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_29.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_30.png" alt=""> <pre class="codeinput"><span class="keyword">end</span>
338 351  
339   -<span class="comment">%add current calibration phases</span>
  352 +<span class="keyword">end</span>
  353 +
  354 +<span class="comment">% add current calibration phases</span>
340 355 <span class="keyword">if</span> strcmp(site,<span class="string">'PAO'</span>) == 1
341 356 load <span class="string">calibration_pao.wera</span>
342 357 figure(14)
... ... @@ -414,7 +429,8 @@ fig_hdl = figure(14);
414 429 print(fig_hdl,<span class="string">'-dpng'</span>,<span class="string">'-r500'</span>,fname)
415 430 <span class="comment">% crop([fname '.png'])</span>
416 431 <span class="keyword">end</span>
417   -</pre><p class="footer"><br>
  432 +</pre><pre class="codeoutput">calibrations started at 29-Oct-2014 17:48:54 and ended at 29-Oct-2014 18:25:54
  433 +</pre><img vspace="5" hspace="5" src="WERAcalibration_01.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_02.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_03.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_04.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_05.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_06.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_07.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_08.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_09.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_10.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_11.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_12.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_13.png" alt=""> <img vspace="5" hspace="5" src="WERAcalibration_14.png" alt=""> <p class="footer"><br>
418 434 Published with MATLAB&reg; 7.14<br></p></div><!--
419 435 ##### SOURCE BEGIN #####
420 436 clear all
... ... @@ -423,8 +439,9 @@ close all
423 439 %cd('/home/bandma01/Documents/WERAs/')
424 440 cd('/home/bandma01/hfr-data-processing/calibrations/')
425 441  
426   -% site = 'PAO';
427   -site = 'PAB';
  442 +site = 'PAO';
  443 +% site = 'PAB';
  444 +
428 445 cd(site)
429 446  
430 447 % system
... ... @@ -455,7 +472,7 @@ for sys = 1:2
455 472 end
456 473  
457 474  
458   - %% phases and time of calibration
  475 + % phases and time of calibration
459 476 if strcmp(site,'PAO') == 1
460 477 cd('2014302')
461 478 elseif strcmp(site,'PAB') == 1
... ... @@ -481,6 +498,8 @@ for sys = 1:2
481 498 end
482 499 clear phase
483 500  
  501 + % for each time stamp (column) compute phase difference between 2
  502 + % consecutive antennas (row)
484 503 delta_phase_obs = modsym(diff(Phase,1),180);
485 504  
486 505 To = datenum(2014,01,01,0,0,0);
... ... @@ -489,35 +508,37 @@ for sys = 1:2
489 508 calib_time = To + calib_jd -1 + calib_hr./24 + calib_min/24/60 + 20/24/3600;
490 509 % datevec(calib_time)
491 510  
  511 + clear calib_year calib_jd calib_hr calib_min
  512 +
492 513 if strcmp(site,'PAO') == 1
493 514 % calib time 2min34sec late wrt GPS
494 515 calib_time = calib_time + 2/24/60 + 34/24/3600;
495 516 end
  517 + disp(['calibrations started at ' datestr(calib_time(1)) ' and ended at ' datestr(calib_time(end))])
496 518  
497   - %% separate between calibrations made with Helzel's system and Pierre's system
  519 + % separate between calibrations made with Helzel's system and Flament's system
498 520 % date of new calibration files
499 521 % Helzel 2014-10-30 16:52 to 17:54
500 522 % Flament 2014-10-30 17:55 to 18:20
501   - T1 = datenum(2014,10,30,17,59,00);
  523 + T1 = datenum(2014,10,30,17,55,00);
502 524  
  525 + clear J
  526 + if strcmp(systm,'Helzel') == 1
  527 + J = find(calib_time > T1);
  528 + elseif strcmp(systm,'Flament') == 1
  529 + J = find(calib_time < T1);
  530 + end
503 531  
504   - clear J
505   - if strcmp(systm,'Helzel') == 1
506   - J = find(calib_time > T1);
507   - elseif strcmp(systm,'Flament') == 1
508   - J = find(calib_time > T1);
509   - end
510   -
511   - %remove unwanted data to avoid confusion!
512   - calib_year(J)=[];
513   - calib_jd(J)=[];
514   - calib_hr(J)=[];
515   - calib_min(J)=[];
516   - Phase(:,J)=[];
517   - delta_phase_obs(:,J)=[];
518   - calib_time(J)=[];
519   - d(J)=[];
520   - clear c Nfiles
  532 + % remove unwanted data to avoid confusion!
  533 +% calib_year(J)=[];
  534 +% calib_jd(J)=[];
  535 +% calib_hr(J)=[];
  536 +% calib_min(J)=[];
  537 + Phase(:,J)=[];
  538 + delta_phase_obs(:,J)=[];
  539 + calib_time(J)=[];
  540 + d(J)=[];
  541 + clear c Nfiles
521 542  
522 543 %+++++++++++++++++++++++++++++++++++++++++
523 544 % for each calibration file (series# 1)
... ... @@ -561,10 +582,14 @@ for sys = 1:2
561 582 hold on
562 583 plot(lon,lat,'x')
563 584 grid on
  585 + legend('gps','antennas','location','southeast')
  586 + xlabel('longitude')
  587 + ylabel('latitude')
564 588 title(['site = ' site ])
  589 +
565 590  
566 591  
567   - %% theoretical values
  592 + % compute theoretical phase values
568 593 % wavelength
569 594 lambda = 3E8/16.15E6;
570 595  
... ... @@ -583,9 +608,8 @@ for sys = 1:2
583 608 r(ant) = ddist*1000; %meters
584 609 end
585 610  
586   -
587   -
588   - %for each pair of antennas
  611 +
  612 + % for each pair of antennas
589 613 for p = 1:11
590 614  
591 615 %delta_r
... ... @@ -623,14 +647,8 @@ for sys = 1:2
623 647 theta_pairs(p,pt) = modsym(zeta(p)-phi(p)+90,180);
624 648 %rad2deg(theta_pairs)
625 649 end
626   -
627   -
628 650 end
629 651  
630   - figure(1)
631   - legend('gps','antennas','location','southeast')
632   - xlabel('longitude')
633   - ylabel('latitude')
634 652  
635 653  
636 654 % for each pair of antenna
... ... @@ -643,8 +661,6 @@ for sys = 1:2
643 661 xlabel('incidence (degrees)')
644 662 ylabel('phase difference (degrees)')
645 663 legend('theoretical','observed','location','best')
646   - title([site ' - antennas ' num2str(p) '-' num2str(p+1)])
647   -
648 664 if strcmp(systm,'Helzel') == 1
649 665 set(h2,'color','b')
650 666 set(h3,'color','r','marker','o')
... ... @@ -652,7 +668,8 @@ for sys = 1:2
652 668 set(h2,'color','k','marker','*')
653 669 set(h3,'color','r','marker','sq','markerfacecolor','r')
654 670 end
655   -
  671 + title([site ' - antennas ' num2str(p) '-' num2str(p+1)])
  672 +
656 673  
657 674 % compute calibration value for antennas pairs
658 675 calib_value(p) = nanmean(modsym(-delta_phase_obs(p,:)-delta_phase_th(p,:),180),2);
... ... @@ -663,9 +680,9 @@ for sys = 1:2
663 680  
664 681 end
665 682  
666   - P=ones(16,3);
667   - P(13:16,:)=0;
668   - P(1:12,3)=phi_ant(:,1);
  683 + P = ones(16,3);
  684 + P(13:16,:) = 0;
  685 + P(1:12,3) = phi_ant(:,1);
669 686  
670 687  
671 688 % save file (we only keep results from Helzel
... ... @@ -678,30 +695,30 @@ for sys = 1:2
678 695  
679 696 figure(13)
680 697 hold on
681   - h4=plot(calib_value,'.-');
  698 + h4 = plot(calib_value,'.-');
682 699 grid on
683 700 xlabel('pairs of antennas')
684   - ylabel({'phase difference between','theoretical and observed values (degrees)'})
685   - title(site)
686   -
  701 + ylabel(['phase differences between theoretical',10 ,'and observed values (degrees)'])
687 702 if strcmp(systm,'Helzel') == 1
688 703 set(h4,'color','b')
  704 + legend('Helzel')
689 705 elseif strcmp(systm,'Flament') == 1
690 706 set(h4,'color','k')
691 707 legend('Helzel','Flament','location','southeast')
692 708 end
  709 + title(site)
693 710  
694   -
695   - figure(14)
696   - hold on
697   - h5=plot([1:12],phi_ant(:,sys),'r');
698   - grid on
699   - if strcmp(systm,'Helzel') == 1
700   - set(h5,'marker','o')
701   - elseif strcmp(systm,'Flament') == 1
702   - set(h5,'marker','sq','markerfacecolor','r')
703   - legend('Helzel','Flament','location','southeast')
704   - end
  711 + figure(14)
  712 + hold on
  713 + h5 = plot([1:12],phi_ant(:,sys),'r');
  714 + grid on
  715 + if strcmp(systm,'Helzel') == 1
  716 + set(h5,'marker','o')
  717 + legend('Helzel')
  718 + elseif strcmp(systm,'Flament') == 1
  719 + set(h5,'marker','sq','markerfacecolor','r')
  720 + legend('Helzel','Flament','location','southeast')
  721 + end
705 722 xlabel('antenna number')
706 723 ylabel('calibration phase to be applied (degrees)')
707 724 title(site)
... ... @@ -712,7 +729,7 @@ for sys = 1:2
712 729  
713 730 end
714 731  
715   -%add current calibration phases
  732 +% add current calibration phases
716 733 if strcmp(site,'PAO') == 1
717 734 load calibration_pao.wera
718 735 figure(14)
... ...
build/html/WERAcalibration.png deleted

3.44 KB

build/html/WERAcalibration_01.png deleted

8.59 KB

build/html/WERAcalibration_02.png deleted

9.03 KB

build/html/WERAcalibration_03.png deleted

9.47 KB

build/html/WERAcalibration_04.png deleted

9.4 KB

build/html/WERAcalibration_05.png deleted

9.67 KB

build/html/WERAcalibration_06.png deleted

9.53 KB

build/html/WERAcalibration_07.png deleted

9.45 KB

build/html/WERAcalibration_08.png deleted

9.4 KB

build/html/WERAcalibration_09.png deleted

9.36 KB

build/html/WERAcalibration_10.png deleted

9.25 KB

build/html/WERAcalibration_11.png deleted

9.24 KB

build/html/WERAcalibration_12.png deleted

9.22 KB

build/html/WERAcalibration_13.png deleted

9 KB

build/html/WERAcalibration_14.png deleted

9.03 KB

build/html/WERAcalibration_15.png deleted

7.22 KB

build/html/WERAcalibration_16.png deleted

9.07 KB

build/html/WERAcalibration_17.png deleted

9.07 KB

build/html/WERAcalibration_18.png deleted

9.59 KB

build/html/WERAcalibration_19.png deleted

9.53 KB

build/html/WERAcalibration_20.png deleted

9.59 KB

build/html/WERAcalibration_21.png deleted

9.75 KB

build/html/WERAcalibration_22.png deleted

9.5 KB

build/html/WERAcalibration_23.png deleted

9.47 KB

build/html/WERAcalibration_24.png deleted

9.46 KB

build/html/WERAcalibration_25.png deleted

9.36 KB

build/html/WERAcalibration_26.png deleted