function mapmaskmanuel(x,y,a,ubar,vbar,vmin,vmax,lc,Nlc)

%------------------------------------
% display a filled contour plot of the matrix a
% with streamlines based upon ubar and vbar
% ------------------------------------
% MARS 2008 : deboggage sur la base du travail de Julie

% DECEMBRE 2016: ajout de l'option 'lc' (ligne de courant)
% si lc == 1, on trace les lignes de courant
% si lc == 0, on ne trace pas les lignes de courant


colormap(jet(64))
v=linspace(vmin,vmax);
%  v=linspace(0,0.07);

contourf(x,y,double(a),v,'LineStyle','none');
% colorbar;
% xlabel('x/H');
% ylabel('Z/H');
% cm=gray;            % pour des niveaux en gris
cm=colormap(jet);   % pour des niveaux en couleur
colormap(cm);
% cl=colorbar;
% ylabel(cl, 'k/(U_{deb})^{2}')
% ylabel(cl, 'U/U_{deb}')

axis equal tight


if lc == 0
    
elseif lc == 1 
    h=streamslice(x,y,ubar,vbar,Nlc);
    set(h,'Color','black')
    % xlabel('x (cm)')
    % ylabel('y (cm)')
end


% colormap(jet(64));
% pcolor(x,y,a);
% shading interp;
% colorbar;

% xlabel('cm');

% ylabel('cm');

%figure;
%pcolor(x, y, a); shading interp; 
% Conserve le contenu de la fenetre graphique :
%hold on; 
% Trace les isocontours du champ z2 en trait noir
%contour(x,y,ubar,10,'k');

