*+ READMAP - read a map from a big image SUBROUTINE READMAP ( STATUS ) * Description : * Read and print a rectangular area from a SuperCOSMOS mapping data set. * Invocation : * CALL READMAP ( STATUS ) * Parameters : * STATUS=INTEGER (given and returned) * global status * Method : * Open the bigimage, extract a rectangular area and list its contents. * Authors : * B.D.Kelly (ROE) * History : * 24.01.1995: original (BDK) * endhistory * Type Definitions : IMPLICIT NONE * Global constants : INCLUDE 'SAE_PAR' * Status : INTEGER STATUS * Local variables : CHARACTER*132 MAPNAME ! name of directory to hold bigimage INTEGER BD ! bigimage descriptor INTEGER*2 VALUES(200) ! values read from the bigimage INTEGER J ! loop counter INTEGER I ! loop counter CHARACTER*80 OUTBUFF ! output print buffer CHARACTER*132 ERRMSG ! buffer for error messages *- IF ( STATUS .NE. SAI__OK ) RETURN * * Initialise the input array * DO J = 1, 200 VALUES(J) = 0 ENDDO * * Get the mapname, which is the name of the directory to hold the * bigimage * CALL PAR_GET0C ( 'MAPNAME', MAPNAME, STATUS ) * * Open the bigimage, extract a 2-D rectangular map from it, and print * out the map * CALL COSMAP_OPEN ( MAPNAME, 'update', BD, STATUS ) IF ( STATUS .NE. SAI__OK ) THEN CALL COSMAP_ERRTRAN ( ERRMSG, STATUS ) CALL ERR_REP ( ' ', ERRMSG, STATUS ) ELSE CALL COSMAP_GETMAP ( BD, 5, 5, 14, 24, VALUES, STATUS ) DO J = 1, 20 WRITE ( OUTBUFF, '(10I3)' ) (VALUES((J-1)*10+I),I=1,10) CALL MSG_OUT ( ' ', OUTBUFF, STATUS ) ENDDO IF ( STATUS .NE. SAI__OK ) THEN CALL COSMAP_ERRTRAN ( ERRMSG, STATUS ) CALL ERR_REP ( ' ', ERRMSG, STATUS ) ENDIF ENDIF END