function [] = femPostPlot(femsol,filename)

N=size(femsol,2);

for n=1:N
    fem=femsol(n);
    postplot(fem, ...
         'tridata',{'Hz','cont','internal','unit','A/m'}, ...
         'trimap','jet(1024)', ...
         'title',['Surface: mag. field Hz [A/m] iter=' num2str(n)], ...
         'axis',[-2E-6,2E-6,-2E-6,2E-6]); 
     if (~isempty(filename))
         saveas(gcf,[filename num2str(n) '.png']);
     end
end