%ENG_ELEV uses ordinary kriging to approximate England elevations.


%NOTE: Requires that 'eng_elev_dat.txt' be loaded in the workspace.

clear('X0','L0','X1','L1');

clear('x','y','x1','y1','x2','y2','z','L','X','XX') ;

z = eng_elev_dat(:,3); %elevation

%Transform Sample Locations


L = eng_elev_dat(:,1:2);



%SAMPLE DATA

L0 = L ;


%Now form grid locations

clear('x','y','X','U','x1','y1','x2','y2','z2','s1','u1','U1',...
      's2','u2','U2','S'); 

s1 = [160:20:580]';
u1 = length(s1);
U1 = ones(u1,1);

s2 = [120:20:580]';
u2 = length(s2);
U2 = ones(u2,1);

S(:,1) = kron(s1,U2);
S(:,2) = kron(U1,s2);



%GRID DATA

L1 = S ;