SUBROUTINE BDI_GETSHP( ID, MXNDIM, DIMS, NDIM, STATUS ) *+ * Name: * BDI_GETSHP * Purpose: * Get dimensions of data model object supported by BDI * Language: * Starlink Fortran * Invocation: * CALL BDI_GETSHP( ID, MXNDIM, DIMS, NDIM, STATUS ) * Description: * Extracts the dimensions of a data model object for subsequent * manipulation by BDI. * Arguments: * ID = INTEGER (given) * ADI identifier of BinDS, Array or Scalar object, or derivatives * thereof * MXNDIM = INTEGER (given) * Maximum number of dimensions to export * DIMS[NDIM] = INTEGER (returned) * Sizes of each of the NDIM dimensions of the object * NDIM = INTEGER (returned) * Dimensionality of object, zero for scalar * 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:public * 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' INCLUDE 'AST_PKG' * Arguments Given: INTEGER ID, MXNDIM * Arguments Returned: INTEGER NDIM, DIMS(*) * Status: INTEGER STATUS ! Global status * External References: LOGICAL AST_QPKGI EXTERNAL AST_QPKGI * Local Variables: INTEGER IDIM ! Loop over dimensions LOGICAL SCALAR ! Input is scalar? *. * Check inherited global status. IF ( STATUS .NE. SAI__OK ) RETURN * Check initialised IF ( .NOT. AST_QPKGI( BDI__PKG ) ) CALL BDI0_INIT( STATUS ) * Object is derived from Scalar? CALL ADI_DERVD( ID, 'Scalar', SCALAR, STATUS ) IF ( SCALAR ) THEN * Set dimensions NDIM = 0 DO IDIM = 1, MXNDIM DIMS(IDIM) = 0 END DO ELSE * How big is dimensions array CALL ADI_CSIZE( ID, 'SHAPE', NDIM, STATUS ) IF ( NDIM .GT. MXNDIM ) THEN CALL BDI0_DESCID( ID, 'ID', STATUS ) CALL MSG_SETI( 'ND', NDIM ) CALL MSG_SETI( 'NA', MXNDIM ) STATUS = SAI__ERROR CALL ERR_REP( 'BDI_GETSHP_0', 'Input object ^ID has '/ : /'^ND dimensions, whereas the calling '/ : /'program can only cope with ^NA', STATUS ) ELSE * Read dimensions CALL ADI_CGET1I( ID, 'SHAPE', MXNDIM, DIMS, NDIM, STATUS ) IF ( STATUS .NE. SAI__OK ) THEN STATUS = SAI__ERROR CALL ERR_REP( 'BDI_GETSHP_1', 'Unable to read data object'/ : /' dimensions', STATUS ) END IF END IF END IF * Report any errors IF ( STATUS .NE. SAI__OK ) CALL AST_REXIT( 'BDI_GETSHP', STATUS ) END