figure(1);%hold on
subplot(2,2,1);
pcolor(c');shading flat;caxis([0 10]);
for i=1:1;
    t=num2str(i*1000+100000);
    fname=['h.',t(2:end),'.dat'];
    h=load(fname);
subplot(2,2,1);
pcolor(c');shading flat;caxis([0 10]);
    pause
end