function DoFitWindData(DataFolder,stations)
% %station=[{'Saskatoon'},{'Pykkvibaer'}];
% %% Fit tides and 2 day wave to data
 for i=1:length(stations)
     disp(['Fit ',stations{i}])
 FitData3(DataFolder,stations{i},0) % incl. 2 day component and 8 h component
 end
% %% Monthly Mean 
% %yearStart=[2001,2000,2000,1993,2007,1993,1995,1994,1996,1995];
% %yearEnd=[2009,2009,2009,2008,2008,2009,2008,2008,2008,2009];
 for i=1:length(stations)    
 ProcessedMonthlyWeightedMeanLST( DataFolder, stations, stations{i},3,1);%%with 2 day and 8 hour
 end
% 
%% plot data coverage mean dayly wind Fit3
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]; 
figure4=figure('Name','DatacoverageLon2','units','centimeters','Position',[0 0 50 10],...                    % setter str på figuren du vil se på skjermen. VIKIT Gat units kommer FØR Position
    'PaperUnits','centimeters','PaperPosition', [0 0 50 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(stations)
 
  %load([station{i},'_Fit2'])
  load(fullfile(DataFolder,stations{i},[stations{i},'_Fit3']))
  Data.Fit.Lon=ones(size(Data.Fit.DaysWithFit))*Data.Lon(1);
  plot(Data.Fit.DaysWithFit,Data.Fit.Lon,'.','Color',C(i,:))
  hold on
  Lon=ones(size(Data.Fit.DaysWithoutFit))*Data.Lon(1);
 plot(Data.Fit.DaysWithoutFit,Lon,'.','Color',[1 1 1])
  clear Lon
  %save([station{i},'_Fit2'],'Data')
  save(fullfile(DataFolder,stations{i},[stations{i},'_Fit3']),'Data')
hold on
end
ylabel('longitude')
xlabel('date [year/month]')
datetick('x','yy/mm');  
legend(stations{1},'',stations{2},'',stations{3},'',stations{4},'',stations{5},'',stations{6},'',stations{7},'',stations{8},'',stations{9},'',stations{10},'','location','EastOutside')
%plot([datenum(2004,8,8),datenum(2004,8,8)],[-180,180],'linewidth',2)
ylim([-181,180])

saveas(figure4,fullfile(DataFolder,['LatitudeCoverage3','.fig']))
makePDF(figure4,[50 10],'LatitudeCoverage3',DataFolder)
