SUBROUTINE TCI2_WRITE( NARG, ARGS, OARG, STATUS )
*+
* Name:
* TCI2_WRITE
* Purpose:
* Write hardware description to a FITS file
* Language:
* Starlink Fortran
* Invocation:
* CALL TCI2_WRITE( NARG, ARGS, OARG, STATUS )
* Description:
* Write data to a FITS describing timing information
* Arguments:
* NARG = INTEGER (given)
* Number of method arguments
* ARGS(*) = INTEGER (given)
* ADI identifier of method arguments
* OARG = INTEGER (returned)
* Output data
* 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:
* TCI Subroutine Guide : http://www.sr.bham.ac.uk/asterix-docs/Programmer/Guides/tci.html
* Keywords:
* package:tci, usage:private
* Copyright:
* Copyright (C) University of Birmingham, 1995
* Authors:
* DJA: David J. Allan (Jet-X, University of Birmingham)
* {enter_new_authors_here}
* History:
* 6 Mar 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'
* Arguments Given:
INTEGER NARG ! # of input arguments
INTEGER ARGS(NARG) ! Input arguments
* Arguments Returned:
INTEGER OARG ! Output structure
* Status:
INTEGER STATUS ! Global status
* Local Variables:
CHARACTER*8 CONTNT ! Content of file
*.
* Check inherited global status.
IF ( STATUS .NE. SAI__OK ) RETURN
* Initialise
OARG = ADI__NULLID
* Write to SPECTRUM too
CALL ADI2_GKEY0C( ARGS(1), ' ', 'CONTENT', .FALSE., .FALSE.,
: CONTNT, ' ', STATUS )
IF ( STATUS .EQ. SAI__OK ) THEN
* OGIP spectrum?
IF ( CONTNT .EQ. 'SPECTRUM' ) THEN
* Write data to primary HDU and SPECTRUM extension
CALL TCI2_WRITE1( NARG, ARGS, ' ', STATUS )
CALL TCI2_WRITE1( NARG, ARGS, 'SPECTRUM', STATUS )
* OGIP response matrix?
ELSE IF ( CONTNT .EQ. 'MATRIX' ) THEN
CALL TCI2_WRITE1( NARG, ARGS, 'MATRIX', STATUS )
CALL TCI2_WRITE1( NARG, ARGS, 'EBOUNDS', STATUS )
END IF
* Write data to primary HDU if CONTENT not defined
ELSE
CALL ERR_ANNUL( STATUS )
END IF
* Report any errors
IF ( STATUS .NE. SAI__OK ) CALL AST_REXIT( 'TCI2_WRITE', STATUS )
END
SUBROUTINE TCI2_WRITE1( NARG, ARGS, HDUNAM, STATUS )
*+
* Name:
* TCI2_WRITE1
* Purpose:
* Write timing description to a specific FITS file HDU
* Language:
* Starlink Fortran
* Invocation:
* CALL TCI2_WRITE1( NARG, ARGS, HDUNAM, STATUS )
* Description:
* Write the data to a FITS file describing timing
* Arguments:
* NARG = INTEGER (given)
* Number of method arguments
* ARGS(*) = INTEGER (given)
* ADI identifier of method arguments
* HDUNAM = CHARACTER*(*) (given)
* HDU to write to
* 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:
* TCI Subroutine Guide : http://www.sr.bham.ac.uk/asterix-docs/Programmer/Guides/tci.html
* Keywords:
* package:tci, usage:private
* Copyright:
* Copyright (C) University of Birmingham, 1995
* Authors:
* DJA: David J. Allan (Jet-X, University of Birmingham)
* {enter_new_authors_here}
* History:
* 6 Mar 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'
* Arguments Given:
INTEGER NARG ! # of input arguments
INTEGER ARGS(NARG) ! Input arguments
CHARACTER*(*) HDUNAM ! HDU name
* Status:
INTEGER STATUS ! Global status
* Local Variables:
CHARACTER*8 DSTR, TSTR ! Date and time strings
DOUBLE PRECISION MJD !
DOUBLE PRECISION OBLEN !
INTEGER FID, TIMID ! Extracted arguments
LOGICAL LOK, MOK ! Things present?
*.
* Check inherited global status.
IF ( STATUS .NE. SAI__OK ) RETURN
* Initialise
FID = ARGS(1)
TIMID = ARGS(2)
* Start date/time of observation
CALL ADI_THERE( TIMID, 'MJDObs', MOK, STATUS )
IF ( MOK ) THEN
CALL ADI_CGET0D( TIMID, 'MJDObs', MJD, STATUS )
CALL TCI_MJD2DT( MJD, DSTR, TSTR, STATUS )
CALL ADI2_PKEY0D( FID, HDUNAM, 'MJD-OBS', MJD,
: 'MJD of the data start', STATUS )
CALL ADI2_PKEY0C( FID, HDUNAM, 'DATE-OBS', DSTR,
: 'Date of data start', STATUS )
CALL ADI2_PKEY0C( FID, HDUNAM, 'TIME-OBS', TSTR,
: 'Time of data start', STATUS )
END IF
* Observation length ok?
CALL ADI_THERE( TIMID, 'ObsLength', LOK, STATUS )
IF ( MOK .AND. LOK ) THEN
* Get length of observation
CALL ADI_CGET0D( TIMID, 'ObsLength', OBLEN, STATUS )
* Add to start time
MJD = MJD + OBLEN / 86400D0
CALL TCI_MJD2DT( MJD, DSTR, TSTR, STATUS )
CALL ADI2_PKEY0D( FID, HDUNAM, 'MJD-END', MJD,
: 'MJD of the data end', STATUS )
CALL ADI2_PKEY0C( FID, HDUNAM, 'DATE-END', DSTR,
: 'Date of data end', STATUS )
CALL ADI2_PKEY0C( FID, HDUNAM, 'TIME-END', TSTR,
: 'Time of data end', STATUS )
END IF
* Clock reference date/time
CALL ADI_THERE( TIMID, 'MJDRef', MOK, STATUS )
IF ( MOK ) THEN
CALL ADI_CGET0D( TIMID, 'MJDRef', MJD, STATUS )
CALL TCI_MJD2DT( MJD, DSTR, TSTR, STATUS )
CALL ADI2_PKEY0D( FID, HDUNAM, 'MJDREF', MJD,
: 'MJD SC clock start', STATUS )
CALL ADI2_PKEY0C( FID, HDUNAM, 'ZERODATE', DSTR,
: 'Date of data start', STATUS )
CALL ADI2_PKEY0C( FID, HDUNAM, 'ZEROTIME', TSTR,
: 'Time of data start', STATUS )
END IF
* Report any errors
IF ( STATUS .NE. SAI__OK ) CALL AST_REXIT( 'TCI2_WRITE1', STATUS )
END