adding contour lines to pcolor plot

Hi everyone,

I'm trying to add contour lines at 13 degrees to my pcolor plot, but I get the following error:

Error using contour

Input arguments must be numeric or objects which can be converted to double.

Here's my code:

t1 = datetime(2021, 5, 1)
t2 = datetime(2021, 9, 30)
dateVector = t1:t2

combinedaverages = [dailyaverages; dailyaverages06; dailyaverages07; dailyaverages08; dailyaverages09];

figure
pcolor(dateVector, depth, combinedaverages')
set(gca, 'YDir','reverse')
shading flat
ylim([0 70])
colorbar
hold on
contour(dateVector,depth,combinedaverages',[13 13],'k');

Variables in the code are:

combine averages: 153x33 double

depth: 33x1 single

dateVector: 1x153 datetime

[Mathieu NOE]

hello

it's good practice to share the data (as mat file) along the code (if you want someone to help you)

all the best

[Giulia]

Hi, the data are about 150 netcd files..how do I share it all?

[Adam Danz]

The data are the variables you mentioned in your question. If you save them to a mat file and attach them we can run the code in your question.

[Giulia]

Hi Adam, thank you.

Here are my mat files:

Answer 1

[Mathieu NOE]

hello again

maybe this ?

figure
pcolor(dateVector, depth, combinedaverages')
set(gca, 'YDir','reverse')
shading flat
ylim([0 70])
colorbar
hold on
contour(1:numel(dateVector),depth,combinedaverages','k');

[Giulia]

Hi Mathieu, thanks a lot for your response!

How do I apply the code to set contour lines at 10 and 20 m?

[Mathieu NOE]

hello Giulia

see code modification below

clc
figure
pcolor(dateVector, depth, combinedaverages')
set(gca, 'YDir','reverse')
shading flat
ylim([0 70])
colorbar
hold on
V = [10 20]; % set contour lines at 10 and 20 m
contour(1:numel(dateVector),depth,combinedaverages',V,'k');

but the contour line at 20 m will not show up as your Z data (combinedaverages) range spans from 9.0550 to 17.2330

[Giulia]

Great! Thank you so much again for your help!

All the best

[Mathieu NOE]

My pleasure !