ellipt_compliant_plot_sim_data.m

function phs = ellipt_compliant_plot_sim_data(model,model_data,sim_data,params)
% function plots the results of the detailed simulation
%
%input:                 model: 
%                       model_data: grid information (p,e,t)
%                       sim_data:
%output:                plot

%shortcuts
p=model_data.grid.p;
e=model_data.grid.e;
t=model_data.grid.t;
u=sim_data.u;
%sim_data.plot_title=
ti = 'Plot of a detailed simulation';

phs = [];
phs2 = [];

if model.plot_isolines == 1
   phs = pdeplot(p,e,t,'xydata',u,'colormap','jet', 'contour','on', 'title',ti);
%   disp('plotted solution isolines');
%   keyboard;
   % das handle succhen, das eine typ 'line' ist und isempty(get(ph),'Zdata');
                
        htypes=get(phs(:),'Type');
        checkline=strcmp('line',htypes);
        lineindex=find(checkline==1);

        hline=phs(lineindex);
        for i=1:length(hline)
                checkcont(i)=isempty(get(hline(i),'ZData'))';
        end
        
        contindex=find(checkcont == 0);
        
        hcont=hline(contindex);

        %%%%definition of contour lines
        set(hcont,'LineStyle','-')
        set(hline(find(checkcont==1)),'LineStyle',':','Linewidth',1)
else
   phs = pdeplot(p,e,t,'xydata',u,'colormap','jet', 'title',ti);
   disp('plotted solution, no isolines');
end

if model.plot_mesh == 1
        figure;
   phs2 = pdemesh(p,e,t);
   disp('plotted mesh');
end

phs = [phs(:); phs2(:)];

%keyboard;