function meshmap (data, xblocks, yblocks, interpmethod)
 %meshamp
 %
 %Use delaunay function to compute the triangles, then use trimesh to plot
 %them
    %
    %delaunay function; TRI = delaunay(x,y,options)
    %Trimesh function; trimesh(Tri,X,Y,Z) 
    %
    %Example:
    %Create vertex vectors and a face matrix, then create a triangular mesh
    %plot.
    %x = rand(1,50);
    %y = rand(1,50);
    %z = peaks(6*x-3,6*x-3);
    %tri = delaunay(x,y);
    %trimesh(tri,x,y,z)
    
    %triangles for the linear interpolation
    %
 %Wendy Bohon, March 9, 2009
    
interpmethod = 'linear'

%Clears the previous figure and makes the new one
figure(5)
clf

%delaunay function computes the triangles
tri = delaunay(e,n);

%trimesh function plots the triangles computed by delaunay
trimesh(tri,e,n,h)
colormap([0 0 0])
view(2)
hold on

%show data:
plot(e,n,'ko')

%labels the axis and names the image
xlabel('easting')    
ylabel('northing')   
title('SAF Mesh Plot')
axis equal
print('meshmap_SAF.jpg','-djpeg');

end