%ENG_KRIGE computes relevant inputs for kriging the england 
%          temperature data.


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

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

clear('x','y','z','L','X') ;

z = eng_temp(:,3); %temperatures

%Transform Sample Locations

x = eng_temp(:,1);
y = eng_temp(:,2);

L = eng_temp(:,1:2);

X(:,1:2) = L;
X(:,3) = x.^2 ;

%SAMPLE DATA

X0 = X ;
L0 = L ;

%Now form grid locations

clear('x','y','X','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);

x = S(:,1) ;
y = S(:,2) ;

X(:,1:2) = S;
X(:,3) = x.^2 ;


%GRID DATA

X1 = X ;
L1 = S ;