function Havoc
% HAVOC Horizontal Adjustment by Variation Of Coordinates

%
% Set up global structures to store information
global STATION OBS JOB X Z VFACT Qxx               

STATION = struct('serial',0,'east',0.0,'north',0.0,'code',0);
OBS     = struct('at',0,'to',0,'left',0,'right',0,'val',0.0,'sigma',0.0,'code',0,'status',0);
JOB     = struct('infile',' ','outfile',' ','fixed',0,'floating',0,'dirns',0,'dists',0,...
                 'cbrgs',0,'cdists',0,'cangles',0,'orients',0,'n',0,'u',0,'c',0,'iter',0);

%
% read in the data
readdata;
%
% determine status of each observation
status;
%
% solve the system
solve;
%
% write out the results
results;
%
% plot the network
HavocPlot;


