SUBROUTINE ADI2_PARSE( FSPEC, LFILEC, HDU, FPATH, FPLEN, STATUS )
*+
* Name:
* ADI2_PARSE
* Purpose:
* Parse a FITS file specification.
* Language:
* Fortran
* Invocation:
* CALL ADI2_PARSE( FSPEC, LFILEC, HDU, FPATH, STATUS )
* Description:
* Attempts to open the named file object as an HDS file. If successful
* the routine stores the logical unit and HDU number on the property
* list of the ID object.
* Arguments:
* FSPEC = CHAR (Given)
* Name of the object on which FITS access to be attempted
* LFILEC = INTEGER (Returned)
* Index of last character in filemname proper
* HDU = INTEGER (Returned)
* HDU number, zero if none
* FPATH = CHARACTER*(*) (Returned)
* Sub-file path
* FPLEN = INTEGER
* Length of FPATH, or zero if none (Returned)
* STATUS = INTEGER (Given and returned)
* The global status.
* Examples:
* {routine_example_text}
* {routine_example_description}
* Notes:
* {routine_notes}...
* Prior Requirements:
* {routine_prior_requirements}...
* Side Effects:
* {routine_side_effects}...
* Algorithm:
* {algorithm_description}...
* External Routines Used:
* {name_of_facility_or_package}:
* {routine_used}...
* Implementation Deficiencies:
* {routine_deficiencies}...
* References:
* ADI Subroutine Guide : http://www.sr.bham.ac.uk/asterix-docs/Programmer/Guides/adi.html
* Keywords:
* package:adi, usage:private
* Copyright:
* Copyright (C) University of Birmingham, 1995
* Authors:
* DJA: David J. Allan (JET-X,University of Birmingham)
* {enter_new_authors_here}
* History:
* 15 Jul 1994 (DJA):
* Original version.
* {enter_changes_here}
* Bugs:
* {note_any_bugs_here}
*-
* Type Definitions:
IMPLICIT NONE ! No implicit typing
* Global Constants:
INCLUDE 'SAE_PAR' ! Standard SAE constants
* Arguments Given:
CHARACTER*(*) FSPEC
* Arguments Returned:
INTEGER LFILEC, HDU, FPLEN
CHARACTER*(*) FPATH
* Status:
INTEGER STATUS ! Global status
* External references:
EXTERNAL CHR_LEN
INTEGER CHR_LEN
EXTERNAL CHR_ISDIG
LOGICAL CHR_ISDIG
* Local Variables:
INTEGER FLEN ! Length of FSPEC
INTEGER ICP ! Character index
INTEGER NCP ! Character index
*.
* Check inherited global status.
IF ( STATUS .NE. SAI__OK ) RETURN
* Length of FSPEC
FLEN = CHR_LEN( FSPEC )
* Look for % delimiter, so that trailing representation can be ignored
ICP = FLEN
DO WHILE ( (ICP.GT.0) .AND. (FSPEC(ICP:ICP).NE.'%') )
ICP = ICP - 1
END DO
IF ( ICP .NE. 0 ) FLEN = ICP - 1
* Default values
LFILEC = FLEN
HDU = 0
FPLEN = 0
* HDU delimiter present?
ICP = INDEX( FSPEC, '[' )
IF ( ICP .EQ. 0 ) THEN
ICP = INDEX( FSPEC, '+' )
IF ( ICP .GT. 0 ) THEN
* Skip over digits of HDU number
ICP = ICP + 1
NCP = ICP
DO WHILE ( CHR_ISDIG( FSPEC(NCP:NCP) ) )
NCP = NCP + 1
END DO
NCP = NCP - 1
* If no digits signal error
IF ( NCP .LT. ICP ) THEN
STATUS = SAI__ERROR
CALL ERR_REP( ' ', 'Missing HDU number in extension '/
: /'specification', STATUS )
END IF
END IF
ELSE
ICP = ICP + 1
NCP = INDEX( FSPEC, ']' )
IF ( NCP .EQ. 0 ) THEN
STATUS = SAI__ERROR
CALL ERR_REP( ' ', 'Missing "]" in HDU extension '/
: /'specification', STATUS )
ELSE IF ( NCP .LE. ICP ) THEN
STATUS = SAI__ERROR
CALL ERR_REP( ' ', 'Misplaced "]" in HDU extension '/
: /'specification', STATUS )
ELSE
NCP = NCP - 1
END IF
END IF
IF ( STATUS .NE. SAI__OK ) GOTO 99
* Sub-HDU structure may only appear after an HDU specification. NCP marks
* the last character of the HDU specification
IF ( ICP .GT. 0 ) THEN
* Last character in the filename proper
LFILEC = ICP - 2
* Extract HDU number
CALL CHR_CTOI( FSPEC(ICP:NCP), HDU, STATUS )
IF ( (STATUS .NE. SAI__OK) .OR. (HDU.LT.1) ) THEN
CALL MSG_SETC( 'NUM', FSPEC(ICP+1:NCP) )
CALL ERR_REP( ' ', 'Illegal HDU number /^NUM/', STATUS )
END IF
* Encode the HDU number into FPATH
CALL CHR_ITOC( HDU, FPATH(2:), FPLEN )
FPATH(1:1) = '+'
FPLEN = FPLEN + 1
* Remaining text after end of HDU spec?
IF ( FSPEC(NCP+1:NCP+1) .EQ. ']' ) THEN
FPATH(FPLEN+1:) = FSPEC(NCP+2:FLEN)
FPLEN = FPLEN + FLEN - NCP - 1
ELSE
FPATH(FPLEN+1:) = FSPEC(NCP+1:FLEN)
FPLEN = FPLEN + FLEN - NCP
END IF
CALL CHR_UCASE( FPATH )
END IF
IF ( STATUS .NE. SAI__OK ) GOTO 99
* Look for sub-file specification
99 CONTINUE
END