Problem About Legend And Axis Of 3D Map Figure

Hi everyone,

I am trying to plot 3D Map by using the code below. However, I cannot display the axis as latitude and longitude values, and also the legend command employed in the code below doesn' t work for showing the color defined below in the legend of figure. Here is my code;

Latitude=[38.42 38.7043 38.843 38.8335];
Longitude=[43.5710 43.4535 43.4585 43.5795];
Depth=[21.9 21 25.6 21.8];
D={'20<=D<30'};
Co3={'b'};
ax = axesm('mercator','MapLatLimit',[37.400 39.200],'MapLonLimit',...
               [42.000 44.400]); view(3);
land = shaperead('landareas', 'UseGeoCoords', true);
geoshow(ax,land, 'FaceColor', 'yellow');
stem3m(Latitude,Longitude,Depth*-0.001,...
       'o', 'MarkerSize', 5,'MarkerFaceColor', [0 0 1],'LineStyle','None');
mlabel('on');
title('Region');
h=legend(D,'Location','SouthEast');
v = get(h,'title');
set(v,'string','Depth Range');

How can I solve this problem?