
[KPM, Datumk]=VolPeriM(DAXT, 1, 300, DATUMT, DATUMJ, 1)
clear F1
for i=1:300
    isokorrel(KPM(:,:,1:250), i/1000, Datumk, 0, jet, .05, 0, 0, 10)
    daspect([50 50 1])
    camlight
    view(100+i*2,(i/100+1)*45)
    F1(i)=getframe(gcf);
    i
end
