function  PlotLonModes(LonMod,S,y1,y2,RadarsWorking)
%plots a Hovmöller-Diagram

%lon=-180:1:180;
d=max([abs(round(min(min(LonMod.(S))))),abs(round(max(max(LonMod.(S)))))]);
if d>5
contourf(LonMod.lon(1:y1:end),LonMod.Time,LonMod.(S)(:,1:y1:end),'LevelStep',1,'LineStyle','none')
else
    contourf(LonMod.lon(1:y1:end),LonMod.Time,LonMod.(S)(:,1:y1:end),'LevelStep',0.1,'LineStyle','none')
end

title([S,', at least ',int2str(RadarsWorking),' radars working']);
if length(RadarsWorking)==1
if RadarsWorking>10
title([S,' WACCM, pressurelevel ',LonMod.pressure(RadarsWorking)]);
end
end
% if y1==1
% ylabel('date [dd/mm/yy]')
% set(gca,'yTick',LonMod.Time(1:20:end))
% datetick('y','dd/mm/yy','keepticks')
% end
xlabel('longitude')
set(gca,'xTick',LonMod.lon(1):180:LonMod.lon(end));


%ylim([LonMod.Time(y1),LonMod.Time(y2)])
if y2>0
    c=colorbar('horiz');
    d=y2;
    caxis([-d,d]); 
    set(c, 'Position',[0.2 0.1 0.7 0.05]);
elseif y2==0 %if y2==0
    c=colorbar('horiz');
    caxis([-d,d]); 
else %y2<0
    c=colorbar('horiz');
    caxis([y2,-y2])
end
if y1==2
xlabel(c,'meridional wind [m/s]')
if length(RadarsWorking)>1
   xlabel(c,['meridional wind [m/s] ',RadarsWorking],'Interpreter','none') 
end

colormap(redblue)
end

