%Programa que me paso Julliette para hacer espectros con multitaper 
addpath /Francia/Velocidades
load S1v
sst = S1v;%sst_nino3;
nfen=4;
nfft=floor(length(sst)/nfen);
Fs=1;
[spect,Pxxc,freq]=pmtm(sst,nfen,nfft,Fs,'adapt',.95);
