NCCOOS Trac Projects: Top | Web | Platforms | Processing | Viz | Sprints | Sandbox | (Wind)

root/gliderproc/trunk/pelagia_DO_hysteresis_plots.m

Revision 516 (checked in by cbc, 10 years ago)

Add green bar as indicator of sample time to glider DO hysteresis plots.

Line 
1 clear;
2 load('GLIDER_DATA_LEVEL1\Pelagia_Deployment1_DO_L1.mat');
3 figure('Name', 'Pelagia Dissolved Oxygen Hysteresis Analysis', 'NumberTitle', 'off');
4
5 start = 250174;
6 middle = 250290;
7 stop = 250401;
8 subplot(3, 4, 1);
9 [ax, ~, h2] = plotyy(ptime_ebd_datenum, -depthi, ptime_ebd_datenum, o2_sat);
10 datetick(ax(1), 'x', 6, 'keeplimits', 'keepticks');
11 datetick(ax(2), 'x', 6, 'keeplimits', 'keepticks');
12 set(h2, 'Color', 'r');
13 set(ax(2), 'YColor', 'r');
14 title({'Pelagia Deployment 1'; 'Depth and O2 Saturation'; '\color[rgb]{0 .5 0}Sample Time as Green Bar'});
15 xlabel('UTC Time');
16 ylabel(ax(1), 'Depth (m)');
17 ylabel(ax(2), 'O2 Saturation (%)');
18 axes(ax(2));
19 yl = get(gca, 'YLim');
20 line([ptime_ebd_datenum(middle) ptime_ebd_datenum(middle)], [yl(1) yl(2)], 'LineStyle', '-', 'Color', 'green', 'LineWidth', 2);
21
22 subplot(3, 4, 2);
23 [ax, ~, h2] = plotyy(ptime_ebd_datenum(start:stop), -depthi(start:stop), ptime_ebd_datenum(start:stop), o2_sat(start:stop));
24 datetick(ax(1), 'x', 15, 'keeplimits', 'keepticks');
25 datetick(ax(2), 'x', 15, 'keeplimits', 'keepticks');
26 set(h2, 'Color', 'r');
27 set(ax(2), 'YColor', 'r');
28 title({'Pelagia Deployment 1'; 'Depth and O2 Saturation'; ['on ', datestr(ptime_ebd_datenum(start), 1)]});
29 xlabel('UTC Time');
30 ylabel(ax(1), 'Depth (m)');
31 ylabel(ax(2), 'O2 Saturation (%)');
32
33 subplot(3, 4, 3);
34 plot(o2_sat(start:middle), -depthi(start:middle), '.');
35 hold on;
36 plot(o2_sat(middle:stop), -depthi(middle:stop), 'r.');
37 title({'Pelagia Deployment 1'; 'Depth and O2 Saturation'; '\color{blue}Downcast = Blue \color{black}/ \color{red}Upcast = Red'});
38 xlabel('O2 Saturation (%)');
39 ylabel('Depth (m)');
40
41 subplot(3, 4, 4);
42 plot(o2_sat(start:middle), tempi(start:middle), '.');
43 hold on;
44 plot(o2_sat(middle:stop), tempi(middle:stop), 'r.');
45 title({'Pelagia Deployment 1'; 'Temp and O2 Saturation'; '\color{blue}Downcast = Blue \color{black}/ \color{red}Upcast = Red'});
46 xlabel('O2 Saturation (%)');
47 ylabel('Temp (Degrees C)');
48
49 clear;
50 load('GLIDER_DATA_LEVEL1\Pelagia_Deployment2_DO_L1.mat');
51
52 start = 303249;
53 middle = 303491;
54 stop = 303678;
55 subplot(3, 4, 5);
56 [ax, ~, h2] = plotyy(ptime_ebd_datenum, -depthi, ptime_ebd_datenum, o2_sat);
57 datetick(ax(1), 'x', 6, 'keeplimits', 'keepticks');
58 datetick(ax(2), 'x', 6, 'keeplimits', 'keepticks');
59 set(h2, 'Color', 'r');
60 set(ax(2), 'YColor', 'r');
61 title({'Pelagia Deployment 2'; 'Depth and O2 Saturation'; '\color[rgb]{0 .5 0}Sample Time as Green Bar'});
62 xlabel('UTC Time');
63 ylabel(ax(1), 'Depth (m)');
64 ylabel(ax(2), 'O2 Saturation (%)');
65 axes(ax(2));
66 yl = get(gca, 'YLim');
67 line([ptime_ebd_datenum(middle) ptime_ebd_datenum(middle)], [yl(1) yl(2)], 'LineStyle', '-', 'Color', 'green', 'LineWidth', 2);
68
69 subplot(3, 4, 6);
70 [ax, ~, h2] = plotyy(ptime_ebd_datenum(start:stop), -depthi(start:stop), ptime_ebd_datenum(start:stop), o2_sat(start:stop));
71 datetick(ax(1), 'x', 15, 'keeplimits', 'keepticks');
72 datetick(ax(2), 'x', 15, 'keeplimits', 'keepticks');
73 set(h2, 'Color', 'r');
74 set(ax(2), 'YColor', 'r');
75 title({'Pelagia Deployment 2'; 'Depth and O2 Saturation'; ['on ', datestr(ptime_ebd_datenum(start), 1)]});
76 xlabel('UTC Time');
77 ylabel(ax(1), 'Depth (m)');
78 ylabel(ax(2), 'O2 Saturation (%)');
79
80 subplot(3, 4, 7);
81 plot(o2_sat(start:middle), -depthi(start:middle), '.');
82 hold on;
83 plot(o2_sat(middle:stop), -depthi(middle:stop), 'r.');
84 title({'Pelagia Deployment 2'; 'Depth and O2 Saturation'; '\color{blue}Downcast = Blue \color{black}/ \color{red}Upcast = Red'});
85 xlabel('O2 Saturation (%)');
86 ylabel('Depth (m)');
87
88 subplot(3, 4, 8);
89 plot(o2_sat(start:middle), tempi(start:middle), '.');
90 hold on;
91 plot(o2_sat(middle:stop), tempi(middle:stop), 'r.');
92 title({'Pelagia Deployment 2'; 'Temp and O2 Saturation'; '\color{blue}Downcast = Blue \color{black}/ \color{red}Upcast = Red'});
93 xlabel('O2 Saturation (%)');
94 ylabel('Temp (Degrees C)');
95
96 clear;
97 load('GLIDER_DATA_LEVEL1\Pelagia_Deployment3_DO_L1.mat');
98
99 start = 134879;
100 middle = 135082;
101 stop = 135263;
102 subplot(3, 4, 9);
103 [ax, ~, h2] = plotyy(ptime_ebd_datenum, -depthi, ptime_ebd_datenum, o2_sat);
104 datetick(ax(1), 'x', 6, 'keeplimits', 'keepticks');
105 datetick(ax(2), 'x', 6, 'keeplimits', 'keepticks');
106 set(h2, 'Color', 'r');
107 set(ax(2), 'YColor', 'r');
108 title({'Pelagia Deployment 3'; 'Depth and O2 Saturation'; '\color[rgb]{0 .5 0}Sample Time as Green Bar'});
109 xlabel('UTC Time');
110 ylabel(ax(1), 'Depth (m)');
111 ylabel(ax(2), 'O2 Saturation (%)');
112 axes(ax(2));
113 yl = get(gca, 'YLim');
114 line([ptime_ebd_datenum(middle) ptime_ebd_datenum(middle)], [yl(1) yl(2)], 'LineStyle', '-', 'Color', 'green', 'LineWidth', 2);
115
116 subplot(3, 4, 10);
117 [ax, ~, h2] = plotyy(ptime_ebd_datenum(start:stop), -depthi(start:stop), ptime_ebd_datenum(start:stop), o2_sat(start:stop));
118 datetick(ax(1), 'x', 15, 'keeplimits', 'keepticks');
119 datetick(ax(2), 'x', 15, 'keeplimits', 'keepticks');
120 set(h2, 'Color', 'r');
121 set(ax(2), 'YColor', 'r');
122 title({'Pelagia Deployment 3'; 'Depth and O2 Saturation'; ['on ', datestr(ptime_ebd_datenum(start), 1)]});
123 xlabel('UTC Time');
124 ylabel(ax(1), 'Depth (m)');
125 ylabel(ax(2), 'O2 Saturation (%)');
126
127 subplot(3, 4, 11);
128 plot(o2_sat(start:middle), -depthi(start:middle), '.');
129 hold on;
130 plot(o2_sat(middle:stop), -depthi(middle:stop), 'r.');
131 title({'Pelagia Deployment 3'; 'Depth and O2 Saturation'; '\color{blue}Downcast = Blue \color{black}/ \color{red}Upcast = Red'});
132 xlabel('O2 Saturation (%)');
133 ylabel('Depth (m)');
134
135 subplot(3, 4, 12);
136 plot(o2_sat(start:middle), tempi(start:middle), '.');
137 hold on;
138 plot(o2_sat(middle:stop), tempi(middle:stop), 'r.');
139 title({'Pelagia Deployment 3'; 'Temp and O2 Saturation'; '\color{blue}Downcast = Blue \color{black}/ \color{red}Upcast = Red'});
140 xlabel('O2 Saturation (%)');
141 ylabel('Temp (Degrees C)');
Note: See TracBrowser for help on using the browser.