nces
%         The University of British Columbia
%         
% Date: 02, 2012
%

%% Camenbert example 
% 
% This is an demonstration of modified gausss-newton algorithm based on Camembert
% model [3], please visit the script <MGNFWI_camenbert.html> for detail 
%
% true model and initial model
z = 0:10:1000;
x = 0:10:1000;
[o,d,n] = grid2odn(z,x);
[zz,xx] = ndgrid(z,x);
m0 = 2500 + 0*xx;

% circular perturbation with radius 250 m and strenth 10%
dv = 0*xx; dv((xx-500).^2 + (zz-500).^2 <= 250^2) = .1*2500;
m 