rdmtbd.f

Go to the documentation of this file.

C> @file
C> @brief Read master Table D information from local file system.
C>
C> @author J. Ator @date 2007-01-19
 
C> This subroutine reads master Table D information from two separate
C> ASCII files (one standard and one local) and then merges the
C> output into a single set of arrays.
C>
C> Each of the two ASCII files must already be individually sorted
C> in ascending order with respect to the FXY numbers.
C>
C> @param[in] LUNSTD -- integer: Fortran logical unit number for
C>                      ASCII file containing standard Table D
C>                      information
C> @param[in] LUNLTD -- integer: Fortran logical unit number for
C>                      ASCII file containing local Table D
C>                      information
C> @param[in] MXMTBD -- integer: Dimensioned size (in integers) of
C>                      merged output arrays; used by the subroutine
C>                      to ensure that it doesn't overflow these
C>                      arrays
C> @param[in] MXELEM -- integer: Maximum number of elements to be
C>                      stored per Table D sequence within merged
C>                      output arrays; used by the subroutine to
C>                      ensure that it doesn't overflow these arrays
C> @param[out] IMT   -- integer: Master table
C>                      - This value is read from both ASCII
C>                        files and must be identical between them.
C> @param[out] IMTV  -- integer: Version number of master table
C>                      - This value is read from the standard ASCII
C>                        file.
C> @param[out] IOGCE -- integer: Originating center
C>                      - This value is read from the local ASCII
C>                        file.
C> @param[out] ILTV  -- integer: Version number of local table
C>                      - This value is read from the local ASCII
C>                        file.
C> @param[out] NMTBD -- integer: Number of entries in merged output
C>                      arrays
C> @param[out] IMFXYN -- integer(*): Merged array containing bit-wise
C>                       representations of FXY numbers
C> @param[out] CMMNEM -- character*8(*): Merged array containing
C>                       mnemonics
C> @param[out] CMDSC -- character*4(*): Merged array containing
C>                      descriptor codes
C> @param[out] CMSEQ -- character*120(*): Merged array containing
C>                      sequence names
C> @param[out] NMELEM -- integer(*): Merged array containing number of
C>                       elements stored for each sequence
C> @param[out] IEFXYN -- integer(*,*): Merged array containing bit-wise
C>                       representations of element FXY numbers
C> @param[out] CEELEM -- character*120(*,*): Merged array containing
C>                       element names
C>
C> @author J. Ator @date 2007-01-19
        SUBROUTINE rdmtbd ( LUNSTD, LUNLTD, MXMTBD, MXELEM,
     .                      IMT, IMTV, IOGCE, ILTV,
     .                      NMTBD, IMFXYN, CMMNEM, CMDSC, CMSEQ,
     .                      NMELEM, IEFXYN, CEELEM )
 
        CHARACTER*200   STLINE, LTLINE
        CHARACTER*128   BORT_STR
        CHARACTER*120   CEELEM(MXMTBD,MXELEM)
        CHARACTER*6     CMATCH, ADN30
        CHARACTER*4     CMDSC(*)
        CHARACTER       CMSEQ(120,*)
        CHARACTER       CMMNEM(8,*)
 
        INTEGER         IMFXYN(*), NMELEM(*),
     .                  iefxyn(mxmtbd,mxelem)
 
C-----------------------------------------------------------------------
C-----------------------------------------------------------------------
 
C       Call WRDLEN to initialize some important information about the
C       local machine, just in case it hasn't already been called.
 
        CALL wrdlen
 
C       Read and parse the header lines of both files.
 
        CALL gettbh ( lunstd, lunltd, 'D', imt, imtv, iogce, iltv )
 
C       Read through the remainder of both files, merging the
C       contents into a unified set of master Table D arrays.
 
        nmtbd = 0
        CALL getntbe ( lunstd, isfxyn, stline, iers )
        CALL getntbe ( lunltd, ilfxyn, ltline, ierl )
        DO WHILE ( ( iers .EQ. 0 ) .OR. ( ierl .EQ. 0 ) )
          IF ( ( iers .EQ. 0 ) .AND. ( ierl .EQ. 0 ) ) THEN
            IF ( isfxyn .EQ. ilfxyn ) THEN
              cmatch = adn30( isfxyn, 6 )
              GOTO 900
            ELSE IF ( isfxyn .LT. ilfxyn ) THEN
              CALL sntbde ( lunstd, isfxyn, stline, mxmtbd, mxelem,
     .                      nmtbd, imfxyn, cmmnem, cmdsc, cmseq,
     .                      nmelem, iefxyn, ceelem )
              CALL getntbe ( lunstd, isfxyn, stline, iers )
            ELSE
              CALL sntbde ( lunltd, ilfxyn, ltline, mxmtbd, mxelem,
     .                      nmtbd, imfxyn, cmmnem, cmdsc, cmseq,
     .                      nmelem, iefxyn, ceelem )
              CALL getntbe ( lunltd, ilfxyn, ltline, ierl )
            ENDIF
          ELSE IF ( iers .EQ. 0 ) THEN
            CALL sntbde ( lunstd, isfxyn, stline, mxmtbd, mxelem,
     .                    nmtbd, imfxyn, cmmnem, cmdsc, cmseq,
     .                    nmelem, iefxyn, ceelem )
            CALL getntbe ( lunstd, isfxyn, stline, iers )
          ELSE IF ( ierl .EQ. 0 ) THEN
            CALL sntbde ( lunltd, ilfxyn, ltline, mxmtbd, mxelem,
     .                    nmtbd, imfxyn, cmmnem, cmdsc, cmseq,
     .                    nmelem, iefxyn, ceelem )
            CALL getntbe ( lunltd, ilfxyn, ltline, ierl )
          ENDIF
        ENDDO
 
        RETURN
 900    WRITE(bort_str,'("BUFRLIB: RDMTBD - STANDARD AND LOCAL'//
     . ' TABLE D FILES BOTH CONTAIN SAME FXY NUMBER: ",5A)')
     .   cmatch(1:1), '-', cmatch(2:3), '-', cmatch(4:6)
        CALL bort(bort_str)
        END