mel(lon) );

         for i=1:numel(lat)
            for j=1:numel(lon)

               n = n+1;

               grid(1,n) = lon(j);
               grid(2,n) = lat(i);

            end
         end

         grid = grid * <a href="cspice_rpd.html">cspice_rpd</a>();

         %
         % Find the surface points corresponding to the grid points.
         %
         % Compute outward normal vectors at the surface points,
         % using both surface representations.
         %
         for i = 1:2
          