function plota_base(parametros,malha,NP,...
    p,A,base,pasta,plota,estagio)
fprintf(2,'--- inicio de plota_base -----------------------\n');
nb = size(base.elementos,2);
np = size(p,1);
%fprintf(2,'%d elementos na base\n',nb);
%fprintf(2,'%d pontos de amostragem\n',np);
%intervalo de elementos da a plotar
ib = escolhe_elemento_central(base,malha.CD);
fb = ib;

if plota.base
    for k = ib:fb
    %for k = 1:nb
        figure(2+plota.janela_base);
        hold off
        if plota.curva_nivel
            plot_contour_funcao(parametros,base,p,A(:,k),malha.CD,malha.LD,NP);
        else
            plota_funcao(p,A(:,k),true,malha.CD,malha.LD,NP,NaN,NaN);
        end
        title(sprintf('elemento %d tipo %d %s',...
            k,parametros.tipobase,estagio),'FontSize',20);
        if plota.salva_eps
            saveas(gcf,sprintf('%s/psi%04d.eps',pasta,k));
        end
        hold off
        
        %hold on
        figure(1 + plota.janela_base);
        plota_suporte(A(:,k),base.elementos(k),parametros.tensorial,p,...
            malha,pasta,plota,estagio);
%         figure(200 + plota.janela_base);
%         desenha_centros(base);
        %hold off
    end
end
