
disps=dhat_series;

j=0;
for k=1:numsites

j=j+1;    

if mod(k,3)==1
    figure
    j=1;
end
    

subplot(3,3,j)
errorbar(InvertSeries{k,1}(:,1)-teq,InvertSeries{k,1}(:,2),InvertSeries{k,1}(:,3),'.'); hold on
plot(InvertSeries{k,1}(:,1)-teq,disps{3*k-2},'r');
axis tight

title(names(k,:))

j=j+1;
subplot(3,3,j)
errorbar(InvertSeries{k,2}(:,1)-teq,InvertSeries{k,2}(:,2),InvertSeries{k,2}(:,3),'.'); hold on
plot(InvertSeries{k,2}(:,1)-teq,disps{3*k-1},'r');
axis tight

j=j+1;
subplot(3,3,j)
errorbar(InvertSeries{k,3}(:,1)-teq,InvertSeries{k,3}(:,2),InvertSeries{k,3}(:,3),'.'); hold on
plot(InvertSeries{k,3}(:,1)-teq,disps{3*k},'r');
axis tight

end

