z,x] = odn2grid(o,d,n);
z = z*1e-3; x = x*1e-3;

%%
scnsize = get(0,'ScreenSize');

% The true and initial model
%

figure('Position',scnsize./[1 1 2 2]);imagesc(x,z,v,[1.5 4.5]);colorbar;set(gca,'plotboxaspectratio',[3 1 1]);
xlabel('x [km]');ylabel('z [km]'); title('true model');

%%
%

figure('Position',scnsize./[1 1 2 2]);imagesc(x,z,v0,[1.5 4.5]);colorbar;set(gca,'plotboxaspectratio',[3 1 1]);
xlabel('x [km]');ylabel('z [km]'); title('initial model');

%%
% reconstruction with batching after 17th. freq