, j) );

               fprintf( [ '\n%s\n'                                       ...
                              'Surface point for grid point %d:\n'       ...
                              '  Latitudinal Coordinates:\n'             ...
                              '    Longitude           (deg): %12.6f\n'  ...
                              '    Latitude            (deg): %12.6f\n'  ...
                              '    Radius              (km): %12.6f\n' ],...
                              char( metho