SUBROUTINE ERI_GETIDS( ID, INDEX, RMFID, ARFID, STATUS )
*+
* Name:
* ERI_GETIDS
* Purpose:
* Construct the ERI data objects given a dataset identifier
* Language:
* Starlink Fortran
* Invocation:
* CALL ERI_GETIDS( ID, INDEX, RMFID, ARFID, STATUS )
* Description:
* Reads in the energy response from the supplied dataset. The response
* consists of at least a redistribution matrix. If that matrix is
* normalised then an area response is also required. The returned
* objects have the classes RedistributionMatrix and AreaResponse, or
* are derived from those classes.
* Arguments:
* ID = INTEGER (given)
* ADI identifier of the dataset object
* INDEX = INTEGER (given)
* Index of spectrum in dataset
* RMFID = INTEGER (returned)
* ADI identifier of the RestributionMatrix derived object
* ARFID = INTEGER (returned)
* ADI identifier of the AreaResponse derived object
* STATUS = INTEGER (given and returned)
* The global status.
* Examples:
* {routine_example_text}
* {routine_example_description}
* Pitfalls:
* {pitfall_description}...
* Notes:
* {routine_notes}...
* Prior Requirements:
* {routine_prior_requirements}...
* Side Effects:
* {routine_side_effects}...
* Algorithm:
* {algorithm_description}...
* Accuracy:
* {routine_accuracy}
* Timing:
* {routine_timing}
* External Routines Used:
* {name_of_facility_or_package}:
* {routine_used}...
* Implementation Deficiencies:
* {routine_deficiencies}...
* References:
* ERI Subroutine Guide : http://www.sr.bham.ac.uk/asterix-docs/Programmer/Guides/eri.html
* Keywords:
* package:eri, usage:public
* Copyright:
* Copyright (C) University of Birmingham, 1995
* Authors:
* DJA: David J. Allan (Jet-X, University of Birmingham)
* {enter_new_authors_here}
* History:
* 27 Feb 1995 (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
INCLUDE 'ADI_PAR'
INCLUDE 'AST_PKG'
* Arguments Given:
INTEGER ID, INDEX
* Arguments Returned:
INTEGER ARFID, RMFID
* Status:
INTEGER STATUS ! Global status
* External references:
LOGICAL AST_QPKGI
EXTERNAL AST_QPKGI
* Local Variables:
CHARACTER*12 PNN ! Property name & no.
INTEGER FILID ! Base file object
INTEGER INDID ! ADI copy of INDEX
INTEGER PNC ! Length of PNN
INTEGER RESID ! Method output data
LOGICAL RTHERE ! .RESPnn exists?
LOGICAL THERE ! Component exists?
*.
* Check inherited global status.
IF ( STATUS .NE. SAI__OK ) RETURN
* Check initialised
IF ( .NOT. AST_QPKGI( ERI__PKG ) ) CALL ERI0_INIT( STATUS )
* Initialise return values
RMFID = ADI__NULLID
ARFID = ADI__NULLID
* Construct the property name number
PNN = '.RESP'
IF ( INDEX .EQ. 0 ) THEN
PNC = 5
ELSE
CALL CHR_ITOC( INDEX, PNN(6:), PNC )
PNC = PNC + 5
END IF
* Is response already computed?
CALL ADI_THERE( ID, PNN(:PNC), RTHERE, STATUS )
IF ( RTHERE ) THEN
CALL ADI_FIND( ID, PNN(:PNC), RESID, STATUS )
ELSE
* Get base file object
CALL ADI_GETFILE( ID, FILID, STATUS )
* Simply invoke the ReadRMF method
IF ( INDEX .EQ. 0 ) THEN
CALL ADI_EXEC( 'ReadRMF', 1, FILID, RESID, STATUS )
ELSE
CALL ADI_NEWV0I( INDEX, INDID, STATUS )
CALL ADI_EXEC2( 'ReadRMF', FILID, INDID, RESID, STATUS )
CALL ADI_ERASE( INDID, STATUS )
END IF
* If ok, extract results
IF ( (STATUS .EQ. SAI__OK) .AND. (RESID.NE.ADI__NULLID) ) THEN
* Store the returned object on the property list of the object
CALL ADI_CPUTID( ID, PNN(:PNC), RESID, STATUS )
END IF
END IF
* Locate sub-components for callee
CALL ADI_THERE( RESID, 'RMF', THERE, STATUS )
IF ( THERE ) THEN
CALL ADI_FIND( RESID, 'RMF', RMFID, STATUS )
END IF
CALL ADI_THERE( RESID, 'ARF', THERE, STATUS )
IF ( THERE ) THEN
CALL ADI_FIND( RESID, 'ARF', ARFID, STATUS )
END IF
* Release response structure if found
IF ( RTHERE ) THEN
CALL ADI_ERASE( RESID, STATUS )
END IF
* Report any errors
IF ( STATUS .NE. SAI__OK ) CALL AST_REXIT( 'ERI_GETIDS', STATUS )
END