SUBROUTINE ADI2_NEWHDU( FID, HDUNAME, STATUS ) *+ * Name: * ADI2_NEWHDU * Purpose: * Create new HDU description block in FITSfile object * Language: * Starlink Fortran * Invocation: * CALL ADI2_NEWHDU( ID, HDUNAME, STATUS ) * Description: * {routine_description} * Arguments: * FID = INTEGER (given) * ADI identifier of the FITSfile object * HDUNAME = CHARACTER*(*) (given) * HDU name * 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: * 14 Jun 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 FID CHARACTER*(*) HDUNAME * Status: INTEGER STATUS ! Global status * External references: EXTERNAL CHR_LEN INTEGER CHR_LEN * Local Variables: CHARACTER*20 LHDU CHARACTER*4 STR INTEGER HCID ! HDU container INTEGER HID ! New HDU INTEGER HLEN INTEGER I ! Loop over LHDU INTEGER NDIG ! Digits used in STR INTEGER NHDU ! # HDUs in file *. * Check inherited global status. IF ( STATUS .NE. SAI__OK ) RETURN * Make local copy of name, and translate embedded spaces to underscores LHDU = HDUNAME HLEN = CHR_LEN(LHDU) DO I = 1, HLEN IF ( LHDU(I:I) .EQ. ' ' ) LHDU(I:I) = '_' END DO * Locate file's HDU container CALL ADI_FIND( FID, 'Hdus', HCID, STATUS ) * Create named component & locate it CALL ADI_CNEW0( HCID, LHDU(:HLEN), 'FITShdu', STATUS ) CALL ADI_FIND( HCID, LHDU(:HLEN), HID, STATUS ) * Set the HDU counter CALL ADI_CGET0I( FID, 'Nhdu', NHDU, STATUS ) NHDU = NHDU + 1 CALL ADI_CPUT0I( FID, 'Nhdu', NHDU, STATUS ) * Store the HDUs self reference number CALL ADI_CPUT0I( HID, 'Ihdu', NHDU, STATUS ) * Write property name so that ADI can do number -> HDU name mapping CALL CHR_ITOC( NHDU, STR, NDIG ) CALL ADI_CPUT0C( FID, '.HDU_'//STR(:NDIG), LHDU(:HLEN), STATUS ) * Write true name CALL ADI_CPUT0C( HID, 'Name', HDUNAME, STATUS ) * Relase new HDU CALL ADI_ERASE( HID, STATUS ) * Release container CALL ADI_ERASE( HCID, STATUS ) * Report any errors IF ( STATUS .NE. SAI__OK ) CALL AST_REXIT( 'ADI2_NEWHDU', STATUS ) END