bang=round(it/20);if bang==it/20; figure(2);subplot(111);imagesc([1:nx]*dx,[1:nz]*dx,p');
title(['Pressure Snapshot at Time Step ',num2str(it*dt),' secs']);
xlabel('X (ft)');ylabel('Z (ft)');colorbar;
pause(.015);end