getgb1r.f

Go to the documentation of this file.

C> @file
C> @brief Reads and unpacks a grib message.
C> @author Mark Iredell @date 1995-10-31
 
C> Program history log:
C> - Mark Iredell 1995-10-31
C> - Chuang 2004-07-22 Add packing bit number nbitss in the argument
C> list because eta grib files need it to repack grib file.
C> @param[in] LUGB Integer unit of the unblocked grib data file.
C> @param[in] LSKIP Integer number of bytes to skip.
C> @param[in] LGRIB Integer number of bytes to read.
C> @param[out] KF Integer number of data points unpacked.
C> @param[out] KPDS Integer (200) unpacked pds parameters.
C> @param[out] KGDS Integer (200) unpacked gds parameters.
C> @param[out] KENS Integer (200) unpacked ensemble pds parms.
C> @param[out] LB Logical*1 (kf) unpacked bitmap if present.
C> @param[out] F Real (kf) unpacked data.
C> @param[out] NBITSS Packaging bit number. Used by GRIB file to repack.
C> @param[out] IRET Integer return code.
C> - 0 All ok.
C> - 97 Error reading grib file.
C> - other w3fi63 grib unpacker return code.
C>
C> @note There is no protection against unpacking too much data.
C> Subprogram can be called from a multiprocessing environment.
C> Do not engage the same logical unit from more than one processor.
C> This subprogram is intended for private use by getgb routines only.
C>
C> @author Mark Iredell @date 1995-10-31
 
C-----------------------------------------------------------------------
      SUBROUTINE getgb1r(LUGB,LSKIP,LGRIB,KF,KPDS,KGDS,KENS,LB,F,NBITSS
     +      ,IRET)
      INTEGER KPDS(200),KGDS(200),KENS(200)
      LOGICAL*1 LB(*)
      REAL F(*)
      INTEGER KPTR(200)
      CHARACTER GRIB(LGRIB)*1
C - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C  READ GRIB RECORD
      CALL baread(lugb,lskip,lgrib,lread,grib)
C - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C  UNPACK GRIB RECORD
      IF(lread.EQ.lgrib) THEN
        CALL w3fi63(grib,kpds,kgds,lb,f,kptr,iret)
        IF(iret.EQ.0.AND.kpds(23).EQ.2) THEN
          CALL pdseup(kens,kprob,xprob,kclust,kmembr,45,grib(9))
        ENDIF
      ELSE
        iret=97
      ENDIF
      nbitss=kptr(20)
C - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C  RETURN NUMBER OF POINTS
      IF(iret.EQ.0) THEN
        kf=kptr(10)
      ELSE
        kf=0
      ENDIF
C - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
      RETURN
      END