function geoplot( minlon, maxlon, minlat, maxlat, coastfile, bathyfile, dx, dy)
%
%    The function geoplot(minlon, maxlon, minlat, maxlat, 
%                         coastfile,bathyfile)
%    makes a plot of black (filled-in) land and isobaths between the
%    coordinates minlon < x < maxlon, minlat < y < maxlat. The variables
%    dx and dy define how much the bathymetry should be shifted to
%    fit the coastline, as the files are independent.  The coast
%    line data should be read from the file 'coastfile', while the
%    bathymetric data (isobaths) from the file 'bathyfile', which 
%    if it doesn't exist, can be created by running my MATLAB sub-
%    routine 'geomake'.  This subroutine is MUCH faster than geosurf,
%    or geocont, and (on the contrary to the forementioned) can be
%    superimposed by float tracks.
%
%    See also geocont, geosurf, geomake
%
%---------------------------------------------------------------- 
%    Written by Vassilis Zervakis, in January 1996, while at
%               AINCO - Interocean, Madrid, Spain.
%----------------------------------------------------------------
%
%----read and plot the coastline file ----------
%    figure
     coastfile;
     eval('bwmap2(minlon, maxlon, minlat, maxlat, coastfile)');
     hold on;
%----read the bathymetry file ----------
filenm = num2str(bathyfile);
n = findstr('.',filenm);
dirs = findstr('\',filenm);
nn = length(dirs);
m = dirs(nn);
filebody = filenm(m+1:n-1);
eval(['load ',filenm]);
data = eval(filebody);
depth = data(1,1);
points= data(1,2);
x = data(2:points+1,1) + dx;
y = data(2:points+1,2) + dy;
plot(x,y);
hold on;
clear x;
clear y;
points_read = points + 1;
while points_read < length(data(:,1)),
 m = points_read+1;
 depth = data(m,1);
 points = data(m,2);
 points_read = points_read + points + 1;
 x = data(m+1:points_read,1) + dx;
 y = data(m+1:points_read,2) + dy;
 plot(x,y);
 hold on;
 clear x;
 clear y;
end
   
