c=load('h_init.dat');
figure(1);
subplot(5,2,1);
pcolor(c');shading flat;caxis([0 10]);axis equal
for i=1:9;
    t=num2str(i*1000+100000);
    fname=['h.',t(2:end),'.dat'];
    h=load(fname);
    subplot(5,2,i+1);
    pcolor(h);shading flat;caxis([0 10]);axis equal
end
