how do I plot contour lines on orthographic projection?

      H= 163842×4 double - solution;
      he=H(:,4);
      lam = linspace(-pi,pi,3901);
      theta = linspace(-pi/2,pi/2,42);
      [L,T] = meshgrid(lam, theta);
      Z = reshape(he,3901,42);
      figure;
      axesm ortho
      geoshow(L*180/pi,T*180/*pi,Z,'DisplayType','contour')
      demcmap(Z)
      daspectm('m',1)

My result is a blank figure. Please help me out. I want my solution to be exactly like figure below. Your prompt response will be highly appreciated.