  srfpts = <a href="cspice_latsrf.html">cspice_latsrf</a>( method(i), target, et, fixref, grid);

            normls = <b>cspice_srfnrm</b>( method(i), target, et, fixref, srfpts);

            for j=1:n

               %
               % Use <a href="cspice_recrad.html">cspice_recrad</a> rather than <a href="cspice_reclat.html">cspice_reclat</a> to produce
               % non-negative longitudes.
               %
               [ xr, xlon, xlat] = <a href="cspice_recrad.html">cspice_recrad</a>( srfpts(1:3