function Data=DoReadWindData(DataFolder,station,stationlist)
smoothing=30;

%% readData
%readNewSuperDARNdata(station,stationlist,DataFolder);
 for i=1:length(station)
 disp(station{i})

readData(DataFolder,station{i});
figure
load(fullfile(DataFolder,station{i},[station{i},'Clean']))
plot(Data.DateNum)
title(station{i})
clear Data
 
 end




%% plot data coverage Lat
C=[0 0 0;0 0.5 0; 0 1 0; 1 0 0; 0 0 1; 1 0 1; 0 1 1;1 1 0; 1 0.5 0;0.5 0 0.5]; 
figure2=figure('Name','DatacoverageLat','units','centimeters','Position',[0 0 20 10],...                    % setter str på figuren du vil se på skjermen. VIKIT Gat units kommer FØR Position
    'PaperUnits','centimeters','PaperPosition', [0 0 20 10],...% setter str på figuren du vil lagre (uavhengig av format) VIKITG at units kommer FØR Position
    'DefaultAxesFontSize',12,'DefaultAxesFontWeight','b');               %setter str på font i figure 

for i=1:length(station)
  load(fullfile(DataFolder,station{i},[station{i},'Clean']))
  
  plot(Data.DateNum,Data.Lat,'line','.','Color',C(i,:))
  hold on
end
ylabel('latitude')
xlabel('date [year/month]')
datetick('x','yy/mm');  
legend(station{1},station{2},station{3},station{4},station{5},station{6},station{7},station{8},station{9},station{10},'location','EastOutside')
plot([datenum(2004,8,8),datenum(2004,8,8)],[50,65],'linewidth',2)

%% plot data coverage Lon
C=[0 0 0;0 0.5 0; 0 1 0; 1 0 0; 0 0 1; 1 0 1; 0 1 1;1 1 0; 1 0.5 0;0.5 0 0.5]; 
figure3=figure('Name','DatacoverageLon','units','centimeters','Position',[0 0 20 10],...                    % setter str på figuren du vil se på skjermen. VIKIT Gat units kommer FØR Position
    'PaperUnits','centimeters','PaperPosition', [0 0 20 10],...% setter str på figuren du vil lagre (uavhengig av format) VIKITG at units kommer FØR Position
    'DefaultAxesFontSize',12,'DefaultAxesFontWeight','b');               %setter str på font i figure 

for i=1:length(station)

  load(fullfile(DataFolder,station{i},[station{i},'Clean']))
  
  plot(Data.DateNum,Data.Lon,'line','.','Color',C(i,:))
  hold on
end
ylabel('longitude')
xlabel('date [year/month]')
datetick('x','yy/mm');  
ylim([-180,180])
legend(station{1},station{2},station{3},station{4},station{5},station{6},station{7},station{8},station{9},station{10},'location','EastOutside')
plot([datenum(2004,8,8),datenum(2004,8,8)],[-180,180],'linewidth',2)
%%
end