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 |
|
---|