function h=plotLonModesYearFit5_3(DataFolder,fs,FS,FW,spacing,startyear,y2)
%y2=min, max colorbar y2<=0 min max calculated for each plot
location0='LonMod';
location=location0;
load(fullfile(DataFolder,location0,'LonModWeightedAtLeast6RadarsWorking_long_Fit5_3'))
i=startyear;
RadarsWorking=6;

h=get_figure(fs,FS,FW);
ind=find(LonMod.Time==datenum(i,1,1));
ind2=find(LonMod.Time==datenum(i,12,31));
Plot.lon=LonMod.lon;
Plot.Time=LonMod.Time(ind:ind2);
Plot.S0=LonMod.S0(ind:ind2,:);
Plot.S1=LonMod.S1(ind:ind2,:);
Plot.S2=LonMod.S2(ind:ind2,:);
Plot.RadarsWorking=LonMod.RadarsWorking(ind:ind2);
if (sum(~isnan(Plot.S0(:,1)))>5)
y1=1;
y2=-1;
axes1=subplot(1,4,1);

PlotLonModes(Plot,'S0',y1,y2,RadarsWorking)
title('a) S0')
set(axes1,'Position',[0.1 0.3 0.085 0.6]);

ylim([datenum(i,1,1) datenum(i,12,31)])
datetick('y','mmm','keeplimits');
ylim([datenum(i,1,1) datenum(i,12,31)])
ylabel(['month ',int2str(i)])
y1=2;
y2=-15;
axes2=subplot(1,4,2);

PlotLonModes(Plot,'S1',y1,y2,RadarsWorking)
title('b) S1')
ylim([datenum(i,1,1) datenum(i,12,31)])
datetick('y','mmm','keeplimits');
set(gca,'yTicklabel',{})
ylim([datenum(i,1,1) datenum(i,12,31)])
set(axes2, 'Position',[0.25 0.3 0.25 0.6]);

y1=3;
y2=-15;
axes3=subplot(1,4,3);

PlotLonModes(Plot,'S2',y1,y2,RadarsWorking)
title('c) S2')
ylim([datenum(i,1,1) datenum(i,12,31)])
datetick('y','mmm','keeplimits');
set(gca,'yTicklabel',{})
ylim([datenum(i,1,1) datenum(i,12,31)])



axes4=subplot(1,4,4);
plot(Plot.RadarsWorking,Plot.Time,'.')
title('d)')
ylim([datenum(i,1,1) datenum(i,12,31)])
datetick('y','mmm','keeplimits');
set(gca,'yTicklabel',{})
ylim([datenum(i,1,1) datenum(i,12,31)])
xlim([5.5,10.5])
set(axes4,'Position',[0.88 0.3 0.1 0.6]);
xlabel('#Radars')
set(axes3,'Position',[0.57 0.3 0.25 0.6]);
%  annotation(h,'textbox',...
%     [0.4 0.88 0.1 0.1],...
%     'String',{'Fit5_3'},...
%     'HorizontalAlignment','center',...
%      'FontSize',12,...
%     'FitBoxToText','off',...
%     'LineStyle','none');

figure(h)
t=['Hovmoeller',int2str(startyear),'_Fit5_3'];
    saveas(gca,fullfile(DataFolder,location,[t,'.fig']),'fig') 
    makePDF(gca,fs,t,fullfile(DataFolder,location))
end
end



