load ubar;
time=tm;
load discharge;
Q=interp1(tm,Q,time);
load eta_BP;
amp1=interp1(tm,amp,time);
hlow1=interp1(tm,hlow,time);
load eta_bat;
amp2=interp1(tm,amp,time);
hlow2=interp1(tm,hlow,time);
load winds_jfk_low;
ue=interp1(tm,ue,time);
un=interp1(tm,un,time);
X(1,:)=v;
X(2,:)=Q;
X(3,:)=amp1;
X(4,:)=hlow1;
X(5,:)=amp2;
X(6,:)=hlow2;
X(7,:)=ue;
X(8,:)=un;
X=X';

corrcoef(X)