How can I create a contourf with costline?

Dear All, I have a 2D matrix about precipitation values in Europe. I would like to create a contourf from this matrix with costline. I tried to the followings:

h = worldmap('Europe');
geoshow('landareas.shp', 'FaceColor', [1 1 1])
hold on
axesm('eqdconic','maplatlimit',[26.5861 71.8699],'maplonlimit',[-34.748 60.8422]);
>> contourf(lon,lat,IU2);

But unfortunately it does not work. Can someone suggest me a solution for creating contourf with costline? I would appreciate if someone helped me.

Answer 1

[Szabó-Takács Beáta]

Meantime I found a solution.

worldmap('Europe');
load coast;
plotm(lat, long, 'k');
clear lat long;
lat=load(lat.mat)
lon=load(lon.mat)
contourfm(lat,lon,IU2);