function SAFmap = orgmap (x, y)
%This code is about gridding.
%If we we do not have evenly spaced observations, than
%We need to grid our unevenly spaced data by determining 
%some set of evenly spaced estimates that approximate the observations.
%EXAMPLE;  Here is the map of points:
    %figure(1)
    %clf
    %plot(e, n, 'k+')
    %xlabel('easting')    
    %ylabel('northing')   
    %title('Observation points')
    %axis equal
    %axis([0 6.5 0 6.5])
    %print -depsc topopts.eps
    
    load ('e.txt')
    load ('n.txt')
    load ('h.txt')
    
    figure (1)
    clf
    plot (e, n, '+k')
    
    %do i need to use griddata here?
    help griddata
    [xi,yi,zi] = griddata(e, n, h, e, n');
    
    xlabel ('easting')
    ylabel ('northing')
    title ('Orginial Data Points')
    axis equal
    
    print('SAFdatapts.jpg','-djpeg');
    
   
end