% Read data
close all

load FIG'URE 5 (x, in meters).txt';
load FIG'URE 5 (y, in meters).txt';
x = FIGURE_5__x__in_meters_;
y = FIGURE_5__y__in_meters_;

Nx = length(x);

% remove first point
% x = x(2:end);
% y = y(2:end);

% Interpolate once. Gives a little more than 12p/wavelength in r-direction
% I_C2F = Interp_1to2_4th_order(Nx);
% y = I_C2F*y;
% x = 0:0.5:1999;

% Remove first two points
% x = x(3:end);
% y = y(3:end);
% Nx = length(x);


[x,y] = manipulate_geometry(x,y);
Nx = length(x);

Nx = length(x);
g = y;
z1 = 200;
Nz = 301; %Nz = 301 => 10p/wavelength i z-dir.
Nr = length(x);
x0 = x(1); x1 = x(end);
[R,Z,dr,dz,Rxi1,Rxi2,Zxi1,Zxi2,CHI] = grid_Kurt(Nx,Nz,x0,x1,g,z1);


