MATLAB: Pcolor

Question

Hi dear MATLAB user,

Im making a FDTD (acoustic preassure) model and at the end of each calculation representing the preassure. The problem is by wieving the contours, it only represents the actual preassure but if the wave hasnt arrived yet you wont see nothing from the contours.

for ...
for ...
end
...
    pcolor((1:nx)*dh,(1:ny)*dh,10*log10(min(max(pp'.^2,1d-6),1)) );
    shading flat,axis equal; colorbar; title(['time=' num2str(round((tt+maxtt*(ttt-1))*1000*dt)) ' ms']);drawnow
    hold on
    [x,y]=find(mapaz == 0);
    contorno=plot(x,y,'+k');
    set(contorno,'MarkerSize',10);
    hold off
    drawnow
end

Im trying to draw the contours stored in matrix 'mapaz': 0 = contour; 1 nothing.

nx and ny are the sizes of the matrix (1600×1200)

I really don't know to mix this together. Could you help me?

Thanks, Daniel Serra

Answer 1

Perhaps

contour((1:nx)*dh, (1:ny)*dh, 0 + (mapaz ~= 0), [1 1])