%figura 3 articulo
clear
load waveun2
load dilatun2
load ridgeun2
load yun2
load period
load wavevn2
load dilatvn2
load ridgevn2
load yvn2


subplot(2,2,3);
plot(yun2,'r')
hold on
plot(ridgeun2,'k','LineWidth',2)
ylabel('\delta u (cm/s)')

subplot(2,2,1);
contourf(abs(waveun2),[0:0.05:1])
caxis([0 1])
view(0,-90)
colorbar('horiz')
hold on
%plot([1 200],[75 75],'w')
ylim([0 200])
r=30:30:200;
%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)))
xlabel('days')
ylabel('period(days)')
plot(dilatun2,'w','LineWidth',2)

subplot(2,2,4);
plot(yvn2,'r')
hold on
plot(ridgevn2,'k','LineWidth',2)

subplot(2,2,2);
contourf(abs(wavevn2),[0:0.05:1])
caxis([0 1])
view(0,-90)
colorbar('horiz')
%set(gca,'position',[0.6 0.1 0.1 0.8])
hold on
%plot([1 200],[75 75],'w')
ylim([0 200])
r=30:30:200;
%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)))%, 'position',[0.6 0.1 0.3 0.4])
xlabel('days')

plot(dilatvn2,'w','LineWidth',2)

