if ( found )
[r, pclon, pclat] = cspice_reclat( point );
re = radii(1);
rp = radii(3);
f = ( re - rp ) / re;
[pdlon, pdlat, alt] = cspice_recgeo( point, re, f );
[phase, solar, emissn] = cspice_illum( satnm, et, ABCORR, scnm, point );
fprintf( 'Intercept planetocentric longitude (deg): %11.6f\n',...
fprintf( 'Intercept planetocentric latitude (deg): %11.6f\n',...
fprintf( 'Intercept planetodetic longitude (deg): %11.6f\n',...
fprintf( 'Intercept planetodetic latitude (deg): %11.6f\n',...
fprintf( 'Range from spacecraft to intercept point (km): %11.6f\n',...
fprintf( 'Range from spacecraft to target center (km): %11.6f\n',...
fprintf( 'Intercept phase angle (deg): %11.6f\n',...
fprintf( 'Intercept solar incidence angle (deg): %11.6f\n',...
fprintf( 'Intercept emission angle (deg): %11.6f\n',...
else
disp( ['No intercept point found at ' time ] )
end