function plotmap(LL,meshele,meshnod,bathy,arg5);
%  plotmap(L,meshele,meshnod,scalarnod,'interp');
%  Makes color contour map of scalar
% plotmap(L,meshele,meshnod,bathy(:,2),'interp');
% or  bit faster don't do the shading interp:
% plotmap(L,meshele,meshnod,bathy(:,2));

L(1) = LL(1) - .1*(LL(2) -LL(1));
L(2) = LL(2) + .1*(LL(2) -LL(1));
L(3) = LL(3) - .1*(LL(4) -LL(3));
L(4) = LL(4) + .1*(LL(4) -LL(3));
clf
tic
%x = meshnod(:,2);
%y = meshnod(:,3);
x3 = meshele(:,2:4);
x3(:)=meshnod(meshele(:,2:4),2);
y3 = meshele(:,2:4);
y3(:)=meshnod(meshele(:,2:4),3);
field3 = meshele(:,2:4);
field3(:)=bathy(meshele(:,2:4));
toc

disp('find the inside numbers')
tic
lplot = find(x3(:,1)>L(1) & x3(:,1)<L(2));
lplot = lplot(find(y3(lplot,1)>L(3) & y3(lplot,1)<L(4) ));
%x3=x3(lplot,:)';
%y3=y3(lplot,:)';
%field3=field3(lplot,:)';
toc
disp(['Do the plot'])
size(x3(lplot,:))
fill(x3(lplot,:)',y3(lplot,:)',field3(lplot,:)')
  if nargin==5
    shading interp
  end
axis('equal')
axis(LL)
caxis([min(field3(lplot,1)) max(field3(lplot,1))]);

Hcolorbar = colorbar('horiz');

