%Figuras Christine
load Pente_vent.mat
fet=date_windTP; feti=fet(1):fet(end);
xxx=feti(1):300:feti(end);
fechaA=datestr(xxx,12);

load vientochristinescasa
load topexscasa

load vientochristinescasaf
load topexscasaf

load catopex %ciclo anual
load cavientochristine

load savientochristine
load satopex

figure(1)
plot(feti(51:end-50),satopex(51:end-50),'g')
hold on
plot(feti(51:end-50),catopex(51:end-50))
hold on
plot(feti(51:end-50),topexscasaf,'r')
set(gca,'xlim',[min(feti)-3 max(feti)+3],'xtick',xxx,'xticklabel',fechaA)

figure(2)
plot(feti(51:end-50),savientochristine(51:end-50),'g')
hold on
plot(feti(51:end-50),cavientochristine(51:end-50))
hold on
plot(feti(51:end-50),vientochristinescasaf,'r')
set(gca,'xlim',[min(feti)-3 max(feti)+3],'xtick',xxx,'xticklabel',fechaA)

figure (3)

topexn=(topexscasaf-mean(topexscasaf))/std(topexscasaf);
vientochristinen=(vientochristinescasaf-mean(vientochristinescasaf))/std(vientochristinescasaf);
plot(feti(51:end-50),topexn)
hold on
plot(feti(51:end-50),vientochristinen,'r')
set(gca,'xlim',[min(feti)-3 max(feti)+3],'xtick',xxx,'xticklabel',fechaA)

