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

root/adcp/trunk/adcp/adcp_matlab/Winfo_plot.m

Revision 168 (checked in by cbc, 16 years ago)

Adding diwasp customizations.

Line 
1 function [fig_sigh,fig_pp,fig_dp,fig_dd]=Winfo_plot(Winfo);
2
3 %This function takes the wave info (sigH, peak period, direction of peak and dominant direction
4 %for a series of bursts and plots it as a function of time
5
6
7 %plot sig wave height
8 scrsz = get(0,'ScreenSize');
9 fig_sigh=figure('Position',[scrsz]);
10 h1=plot(Winfo.time,Winfo.hsig(1,:),'b.-');
11 hold on
12 h2=plot(Winfo.time,Winfo.hsig(2,:),'r.-');
13 h3=plot(Winfo.time,Winfo.hsig(3,:),'g.-');
14 h4=plot(Winfo.time,Winfo.hsig(4,:),'c.-');
15 h5=plot(Winfo.time,Winfo.hsig(5,:),'m.-');
16 h6=plot(Winfo.time,Winfo.hsig(7,:),'k.-');
17
18 datetick('x',15,'keepticks');
19 legend([h1, h2, h3, h4, h5, h6],'EMEP uvw','IMLM uvw','EMEP range','IMLM range','EMEP radial','wavesmon','location','best');
20 title(['Significant Wave Height starting on ',datestr(Winfo.time(:,1),0)]);
21 xlabel('time');
22 ylabel('Wave Height (m)');
23
24 %plot peak period
25 fig_pp=figure('Position',[scrsz]);
26 h1=plot(Winfo.time,Winfo.peakP(1,:),'b.-');
27 hold on
28 h2=plot(Winfo.time,Winfo.peakP(2,:),'r.-');
29 h3=plot(Winfo.time,Winfo.peakP(3,:),'g.-');
30 h4=plot(Winfo.time,Winfo.peakP(4,:),'c.-');
31 h5=plot(Winfo.time,Winfo.peakP(5,:),'m.-');
32 h6=plot(Winfo.time,Winfo.peakP(7,:),'k.-');
33
34 datetick('x',15,'keepticks');
35 legend([h1, h2, h3, h4, h5, h6],'EMEP uvw','IMLM uvw','EMEP range','IMLM range','EMEP radial','wavesmon','location','best');
36 title(['Peak Period starting on ',datestr(Winfo.time(:,1),0)]);
37 xlabel('time');
38 ylabel('Period (seconds)');
39
40 %plot dir of peak period
41 fig_dp=figure('Position',[scrsz]);
42 h1=plot(Winfo.time,Winfo.dirP(1,:),'b.-');
43 hold on
44 h2=plot(Winfo.time,Winfo.dirP(2,:),'r.-');
45 h3=plot(Winfo.time,Winfo.dirP(3,:),'g.-');
46 h4=plot(Winfo.time,Winfo.dirP(4,:),'c.-');
47 h5=plot(Winfo.time,Winfo.dirP(5,:),'m.-');
48 h6=plot(Winfo.time,Winfo.dirP(7,:),'k.-');
49
50 datetick('x',15,'keepticks');
51 legend([h1, h2, h3, h4, h5, h6],'EMEP uvw','IMLM uvw','EMEP range','IMLM range','EMEP radial','wavesmon','location','best');
52 title(['Direction of Peak Period starting on ',datestr(Winfo.time(:,1),0)]);
53 xlabel('time');
54 ylabel('Direction (Degrees True)');
55
56 %plot the dominant direction
57
58 fig_dd=figure('Position',[scrsz]);
59 h1=plot(Winfo.time,Winfo.Ddir(1,:),'b.-');
60 hold on
61 h2=plot(Winfo.time,Winfo.Ddir(2,:),'r.-');
62 h3=plot(Winfo.time,Winfo.Ddir(3,:),'g.-');
63 h4=plot(Winfo.time,Winfo.Ddir(4,:),'c.-');
64 h5=plot(Winfo.time,Winfo.Ddir(5,:),'m.-');
65 h6=plot(Winfo.time,Winfo.Ddir(7,:),'k.-');
66
67 datetick('x',15,'keepticks');
68 legend([h1, h2, h3, h4, h5, h6],'EMEP uvw','IMLM uvw','EMEP range','IMLM range','EMEP radial','wavesmon','location','best');
69 title(['Dominant Direction starting on ',datestr(Winfo.time(:,1),0)] );
70 xlabel('time');
71 ylabel('Direction (Degrees True)');
72
Note: See TracBrowser for help on using the browser.