dxtable.F90

Go to the documentation of this file.

 
 
subroutine readdx(lunit,lun,lundx)
 
  use modv_vars, only: iprt
 
  implicit none
 
  integer, intent(in) :: lunit, lun, lundx
  integer lud, ildx, imdx
 
  character*128 errstr
 
  ! Get the status of unit lundx
 
  call status(lundx,lud,ildx,imdx)
 
  ! Read a dictionary table from the indicated source
 
  if (lunit==lundx) then
    ! Source is input BUFR file in lunit
    if(iprt>=2) then
      call errwrt('++++++++++++BUFR ARCHIVE LIBRARY+++++++++++++++')
      write ( unit=errstr, fmt='(A,A,I3,A)' ) 'BUFRLIB: READDX - READING BUFR DICTIONARY TABLE FROM ', &
        'INPUT BUFR FILE IN UNIT ', lundx, ' INTO INTERNAL ARRAYS'
      call errwrt(errstr)
      call errwrt('++++++++++++BUFR ARCHIVE LIBRARY+++++++++++++++')
      call errwrt(' ')
    endif
    rewind lunit
    call rdbfdx(lunit,lun)
  elseif(ildx==-1) then
    ! Source is input BUFR file in lundx; BUFR file in lunit may be input or output
    if(iprt>=2) then
      call errwrt('++++++++++++BUFR ARCHIVE LIBRARY+++++++++++++++')
      write ( unit=errstr, fmt='(A,A,I3,A,A,I3)' ) 'BUFRLIB: READDX - COPYING BUFR DCTY TBL FROM INTERNAL ', &
        'ARRAYS ASSOC. W/ INPUT UNIT ', lundx, ' TO THOSE ASSOC. W/ UNIT ', lunit
      call errwrt(errstr)
      call errwrt('++++++++++++BUFR ARCHIVE LIBRARY+++++++++++++++')
      call errwrt(' ')
    endif
    call cpbfdx(lud,lun)
    call makestab
  elseif(ildx==1) then
    ! Source is output BUFR file in lundx; BUFR file in lunit may be input or output
    if(iprt>=2) then
      call errwrt('++++++++++++BUFR ARCHIVE LIBRARY+++++++++++++++')
      write ( unit=errstr, fmt='(A,A,I3,A,A,I3)' ) 'BUFRLIB: READDX - COPYING BUFR DCTY TBL FROM INTERNAL ', &
        'ARRAYS ASSOC. W/ OUTPUT UNIT ', lundx, ' TO THOSE ASSOC. W/ UNIT ', lunit
      call errwrt(errstr)
      call errwrt('++++++++++++BUFR ARCHIVE LIBRARY+++++++++++++++')
      call errwrt(' ')
    endif
    call cpbfdx(lud,lun)
    call makestab
  elseif(ildx==0) then
    ! Source is user-supplied character table in lundx; BUFR file in lunit may be input or output
    if(iprt>=2) then
      call errwrt('++++++++++++BUFR ARCHIVE LIBRARY+++++++++++++++')
      write ( unit=errstr, fmt='(A,A,I3,A)' ) 'BUFRLIB: READDX - READING BUFR DICTIONARY TABLE FROM ', &
        'USER-SUPPLIED TEXT FILE IN UNIT ', lundx, ' INTO INTERNAL ARRAYS'
      call errwrt(errstr)
      call errwrt('++++++++++++BUFR ARCHIVE LIBRARY+++++++++++++++')
      call errwrt(' ')
    endif
    rewind lundx
    call rdusdx(lundx,lun)
  else
    call bort('BUFRLIB: READDX - CANNOT DETERMINE SOURCE OF INPUT DICTIONARY TABLE')
  endif
 
  return
end subroutine readdx
 
subroutine rdbfdx(lunit,lun)
 
  use bufrlib
 
  use modv_vars, only: iprt
 
  use moda_mgwa
 
  implicit none
 
  integer, intent(in) :: lunit, lun
  integer ict, ier, idxmsg, iupbs3
 
  character*128 errstr
 
  logical done
 
  call dxinit(lun,0)
 
  ict = 0
  done = .false.
 
  ! Read a complete dictionary table from lunit, as a set of one or more DX dictionary messages.
 
  do while ( .not. done )
    call rdmsgw ( lunit, mgwa, ier )
    if ( ier == -1 ) then
      ! Don't abort for an end-of-file condition, since it may be possible for a file to end with dictionary messages.
      ! Instead, backspace the file pointer and let the calling routine diagnose the end-of-file condition and deal with
      ! it as it sees fit.
      call backbufr_c(lun)
      done = .true.
    else if ( ier == -2 ) then
      call bort('BUFRLIB: RDBFDX - ERROR READING A BUFR DICTIONARY MESSAGE')
    else if ( idxmsg(mgwa) /= 1 ) then
      ! This is a non-DX dictionary message.  Assume we've reached the end of the dictionary table, and backspace lunit
      ! so that the next read (e.g. in the calling routine) will get this same message.
      call backbufr_c(lun)
      done = .true.
    else if ( iupbs3(mgwa,'NSUB') == 0 ) then
      ! This is a DX dictionary message, but it doesn't contain any actual dictionary information.  Assume we've reached
      ! the end of the dictionary table.
      done = .true.
    else
      ! Store this message into module @ref moda_tababd.
      ict = ict + 1
      call stbfdx(lun,mgwa)
    endif
  enddo
 
  if ( iprt >= 2 ) then
    call errwrt('+++++++++++++++++++++++++++++++++++++++++++++++++')
    write ( unit=errstr, fmt='(A,I3,A)' ) 'BUFRLIB: RDBFDX - STORED NEW DX TABLE CONSISTING OF (', ict, ') MESSAGES;'
    call errwrt(errstr)
    errstr = 'WILL APPLY THIS TABLE TO ALL SUBSEQUENT DATA IN FILE UNTIL NEXT DX TABLE IS FOUND'
    call errwrt(errstr)
    call errwrt('+++++++++++++++++++++++++++++++++++++++++++++++++')
    call errwrt(' ')
  endif
 
  call makestab
 
  return
end subroutine rdbfdx
 
subroutine rdusdx(lundx,lun)
 
  use moda_tababd
 
  implicit none
 
  integer, intent(in) :: lundx, lun
  integer ios, iret, n, numbck, nemock, igetntbi
 
  character*128 bort_str1
  character*156 bort_str2
  character*80 card
  character*8 nemo
  character*6 numb, nmb2
 
  ! Initialize the dictionary table control word partition arrays with apriori Table B and D entries
 
  call dxinit(lun,1)
  rewind lundx
 
  ! Read user cards until there are no more
 
  do while (.true.)
 
    read(lundx, '(A80)', iostat = ios) card
    if (ios/=0) then
      call makestab
      return
    endif
 
    if(card(1: 1)==       '*') cycle  ! comment line
    if(card(3:10)=='--------') cycle  ! separation line
    if(card(3:10)=='        ') cycle  ! blank line
    if(card(3:10)=='MNEMONIC') cycle  ! header line
    if(card(3:10)=='TABLE  D') cycle  ! header line
    if(card(3:10)=='TABLE  B') cycle  ! header line
 
    if(card(12:12)=='|' .and. card(21:21)=='|') then
 
      ! Parse a descriptor definition card
      nemo = card(3:10)  ! nemo is the (up to) 8-character mnemonic
      iret=nemock(nemo)
      if(iret==-2) then
        write(bort_str1,'("BUFRLIB: RDUSDX - CARD READ IN IS: ",A)') card
        write(bort_str2,'(18X,"MNEMONIC ",A," IN USER DICTIONARY HAS INVALID CHARACTERS")') nemo
        call bort2(bort_str1,bort_str2)
      endif
      numb = card(14:19) ! numb is the 6-character FXY value corresponding to nemo
      nmb2 = numb
      if(nmb2(1:1)=='A') nmb2(1:1) = '3'
      iret=numbck(nmb2)
      if(iret==-1) then
        write(bort_str1,'("BUFRLIB: RDUSDX - CARD READ IN IS: ",A)') card
        write(bort_str2,'(18X,"DESCRIPTOR NUMBER ",A," IN USER '// &
          'DICTIONARY HAS AN INVALID FIRST CHARACTER (F VALUE) - MUST BE A, 0 OR 3")') numb
        call bort2(bort_str1,bort_str2)
      endif
      if(iret==-2) then
        write(bort_str1,'("BUFRLIB: RDUSDX - CARD READ IN IS: ",A)') card
        write(bort_str2,'(18X,"DESCRIPTOR NUMBER ",A," IN USER '// &
          'DICTIONARY HAS NON-NUMERIC VALUES IN CHARACTERS 2-6 (X AND Y VALUES)")') numb
        call bort2(bort_str1,bort_str2)
      endif
      if(iret==-3) then
        write(bort_str1,'("BUFRLIB: RDUSDX - CARD READ IN IS: ",A)') card
        write(bort_str2,'(18X,"DESCRIPTOR NUMBER ",A," IN USER '// &
          'DICTIONARY HAS INVALID NUMBER IN CHARACTERS 2-3 (X VALUE) - MUST BE BETWEEN 00 AND 63")') numb
        call bort2(bort_str1,bort_str2)
      endif
      if(iret==-4) then
        write(bort_str1,'("BUFRLIB: RDUSDX - CARD READ IN IS: ",A)') card
        write(bort_str2,'(18X,"DESCRIPTOR NUMBER ",A," IN USER '// &
          'DICTIONARY HAS INVALID NUMBER IN CHARACTERS 4-6 (Y VALUE) - MUST BE BETWEEN 000 AND 255")') numb
        call bort2(bort_str1,bort_str2)
      endif
 
      if(numb(1:1)=='A') then
        ! Table A descriptor found
        n = igetntbi( lun, 'A' )
        call stntbia ( n, lun, numb, nemo, card(23:) )
        if ( idna(n,lun,1) == 11 ) then
          write(bort_str1,'("BUFRLIB: RDUSDX - CARD READ IN IS: ",A)') card
          write(bort_str2,'(18X,"USER-DEFINED MESSAGE TYPE ""011"" IS RESERVED FOR DICTIONARY MESSAGES")')
          call bort2(bort_str1,bort_str2)
        endif
        ! Replace "A" with "3" so Table D descriptor will be found in card as well (see below).
        numb(1:1) = '3'
      endif
 
      if(numb(1:1)=='0') then
        ! Table B descriptor found
        call stntbi ( igetntbi(lun,'B'), lun, numb, nemo, card(23:) )
        cycle
      endif
 
      if(numb(1:1)=='3') then
        ! Table D descriptor found
        call stntbi ( igetntbi(lun,'D'), lun, numb, nemo, card(23:) )
        cycle
      endif
 
      write(bort_str1,'("BUFRLIB: RDUSDX - CARD READ IN IS: ",A)') card
      write(bort_str2,'(18X,"DESCRIPTOR NUMBER ",A," IN USER '// &
        'DICTIONARY HAS AN INVALID FIRST CHARACTER (F VALUE) - MUST BE A, 0 OR 3")') numb
      call bort2(bort_str1,bort_str2)
 
    endif
 
    if(card(12:12)=='|' .and. card(19:19)/='|') then
      ! Parse a sequence definition card
      call seqsdx(card,lun)
      cycle
    endif
 
    if(card(12:12)=='|' .and. card(19:19)=='|') then
      ! Parse an element definition card
      call elemdx(card,lun)
      cycle
    endif
 
    ! Can't figure out what kind of card it is
    write(bort_str1,'("BUFRLIB: RDUSDX - CARD READ IN IS: ",A)') card
    write(bort_str2,'(18X,"THIS CARD HAS A BAD FORMAT - IT IS NOT RECOGNIZED BY THIS SUBROUTINE")')
    call bort2(bort_str1,bort_str2)
  enddo
 
end subroutine rdusdx
 
subroutine seqsdx(card,lun)
 
  use modv_vars, only: reps, idnr
 
  implicit none
 
  integer, intent(in) :: lun
  integer ntag, idn, jdn, iseq, irep, i, j, n, itab, iret, ier, numr, nemock
  integer, parameter :: maxtgs = 250, maxtag = 13
 
  character*128 bort_str1, bort_str2
  character*80 seqs
  character*80, intent(in) :: card
  character*(maxtag) atag, tags(maxtgs)
  character*8 nemo, nema, nemb
  character*6 adn30, clemon
  character tab
 
  ! Find the sequence tag in Table D and parse the sequence string
 
  nemo = card( 3:10)
  seqs = card(14:78)
 
  ! Note that an entry for this mnemonic should already exist within the internal BUFR Table D array tabd(*,LUN); this entry
  ! should have been created by subroutine rdusdx() when the mnemonic and its associated FXY value and description were
  ! initially defined within a card read from the "Descriptor Definition" section at the top of the user-supplied DX BUFR
  ! table in character format.  Now, we need to retrieve the positional index for that  entry within tabd(*,lun) so that we
  ! can access the entry and then add the decoded sequence information to it.
 
  call nemtab(lun,nemo,idn,tab,iseq)
  if(tab/='D') then
    write(bort_str1,'("BUFRLIB: SEQSDX - CARD READ IN IS: ",A)') card
    write(bort_str2,'(18X,"MNEMONIC ",A," IS NOT A TABLE D ENTRY (UNDEFINED, TAB=",A,")")') nemo,tab
    call bort2(bort_str1,bort_str2)
  endif
  call parstr(seqs,tags,maxtgs,ntag,' ',.true.)
  if(ntag==0) then
    write(bort_str1,'("BUFRLIB: SEQSDX - CARD READ IN IS: ",A)') card
    write(bort_str2,'(18X,"TABLE D SEQUENCE (PARENT) MNEMONIC ",A," DOES NOT CONTAIN ANY CHILD MNEMONICS")') nemo
    call bort2(bort_str1,bort_str2)
  endif
 
  do n=1,ntag
    atag = tags(n)
    irep = 0
 
    ! Check for a replicator
 
    outer: do i=1,5
      if(atag(1:1)==reps(i)) then
        ! Note that reps(*), which contains all of the symbols used to denote all of the various replication schemes that
        ! are possible within a user-supplied BUFR dictionary table in character format, was previously defined within
        ! subroutine bfrini().
        do j=2,maxtag
          if(atag(j:j)==reps(i+5)) then
            ! Note that subroutine strnum() will return numr = 0 if the string passed to it contains all blanks
            ! (as *should* be the case whenever i = 2 '(' ')', 3 '{' '}', 4 '[' ']', or 5 '<' '>').
            ! However, when i = 1 '"' '"', then subroutine strnum() will return numr = (the number of replications for
            ! the mnemonic using F=1 "regular" (i.e. non-delayed) replication).
            call strnum(atag(j+1:maxtag),numr,ier)
            if(i==1 .and. numr<=0) then
              write(bort_str1,'("BUFRLIB: SEQSDX - CARD READ IN IS: ",A)') card
              write(bort_str2,'(9X,"TBL D MNEM. ",A," CONTAINS REG. REPL. '// &
                'CHILD MNEM. ",A," W/ INVALID # OF REPLICATIONS (",I3,") AFTER 2ND QUOTE")') nemo,tags(n),numr
              call bort2(bort_str1,bort_str2)
            endif
            if(i==1 .and. numr>255) then
              write(bort_str1,'("BUFRLIB: SEQSDX - CARD READ IN IS: ",A)') card
              write(bort_str2,'(18X,"TBL D MNEM. ",A," CONTAINS REG. REPL. '// &
                'CHILD MNEM. ",A," W/ # OF REPLICATIONS (",I3,") > LIMIT OF 255")') nemo,tags(n),numr
              call bort2(bort_str1,bort_str2)
            endif
            if(i/=1 .and. numr/=0) then
              write(bort_str1,'("BUFRLIB: SEQSDX - CARD READ IN IS: ",A)') card
              write(bort_str2,'(18X,"TBL D MNEM. ",A," CONTAINS DELAYED REPL. '// &
                'CHILD MNEM. ",A," W/ # OF REPL. (",I3,") SPECIFIED - A NO-NO")') nemo,tags(n),numr
              call bort2(bort_str1,bort_str2)
            endif
            atag = atag(2:j-1)
            irep = i
            exit outer
          endif
        enddo
        write(bort_str1,'("BUFRLIB: SEQSDX - CARD READ IN IS: ",A)') card
        write(bort_str2,'(18X,"TABLE D SEQUENCE (PARENT) MNEMONIC ",A,'// &
          '" CONTAINS A BADLY FORMED CHILD MNEMONIC ",A)') nemo,tags(n)
        call bort2(bort_str1,bort_str2)
      endif
    enddo outer
 
    ! Check for a valid tag
 
    iret=nemock(atag)
    if(iret==-1) then
      write(bort_str1,'("BUFRLIB: SEQSDX - CARD READ IN IS: ",A)') card
      write(bort_str2,'(18X,"TABLE D (PARENT) MNEMONIC ",A," CONTAINS'// &
        ' A CHILD MNEMONIC ",A," NOT BETWEEN 1 & 8 CHARACTERS")') nemo,tags(n)
      call bort2(bort_str1,bort_str2)
    endif
    if(iret==-2) then
      write(bort_str1,'("BUFRLIB: SEQSDX - CARD READ IN IS: ",A)') card
      write(bort_str2,'(18X,"TABLE D (PARENT) MNEMONIC ",A," CONTAINS'// &
        ' A CHILD MNEMONIC ",A," WITH INVALID CHARACTERS")') nemo,tags(n)
      call bort2(bort_str1,bort_str2)
    endif
    call nemtab(lun,atag,idn,tab,iret)
    if(iret>0) then
      ! Note that the next code line checks that we are not trying to replicate a Table B mnemonic (which is currently not
      ! allowed).  The logic works because, for replicated mnemonics, irep = i = (the index within reps(*) of the symbol
      ! associated with the type of replication in question (e.g. "{, "<", etc.))
      if(tab=='B' .and. irep/=0) then
        write(bort_str1,'("BUFRLIB: SEQSDX - CARD READ IN IS: ",A)') card
        write(bort_str2,'(18X,"TABLE D (PARENT) MNEMONIC ",A," CONTAINS'// &
          ' A REPLICATED CHILD TABLE B MNEMONIC ",A," - A NO-NO")') nemo,tags(n)
        call bort2(bort_str1,bort_str2)
      endif
      if(atag(1:1)=='.') then
        ! This mnemonic is a "following value" mnemonic (i.e. it relates to the mnemonic that immediately follows it within
        ! the user-supplied character-format BUFR dictionary table sequence), so confirm that it contains, as a substring,
        ! this mnemonic that immediately follows it.
        if(n==ntag) then
          write(bort_str1,'("BUFRLIB: SEQSDX - CARD READ IN IS: ",A)') card
          write(bort_str2,'(18X,"TBL D (PARENT) MNEM. ",A," CONTAINS A '// &
            '''FOLLOWING VALUE'' MNEMONIC WHICH IS LAST IN THE STRING")') nemo
          call bort2(bort_str1,bort_str2)
        endif
        nemb = tags(n+1)(1:8)
        call numtab(lun,idn,nema,tab,itab)
        call nemtab(lun,nemb,jdn,tab,iret)
        call rsvfvm(nema,nemb)
        if(nema/=atag) then
          write(bort_str1,'("BUFRLIB: SEQSDX - CARD READ IN IS: ",A)') card
          write(bort_str2,'(18X,"TBL D (PARENT) MNEM. ",A," CONTAINS AN '// &
            'INVALID ''FOLLOWING VALUE'' MNEMONIC ",A,"(SHOULD BE ",A,")")') nemo,tags(n),nema
          call bort2(bort_str1,bort_str2)
        endif
        if(tab/='B') then
          write(bort_str1,'("BUFRLIB: SEQSDX - CARD READ IN IS: ",A)') card
          write(bort_str2,'(18X,"TBL D (PARENT) MNEM. ",A,", THE MNEM. ",'// &
            'A," FOLLOWING A ''FOLLOWING VALUE'' MNEM. IS NOT A TBL B ENTRY")') nemo,nemb
          call bort2(bort_str1,bort_str2)
        endif
      endif
    else
      write(bort_str1,'("BUFRLIB: SEQSDX - CARD READ IN IS: ",A)') card
      write(bort_str2,'(18X,"TABLE D SEQUENCE (PARENT) MNEMONIC ",A,'// &
        '" CONTAINS A CHILD MNEMONIC ",A," NOT FOUND IN ANY TABLE")') nemo,tags(n)
      call bort2(bort_str1,bort_str2)
    endif
 
    ! Write the descriptor string into the tabd array, but first look for a replication descriptor
    if(irep>0) call pktdd(iseq,lun,idnr(irep)+numr,iret)
    if(iret<0) then
      clemon = adn30(idnr(irep)+numr,6)
      write(bort_str1,'("BUFRLIB: SEQSDX - CARD READ IN IS: ",A)') card
      write(bort_str2,'(9X,"TBL D (PARENT) MNEM. ",A," - BAD RETURN '// &
        'FROM PKTDD TRYING TO STORE REPL. DESC. ",A,", SEE PREV. WARNING MSG")') nemo,clemon
      call bort2(bort_str1,bort_str2)
    endif
    call pktdd(iseq,lun,idn,iret)
    if(iret<0) then
      write(bort_str1,'("BUFRLIB: SEQSDX - CARD READ IN IS: ",A)') card
      write(bort_str2,'(9X,"TBL D (PARENT) MNEM. ",A," - BAD RETURN '// &
        'FROM PKTDD TRYING TO STORE CHILD MNEM. ",A,", SEE PREV. WARNING MSG")') nemo,tags(n)
      call bort2(bort_str1,bort_str2)
    endif
 
  enddo
 
  return
end subroutine seqsdx
 
subroutine elemdx(card,lun)
 
  use moda_tababd
 
  implicit none
 
  integer, intent(in) :: lun
  integer idsn, iele, iret
 
  character*128 bort_str1, bort_str2
  character*80, intent(in) :: card
  character*24 unit
  character*11 refr, refr_orig
  character*8 nemo
  character*4 scal, scal_orig
  character*3 bitw, bitw_orig
  character sign, tab
 
  ! Capture the various elements characteristics
 
  nemo = card( 3:10)
  scal = card(14:17)
  refr = card(21:31)
  bitw = card(35:37)
  unit = card(41:64)
  ! Make sure the units are all capitalized
  call capit(unit)
 
  ! Find the element tag in Table B.  Note that an entry for this mnemonic should already exist within the internal
  ! BUFR Table B array tabb(*,lun).  We now need to retrieve the positional index for that entry within tabb(*,lun)
  ! so that we can access the entry and then add the scale factor, reference value, bit width, and units to it.
 
  call nemtab(lun,nemo,idsn,tab,iele)
  if(tab/='B') then
    write(bort_str1,'("BUFRLIB: ELEMDX - CARD READ IN IS: ",A)') card
    write(bort_str2,'(18X,"MNEMONIC ",A," IS NOT A TABLE B ENTRY (UNDEFINED, TAB=",A,")")') nemo,tab
    call bort2(bort_str1,bort_str2)
  endif
 
  ! Left justify and store characteristics
 
  unit = adjustl(unit)
  if(unit==' ') then
    write(bort_str1,'("BUFRLIB: ELEMDX - CARD READ IN IS: ",A)') card
    write(bort_str2,'(18X,"UNITS FIELD IS EMPTY")')
    call bort2(bort_str1,bort_str2)
  endif
  tabb(iele,lun)(71:94) = unit
 
  scal_orig=scal
  call jstnum(scal,sign,iret)
  if(iret/=0) then
    write(bort_str1,'("BUFRLIB: ELEMDX - CARD READ IN IS: ",A)') card
    write(bort_str2,'(18X,"PARSED SCALE VALUE (=",A,") IS NOT NUMERIC")') scal_orig
    call bort2(bort_str1,bort_str2)
  endif
  tabb(iele,lun)(95:95) = sign
  tabb(iele,lun)(96:98) = scal(1:3)
 
  refr_orig=refr
  call jstnum(refr,sign,iret)
  if(iret/=0) then
    write(bort_str1,'("BUFRLIB: ELEMDX - CARD READ IN IS: ",A)') card
    write(bort_str2,'(18X,"PARSED REFERENCE VALUE (=",A,") IS NOT NUMERIC")') refr_orig
    call bort2(bort_str1,bort_str2)
  endif
  tabb(iele,lun)( 99: 99) = sign
  tabb(iele,lun)(100:109) = refr(1:10)
 
  bitw_orig=bitw
  call jstnum(bitw,sign,iret)
  if(iret/=0 .or. sign=='-') then
    write(bort_str1,'("BUFRLIB: ELEMDX - CARD READ IN IS: ",A)') card
    write(bort_str2,'(18X,"PARSED BIT WIDTH VALUE (=",A,") IS NOT NUMERIC")') bitw_orig
    call bort2(bort_str1,bort_str2)
  endif
  tabb(iele,lun)(110:112) = bitw
 
  return
end subroutine elemdx
 
subroutine dxinit(lun,ioi)
 
  use modv_vars, only: idnr, fxy_fbit, fxy_sbyct, fxy_drf16, fxy_drf8, fxy_drf1
 
  use moda_tababd
 
  implicit none
 
  integer, intent(in) :: lun, ioi
  integer ninib, ninid, n, i, iret, ifxy
 
  character*8 inib(6,5),inid(5)
  character*6 adn30
 
  data inib  / '------','BYTCNT  ','BYTES  ','+0','+0','16', &
               '------','BITPAD  ','NONE   ','+0','+0','1 ', &
               fxy_drf1,'DRF1BIT ','NUMERIC','+0','+0','1 ', &
               fxy_drf8,'DRF8BIT ','NUMERIC','+0','+0','8 ', &
              fxy_drf16,'DRF16BIT','NUMERIC','+0','+0','16'/
  data ninib /5/
 
  data inid  /'        ', &
              'DRP16BIT', &
              'DRP8BIT ', &
              'DRPSTAK ', &
              'DRP1BIT '/
  data ninid /5/
 
  ! Clear out a table partition
 
  ntba(lun) = 0
  do i=1,ntba(0)
    taba(i,lun) = ' '
    mtab(i,lun) = 0
  enddo
 
  ntbb(lun) = 0
  tabb(1:ntbb(0),lun) = ' '
 
  ntbd(lun) = 0
  do i=1,ntbd(0)
    tabd(i,lun) = ' '
    call pktdd(i,lun,0,iret)
  enddo
 
  if(ioi==0) return
 
  ! Initialize table with apriori Table B and D entries
 
  inib(1,1) = fxy_sbyct
  inib(1,2) = fxy_fbit
 
  do i=1,ninib
    ntbb(lun) = ntbb(lun)+1
    idnb(i,lun) = ifxy(inib(1,i))
    tabb(i,lun)(  1:  6) = inib(1,i)(1:6)
    tabb(i,lun)(  7: 70) = inib(2,i)
    tabb(i,lun)( 71: 94) = inib(3,i)
    tabb(i,lun)( 95: 98) = inib(4,i)(1:4)
    tabb(i,lun)( 99:109) = inib(5,i)
    tabb(i,lun)(110:112) = inib(6,i)(1:3)
  enddo
 
  do i=2,ninid
    n = ntbd(lun)+1
    idnd(n,lun) = idnr(i)
    tabd(n,lun)(1: 6) = adn30(idnr(i),6)
    tabd(n,lun)(7:70) = inid(i)
    call pktdd(n,lun,idnr(1),iret)
    call pktdd(n,lun,idnr(i+5),iret)
    ntbd(lun) = n
  enddo
 
  return
end subroutine dxinit
 
subroutine dxmini(mbay,mbyt,mb4,mba,mbb,mbd)
 
  use modv_vars, only: mxmsgld4, mtv, nby0, nby1, nby2, nby5, bmostr, idxv
 
  implicit none
 
  integer, intent(out) :: mbay(*), mbyt, mb4, mba, mbb, mbd
  integer nxstr, ldxa, ldxb, ldxd, ld30, mtyp, msbt, mbit, ih, id, im, iy, i, nsub, idxs, ldxs, &
    len3, nby4, iupm
 
  character*128 bort_str
  character*56 dxstr
 
  common /dxtab/ nxstr(10), ldxa(10), ldxb(10), ldxd(10), ld30(10), dxstr(10)
 
  msbt = idxv
 
  ! Initialize the message
 
  mbit = 0
  mbay(1:mxmsgld4) = 0
 
  ! For DX table messages, the Section 1 date is simply zeroed out.  Note that there is logic in function idxmsg()
  ! which relies on this.
  ih = 0
  id = 0
  im = 0
  iy = 0
 
  mtyp = 11  ! DX table messages are always type 11, per WMO BUFR Table A
  nsub = 1
 
  idxs = idxv+1
  ldxs = nxstr(idxs)
 
  len3 = 7 + nxstr(idxs) + 1
  nby4 = 7
  mbyt = nby0+nby1+nby2+len3+nby4+nby5
 
  if(mod(len3,2)/=0) call bort ('BUFRLIB: DXMINI - LENGTH OF SECTION 3 IS NOT A MULTIPLE OF 2')
 
  ! Section 0
 
  call pkc(bmostr ,  4 , mbay,mbit)
  call pkb(  mbyt , 24 , mbay,mbit)
  call pkb(     3 ,  8 , mbay,mbit)
 
  ! Section 1
 
  call pkb(  nby1 , 24 , mbay,mbit)
  call pkb(     0 ,  8 , mbay,mbit)
  call pkb(     3 ,  8 , mbay,mbit)
  call pkb(     7 ,  8 , mbay,mbit)
  call pkb(     0 ,  8 , mbay,mbit)
  call pkb(     0 ,  8 , mbay,mbit)
  call pkb(  mtyp ,  8 , mbay,mbit)
  call pkb(  msbt ,  8 , mbay,mbit)
  call pkb(   mtv ,  8 , mbay,mbit)
  call pkb(  idxv ,  8 , mbay,mbit)
  call pkb(    iy ,  8 , mbay,mbit)
  call pkb(    im ,  8 , mbay,mbit)
  call pkb(    id ,  8 , mbay,mbit)
  call pkb(    ih ,  8 , mbay,mbit)
  call pkb(     0 ,  8 , mbay,mbit)
  call pkb(     0 ,  8 , mbay,mbit)
 
  ! Section 3
 
  call pkb(  len3 , 24 , mbay,mbit)
  call pkb(     0 ,  8 , mbay,mbit)
  call pkb(     1 , 16 , mbay,mbit)
  call pkb(  2**7 ,  8 , mbay,mbit)
  do i=1,ldxs
    call pkb(iupm(dxstr(idxs)(i:i),8),8,mbay,mbit)
  enddo
  call pkb(     0 ,  8 , mbay,mbit)
 
  ! Section 4
 
  mb4 = mbit/8+1
  call pkb(  nby4 , 24 , mbay,mbit)
  call pkb(     0 ,  8 , mbay,mbit)
  mba = mbit/8+1
  call pkb(     0 ,  8 , mbay,mbit)
  mbb = mbit/8+1
  call pkb(     0 ,  8 , mbay,mbit)
  mbd = mbit/8+1
  call pkb(     0 ,  8 , mbay,mbit)
 
  if(mbit/8+nby5/=mbyt) then
    write(bort_str,'("BUFRLIB: DXMINI - NUMBER OF BYTES STORED FOR '// &
      'A MESSAGE (",I6,") IS NOT THE SAME AS FIRST CALCULATED, MBYT (",I6)') mbit/8+nby5,mbyt
    call bort(bort_str)
  endif
 
  return
end subroutine dxmini
 
subroutine writdx(lunit,lun,lundx)
 
  implicit none
 
  integer, intent(in) :: lunit, lun, lundx
 
  character*128 bort_str
 
  ! The table must be coming from an input file
 
  if(lunit==lundx) then
    write(bort_str,'("BUFRLIB: WRITDX - FILES CONTAINING BUFR DATA '// &
      'AND DICTIONARY TABLE CANNOT BE THE SAME (HERE BOTH SHARE FORTRAN UNIT NUMBER ",I3,")")') lunit
    call bort(bort_str)
  endif
 
  ! Must first call readdx() to generate internal dictionary table arrays, before calling wrdxtb()
 
  call readdx(lunit,lun,lundx)
  call wrdxtb(lunit,lunit)
 
  return
end subroutine writdx
 
recursive subroutine wrdxtb(lundx,lunot)
 
  use modv_vars, only: im8b, idxv
 
  use moda_tababd
  use moda_mgwa
  use moda_bitbuf, only: maxbyt
 
  implicit none
 
  integer, intent(in) :: lundx, lunot
  integer nxstr, ldxa, ldxb, ldxd, ld30, my_lundx, my_lunot, ldx, lot, il, im, lda, ldb, ldd, l30, nseq, &
    mbit, mbyt, mby4, mbya, mbyb, mbyd, i, j, jj, idn, lend, len0, len1, len2, l3, l4, l5, iupb, iupm
 
  common /dxtab/ nxstr(10), ldxa(10), ldxb(10), ldxd(10), ld30(10), dxstr(10)
 
  character*56 dxstr
  character*6 adn30
 
  logical msgfull
 
  ! Check for I8 integers
 
  if(im8b) then
    im8b=.false.
 
    call x84(lundx,my_lundx,1)
    call x84(lunot,my_lunot,1)
    call wrdxtb(my_lundx,my_lunot)
 
    im8b=.true.
    return
  endif
 
  ! Check file statuses
 
  call status(lunot,lot,il,im)
  if(il==0) call bort('BUFRLIB: WRDXTB - OUTPUT BUFR FILE IS CLOSED, IT MUST BE OPEN FOR OUTPUT')
  if(il<0) call bort('BUFRLIB: WRDXTB - OUTPUT BUFR FILE IS OPEN FOR INPUT, IT MUST BE OPEN FOR OUTPUT')
 
  call status(lundx,ldx,il,im)
  if(il==0) call bort('BUFRLIB: WRDXTB - DX TABLE FILE IS CLOSED, IT MUST BE OPEN')
 
  ! If files are different, copy internal table information from lundx to lunot
 
  if(lundx/=lunot) call cpbfdx(ldx,lot)
 
  ! Generate and write out BUFR dictionary messages to lunot
 
  call dxmini(mgwa,mbyt,mby4,mbya,mbyb,mbyd)
 
  lda = ldxa(idxv+1)
  ldb = ldxb(idxv+1)
  ldd = ldxd(idxv+1)
  l30 = ld30(idxv+1)
 
  ! Table A information
 
  do i=1,ntba(lot)
    if(msgfull(mbyt,lda,maxbyt).or.(iupb(mgwa,mbya,8)==255)) then
      call msgwrt(lunot,mgwa,mbyt)
      call dxmini(mgwa,mbyt,mby4,mbya,mbyb,mbyd)
    endif
    mbit = 8*(mby4-1)
    call pkb(iupb(mgwa,mby4,24)+lda,24,mgwa,mbit)
    mbit = 8*(mbya-1)
    call pkb(iupb(mgwa,mbya,8)+1,8,mgwa,mbit)
    mbit = 8*(mbyb-1)
    call pkc(taba(i,lot),lda,mgwa,mbit)
    call pkb(0,8,mgwa,mbit)
    call pkb(0,8,mgwa,mbit)
    mbyt = mbyt+lda
    mbyb = mbyb+lda
    mbyd = mbyd+lda
  enddo
 
  ! Table B information
 
  do i=1,ntbb(lot)
    if(msgfull(mbyt,ldb,maxbyt).or.(iupb(mgwa,mbyb,8)==255)) then
      call msgwrt(lunot,mgwa,mbyt)
      call dxmini(mgwa,mbyt,mby4,mbya,mbyb,mbyd)
    endif
    mbit = 8*(mby4-1)
    call pkb(iupb(mgwa,mby4,24)+ldb,24,mgwa,mbit)
    mbit = 8*(mbyb-1)
    call pkb(iupb(mgwa,mbyb,8)+1,8,mgwa,mbit)
    mbit = 8*(mbyd-1)
    call pkc(tabb(i,lot),ldb,mgwa,mbit)
    call pkb(0,8,mgwa,mbit)
    mbyt = mbyt+ldb
    mbyd = mbyd+ldb
  enddo
 
  ! Table D information
 
  do i=1,ntbd(lot)
    nseq = iupm(tabd(i,lot)(ldd+1:ldd+1),8)
    lend = ldd+1 + l30*nseq
    if(msgfull(mbyt,lend,maxbyt).or.(iupb(mgwa,mbyd,8)==255)) then
      call msgwrt(lunot,mgwa,mbyt)
      call dxmini(mgwa,mbyt,mby4,mbya,mbyb,mbyd)
    endif
    mbit = 8*(mby4-1)
    call pkb(iupb(mgwa,mby4,24)+lend,24,mgwa,mbit)
    mbit = 8*(mbyd-1)
    call pkb(iupb(mgwa,mbyd,8)+1,8,mgwa,mbit)
    mbit = 8*(mbyt-4)
    call pkc(tabd(i,lot),ldd,mgwa,mbit)
    call pkb(nseq,8,mgwa,mbit)
    do j=1,nseq
      jj  = ldd+2 + (j-1)*2
      idn = iupm(tabd(i,lot)(jj:jj+1),16)
      call pkc(adn30(idn,l30),l30,mgwa,mbit)
    enddo
    mbyt = mbyt+lend
  enddo
 
  ! Write the unwritten (leftover) message.
 
  call msgwrt(lunot,mgwa,mbyt)
 
  ! Write out one additional (dummy) DX message containing zero subsets.  This will serve as a delimiter for this set of
  ! table messages within output unit lunot, just in case the next thing written to lunot ends up being another set of
  ! table messages.
 
  call dxmini(mgwa,mbyt,mby4,mbya,mbyb,mbyd)
  call getlens(mgwa,2,len0,len1,len2,l3,l4,l5)
  mbit = (len0+len1+len2+4)*8
  call pkb(0,16,mgwa,mbit)
  call msgwrt(lunot,mgwa,mbyt)
 
  return
end subroutine wrdxtb
 
subroutine stbfdx(lun,mesg)
 
  use modv_vars, only: maxcd, idxv
 
  use moda_tababd
 
  implicit none
 
  integer, intent(in) :: lun, mesg(*)
  integer nxstr, ldxa, ldxb, ldxd, ld30, ldxbd(10), ldxbe(10), ja, jb, idxs, i3, i, j, n, nd, ndd, idn, &
    jbit, len0, len1, len2, len3, l4, l5, lda, ldb, ldd, ldbd, ldbe, l30, ia, la, ib, lb, id, ld, iret, &
    ifxy, iupb, iupbs01, igetntbi, idn30
 
  character*128 bort_str
  character*128 tabb1, tabb2
  character*56 dxstr
  character*55 cseq
  character*50 dxcmp
  character*24 unit
  character*8 nemo
  character*6 numb, cidn
 
  common /dxtab/ nxstr(10), ldxa(10), ldxb(10), ldxd(10), ld30(10), dxstr(10)
 
  data ldxbd /38, 70, 8*0/
  data ldxbe /42, 42, 8*0/
 
  ! Statement functions
  ja(i) = ia+1+lda*(i-1)
  jb(i) = ib+1+ldb*(i-1)
 
  ! Get some preliminary information from the message
 
  idxs = iupbs01(mesg,'MSBT')+1
  if(idxs>idxv+1) idxs = iupbs01(mesg,'MTVL')+1
  if(ldxa(idxs)==0 .or. ldxb(idxs)==0 .or. ldxd(idxs)==0) call bort('BUFRLIB: STBFDX - UNEXPECTED DICTIONARY '// &
    'MESSAGE SUBTYPE OR LOCAL VERSION NUMBER (E.G., L.V.N. HIGHER THAN KNOWN)')
 
  call getlens(mesg,3,len0,len1,len2,len3,l4,l5)
  i3 = len0+len1+len2
  dxcmp = ' '
  jbit = 8*(i3+7)
  call upc(dxcmp,nxstr(idxs),mesg,jbit,.false.)
  if(dxcmp/=dxstr(idxs)) call bort('BUFRLIB: STBFDX - UNEXPECTED DICTIONARY MESSAGE CONTENTS')
 
  ! Section 4 - read definitions for Tables A, B and D
 
  lda  = ldxa(idxs)
  ldb  = ldxb(idxs)
  ldd  = ldxd(idxs)
  ldbd = ldxbd(idxs)
  ldbe = ldxbe(idxs)
  l30  = ld30(idxs)
 
  ia = i3+len3+5
  la = iupb(mesg,ia,8)
  ib = ja(la+1)
  lb = iupb(mesg,ib,8)
  id = jb(lb+1)
  ld = iupb(mesg,id,8)
 
  ! Table A
 
  do i=1,la
    n = igetntbi(lun,'A')
    jbit = 8*(ja(i)-1)
    call upc(taba(n,lun),lda,mesg,jbit,.true.)
    numb = '   '//taba(n,lun)(1:3)
    nemo = taba(n,lun)(4:11)
    cseq = taba(n,lun)(13:67)
    call stntbia(n,lun,numb,nemo,cseq)
  enddo
 
  ! Table B
 
  do i=1,lb
    n = igetntbi(lun,'B')
    jbit = 8*(jb(i)-1)
    call upc(tabb1,ldbd,mesg,jbit,.true.)
    jbit = 8*(jb(i)+ldbd-1)
    call upc(tabb2,ldbe,mesg,jbit,.true.)
    tabb(n,lun) = tabb1(1:ldxbd(idxv+1))//tabb2(1:ldxbe(idxv+1))
    numb = tabb(n,lun)(1:6)
    nemo = tabb(n,lun)(7:14)
    call nenubd(nemo,numb,lun)
    idnb(n,lun) = ifxy(numb)
    unit = tabb(n,lun)(71:94)
    call capit(unit)
    tabb(n,lun)(71:94) = unit
    ntbb(lun) = n
  enddo
 
  ! Table D
 
  do i=1,ld
    n = igetntbi(lun,'D')
    jbit = 8*id
    call upc(tabd(n,lun),ldd,mesg,jbit,.true.)
    numb = tabd(n,lun)(1:6)
    nemo = tabd(n,lun)(7:14)
    call nenubd(nemo,numb,lun)
    idnd(n,lun) = ifxy(numb)
    nd = iupb(mesg,id+ldd+1,8)
    if(nd>maxcd) then
      write(bort_str,'("BUFRLIB: STBFDX - NUMBER OF DESCRIPTORS IN '// &
        'TABLE D ENTRY ",A," IN BUFR TABLE (",I4,") EXCEEDS THE LIMIT (",I4,")")') nemo,nd,maxcd
      call bort(bort_str)
    endif
    do j=1,nd
      ndd = id+ldd+2 + (j-1)*l30
      jbit = 8*(ndd-1)
      call upc(cidn,l30,mesg,jbit,.true.)
      idn = idn30(cidn,l30)
      call pktdd(n,lun,idn,iret)
      if(iret<0) call bort('BUFRLIB: STBFDX - BAD RETURN FROM BUFRLIB ROUTINE PKTDD, SEE PREVIOUS WARNING MESSAGE')
    enddo
    id = id+ldd+1 + nd*l30
    if(iupb(mesg,id+1,8)==0) id = id+1
    ntbd(lun) = n
  enddo
 
  return
end subroutine stbfdx
 
integer function idxmsg( mesg ) result( iret )
 
  implicit none
 
  integer, intent(in) :: mesg(*)
  integer iupbs01
 
  ! Note that the following test relies upon logic within subroutine dxmini() which zeroes out the Section 1 date of
  ! all DX dictionary messages.
 
  if ( (iupbs01(mesg,'MTYP')==11) .and. &
         (iupbs01(mesg,'MNTH')==0) .and. (iupbs01(mesg,'DAYS')==0) .and. (iupbs01(mesg,'HOUR')==0) ) then
    iret = 1
  else
    iret = 0
  end if
 
  return
end function idxmsg
 
integer function igetntbi ( lun, ctb ) result(iret)
 
  use moda_tababd
 
  implicit none
 
  integer, intent(in) :: lun
  integer imax
 
  character, intent(in) :: ctb
  character*128 bort_str
 
  if ( ctb == 'A' ) then
    iret = ntba(lun) + 1
    imax = ntba(0)
  else if ( ctb == 'B' ) then
    iret = ntbb(lun) + 1
    imax = ntbb(0)
  else  ! ctb == 'D'
    iret = ntbd(lun) + 1
    imax = ntbd(0)
  endif
  if ( iret > imax ) then
    write(bort_str,'("BUFRLIB: IGETNTBI - NUMBER OF INTERNAL TABLE",A1," ENTRIES EXCEEDS THE LIMIT (",I4,")")') ctb, imax
    call bort(bort_str)
  endif
 
  return
end function igetntbi
 
subroutine nemtbax(lun,nemo,mtyp,msbt,inod)
 
  use moda_tababd
 
  implicit none
 
  integer, intent(in) :: lun
  integer, intent(out) :: mtyp, msbt, inod
  integer i
 
  character*(*), intent(in) :: nemo
  character*128 bort_str
 
  inod = 0
 
  ! Look for nemo in Table A
 
  do i=1,ntba(lun)
    if(taba(i,lun)(4:11)==nemo) then
      mtyp = idna(i,lun,1)
      msbt = idna(i,lun,2)
      inod = mtab(i,lun)
      if(mtyp<0 .or. mtyp>255) then
        write(bort_str,'("BUFRLIB: NEMTBAX - INVALID MESSAGE TYPE (",I4,") RETURNED FOR MENMONIC ",A)') mtyp, nemo
        call bort(bort_str)
      endif
      if(msbt<0 .or. msbt>255) then
        write(bort_str,'("BUFRLIB: NEMTBAX - INVALID MESSAGE SUBTYPE (",I4,") RETURNED FOR MENMONIC ",A)') msbt, nemo
        call bort(bort_str)
      endif
      exit
    endif
  enddo
 
  return
end subroutine nemtbax
 
subroutine nemtba(lun,nemo,mtyp,msbt,inod)
 
  implicit none
 
  integer, intent(in) :: lun
  integer, intent(out) :: mtyp, msbt, inod
 
  character*(*), intent(in) :: nemo
  character*128 bort_str
 
  ! Look for nemo in Table A
 
  call nemtbax(lun,nemo,mtyp,msbt,inod)
  if(inod==0) then
    write(bort_str,'("BUFRLIB: NEMTBA - CAN''T FIND MNEMONIC ",A)') nemo
    call bort(bort_str)
  endif
 
  return
end subroutine nemtba
 
subroutine nemtbb(lun,itab,unit,iscl,iref,ibit)
 
  use moda_tababd
 
  implicit none
 
  integer, intent(in) :: lun, itab
  integer, intent(out) :: iscl, iref, ibit
  integer idn, ierns
 
  character*128 bort_str
  character*24, intent(out) :: unit
  character*8 nemo
 
  if(itab<=0 .or. itab>ntbb(lun)) then
    write(bort_str,'("BUFRLIB: NEMTBB - ITAB (",I7,") NOT FOUND IN TABLE B")') itab
    call bort(bort_str)
  endif
 
  ! Pull out Table B information
 
  idn  = idnb(itab,lun)
  nemo = tabb(itab,lun)( 7:14)
  unit = tabb(itab,lun)(71:94)
  call strnum(tabb(itab,lun)( 95: 98),iscl,ierns)
  call strnum(tabb(itab,lun)( 99:109),iref,ierns)
  call strnum(tabb(itab,lun)(110:112),ibit,ierns)
 
  ! Check Table B contents
 
  if(unit(1:5)/='CCITT' .and. ibit>32) then
    write(bort_str,'("BUFRLIB: NEMTBB - BIT WIDTH FOR NON-CHARACTER TABLE B MNEMONIC ",A," (",I7,") IS > 32")') nemo,ibit
    call bort(bort_str)
  endif
  if(unit(1:5)=='CCITT' .and. mod(ibit,8)/=0) then
    write(bort_str,'("BUFRLIB: NEMTBB - BIT WIDTH FOR CHARACTER TABLE B MNEMONIC ",A," (",I7,") IS NOT A MULTIPLE OF 8")') &
      nemo,ibit
    call bort(bort_str)
  endif
 
  return
end subroutine nemtbb
 
subroutine nemtbd(lun,itab,nseq,nems,irps,knts)
 
  use modv_vars, only: maxcd
 
  use moda_tababd
 
  implicit none
 
  integer, intent(in) :: lun, itab
  integer, intent(out) :: nseq, irps(*), knts(*)
  integer i, j, ndsc, idsc, iret
 
  character*128 bort_str
  character*8, intent(out) :: nems(*)
  character*8 nemo, nemt, nemf
  character tab
 
  if(itab<=0 .or. itab>ntbd(lun)) then
    write(bort_str,'("BUFRLIB: NEMTBD - ITAB (",I7,") NOT FOUND IN TABLE D")') itab
    call bort(bort_str)
  endif
 
  ! Clear the return values
 
  nseq = 0
 
  do i=1,maxcd
    nems(i) = ' '
    irps(i) = 0
    knts(i) = 0
  enddo
 
  ! Parse the Table D entry
 
  nemo = tabd(itab,lun)(7:14)
  idsc = idnd(itab,lun)
  call uptdd(itab,lun,0,ndsc)
 
  ! Loop through each child mnemonic
 
  do j=1,ndsc
    if(nseq+1>maxcd) then
      write(bort_str,'("BUFRLIB: NEMTBD - THERE ARE MORE THAN '// &
        '(",I4,") DESCRIPTORS (THE LIMIT) IN TABLE D SEQUENCE MNEMONIC ",A)') maxcd, nemo
      call bort(bort_str)
    endif
    call uptdd(itab,lun,j,idsc)
    call numtab(lun,idsc,nemt,tab,iret)
    if(tab=='R') then
      if(iret<0) then
        ! Regular (i.e. non-delayed) replication
        irps(nseq+1) = 1
        knts(nseq+1) = abs(iret)
      elseif(iret>0) then
        ! Delayed replication
        irps(nseq+1) = iret
      endif
    elseif(tab=='F') then
      ! Replication factor
      irps(nseq+1) = iret
    elseif(tab=='D'.or.tab=='C') then
      nseq = nseq+1
      nems(nseq) = nemt
    elseif(tab=='B') then
      nseq = nseq+1
      if((nemt(1:1)=='.').and.(j<ndsc)) then
        ! This is a "following value" mnemonic
        call uptdd(itab,lun,j+1,idsc)
        call numtab(lun,idsc,nemf,tab,iret)
        call rsvfvm(nemt,nemf)
      endif
      nems(nseq) = nemt
    endif
  enddo
 
  return
end subroutine nemtbd
 
recursive subroutine nemdefs ( lunit, nemo, celem, cunit, iret )
 
  use modv_vars, only: im8b
 
  use moda_tababd
 
  implicit none
 
  integer, intent(in) :: lunit
  integer, intent(out) :: iret
  integer my_lunit, lun, il, im, idn, iloc, ls
 
  character*(*), intent(in) :: nemo
  character*(*), intent(out) :: celem, cunit
  character tab
 
  ! Check for I8 integers.
 
  if(im8b) then
    im8b=.false.
    call x84 ( lunit, my_lunit, 1 )
    call nemdefs ( my_lunit, nemo, celem, cunit, iret )
    call x48 ( iret, iret, 1 )
    im8b=.true.
    return
  endif
 
  iret = -1
 
  ! Get lun from lunit.
 
  call status( lunit, lun, il, im )
  if ( il == 0 ) return
 
  ! Find the requested mnemonic in the internal Table B arrays.
 
  call nemtab( lun, nemo, idn, tab, iloc )
  if ( ( iloc == 0 ) .or. ( tab /= 'B' ) ) return
 
  ! Get the element name and units of the requested mnemonic.
 
  celem = ' '
  ls = min(len(celem),55)
  celem(1:ls) = tabb(iloc,lun)(16:15+ls)
 
  cunit = ' '
  ls = min(len(cunit),24)
  cunit(1:ls) = tabb(iloc,lun)(71:70+ls)
 
  iret = 0
 
  return
end subroutine nemdefs
 
subroutine nenubd(nemo,numb,lun)
 
  use moda_tababd
 
  implicit none
 
  character, intent(in) :: nemo*8, numb*6
  character*128 bort_str
 
  integer, intent(in) :: lun
  integer n
 
  do n=1,ntbb(lun)
    if(numb==tabb(n,lun)(1:6)) then
      write(bort_str,'("BUFRLIB: NENUBD - TABLE B FXY VALUE (",A,") HAS ALREADY BEEN DEFINED (DUPLICATE)")') numb
      call bort(bort_str)
    endif
    if(nemo==tabb(n,lun)(7:14)) then
      write(bort_str,'("BUFRLIB: NENUBD - TABLE B MNEMONIC (",A,") HAS ALREADY BEEN DEFINED (DUPLICATE)")') nemo
      call bort(bort_str)
    endif
  enddo
 
  do n=1,ntbd(lun)
    if(numb==tabd(n,lun)(1:6)) then
      write(bort_str,'("BUFRLIB: NENUBD - TABLE D FXY VALUE (",A,") HAS ALREADY BEEN DEFINED (DUPLICATE)")') numb
      call bort(bort_str)
    endif
    if(nemo==tabd(n,lun)(7:14)) then
      write(bort_str,'("BUFRLIB: NENUBD - TABLE D MNEMONIC (",A,") HAS ALREADY BEEN DEFINED (DUPLICATE)")') nemo
      call bort(bort_str)
    endif
  enddo
 
  return
end subroutine nenubd
 
subroutine stntbia ( n, lun, numb, nemo, celsq )
 
  use moda_tababd
 
  implicit none
 
  integer, intent(in) :: n, lun
  integer i, mtyp, msbt
 
  character*(*), intent(in) :: numb, nemo, celsq
  character*128 bort_str
 
  ! Confirm that neither nemo nor numb has already been defined within the internal BUFR Table A for the given lun.
 
  do i=1,ntba(lun)
    if(numb(4:6)==taba(i,lun)(1:3)) then
      write(bort_str,'("BUFRLIB: STNTBIA - TABLE A FXY VALUE (",A,") HAS ALREADY BEEN DEFINED (DUPLICATE)")') numb
      call bort(bort_str)
    endif
    if(nemo(1:8)==taba(i,lun)(4:11)) then
      write(bort_str,'("BUFRLIB: STNTBIA - TABLE A MNEMONIC (",A,") HAS ALREADY BEEN DEFINED (DUPLICATE)")') nemo
      call bort(bort_str)
    endif
  enddo
 
  ! Store the values within the internal BUFR Table A.
 
  taba(n,lun)(1:3) = numb(4:6)
  taba(n,lun)(4:11) = nemo(1:8)
  taba(n,lun)(13:67) = celsq(1:55)
 
  ! Decode and store the message type and subtype.
 
  if ( verify( nemo(3:8), '1234567890' ) == 0 ) then
    ! Message type & subtype obtained directly from Table A mnemonic
    read ( nemo,'(2X,2I3)') mtyp, msbt
    idna(n,lun,1) = mtyp
    idna(n,lun,2) = msbt
  else
    ! Message type obtained from Y value of Table A seq. descriptor
    read ( numb(4:6),'(I3)') idna(n,lun,1)
    ! Message subtype hardwired to zero
    idna(n,lun,2) = 0
  endif
 
  ! Update the count of internal Table A entries.
 
  ntba(lun) = n
 
  return
end subroutine stntbia
 
subroutine stntbi ( n, lun, numb, nemo, celsq )
 
  use moda_tababd
 
  implicit none
 
  integer, intent(in) :: n, lun
  integer ifxy
 
  character*(*), intent(in) :: numb, nemo, celsq
 
  call nenubd ( nemo, numb, lun )
 
  if ( numb(1:1) == '0') then
    idnb(n,lun) = ifxy(numb)
    tabb(n,lun)(1:6) = numb(1:6)
    tabb(n,lun)(7:14) = nemo(1:8)
    tabb(n,lun)(16:70) = celsq(1:55)
    ntbb(lun) = n
  else  ! numb(1:1) == '3'
    idnd(n,lun) = ifxy(numb)
    tabd(n,lun)(1:6) = numb(1:6)
    tabd(n,lun)(7:14) = nemo(1:8)
    tabd(n,lun)(16:70) = celsq(1:55)
    ntbd(lun) = n
  endif
 
  return
end subroutine stntbi
 
subroutine pktdd(id,lun,idn,iret)
 
  use modv_vars, only: maxcd, iprt, idxv
 
  use moda_tababd
 
  implicit none
 
  integer, intent(in) :: id, lun, idn
  integer, intent(out) :: iret
  integer nxstr, ldxa, ldxb, ldxd, ld30, ldd, nd, idm, iupm
 
  character*128 errstr
  character*56 dxstr
 
  common /dxtab/ nxstr(10), ldxa(10), ldxb(10), ldxd(10), ld30(10), dxstr(10)
 
  ! ldd points to the byte within tabd(id,lun) which contains (in packed integer format) a count of the number of child
  ! mnemonics stored thus far for this parent mnemonic.
  ldd = ldxd(idxv+1)+1
 
  ! Zero the counter if idn is zero
  if(idn==0) then
    call ipkm(tabd(id,lun)(ldd:ldd),1,0)
    iret = 0
    return
  endif
 
  ! Update the stored descriptor count for this Table D entry.  nd is the (unpacked) count of the number of child mnemonics
  ! stored thus far for this parent mnemonic.
  nd = iupm(tabd(id,lun)(ldd:ldd),8)
 
  if(nd<0 .or. nd==maxcd) then
    if(iprt>=0) then
      call errwrt('+++++++++++++++++++++WARNING+++++++++++++++++++++++')
      if(nd<0) then
        write ( unit=errstr, fmt='(A,I4,A)' ) 'BUFRLIB: PKTDD - BAD COUNTER VALUE (=', nd, ') - RETURN WITH IRET = -1'
      else
        write ( unit=errstr, fmt='(A,I4,A,A)' ) 'BUFRLIB: PKTDD - MAXIMUM NUMBER OF CHILD MNEMONICS (=', &
          maxcd, ') ALREADY STORED FOR THIS PARENT - RETURN WITH IRET = -1'
      endif
      call errwrt(errstr)
      call errwrt('+++++++++++++++++++++WARNING+++++++++++++++++++++++')
      call errwrt(' ')
    endif
    iret = -1
    return
  else
    nd = nd+1
    call ipkm(tabd(id,lun)(ldd:ldd),1,nd)
    iret = nd
  endif
 
  ! Pack and store the descriptor.  idm points to the starting byte within tabd(id,lun) at which the idn value for this
  ! child mnemonic will be stored (as a packed integer of width = 2 bytes).
  idm = ldd+1 + (nd-1)*2
  call ipkm(tabd(id,lun)(idm:idm+1),2,idn)
 
  return
end subroutine pktdd
 
subroutine uptdd(id,lun,ient,iret)
 
  use modv_vars, only: idxv
 
  use moda_tababd
 
  implicit none
 
  integer, intent(in) :: id, lun, ient
  integer, intent(out) :: iret
  integer nxstr, ldxa, ldxb, ldxd, ld30, ldd, ndsc, idsc, iupm
 
  character*128 bort_str
  character*56 dxstr
 
  common /dxtab/ nxstr(10), ldxa(10), ldxb(10), ldxd(10), ld30(10), dxstr(10)
 
  ! Check if ient is in bounds
 
  ldd = ldxd(idxv+1)+1
  ndsc = iupm(tabd(id,lun)(ldd:ldd),8)
  if(ient==0) then
    iret = ndsc
    return
  elseif(ient<0 .or. ient>ndsc) then
    write(bort_str,'("BUFRLIB: UPTDD - VALUE OF THIRD ARGUMENT IENT (INPUT) IS OUT OF RANGE (IENT =",I4,")")') ient
    call bort(bort_str)
  endif
 
  ! Return the descriptor indicated by ient
 
  idsc = ldd+1 + (ient-1)*2
  iret = iupm(tabd(id,lun)(idsc:idsc+1),16)
 
  return
end subroutine uptdd
 
subroutine rsvfvm(nem1,nem2)
 
  implicit none
 
  character*8, intent(inout) :: nem1
  character*8, intent(in) :: nem2
 
  integer i, j
 
  do i=1,len(nem1)
    if(i==1) then
      ! Skip the initial ".", and initialize J.
      j = 1
    else
      if(nem1(i:i)=='.') then
        nem1(i:i) = nem2(j:j)
        j = j+1
      endif
    endif
  enddo
 
  return
end subroutine rsvfvm