clear all
ngx = 125;
ngy = 125;
load solution.out
k = 1;
for iy=1:ngx
for ix=1:ngy
x(ix,iy) = solution(k,1);
y(ix,iy) = solution(k,2);
w(ix,iy) = solution(k,3);
k = k+1;
end
end


figure(2)
imagesc(solution(:,1),solution(:,2),w')
axis square
caxis([-0.1 0.1])
colormap gray
