stbfdx.f

Go to the documentation of this file.

C> @file
C> @brief Store a DX BUFR tables message into internal arrays.
C>
C> @author J Ator @date 2009-03-23
 
C> This subroutine copies a DX BUFR tables message
C> from the input array mesg into the internal memory arrays in
C> module @ref moda_tababd.
C>
C> @param[in] LUN - integer: I/O stream index into internal memory arrays.
C> @param[in] MESG - integer(*): DX BUFR tables message
C>
C> @author J. Ator @date 2009-03-23
      SUBROUTINE stbfdx(LUN,MESG)
 
      USE modv_maxcd
      USE moda_tababd
 
      COMMON /dxtab / maxdx,idxv,nxstr(10),ldxa(10),ldxb(10),ldxd(10),
     .                ld30(10),dxstr(10)
 
      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
      dimension     ldxbd(10),ldxbe(10)
 
      dimension     mesg(*)
 
      DATA ldxbd /38,70,8*0/
      DATA ldxbe /42,42,8*0/
 
C-----------------------------------------------------------------------
      ja(i) = ia+1+lda*(i-1)
      jb(i) = ib+1+ldb*(i-1)
C-----------------------------------------------------------------------
 
C  GET SOME PRELIMINARY INFORMATION FROM THE MESSAGE
C  -------------------------------------------------
 
      idxs = iupbs01(mesg,'MSBT')+1
      IF(idxs.GT.idxv+1) idxs = iupbs01(mesg,'MTVL')+1
      IF(ldxa(idxs).EQ.0) GOTO 901
      IF(ldxb(idxs).EQ.0) GOTO 901
      IF(ldxd(idxs).EQ.0) GOTO 901
 
      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.NE.dxstr(idxs)) GOTO 902
 
C  SECTION 4 - READ DEFINITIONS FOR TABLES A, B AND D
C  --------------------------------------------------
 
      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)
 
C  TABLE A
C  -------
 
      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
 
C  TABLE B
C  -------
 
      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
 
C  TABLE D
C  -------
 
      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.GT.maxcd) GOTO 903
        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.LT.0) GOTO 904
        ENDDO
        id = id+ldd+1 + nd*l30
        IF(iupb(mesg,id+1,8).EQ.0) id = id+1
        ntbd(lun) = n
      ENDDO
 
C  EXITS
C  -----
 
      RETURN
901   CALL bort('BUFRLIB: STBFDX - UNEXPECTED DICTIONARY MESSAGE '//
     . 'SUBTYPE OR LOCAL VERSION NUMBER (E.G., L.V.N. HIGHER THAN '//
     . 'KNOWN)')
902   CALL bort('BUFRLIB: STBFDX - UNEXPECTED DICTIONARY MESSAGE '//
     . 'CONTENTS')
903   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)
904   CALL bort('BUFRLIB: STBFDX - BAD RETURN FROM BUFRLIB ROUTINE '//
     . 'PKTDD, SEE PREVIOUS WARNING MESSAGE')
      END