if(mod(t,100)==0)
figure(1), cla, box
set(gcf, 'Position',  [100, 100, 1000, 800])

    imagesc([min(min(X)) max(max(X))],[min(min(Y)) max(max(Y))], T(iny,inx)), axis equal, ...
        axis tight,
        ylabel('y [m]')
    xlabel('x [m]'), colormap jet, drawnow, colorbar, pause(0.1)

%spreadfigures

end