4 % |
4 % |
5 |
5 |
6 %====================================================================== |
6 %====================================================================== |
7 % P L O T R A W . M |
7 % P L O T R A W . M |
8 % doc: Fri Jan 5 15:38:43 2007 |
8 % doc: Fri Jan 5 15:38:43 2007 |
9 % dlm: Fri Jan 5 15:46:39 2007 |
9 % dlm: Mon Jun 3 14:09:21 2013 |
10 % (c) 2007 M. Visbeck with contribs from A. Thurnherr |
10 % (c) 2007 M. Visbeck with contribs from A. Thurnherr |
11 % uE-Info: 15 37 NIL 0 0 72 0 2 4 NIL ofnI |
11 % uE-Info: 79 0 NIL 0 0 72 0 2 4 NIL ofnI |
12 %====================================================================== |
12 %====================================================================== |
13 |
13 |
14 % MODIFICATIONS BY ANT: |
14 % MODIFICATIONS BY ANT: |
15 % Jan 5, 2007: - fixed checkbeam() as suggested by B. Huber |
15 % Jan 5, 2007: - fixed checkbeam() as suggested by B. Huber |
|
16 % Jun 3, 2013: - BUG: top panel of Fig. 2 was wrong for dual-headed |
|
17 % LADCPs with different UL/DL bin sizes |
16 |
18 |
17 pmax=200; |
19 pmax=200; |
18 orient tall |
20 orient tall |
19 |
21 |
20 if length(d.time_jul)>pmax |
22 if length(d.time_jul)>pmax |
54 ie=sum(~isnan(d.re(iz,ii))); |
56 ie=sum(~isnan(d.re(iz,ii))); |
55 i3bd=find(iw>(2*ie)); |
57 i3bd=find(iw>(2*ie)); |
56 n3bd=length(i3bd)/length(ii)*100; |
58 n3bd=length(i3bd)/length(ii)*100; |
57 end |
59 end |
58 |
60 |
59 % contour results |
61 % plot vertical velocities |
60 col=jet(128); |
62 col=jet(128); |
61 col=([[1 1 1]; col]); |
63 col=([[1 1 1]; col]); |
62 colormap(col) |
64 colormap(col) |
63 pcolorn(ii,zz,rw), shading flat |
65 |
|
66 pcolorn(ii,zz,rw) % this is neccessary even for dual-headed data sets |
|
67 |
|
68 if (length(d.izu) > 0) |
|
69 pcolorn(ii,zz(1:length(d.izu)),rw(1:length(d.izu),:)) |
|
70 hold on, shading flat |
|
71 ax = axis; |
|
72 ax(1) = ii(1); |
|
73 ax(2) = ii(end); |
|
74 ax(4) = -zz(end); |
|
75 axis(ax); |
|
76 pcolorn(ii,zz(length(d.izu)+1:end),rw(length(d.izu)+1:end,:)) |
|
77 end |
|
78 |
64 hold on |
79 hold on |
65 ax=axis; |
80 plot([ii(1),ii(end)],[0 0],'-k') |
66 plot(ax(1:2),[0 0],'-k') |
81 colorbar('horiz') |
67 axis(ax), colorbar('horiz') |
82 |
68 % mark 3-beam solutions |
83 % mark 3-beam solutions |
69 if n3bu>10 |
84 if n3bu>10 |
70 l3b=NaN*ii; |
85 l3b=NaN*ii; |
71 l3b(i3bu)=ax(3); |
86 l3b(i3bu)=ax(3); |
72 plot(ii,l3b,'-r','linewidth',8) |
87 plot(ii,l3b,'-r','linewidth',8) |