%Para hacer las figuras del poster de Brasil

clear
load y
load wave
%load uy136
load temp60
%load vs180
load period

%eje de tiempo
fecha=datenum(1999,11,11);
fet=1:375;
fet=fet+fecha;
fety1=fet(121:end);
x=fet(1):60:fet(end);
a=datestr(x,1);

subplot(2,1,2)
plot(fet,y,'r')
hold on
plot(fet,temp60,'k','LineWidth',2)
set(gca,'xlim',[min(fet)-3 max(fet)+3],'xtick',x,'xticklabel',[])
title('56-71 day period','fontsize',8)
subplot(3,1,3)
plot(fety1,y,'r')
hold on
plot(fety1,uy160,'k','LineWidth',2)
%set(gca,'xlim',[min(fet)-3 max(fet)+3],'xtick',x,'xticklabel',[])
title('56-71 day period','fontsize',8)

%subplot(4,1,4)
%plot(fet,y,'r')
%hold on
%plot(fet,vs180,'k','LineWidth',2)
%title('73-84 day period','fontsize',8)

set(gca,'xlim',[min(fet)-3 max(fet)+3],'xtick',x,'xticklabel',a,'fontsize',8)
X=121:375;
Y=1:241;
subplot(2,1,1)
contourf(abs(wave),[0:0.05:1])
caxis([0 1])
view(0,-90)
%colorbar
hold on
%plot([1 200],[75 75],'w')
ylim([0 210])
r=30:30:240;
%set(gca,'YTick',[6 46 70 86 99 109 118 126 133 139] ,'YTickLabel',[10:10:100])
set(gca,'YTick',r,'YTickLabel',floor(period(r)))
set(gca,'xlim',[-3 378],'xtick',x,'xticklabel',[])

ylabel('period (days)','fontsize',8)
title('Equator 80 m, "temperature" wavelet power spectrum')
