Superimposing contour plot on pcolor

I am looking to superimpose a contour plot onto a pcolor plot that I have.

The pcolor plot is attached. I'd like to add a contour plot of topography, such that all points below the surface topography (defined by the y axis of the pcolor plot) are filled black (thereby blocking the background pcolor plot in that area). I'm unsure how to go about this, so any help would be much appreciated.

Answer 1

[Jess Lovering]

To clarify, are you saying that you want to impose a black area over regions that are within a certain value range?

If so, then you can use [C,h] = countour(X,Y,Z,[val val]) where "val" is the number you are using as your cut-off value. If it is zero then [C,h] = countour(X,Y,Z,[0 0]) will return a value in C that is the contour line points. You can use that with the area function (ex: area(C(1,:),C(2,:))) to generate a mask of black for those areas. You will have to modify this to fit your exact work but it could be a good starting point.

[Bruno Rodriguez]

The variable being plotted in the pcolor plot is reflectivity, and it's plotted in the horizontal and vertical 2D plane. However, what I'm wanting is to superimpose a field whose variable is height (as a function of horizontal distance) onto this pcolor plot. Since it's an entirely different variable to that being plotted in the pcolor plot, I don't believe this approach would work...