C$$$  sub  PROGRAM DOCUMENTATION BLOCK
C                .      .    .                                       .
C sub  PROGRAM:  gribst      UNPACK GLOBAL Grib SST
C   PRGMMR: GERALD           ORG: W/NMC21     DATE: 95-08-01
C
C ABSTRACT: DECODES GLOBAL AVN, FNL, OR MRF GRIB GRIDDED FIELDS.
C
C PROGRAM HISTORY LOG:
C   95-08-01  GERALD
C   96-07-23  PETERS - UPDATED FOR USE WITH ETA MODEL, GETS 
C                      REYNOLDS 1X1 OISST 
Cmp
Cmp	Changed 99-02-22 by Pyle to use a normal version of GETGB
C
C	CHANGED 01-02-16 by PYLE to use 0.5 X 0.5 analysis
Cmp
C
C
C USAGE:
C   INPUT FILES:
C     DDNAME1  - GENERIC NAMES & USAGE
C     FTNNF001 - START IN COL 7
C              - BUT TAB CONTINUATIONS TO LINE UP LIKE THIS
C     FXN      - LINE UP DASHES WHEN DESCRIBING FILES
C     TAPENN   - CRAY TYPE DESIGNATION
C     UNITNN   - CRAY TYPE DESIGNATION
C     PARM     - IF PARM FIELD IS USED, DESCRIBE HERE
C
C   OUTPUT FILES:  (INCLUDING SCRATCH FILES)
C     FTMMF001 - NAMES & USAGE AS ABOVE IN THE INPUT SECTION
C     FT06F001 - INCLUDE IF ANY PRINTOUT
C
C   SUBPROGRAMS CALLED: (LIST ALL CALLED FROM ANYWHERE IN CODES)
C     UNIQUE:    - ROUTINES THAT ACCOMPANY SOURCE FOR COMPILE
C     LIBRARY:
C       COMMON   - LIST COMMON LIBRARY ROUTINES, E.G., CONSOL
C       W3LIB    -
C       W4LIB    - DELETE THE CORRESPONDING LINE OR LINES
C       GRAPHICS - IF LIBRARY IS UNNEEDED
C
C   EXIT STATES:
C     COND =   0 - SUCCESSFUL RUN
C          =NNNN - TROUBLE OR SPECIAL FLAG - SPECIFY NATURE
C
C REMARKS: LIST CAVEATS, OTHER HELPFUL HINTS OR INFORMATION
C
C ATTRIBUTES:
C   LANGUAGE: INDICATE EXTENSIONS, COMPILER OPTIONS
C   MACHINE:  CRAY
C
C$$$
C
      subroutine gribst(insst,indxst,gsst,ierr)
c
      Parameter(jf=720*360)
      INTEGER PDS,GDS,GRID
      INTEGER JPDS(200),JGDS(200),IGRD(5,3)
      INTEGER KPDS(200),KGDS(200)
      REAl FLD(jf),out(720,360)
      dimension gsst(721,360)
      logical lb(jf)
      equivalence(out(1,1),fld(1))
c
C
C                sst
      DATA IGRD/ 11, 34,  2, 52, 11,
     *            1,105,102,105,105,
     *            0, 10,  0,  2,  2/
c
C            INPUT UNITS FOR DECODING GRIB FILE
C
          LUGB=INSST
Cmp          LUGI=INDXST
	LUGI=0
          ierr = 0
          j = 0
C
C........    DECODE THE FIELDS
C
        DO 30 GRID = 1,1
C
          DO 10 PDS=1,200
          JPDS(PDS) = -1
   10      CONTINUE
C
             DO 20 GDS = 1,200
           JGDS(GDS) = -1
   20     CONTINUE
C
C........   GET SELECT FIELDS
C
           jPDS(5) = IGRD(GRID,1)
           jPDS(6) = IGRD(GRID,2)
           jPDS(7) = IGRD(GRID,3)
C
C             print *,'call getgb with unit #s ', LUGB,LUGI
	   CALL BAOPEN(LUGB,'hires_sst',IRETBA)
 	write(6,*) 'back from baopen, IRET= ', IRETBA
           CALL GETGB(LUGB,LUGI,JF,j,JPDS,JGDS,
     *                          KF,K,KPDS,KGDS,LB,Fld,IRET)
C
          IF(IRET.NE.0) THEN
            ierr = 1
            WRITE(6,60)IRET
            GO TO 999
          ENDIF
   60  FORMAT(1X,' IRET =',I5)
C
c..   flip grid to PT(1,1) =(0e,-90.0)
c
	write(6,*) 'flipping the data'
              do  jj= 1,360
                do  kk= 1,720
                 gsst(kk,360-jj+1) = out(kk,jj)
                end do
              end do
c
c...   add greenich to right side of grid
c
             do jj = 1,360
               gsst(721,jj) = gsst(1,jj)
             end do
c
Cmp            WRITE(6,61)KPDS,KF,KGDS
Cmp   61  FORMAT(2(/,2X,'PDS=',13I7),2(/,2X,' GDS=',11I7 ))
	write(6,*) 'KPDS= ', (KPDS(I),I=1,12)
	write(6,*) 'KGDS= ', (KGDS(I),I=1,25)
C
   30        CONTINUE
C
  999           CONTINUE
	write(6,*) 'leaving GRIBST'
             return
             END
