1 %====================================================================== |
1 %====================================================================== |
2 % G E T E R R . M |
2 % G E T E R R . M |
3 % doc: Wed Jun 30 23:24:51 2004 |
3 % doc: Wed Jun 30 23:24:51 2004 |
4 % dlm: Wed Jul 6 20:26:14 2011 |
4 % dlm: Wed Jan 28 09:46:12 2015 |
5 % (c) 2004 ladcp@ |
5 % (c) 2004 ladcp@ |
6 % uE-Info: 19 35 NIL 0 0 72 2 2 8 NIL ofnI |
6 % uE-Info: 21 42 NIL 0 0 72 2 2 8 NIL ofnI |
7 %====================================================================== |
7 %====================================================================== |
8 |
8 |
9 % MODIFICATIONS BY ANT: |
9 % MODIFICATIONS BY ANT: |
10 % Jun 30, 2004: - BUG: bin numbering was wrong for asymmetric up/down |
10 % Jun 30, 2004: - BUG: bin numbering was wrong for asymmetric up/down |
11 % bin setup |
11 % bin setup |
15 % Oct 7, 2008: - extensively modified procfig 3 for version IX_6 |
15 % Oct 7, 2008: - extensively modified procfig 3 for version IX_6 |
16 % Jun 29, 2011: - removed skp_figure_3 |
16 % Jun 29, 2011: - removed skp_figure_3 |
17 % - added ps.fig3_colormap, ps.fig3_err_y_axis, ps.fig3_avgerr |
17 % - added ps.fig3_colormap, ps.fig3_err_y_axis, ps.fig3_avgerr |
18 % Jun 30, 2011: - fixed fig.3 middle column plot title for median plot |
18 % Jun 30, 2011: - fixed fig.3 middle column plot title for median plot |
19 % Jul 6, 2001: - fixed plot title |
19 % Jul 6, 2001: - fixed plot title |
|
20 % Jan 25, 2015: - separated uc/dc in bin-averaged residual plots |
|
21 % Jan 28, 2015: - BUG: figure legend typo |
20 |
22 |
21 function l=geterr(ps,dr,d,iplot) |
23 function l=geterr(ps,dr,d,iplot) |
22 % function l=geterr(dr,d,iplot) |
24 % function l=geterr(dr,d,iplot) |
23 % returns predicitons of U_ocean and |
25 % returns predicitons of U_ocean and |
24 % U_ctd on the raw data grid |
26 % U_ctd on the raw data grid |
230 xlabel('Super Ensemble #'); |
235 xlabel('Super Ensemble #'); |
231 title(sprintf('U-err std: %.03f',meannan(stdnan(l.ru_err')))) |
236 title(sprintf('U-err std: %.03f',meannan(stdnan(l.ru_err')))) |
232 |
237 |
233 subplot(232) |
238 subplot(232) |
234 if ps.fig3_avgerr == 2 |
239 if ps.fig3_avgerr == 2 |
235 plot(medianan(l.ru_err')',-ib) |
240 plot(medianan(l.ru_err(:,1:btmi)')',-ib,'r',medianan(l.ru_err(:,btmi:end)')',-ib,'b'); |
236 title('median(U-err)') |
241 title('median(U-err) [r/b: down-/up-cast]') |
237 else |
242 else |
238 plot(meannan(l.ru_err')',-ib) |
243 plot(meanan(l.ru_err(:,1:btmi)')',-ib,'r',meanan(l.ru_err(:,btmi:end)')',-ib,'b'); |
239 title('mean(U-err)') |
244 title('mean(U-err) [r/b: down-/up-cast]') |
240 end |
245 end |
241 set(gca,'XLim',[-0.05 0.05]); |
246 set(gca,'XLim',[-0.05 0.05]); |
242 set(gca,'Ylim',[-ib(end) -ib(1)]); |
247 set(gca,'Ylim',[-ib(end) -ib(1)]); |
243 set(gca,'Xtick',[-0.04:0.02:0.04]); |
248 set(gca,'Xtick',[-0.04:0.02:0.04]); |
244 grid |
249 grid |
296 ylabel('Bin #'); |
301 ylabel('Bin #'); |
297 title(sprintf('V-err std: %.03f',meannan(stdnan(l.rv_err')))) |
302 title(sprintf('V-err std: %.03f',meannan(stdnan(l.rv_err')))) |
298 |
303 |
299 subplot(235) |
304 subplot(235) |
300 if ps.fig3_avgerr == 2 |
305 if ps.fig3_avgerr == 2 |
301 plot(medianan(l.rv_err')',-ib) |
306 plot(medianan(l.rv_err(:,1:btmi)')',-ib,'r',medianan(l.rv_err(:,btmi:end)')',-ib,'b'); |
302 title('median(V-err)') |
307 title('median(V-err) [r/b: down-/up-cast]') |
303 else |
308 else |
304 plot(meannan(l.rv_err')',-ib) |
309 plot(meanan(l.rv_err(:,1:btmi)')',-ib,'r',meanan(l.rv_err(:,btmi:end)')',-ib,'b'); |
305 title('mean(V-err)') |
310 title('mean(V-err) [r/b: down-/up-cast]') |
306 end |
311 end |
307 set(gca,'XLim',[-0.05 0.05]); |
312 set(gca,'XLim',[-0.05 0.05]); |
308 set(gca,'Ylim',[-ib(end) -ib(1)]); |
313 set(gca,'Ylim',[-ib(end) -ib(1)]); |
309 set(gca,'Xtick',[-0.04:0.02:0.04]); |
314 set(gca,'Xtick',[-0.04:0.02:0.04]); |
310 grid |
315 grid |