PROGRAM ssmis
! ---------------------------------------------------------
!
!
! Program to decode the ssmi/su tank
      IMPLICIT none

      REAL bmiss
      PARAMETER (bmiss=10.0E6)
      CHARACTER(8) SUBSET
      REAL(8) hdr(9),ident(3), ssmischn(4,24)
!
      INTEGER ichan, jj, ii, isiid, ntank, lubfi, kt, jt
      INTEGER iorbn, islnm, isaid, idate, iret
      INTEGER mhr, mins, secs, mda, myr, mmo
      INTEGER icount, err, freq
      REAL xlat, xlon

      INTEGER openout, outunit, creturns, textout
!.................................

        CALL W3TAGB('SSMISU_Decode',0095,0333,0077,'NP11   ')

        CALL DATELEN(10)

        lubfi = 11
        icount = 0
        kt = 0
        ntank = 1

      freq = 10 * 1000 !Text output frequency
      textout = 52

      outunit = 51
      creturns = openout(outunit)
      PRINT *,'ssmisubufr creturned = ',creturns


      DO JT = 1, ntank
!-----------------------------------------------------------------------
!                       DECODE BUFR FILE
!-----------------------------------------------------------------------
!
!  OPEN THE BUFR INPUT FILE
!  ------------------------
!
      CALL OPENBF(lubfi,'IN',lubfi)
      CALL READMG(lubfi,SUBSET,IDATE,IRET)
      IF(IRET.NE.0) GOTO 100
      PRINT*,'READING DATA FOR ',IDATE

!  READ A SUBSET - READ A MESSAGE WHEN NO SUBSETS - END WHEN NO MESSAGES
!  ---------------------------------------------------------------------
  150 CONTINUE
      CALL READSB(lubfi,IRET)
!          print*,' readsb iret ',iret
      IF(IRET.NE.0) THEN
         CALL READMG(lubfi,SUBSET,IDATE,IRET)
!          print*,' readms iret ',iret
         IF(IRET.NE.0) GOTO 100
         GOTO 150
      ENDIF

!  -----------------------------------------------------------------
!  CALL UFBINT, UFBREP TO GET THE DATA
!
!     MNEMONICS are embedded in CALLS
!
!  -----------------------------------------------
      !For testing:
      ! IF (icount .GT. 100*freq) GO TO 100

      call ufbint(lubfi,hdr,9,1,iret, &
                  'SAID YEAR MNTH DAYS HOUR MINU SECO CLAT CLON')
!     WRITE(6,*) 'IRET= ',IRET
      call ufbint(lubfi,ident,3,1,iret,'SAID ORBN SLNM')

      isaid = NINT(ident(1))
      iorbn = NINT(ident(2))
      islnm = NINT(ident(3))

      isiid = NINT(hdr(1))
      myr   = NINT(hdr(2))
      mmo   = NINT(hdr(3))
      mda   = NINT(hdr(4))
      mhr   = NINT(hdr(5))
      mins  = NINT(hdr(6))
      secs  = NINT(hdr(7))
      xlat  = hdr(8)
      xlon  = hdr(9)

!    "SSMISCHN"24  =  24 btmps
!    each one is "CHNM  TMBR  WTCA   CTCA "
!      TCA are target calibrations, warm and cold

      call ufbseq(lubfi,ssmischn,4,24,iret,'SSMISCHN')

               icount = icount + 1

      if(mod(icount,freq) .eq. 0) then

        write(textout,50)isiid,myr,mmo,mda,mhr,mins,secs, xlat, &
          xlon,isaid,islnm,iorbn
   50   format(1x,i3,1x,i4,1x,5i2,1x,2f12.5,1x,2i6,1x,i8)

        do ii = 1,4
         do jj = 1,24
          if(ssmischn(ii,jj).ge.bmiss)ssmischn(ii,jj)= -999.99
         enddo
        enddo

        do jj = 1, 24
          if(ssmischn(1,jj).le.-999.99) then
              ichan = -999
            else
             ichan = NINT(ssmischn(1,jj))
          endif 
          write(textout,502)ichan,ssmischn(2,jj)
        enddo
  502   format('chan # ',i4,' tmbr ',f8.2)

      endif

! Output in binary all steps' info:
          CALL ssmisout(hdr, ident, ssmischn)

!
!       get next report
!
      GO TO 150
  100 CONTINUE

         call   closbf(lubfi)

          if(icount.lt.1) then
           WRITE(6,*) 'READMG ERROR W/ BUFR MSG:IRET=',IRET
           err = 23
!            call exit(23)
            STOP
          endif
!
!
        WRITE(6,*) 'AFTER WOOLENS INTERFACE:  icount= ',icount
!
      WRITE(6,617)
  617 FORMAT(1x,'  TIMES AND TOTAL NUMBER OF btmps OVER THE AREA')
!D     WRITE(6,618) IDAT1,IDAT2,icount
!D  618 FORMAT(2X,3I10)
!
      ENDDO
!
!.................
      WRITE(6,692)icount
  692 FORMAT(1x,'END OF JOB','count= ',i8)

        CALL W3TAGE('DCODNET3')

      END