function KNMI_plot(xx,yy1,yy2)
%KNMI_PLOT  plot result from POTWIND_330_2001.TXT (this file is modified by hand)
%
%
%See also: SILTPROFILER_PLOT_FITTED_OBS; SILTPROFILER_PLOT_DEBUG

%% Copyright notice
%   --------------------------------------------------------------------
%   Copyright (C) 2010 Delft University of Technology
%       Gerben J. de Boer
%
%       g.j.deboer@tudelft.nl	
%
%       Fluid Mechanics Section
%       Faculty of Civil Engineering and Geosciences
%       PO Box 5048
%       2600 GA Delft
%       The Netherlands
%
%   This library is free software: you can redistribute it and/or
%   modify it under the terms of the GNU Lesser General Public
%   License as published by the Free Software Foundation, either
%   version 2.1 of the License, or (at your option) any later version.
%
%   This library is distributed in the hope that it will be useful,
%   but WITHOUT ANY WARRANTY; without even the implied warranty of
%   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
%   Lesser General Public License for more details.
%
%   You should have received a copy of the GNU Lesser General Public
%   License along with this library. If not, see <http://www.gnu.org/licenses/>.
%   --------------------------------------------------------------------
%subplot(2,1,1);
x1   = xx(75211:75227);
%x2   = x(75000:75048);
y1  = yy1(75211:75227);
y2  = yy2(75211:75227);
%
%hl1 = line(x1,y1,'Color','r');
%ax1 = gca;
%set(ax1,'XColor','k','YColor','r')
%AXIS([x(529) x(577) 0 15])
%datetick('x')
%ax2 = axes('Position',get(ax1,'Position'),...
%           'XAxisLocation','top',...
%           'YAxisLocation','right',...
%           'Color','none',...
%           'XColor','k','YColor','b');
%
%hl2 = line(x2,y2,'Color','b','Parent',ax2);
%AXIS([x(529) x(577) 0 360])
%datetick('x')
figure
Y1=mod(-270-y2,360);
wind_rose(Y1,y1','ci',[5 10 15 30 40],'legtype',1,'lablegend','Windspeed [m/s]','labtitle',...
    'Wind direction and magnitude for 31st of July 2009')


y2=mod(-90-y2,360);
Y3 = deg2rad(y2);
[U,V]=pol2cart(Y3,y1);

x2   = xx(75211:75217);
%x2   = x(75000:75048);
y3   = yy1(75211:75217);
y4   = yy2(75211:75217);
%

figure
Y3=mod(-270-y4,360);
wind_rose(Y3,y3','ci',[5 10 30 40],'legtype',1,'lablegend','Windspeed [m/s]','labtitle',...
    'Wind direction and magnitude for 30th and 31st of July 2009')



%figure
%compass(U,V)

%figure
%feather(U,V)


