function movie01(colors, isobaths)
%
% the function movie01(colors, isobaths) is used by the ARPRO
% software to plot a base map on which later to overlay float
% tracks.  The input parameters should be as follows:
%        colors :   'yes', if we want a color map, with cyan
%                          sea and yellow land, and (if we 
%                          choose so,) blue isobaths, or
%                   'no',  for black sea and white land and
%                          isobaths (if we print it, it's 
%                          reversed).
%        isobaths:  'yes' for isobath plotting, or
%                   'no'  for isobath omission.
%
% -------------------------------------------------------------
%  Written by Vassilis Zervakis, in February 1996,
%       modified by Mike Sparrow, June 1997
%             for AINCO - Interocean, Madrid, Spain
%--------------------------------------------------------------
%
global moviefigure
global minlon maxlon minlat maxlat
global linelength gnome 
global h1 timehandle
global firstday firstmonth firstyear
global COASTPATH BATHYPATH
global isobathy
%
datafound = 'no ';
% check whether the desired coordinates are within the ----
% existing capabilities -----------------------------------
% ----check whether we'll plot N.E. Atlantic map ---------
if (minlon>=-30)&(maxlon<=0)&(minlat>=20)&(maxlat<=60),
   dx = 0.1;
   dy = 0.1;
   coastfile=[ COASTPATH 'neatlant.map'];
   bathyfile=[ BATHYPATH 'neatlant.iso'];
   datafound = 'yes';
end
% ----check whether we'll plot N.W. Atlantic map ---------
if (minlon>=-100)&(maxlon<=-30)&(minlat>=0)&(maxlat<=65),
   dx = 0.1;
   dy = 0.1;
   coastfile=[ COASTPATH 'nwatlant.map'];
   bathyfile=[ BATHYPATH 'nwatlant.iso'];
   datafound = 'yes';
end
% ----check whether we'll plot S. Atlantic map ---------
if (minlon>=-70)&(maxlon<=20)&(minlat>=-60)&(maxlat<=0.),
   dx = 0.0;
   dy = 0.0;
   coastfile=[ COASTPATH 'satlant.map'];
   bathyfile=[ BATHYPATH 'satlant.iso'];
   datafound = 'yes';
end
% ----check whether we'll plot Mediterranean map--------
if (minlon>=-10)&(maxlon<=45)&(minlat>=30)&(maxlat<=50),
   dx = 0.1;
   dy = 0.1;
   coastfile=[ COASTPATH 'mediterr.map'];
   bathyfile=[ BATHYPATH 'mediterr.iso'];
   datafound = 'yes';
end
% ----check whether we'll plot Hellenic Seas map --------------
if (minlon>=+19)&(maxlon<=+29)&(minlat>=+34)&(maxlat<=+41),
   dx = 0.1;
   dy = 0.1;
   coastfile=[ COASTPATH 'shoregr.bln'];
   bathyfile=[ BATHYPATH 'mediterr.iso'];
   datafound = 'yes';
end
if datafound =='yes',
   if (colors == 'yes') & (isobaths == 'yes'), 
      global moviefigure ;
      moviefigure = figure('name','Movie',...
                         'Color',[ 0 0 0],...
                         'Resize','on',...
                         'Clipping','off',...
                         'NumberTitle','off');
      h = axes('Position',[0.1 0.2 0.8 0.7],...
               'Box','on');
      geoplot2(minlon,maxlon,minlat,maxlat,coastfile,bathyfile,isobathy, dx, dy);
      global seacolor
      seacolor = 'c';
   end
   if (colors == 'yes') & (isobaths == 'no '), 
      global moviefigure ;
      moviefigure = figure('name','Movie',...
                         'Color',[ 0 0 0],...
                         'Resize','on',...
                         'NumberTitle','off');
      h = axes('Position',[0.1 0.2 0.8 0.7],...
                'Box','on') ;
      map2(minlon,maxlon,minlat,maxlat,coastfile);
      global seacolor
      seacolor = 'c';
   end
   if (colors == 'no ') & (isobaths == 'yes'), 
      global moviefigure ;
      moviefigure = figure('name','Movie',...
                         'Resize','on',...
                         'NumberTitle','off');
      h = axes('Position',[0.1 0.2 0.8 0.7],...
               'Box','on');
      bwgeopl2(minlon,maxlon,minlat,maxlat,coastfile,bathyfile,isobathy,dx,dy,gnome);
      global seacolor
      seacolor = 'k';
   end
   if (colors == 'no ') & (isobaths == 'no '), 
      global moviefigure ;
      moviefigure = figure('name','Movie',...
                         'Resize','on',...
                         'NumberTitle','off');
      h = axes('Position',[0.1 0.2 0.8 0.7],...
                'Box','on');
      bwmap2(minlon,maxlon,minlat,maxlat,coastfile,gnome);
      global seacolor
      seacolor = 'k';
   end
   pos = [0.18 0.05 0.18 0.04];
   global h1 h9
   h1 = uicontrol('Style','pushbutton',...
                  'String','Replay',...
                  'Enable','off',...
                  'Units','normalized',...
                  'Position',pos);
   pos = [0.65 0.05 0.18 0.04];
   h9 = uicontrol('Style','pushbutton',...
                  'String','Plot dispersion',...
                  'Enable','off',...
                  'Units','normalized',...
                  'Position',pos);
   timehandle = uicontrol('Style','text',...
                  'Backgroundcolor',[0 0 0],...
                  'ForegroundColor',[1 1 1],...
                  'String',['Date :   ' num2str(firstday),...
                  '/' num2str(firstmonth) '/' num2str(firstyear)],...
                  'Units','normalized',...
                  'Position',[0.7 0.9 0.3 0.05]);
else,
   mbox('OK','Coordinates out of currently',...
             'available range');
end
