      ...
                          '        The emission and phase '    ...
                          'angles should be equal.\n' ] );
            end

            fprintf( [ '\n'                                    ...
                       '      Illumination angles at the '  ...
                       'sub-s/c point:\n' ...
                       '\n'...
                       '      Phase angle            (deg): %15.9f\n'...
                       '      Solar incidence angle  (deg): %15.9f\n'...
     