SUBROUTINE BDI2_IMMAP( NARG, ARGS, OARG, STATUS ) *+ * Name: * BDI2_IMMAP * Purpose: * Service FileItemMap(Spectrum,FITSfile,item,type,mode) requests * Language: * Starlink Fortran * Invocation: * CALL BDI2_IMMAP( NARG, ARGS, OARG, STATUS ) * Description: * Service FileItemGet(Spectrum,FITSfile,item) requests from the BDI system * 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: * BDI Subroutine Guide : http://www.sr.bham.ac.uk/asterix-docs/Programmer/Guides/bdi.html * Keywords: * package:bdi, usage:private * Copyright: * Copyright (C) University of Birmingham, 1995 * Authors: * DJA: David J. Allan (Jet-X, University of Birmingham) * {enter_new_authors_here} * History: * 9 Aug 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 ! # arguments INTEGER ARGS(*) ! Method arguments * Arguments Returned: INTEGER OARG ! Returned data * Status: INTEGER STATUS ! Global status * Local Variables: CHARACTER*20 ITEM CHARACTER*6 MODE ! Mapping mode CHARACTER*8 TYPE ! Mapping type DOUBLE PRECISION BASE ! SGP/38 style base DOUBLE PRECISION CDELT, CRPIX ! Axis keywords INTEGER DIM ! Size of spectrum INTEGER IAX ! An axis number INTEGER PTR ! Mapped data INTEGER PSID ! Private storage INTEGER IMHDU ! SPECTRUM hdu id *. * Check inherited global status. IF ( STATUS .NE. SAI__OK ) RETURN print *,1 call flush(6) * Extract the arguments CALL ADI_GET0C( ARGS(3), ITEM, STATUS ) CALL ADI_GET0C( ARGS(4), TYPE, STATUS ) CALL ADI_GET0C( ARGS(5), MODE, STATUS ) print *,2 call flush(6) * Locate the IMAGE hdu CALL ADI2_FNDHDU( ARGS(2), ' ', .FALSE.,IMHDU, STATUS ) print *,3 call flush(6) * Locate the BDI private storage for the item, creating if required CALL BDI0_LOCPST( ARGS(1), ITEM, .TRUE., PSID, STATUS ) print *,4 call flush(6) * Switch on the various items IF ( ITEM .EQ. 'Data' ) THEN * Simply map primary data CALL ADI2_MAPIMG( IMHDU, TYPE, MODE, PSID, PTR, STATUS ) print *,5 call flush(6) * Axis data ELSE IF ( (ITEM(:5) .EQ. 'Axis_') .AND. : (ITEM(8:).EQ.'Data') ) THEN * Get axis number CALL CHR_CTOI( ITEM(6:6), IAX, STATUS ) * Do standard keywords exist? CALL ADI2_HGKYID( IMHDU, 'CRPIX', IAX, CRPIX, STATUS ) CALL ADI2_HGKYID( IMHDU, 'CDELT', IAX, CDELT, STATUS ) IF ( STATUS .NE. SAI__OK ) THEN CALL ERR_ANNUL( STATUS ) CRPIX = DBLE(DIM)/2.0D0 + 0.5D0 CDELT = 1.0D0 END IF * Convert CRPIX to SGP/38 style base BASE = (1.0D0-CRPIX) * CDELT * Get axis size CALL ADI2_HGKYII( IMHDU, 'NAXIS', IAX, DIM, STATUS ) * Create dynamic array CALL DYN_MAPT( 1, DIM, TYPE, PTR, STATUS ) IF ( TYPE .EQ. 'REAL' ) THEN CALL ARR_REG1R( REAL(BASE), REAL(CDELT), DIM, : %VAL(PTR), STATUS ) ELSE IF ( TYPE .EQ. 'DOUBLE' ) THEN CALL ARR_REG1D( BASE, CDELT, DIM, %VAL(PTR), STATUS ) END IF * Axis widths ELSE IF ( (ITEM(:5) .EQ. 'Axis_') .AND. : (ITEM(8:).EQ.'Width') ) THEN * Get axis number CALL CHR_CTOI( ITEM(6:6), IAX, STATUS ) * Does standard keywords exist? CALL ADI2_HGKYID( IMHDU, 'CDELT', IAX, CDELT, STATUS ) IF ( STATUS .NE. SAI__OK ) THEN CALL ERR_ANNUL( STATUS ) CDELT = 1.0D0 END IF * Get axis size CALL ADI2_HGKYII( IMHDU, 'NAXIS', IAX, DIM, STATUS ) * Create dynamic array CALL DYN_MAPT( 1, DIM, TYPE, PTR, STATUS ) IF ( TYPE .EQ. 'REAL' ) THEN CALL ARR_INIT1R( REAL(CDELT), DIM, %VAL(PTR), STATUS ) ELSE IF ( TYPE .EQ. 'DOUBLE' ) THEN CALL ARR_INIT1D( CDELT, DIM, %VAL(PTR), STATUS ) END IF END IF print *,7 call flush(6) * Store details of mapping c CALL ADI2_STOMAP( PSID, IMHDU, 'I', PTR, TYPE, MODE, STATUS ) print *,8 call flush(6) * Release storage CALL ADI_ERASE( PSID, STATUS ) * Release the IMAGE hdu CALL ADI_ERASE( IMHDU, STATUS ) * Store the pointer CALL ADI_NEWV0I( PTR, OARG, STATUS ) * Report any errors IF ( STATUS .NE. SAI__OK ) CALL AST_REXIT( 'BDI2_IMMAP', STATUS ) END