%VENICE_FUNC uses the Venice coordinate data to construct spatial 
%models of 'draw-down'functions

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

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


V = venice ;
z = V(:,4); %1973 data

%Transform Sample Locations

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

L = V(:,1:2);

%Industry-centered coordinates

x1 = .001*(.873*(x - 418) - .488*(y - 458));
y1 = .001*(.488*(x - 418) + .873*(y - 458));

%Approximation of Elevation Effects

X(:,1) = abs(x1) ;

%Industry draw-down function

X(:,2) = exp((1.5)*(x1.^2) + (y1.^2));

%Venice Island draw-down function

x2 = .001*(x - 545) ;
y2 = .001*(y - 387) ;

X(:,3) = exp(-(x2.^2 + y2.^2)) ;

%SAMPLE DATA

X0 = X ;
L0 = L ;