SUBROUTINE BDI1_GET( NARG, ARGS, OARG, STATUS ) *+ * Name: * BDI1_GET * Purpose: * Service FileItemGet requests from the BDI system for HDS files * Language: * Starlink Fortran * Invocation: * CALL BDI1_GET( NARG, ARGS, OARG, STATUS ) * Description: * Services BDI get requests for HDS files. In get mode we copy the * data from the file into a dynamic ADI object. The BDI top-level * makes the decision about the extraction type and other matters. * 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. * 8 Nov 1995 (DJA): * Handle GET operation for whole axes * 6 Dec 1995 (DJA): * Unified invention scheme * {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 'DAT_PAR' * Arguments Given: INTEGER NARG, ARGS(*) * Arguments Returned: INTEGER OARG * Status: INTEGER STATUS ! Global status * Local Variables: CHARACTER*(DAT__SZLOC) CLOC ! New component CHARACTER*20 ITEM ! Item to get INTEGER ITID ! Invented item id INTEGER NDIM, DIMS(DAT__MXDIM) ! Model object dims INTEGER WBPTR ! Write back address LOGICAL OK ! Data is valid? LOGICAL STRUC ! Object is structure *. * Check inherited global status. IF ( STATUS .NE. SAI__OK ) RETURN * Default return value OARG = ADI__NULLID * Extract the arguments CALL ADI_GET0C( ARGS(3), ITEM, STATUS ) * Locate object to be got CALL BDI1_CFIND( ARGS(1), ARGS(2), ITEM, .FALSE., .FALSE., CLOC, : NDIM, DIMS, STATUS ) * Everything ok? IF ( (STATUS .EQ. SAI__OK) .AND. (CLOC.NE.DAT__NOLOC) ) THEN * Whole axis? IF ( (ITEM(:5) .EQ. 'Axis_') .AND. : (INDEX(ITEM(6:),'_').EQ.0) ) THEN * Create ADI axis description object CALL ADI_NEW0( 'BinDSAxis', OARG, STATUS ) * Copy each component of the axis structure CALL ADI1_CCH2AC( CLOC, 'LABEL', OARG, 'Label', STATUS ) CALL ADI1_CCH2AC( CLOC, 'UNITS', OARG, 'Units', STATUS ) CALL ADI1_CCH2AL( CLOC, 'NORMALISED', OARG, : 'Normalised', STATUS ) CALL ADI1_CCH2AT( CLOC, 'DATA_ARRAY', OARG, 'Data', STATUS ) CALL ADI1_CCH2AT( CLOC, 'WIDTH', OARG, 'Width', STATUS ) CALL ADI1_CCH2AT( CLOC, 'LOWIDTH', OARG, 'LoWidth', STATUS ) CALL ADI1_CCH2AT( CLOC, 'HIWIDTH', OARG, 'HiWidth', STATUS ) ELSE * Is its data valid? CALL DAT_STRUC( CLOC, STRUC, STATUS ) IF ( STRUC ) THEN OK = .TRUE. ELSE CALL DAT_STATE( CLOC, OK, STATUS ) END IF IF ( OK ) THEN * Copy from HDS to ADI CALL ADI1_CCH2AT( CLOC, ' ', OARG, ' ', STATUS ) ELSE STATUS = SAI__ERROR CALL MSG_SETC( 'IT', ITEM ) CALL ERR_REP( 'BDI1_GET_1', 'Item ^IT data is invalid', : STATUS ) END IF * Free the HDS object CALL DAT_ANNUL( CLOC, STATUS ) END IF ELSE * Try to invent data CALL ERR_BEGIN( STATUS ) CALL BDI1_INVNT( ARGS(1), ARGS(2), ITEM, 'REAL', 'READ', : ITID, NDIM, DIMS, WBPTR, STATUS ) CALL ERR_END( STATUS ) * Invented ok? IF ( ITID .NE. ADI__NULLID ) THEN * We should store this and re-use it CALL ERR_ANNUL( STATUS ) * Return to user OARG = ITID ELSE STATUS = SAI__ERROR CALL MSG_SETC( 'IT', ITEM ) CALL ERR_REP( ' ', 'Unable to get item ^IT data from '/ : /'HDS file', STATUS ) END IF END IF * Report any errors IF ( STATUS .NE. SAI__OK ) CALL AST_REXIT( 'BDI1_GET', STATUS ) END