dumpdata.F90

Go to the documentation of this file.

 
 
recursive subroutine ufbdmp(lunin,luprt)
 
  use bufrlib
 
  use modv_vars, only: im8b
 
  use moda_usrint
  use moda_msgcwd
  use moda_tababd
  use moda_tables
 
  implicit none
 
  integer, intent(in) :: lunin, luprt
  integer, parameter :: mxfv = 31
  integer ifv(mxfv), my_lunin, my_luprt, luout, lunit, lun, il, im, nv, nd, it, ib, is, ir, jp, lk, jb, &
    idn, nifv, nchr, n, ii, ipt, isz, isize, ibfms, icbfms, bort_target_set
 
  character lchr2*120, lchr*20, pmiss*20, bits*14, tg*10, tg_rj*10, vc*8, fmtf*7, tp*3, tab, you
 
  real*8 vl
 
  equivalence(vl,vc)
 
  data pmiss /'             MISSING'/
  data you /'Y'/
 
  ! Check for I8 integers
 
  if(im8b) then
    im8b=.false.
    call x84(lunin,my_lunin,1)
    call x84(luprt,my_luprt,1)
    call ufbdmp(my_lunin,my_luprt)
    im8b=.true.
    return
  endif
 
  ! If we're catching bort errors, set a target return location if one doesn't already exist.
 
  if (bort_target_set() == 1) then
    call catch_bort_ufbdmp_c(lunin,luprt)
    call bort_target_unset
    return
  endif
 
  if(luprt==0) then
    luout = 6
  else
    luout = luprt
  endif
 
  ! Check the file status and inode
 
  lunit = abs(lunin)
  call status(lunit,lun,il,im)
  if(il==0) call bort('BUFRLIB: UFBDMP - INPUT BUFR FILE IS CLOSED, IT MUST BE OPEN FOR INPUT')
  if(il>0) call bort('BUFRLIB: UFBDMP - INPUT BUFR FILE IS OPEN FOR OUTPUT, IT MUST BE OPEN FOR INPUT')
  if(im==0) call bort('BUFRLIB: UFBDMP - A MESSAGE MUST BE OPEN IN INPUT BUFR FILE, NONE ARE')
  if(inode(lun)/=inv(1,lun)) call bort('BUFRLIB: UFBDMP - LOCATION OF INTERNAL TABLE FOR '// &
    'INPUT BUFR FILE DOES NOT AGREE WITH EXPECTED LOCATION IN INTERNAL SUBSET ARRAY')
 
  ! Dump the contents of @ref moda_usrint for unit abs(lunin)
 
  do nv=1,nval(lun)
    if(luprt==0 .and. mod(nv,20)==0) then
 
      ! When luprt=0, the output will be scrolled, 20 elements at a time
 
      print*,'(<enter> for MORE, q <enter> to QUIT)'
      read(5,'(A1)') you
 
      ! If the terminal enters "q" followed by "<enter>", then scrolling will end and the subroutine will return to the
      ! calling program
 
      if(you=='q') then
        print*
        print*,'==> You have chosen to stop the dumping of this subset'
        print*
        return
      endif
    endif
    nd = inv(nv,lun)
    vl = val(nv,lun)
    tg = tag(nd)
    tp = typ(nd)
    it = itp(nd)
    ib = ibt(nd)
    is = isc(nd)
    ir = irf(nd)
    jp = jump(nd)
    lk = link(nd)
    jb = jmpb(nd)
    tg_rj = adjustr(tg)
    if(tp/='CHR') then
      bits = '              '
      if(it==2) then
        call nemtab(lun,tg,idn,tab,n)
        if(tabb(n,lun)(71:75)=='FLAG') then
 
          ! Print a listing of the bits corresponding to this value.
 
          call upftbv(lunit,tg,vl,mxfv,ifv,nifv)
          if(nifv>0) then
            bits(1:1) = '('
            ipt = 2
            do ii=1,nifv
              isz = isize(ifv(ii))
              write(fmtf,'(A2,I1,A4)') '(I', isz, ',A1)'
              if((ipt+isz)<=14) then
                write(bits(ipt:ipt+isz),fmtf) ifv(ii), ','
                ipt = ipt + isz + 1
              else
                bits(2:13) = 'MANY BITS ON'
                ipt = 15
              endif
            enddo
            bits(ipt-1:ipt-1) = ')'
          endif
        endif
      endif
      if(ibfms(vl)/=0) then
        write(luout,'(I5,1X,A3,A1,I1,1X,A10,1X,   A20,  14X,7(1X,I5))') nv,tp,'-',it,tg_rj,pmiss,ib,is,ir,nd,jp,lk,jb
      else
        if(lunit==lunin) then
          write(luout,'(I5,1X,A3,A1,I1,1X,A10,5X,G15.6,1X,A14,7(1X,I5))') nv,tp,'-',it,tg_rj,vl,bits,ib,is,ir,nd,jp,lk,jb
        else
          write(luout,'(I5,1X,A3,A1,I1,1X,A10,5X,F15.6,1X,A14,7(1X,I5))') nv,tp,'-',it,tg_rj,vl,bits,ib,is,ir,nd,jp,lk,jb
        endif
      endif
    else
      nchr=ib/8
      if(nchr>8) then
        call readlc(lunit,lchr2,tg_rj)
        if (icbfms(lchr2,nchr)/=0) then
          lchr = pmiss
        else
          lchr = lchr2(1:20)
        endif
      else
        if(ibfms(vl)/=0) then
          lchr = pmiss
        else
          lchr = vc
        endif
      endif
      if ( nchr<=20 .or. lchr==pmiss ) then
        lchr = adjustr(lchr)
        write(luout,'(I5,1X,A3,A1,I1,1X,A10,1X,   A20,  14X,7(1X,I5))') nv,tp,'-',it,tg_rj,lchr,ib,is,ir,nd,jp,lk,jb
      else
        write(luout,'(I5,1X,A3,A1,I1,1X,A10,1X,     A,      7(1X,I5))') nv,tp,'-',it,tg_rj,lchr2(1:nchr),ib,is,ir,nd,jp,lk,jb
      endif
    endif
  enddo
 
  write(luout,'(/A/)') ' >>> END OF SUBSET <<< '
 
  return
end subroutine ufbdmp
 
recursive subroutine ufdump(lunit,luprt)
 
  use bufrlib
 
  use modv_vars, only: im8b
 
  use moda_usrint
  use moda_msgcwd
  use moda_tababd
  use moda_tables
  use moda_tablef
  use moda_nrv203
 
  implicit none
 
  integer, intent(in) :: lunit, luprt
  integer, parameter :: mxfv = 31 , mxcfdp = 5, mxseq = 10, mxls = 10
  integer ifv(mxfv), icfdp(mxcfdp), idxrep(mxseq), numrep(mxseq), lsqnam(mxseq), lsct(mxls), my_lunit, my_luprt, &
    nseq, nls, lcfmeang, luout, lun, il, im, node, lnm2, lnm3, itmp, ityp, ii, jj, nifv, nv, n, nchr, idn, ipt, &
    nrfe, nout, lcfmg, ifvd, iersf, ierbd, ierft, isz, isize, iscl, ireadmt, ibfms, icbfms, imrkopr, bort_target_set
  integer*8 ival
 
  real*8 rval
 
  character cfmeang*120, lchr2*120, fmt*80, desc*64, unit*24, lchr*20, pmiss*20, nemo3*15, nemo*10, nemo2*10, tagrfe*10, &
    seqnam(mxseq)*10, lsnemo(mxls)*10, nemod*8, cval*8, fmtf*7, numb*6, type*3, tab, you
 
  logical track, found, rdrv
 
  equivalence(rval,cval)
 
  data pmiss /'             MISSING'/
  data you /'Y'/
 
  ! Check for I8 integers
 
  if(im8b) then
    im8b=.false.
    call x84(lunit,my_lunit,1)
    call x84(luprt,my_luprt,1)
    call ufdump(my_lunit,my_luprt)
    im8b=.true.
    return
  endif
 
  ! If we're catching bort errors, set a target return location if one doesn't already exist.
 
  if (bort_target_set() == 1) then
    call catch_bort_ufdump_c(lunit,luprt)
    call bort_target_unset
    return
  endif
 
  nseq = 0
  nls = 0
  lcfmeang = len(cfmeang)
 
  if(luprt==0) then
    luout = 6
  else
    luout = luprt
  endif
 
  ! Check the file status and inode
 
  call status(lunit,lun,il,im)
  if(il==0) call bort('BUFRLIB: UFDUMP - INPUT BUFR FILE IS CLOSED, IT MUST BE OPEN FOR INPUT')
  if(il>0) call bort('BUFRLIB: UFDUMP - INPUT BUFR FILE IS OPEN FOR OUTPUT, IT MUST BE OPEN FOR INPUT')
  if(im==0) call bort('BUFRLIB: UFDUMP - A MESSAGE MUST BE OPEN IN INPUT BUFR FILE, NONE ARE')
  if(inode(lun)/=inv(1,lun)) call bort('BUFRLIB: UFDUMP - LOCATION OF INTERNAL TABLE FOR '// &
    'INPUT BUFR FILE DOES NOT AGREE WITH EXPECTED LOCATION IN INTERNAL SUBSET ARRAY')
 
  write(luout,'(/,2A,/)') 'MESSAGE TYPE ',tag(inode(lun))
 
  ! If code/flag table details are being printed, and if this is the first subset of a new message, then
  ! make sure the appropriate master tables have been read in to memory for this message.
 
  if(cdmf=='Y' .and. nsub(lun)==1) itmp = ireadmt(lun)
 
  ! Dump the contents of @ref moda_usrint for unit lunit
 
  do nv=1,nval(lun)
    if(luprt==0 .and. mod(nv,20)==0) then
 
      ! When luprt=0, the output will be scrolled, 20 elements at a time
 
      print*,'(<enter> for MORE, q <enter> to QUIT)'
      read(5,'(A1)') you
 
      ! If the terminal enters "q" followed by "<enter>", then scrolling will end and the subroutine will return to the
      ! calling program
 
      if(you=='q') then
        print*
        print*,'==> You have chosen to stop the dumping of this subset'
        print*
        return
      endif
    endif
 
    node = inv(nv,lun)
    nemo = tag(node)
    ityp = itp(node)
    type = typ(node)
 
    if(ityp>=1.and.ityp<=3) then
      call nemtab(lun,nemo,idn,tab,n)
      if(n>0) then
         numb = tabb(n,lun)(1:6)
         desc = tabb(n,lun)(16:70)
         unit = tabb(n,lun)(71:94)
      endif
      rval = val(nv,lun)
    endif
 
    if((ityp==0).or.(ityp==1)) then
 
      ! Sequence descriptor or delayed descriptor replication factor
 
      if((type=='REP').or.(type=='DRP').or.(type=='DRB').or.(type=='DRS')) then
 
        ! Print the number of replications
        nseq = nseq+1
        if(nseq>mxseq) call bort('BUFRLIB: UFDUMP - MXSEQ OVERFLOW')
        if(type=='REP') then
          numrep(nseq) = irf(node)
        else
          numrep(nseq) = nint(rval)
        endif
        call strsuc(nemo,nemo2,lnm2)
        fmt = '(11X,A,I6,1X,A)'
        write(luout,fmt) nemo2(1:lnm2), numrep(nseq), 'REPLICATIONS'
 
        ! How many times is this sequence replicated?
        if(numrep(nseq)>1) then
          ! Track the sequence
          seqnam(nseq) = nemo2
          lsqnam(nseq) = lnm2
          idxrep(nseq) = 1
        else
          ! Don't bother
          nseq = nseq-1
        endif
      elseif( ((type=='SEQ').or.(type=='RPC').or.(type=='RPS')) .and. (nseq>0) ) then
 
        ! Is this one of the sequences being tracked?
        ii = nseq
        track = .false.
        call strsuc(nemo,nemo2,lnm2)
        do while ((ii>=1).and.(.not.track))
          if(nemo2(1:lnm2)==seqnam(ii)(2:lsqnam(ii)-1)) then
            track = .true.
            ! Mark this level in the output
            fmt = '(4X,A,2X,A,2X,A,I6,2X,A)'
            write(luout,fmt) '++++++', nemo2(1:lnm2), 'REPLICATION #', idxrep(ii), '++++++'
            if(idxrep(ii)<numrep(ii)) then
              ! There are more levels to come
              idxrep(ii) = idxrep(ii)+1
            else
              ! This was the last level for this sequence, so stop tracking it
              nseq = nseq-1
            endif
          else
            ii = ii-1
          endif
        enddo
      endif
    elseif(ityp==2) then
 
      ! Numeric value.
 
      ! First check if this node contains a redefined reference value.  If so, modify the desc field to label it as such.
      jj = 1
      rdrv = .false.
      do while ((jj<=nnrv).and.(.not.rdrv))
        if (node==inodnrv(jj)) then
          rdrv = .true.
          desc = 'New reference value for ' // nemo
          unit = ' '
        else
          jj = jj + 1
        endif
      enddo
 
      ! Next check if this element refers to another element via a bitmap.  If so, modify the desc field to identify the
      ! referred element.
      nrfe = nrfelm(nv,lun)
      if(nrfe>0) then
        tagrfe = tag(inv(nrfe,lun))
        jj = 48
        do while((jj>=1).and.(desc(jj:jj)==' '))
          jj = jj - 1
        enddo
        if(jj<=33) desc(jj+1:jj+15) = ' for ' // tagrfe
      endif
 
      ! Now print the value
      if(ibfms(rval)/=0) then
        ! The value is "missing".
        fmt = '(A6,2X,A10,2X,A20,2X,A24,6X,A48)'
        write(luout,fmt) numb,nemo,pmiss,unit,desc
      else
        fmt = '(A6,2X,A10,2X,      ,2X,A24,6X,A48)'
        ! Based upon the corresponding scale factor, select an appropriate format for the printing of this value.
        if(imrkopr(nemo)==1) then
          iscl = isc(inv(nrfe,lun))
        else
          iscl = isc(node)
        endif
        if(iscl>0) then
          write(fmt(15:20),'(A,I2)') 'F20.', iscl
        else
          write(fmt(18:20),'(A)') 'I20'
        endif
        if(unit(1:4)=='FLAG') then
          ! Print a listing of the bits corresponding to this value.
           call upftbv(lunit,nemo,rval,mxfv,ifv,nifv)
           if(nifv>0) then
             unit(11:11) = '('
             ipt = 12
             do ii=1,nifv
               isz = isize(ifv(ii))
               write(fmtf,'(A2,I1,A4)') '(I', isz, ',A1)'
               if((ipt+isz)<=24) then
                 write(unit(ipt:ipt+isz),fmtf) ifv(ii), ','
                 ipt = ipt + isz + 1
               else
                 unit(12:23) = 'MANY BITS ON'
                 ipt = 25
               endif
             enddo
             unit(ipt-1:ipt-1) = ')'
           endif
        endif
        if(iscl>0) then
          write(luout,fmt) numb,nemo,rval,unit,desc
        else
          ival = nint(rval,8)
          write(luout,fmt) numb,nemo,ival,unit,desc
        endif
        if( (unit(1:4)=='FLAG' .or. unit(1:4)=='CODE') .and. (cdmf=='Y') ) then
          ! Print the meanings of the code and flag values.
          fmt = '(31X,I8,A,A)'
          if(unit(1:4)=='CODE') then
            nifv = 1
            ifv(nifv) = nint(rval)
          endif
          do ii=1,nifv
            icfdp(1) = (-1)
            call srchtbf_c(idn,ifv(ii),icfdp(1),mxcfdp,-1,cfmeang,lcfmeang,lcfmg,iersf)
            if(iersf==0) then
              write(luout,fmt) ifv(ii),' = ',cfmeang(1:lcfmg)
            elseif(iersf<0) then
              write(luout,fmt) ifv(ii),' = ','***THIS IS AN ILLEGAL/UNDEFINED VALUE***'
            else
              ! The meaning of this value is dependent on the value of another mnemonic in the report.  Look for
              ! that other mnemonic within the report and then use it and its associated value to retrieve and print
              ! the proper meaning from the code/flag tables.
              ierft = (-1)
              jj = 0
              do while((jj<iersf).and.(ierft<0))
                jj = jj + 1
                call numtbd(lun,icfdp(jj),nemod,tab,ierbd)
                if((ierbd>0).and.(tab=='B')) call fstag(lun,nemod,-1,nv,nout,ierft)
              enddo
              if(ierft==0) then
                ifvd = nint(val(nout,lun))
                if(jj>1) icfdp(1) = icfdp(jj)
                call srchtbf_c(idn,ifv(ii),icfdp(1),mxcfdp,ifvd,cfmeang,lcfmeang,lcfmg,iersf)
                if(iersf==0) write(luout,fmt) ifv(ii),' = ', cfmeang(1:lcfmg)
              endif
            endif
          enddo
        endif
      endif
    elseif(ityp==3) then
 
      ! Character (CCITT IA5) value
 
      nchr = ibt(node)/8
 
      if(ibfms(rval)/=0) then
        lchr = pmiss
      else if(nchr<=8) then
        lchr = cval
      else
        ! Track the number of occurrences of this long character string, so that we can properly output each one.
        ii = 1
        found = .false.
        do while((ii<=nls).and.(.not.found))
          if(nemo==lsnemo(ii)) then
            found = .true.
          else
            ii = ii + 1
          endif
        enddo
        if(.not.found) then
          nls = nls+1
          if(nls>mxls) call bort('BUFRLIB: UFDUMP - MXLS OVERFLOW')
          lsnemo(nls) = nemo
          lsct(nls) = 1
          nemo3 = nemo
        else
          call strsuc(nemo,nemo3,lnm3)
          lsct(ii) = lsct(ii) + 1
          write(fmtf,'(A,I1,A)') '(2A,I', isize(lsct(ii)), ')'
          write(nemo3,fmtf) nemo(1:lnm3), '#', lsct(ii)
        endif
 
        call readlc(lunit,lchr2,nemo3)
        if (icbfms(lchr2,nchr)/=0) then
          lchr = pmiss
        else
          lchr = lchr2(1:20)
        endif
      endif
 
      if ( nchr<=20 .or. lchr==pmiss ) then
        lchr = adjustr(lchr)
        fmt = '(A6,2X,A10,2X,A20,2X,"(",I2,")",A24,2X,A48)'
        write(luout,fmt) numb,nemo,lchr,nchr,unit,desc
      else
        fmt = '(A6,2X,A10,2X,A,2X,"(",I3,")",A23,2X,A48)'
        write(luout,fmt) numb,nemo,lchr2(1:nchr),nchr,unit,desc
      endif
    endif
 
  enddo
 
  write(luout,'(/A/)') ' >>> END OF SUBSET <<< '
 
  return
end subroutine ufdump
 
recursive subroutine dxdump(lunit,ldxot)
 
  use bufrlib
 
  use modv_vars, only: im8b, reps, fxy_fbit, fxy_sbyct, fxy_drp16, fxy_drp8, fxy_drp8s, fxy_drp1, &
    fxy_drf16, fxy_drf8, fxy_drf1
 
  use moda_tababd
  use moda_nmikrp
 
  implicit none
 
  integer, intent(in) :: lunit, ldxot
  integer my_lunit, my_ldxot, lun, il, im, n, na, nc, nch, ic, icms, nseq, bort_target_set
 
  character card*80, cardi1*80, cardi2*80, cardi3*80, cardi4*80, cmstr*20, wrk3*10, wrk1*8, wrk2*8, adn*6
 
  logical tbskip, tdskip, xtrci1
 
  data cardi1  /'|          |        |                                                          |'/
  data cardi2  /'|          |                                                                   |'/
  data cardi3  /'|          |      |             |     |                          |-------------|'/
  data cardi4  /'|------------------------------------------------------------------------------|'/
 
  ! Statement functions
  tbskip(adn) = ((adn==fxy_sbyct).or.(adn==fxy_fbit).or.(adn==fxy_drf16).or.(adn==fxy_drf8).or.(adn==fxy_drf1))
  tdskip(adn) = ((adn==fxy_drp16).or.(adn==fxy_drp8).or.(adn==fxy_drp8s).or.(adn==fxy_drp1))
 
  ! Check for I8 integers.
 
  if(im8b) then
    im8b=.false.
    call x84(lunit,my_lunit,1)
    call x84(ldxot,my_ldxot,1)
    call dxdump(my_lunit,my_ldxot)
    im8b=.true.
    return
  endif
 
  ! If we're catching bort errors, set a target return location if one doesn't already exist.
 
  if (bort_target_set() == 1) then
    call catch_bort_dxdump_c(lunit,ldxot)
    call bort_target_unset
    return
  endif
 
  ! Determine lun from lunit.
 
  call status(lunit,lun,il,im)
  if(il==0) call bort('BUFRLIB: DXDUMP - BUFR FILE IS CLOSED, IT MUST BE OPEN')
 
  ! Create and write out (to ldxot) the header cards for the descriptor definition section.
 
  card=cardi4
  card( 1: 1)='.'
  card(80:80)='.'
  write (ldxot,'(A)') card
 
  card=cardi4
  card( 2: 2)=' '
  card(79:79)=' '
  card(15:64)='   USER DEFINITIONS FOR TABLE-A TABLE-B TABLE D   '
  write (ldxot,'(A)') card
 
  write (ldxot,'(A)') cardi4
 
  card=cardi1
  card( 3:10)='MNEMONIC'
  card(14:19)='NUMBER'
  card(23:33)='DESCRIPTION'
  write (ldxot,'(A)') card
 
  card=cardi4
  card(12:12)='|'
  card(21:21)='|'
  write (ldxot,'(A)') card
 
  ! Create and write out (to ldxot) the Table D descriptor definition cards.
 
  write (ldxot,'(A)') cardi1
 
  xtrci1=.false.
  do n=1,ntbd(lun)
    if(.not.tdskip(tabd(n,lun)(1:6))) then
      card=cardi1
      card( 3:10)=tabd(n,lun)( 7:14)
      card(14:19)=tabd(n,lun)( 1: 6)
      card(23:77)=tabd(n,lun)(16:70)
      ! Check if this Table D mnemonic is also a Table A mnemonic. If so, then label it as such and also check if it is the
      ! last of the Table A mnemonics, in which case an extra cardi1 line will be written to ldxot in order to separate
      ! the Table A mnemonics from the other Table D mnemonics.
      do na=1,ntba(lun)
        if(taba(na,lun)(4:11)==tabd(n,lun)(7:14)) then
          card(14:14)='A'
          if(na==ntba(lun)) xtrci1=.true.
          exit
        end if
      end do
      write (ldxot,'(A)') card
      if(xtrci1) then
        write (ldxot,'(A)') cardi1
        xtrci1=.false.
      end if
    end if
  end do
 
  ! Create and write out (to ldxot) the Table B descriptor definition cards.
 
  write (ldxot,'(A)') cardi1
 
  do n=1,ntbb(lun)
    if(.not.tbskip(tabb(n,lun)(1:6))) then
      card=cardi1
      card( 3:10)=tabb(n,lun)( 7:14)
      card(14:19)=tabb(n,lun)( 1: 6)
      card(23:77)=tabb(n,lun)(16:70)
      write (ldxot,'(A)') card
    end if
  end do
 
  write (ldxot,'(A)') cardi1
 
  ! Create and write out (to ldxot) the header cards for the sequence definition section.
 
  write (ldxot,'(A)') cardi4
 
  card=cardi2
  card( 3:10)='MNEMONIC'
  card(14:21)='SEQUENCE'
  write (ldxot,'(A)') card
 
  card=cardi4
  card(12:12)='|'
  write (ldxot,'(A)') card
 
  ! Create and write out (to ldxot) the Table D sequence definition cards.
 
  write (ldxot,'(A)') cardi2
 
  do n=1,ntbd(lun)
    if(.not.tdskip(tabd(n,lun)(1:6))) then
      card=cardi2
      card( 3:10)=tabd(n,lun)( 7:14)
      ic = 14
      ! Get the list of child mnemonics for this Table D descriptor, and then add each one (including any replication tags)
      ! to the sequence definition card for this Table D descriptor.
      call nemtbd(lun,n,nseq,nem(1,1),irp(1,1),krp(1,1))
      if(nseq>0) then
        do nc=1,nseq
          cmstr=' '
          icms=0
          call strsuc(nem(nc,1),wrk2,nch)
          if(irp(nc,1)/=0) then
            ! Add the opening replication tag.
            icms=icms+1
            cmstr(icms:icms)=reps(irp(nc,1))
          end if
          cmstr(icms+1:icms+nch)=wrk2(1:nch)
          icms=icms+nch
          if(irp(nc,1)/=0) then
            ! Add the closing replication tag.
            icms=icms+1
            cmstr(icms:icms)=reps(irp(nc,1)+5)
          end if
          if(krp(nc,1)/=0) then
            ! Add the fixed replication count.
            wrk1=' '
            write (wrk1,'(I3)') krp(nc,1)
            call strsuc(wrk1,wrk2,nch)
            cmstr(icms+1:icms+nch)=wrk2(1:nch)
            icms=icms+nch
          end if
          ! Will this child (and its replication tags, if any) fit into the current sequence definition card?  If not, then
          ! write out (to ldxot) the current card and initialize a new one to hold this child.
          if(ic>(79-icms)) then
            write (ldxot,'(A)') card
            card=cardi2
            card( 3:10)=tabd(n,lun)( 7:14)
            ic = 14
          end if
          card(ic:ic+icms-1)=cmstr(1:icms)
          ! Note that we want to leave 2 blank spaces between each child within the sequence definition card, to improve
          ! readability.
          ic=ic+icms+2
        end do
        write (ldxot,'(A)') card
        write (ldxot,'(A)') cardi2
      end if
    end if
  end do
 
  ! Create and write out (to ldxot) the header cards for the element definition section.
 
  write (ldxot,'(A)') cardi4
 
  card=cardi3
  card( 3:10)='MNEMONIC'
  card(14:17)='SCAL'
  card(21:29)='REFERENCE'
  card(35:37)='BIT'
  card(41:45)='UNITS'
  write (ldxot,'(A)') card
 
  card=cardi4
  card(12:12)='|'
  card(19:19)='|'
  card(33:33)='|'
  card(39:39)='|'
  card(66:66)='|'
  write (ldxot,'(A)') card
 
  ! Create and write out (to ldxot) the Table B element definition cards.
 
  write (ldxot,'(A)') cardi3
 
  do n=1,ntbb(lun)
    if(.not.tbskip(tabb(n,lun)(1:6))) then
      card=cardi3
      card( 3:10)=tabb(n,lun)( 7:14)
      card(41:64)=tabb(n,lun)(71:94)
      ! Add the scale factor.
      call strsuc(tabb(n,lun)(96:98),wrk2,nch)
      card(17-nch+1:17)=wrk2
      if(tabb(n,lun)(95:95)=='-') card(17-nch:17-nch)='-'
      ! Add the reference value.
      call strsuc(tabb(n,lun)(100:109),wrk3,nch)
      card(31-nch+1:31)=wrk3
      if(tabb(n,lun)(99:99)=='-') card(31-nch:31-nch)='-'
      ! Add the bit width.
      call strsuc(tabb(n,lun)(110:112),wrk2,nch)
      card(37-nch+1:37)=wrk2
      write (ldxot,'(a)') card
    end if
  end do
 
  write (ldxot,'(A)') cardi3
 
  ! Create and write out (to ldxot) the closing card.
 
  card=cardi4
  card( 1: 1)='`'
  card(80:80)=''''
  write (ldxot,'(A)') card
 
  return
end subroutine dxdump
 
recursive subroutine getabdb(lunit,tabdb,itab,jtab)
 
  use bufrlib
 
  use modv_vars, only: im8b
 
  use moda_tababd
  use moda_nmikrp
 
  implicit none
 
  integer, intent(in) :: lunit, itab
  integer, intent(out) :: jtab
  integer my_lunit, my_itab, lun, il, im, i, j, k, nseq, bort_target_set
 
  character*128, intent(out) :: tabdb(*)
  character, allocatable :: ctabdb(:,:)
  character*8 nemo
 
  ! Check for I8 integers
 
  if(im8b) then
    im8b=.false.
    call x84(lunit,my_lunit,1)
    call x84(itab,my_itab,1)
    call getabdb(my_lunit,tabdb,my_itab,jtab)
    call x48(jtab,jtab,1)
    im8b=.true.
    return
  endif
 
  ! If we're catching bort errors, set a target return location if one doesn't already exist.
 
  if (bort_target_set() == 1) then
    allocate(ctabdb(128,itab))
    call catch_bort_getabdb_c(lunit,itab,ctabdb,jtab)
    do i = 1, jtab
      do j = 1, 128
        tabdb(i)(j:j) = ctabdb(j,i)
      enddo
    enddo
    deallocate(ctabdb)
    call bort_target_unset
    return
  endif
 
  jtab = 0
 
  ! Make sure the file is open
 
  call status(lunit,lun,il,im)
  if(il==0) return
 
  ! Write out the Table D entries for this file
 
  do i=1,ntbd(lun)
    nemo = tabd(i,lun)(7:14)
    call nemtbd(lun,i,nseq,nem(1,1),irp(1,1),krp(1,1))
    do j=1,nseq,10
      jtab = jtab+1
      if(jtab<=itab) then
        write(tabdb(jtab),fmt='(A,A8,10(1X,A10))') 'D ', nemo, (nem(k,1),k=j,min(j+9,nseq))
      endif
    enddo
  enddo
 
  ! Add the Table B entries
 
  do i=1,ntbb(lun)
    jtab = jtab+1
    if(jtab<=itab) then
      write(tabdb(jtab),fmt='(A,A8,1X,A42)') 'B ', tabb(i,lun)(7:14), tabb(i,lun)(71:112)
    endif
  enddo
 
  return
end subroutine getabdb