Inserting black lines on pcolor plot

I have a matrices of values which I want to use to create a pcolor plot.

My data is constrained by water depths where I have set all depth >100m to NaN and all values of 0 to NaN.

This has created a confusing plot where I'm unable to show coastal outlines properly.

How would I change this plot so that I retain my '0 depth' or land values but create a black outline that will clearly show the outline of the coasts.

Thanks.

Answer 1

[KSSV]

You can use contour.

Let X, Y, H be your data of (x,y) coordinates and H be your depth.

figure
hold on
pcolor(X,Y,H)
shading interp
contour(X,Y,H,[0 0],'k')

[Charlie Milford]

Thank you for your answer it works perfectly and will use in future. I ended up using the m_map addon and imported a coastal data zip file in the end.