SUBROUTINE ADI2_CHGIMG( HDUID, NDIF, NDIM, DIMS, STATUS )
*+
* Name:
* ADI2_CREIMG
* Purpose:
* Define a new FITS extension called HDU
* Language:
* Starlink Fortran
* Invocation:
* CALL ADI2_CHGIMG( HDUID, NDIM, DIMS, STATUS )
* Description:
* {routine_description}
* Arguments:
* HDUID = INTEGER (given)
* The ADI identifier of the FITShdu
* NDIF = INTEGER (given)
* Amount by which this extension is larger than the old extension
* NDIM = INTEGER (given)
* Number of dimensions of image extension
* DIMS[NDIM] = INTEGER (given)
* Sizes of the dimensions
* 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:
* 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:
* 28 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
* Arguments Given:
INTEGER HDUID, NDIF, NDIM, DIMS(*)
* Status:
INTEGER STATUS ! Global status
* Local Variables:
CHARACTER*1 STR ! Character of IAX
INTEGER FSTAT ! FITSIO status
INTEGER IAX ! Axis loop
INTEGER LUN ! Logical unit
INTEGER NELM ! New # elements
*.
* Check inherited global status.
IF ( STATUS .NE. SAI__OK ) RETURN
* Get the logical unit
CALL ADI2_HDULUN( HDUID, LUN, STATUS )
* If extension is being shortened fill the remaining data area with zeroes
IF ( NDIF .LT. 0 ) THEN
CALL ARR_SUMDIM( NDIM, DIMS, NELM )
FSTAT = 0
CALL FTPPRU( LUN, 1, NELM + NDIF + 1, -NDIF, FSTAT )
END IF
* Update the dimension keywords
FSTAT = 0
CALL ADI2_HPKYI( HDUID, '@NAXIS', NDIM, ' ', STATUS )
CALL FTMKYJ( LUN, 'NAXIS', NDIM, '&', FSTAT )
DO IAX = 1, NDIM
WRITE( STR, '(I1)' ) IAX
CALL ADI2_HPKYI( HDUID, '@NAXIS'//STR, DIMS(IAX), ' ', STATUS )
CALL FTMKYJ( LUN, 'NAXIS'//STR, DIMS(IAX), '&', FSTAT )
END DO
* Redefine IMAGE extension
IF ( FSTAT .EQ. 0 ) THEN
CALL FTRDEF( LUN, FSTAT )
END IF
IF ( FSTAT .NE. 0 ) CALL ADI2_FITERP( FSTAT, STATUS )
* Mark HDU as changed
CALL ADI_CPUT0L( HDUID, '.Changed', .TRUE., STATUS )
* Report failure
89 IF ( STATUS .NE. SAI__OK ) THEN
CALL ERR_REP( ' ', 'Unable to resize image extension', STATUS )
END IF
* Report any errors
IF ( STATUS .NE. SAI__OK ) THEN
CALL AST_REXIT( 'ADI2_CHGIMG', STATUS )
END IF
END