makgds.f90

 
SUBROUTINE makgds(IOPT,KGDS,GDS,LENGDS,IRET)
  IMPLICIT NONE
  !
  CHARACTER,          INTENT(INOUT) :: GDS(400)
  !
  INTEGER,            INTENT(IN   ) :: IOPT
  INTEGER,            INTENT(INOUT) :: KGDS(200)
  INTEGER,            INTENT(  OUT) :: IRET, LENGDS
  !
  INTEGER                           :: ICOMP, IPDS(200), IGDS(200)
  INTEGER                           :: KPTR(200), KPDS(200), NPTS
  !
  DATA kptr/200*0/, kpds/200*0/
  ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
  !  UNPACK GDS INTO KGDS
  IF(iopt.EQ.-1) THEN
     CALL fi633(gds,kptr,kgds,iret)
     ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
     !  USE KGDS TO PACK GDS
  ELSEIF(iopt.EQ.255) THEN
     CALL r63w72(kpds,kgds,ipds,igds)
     icomp=mod(igds(8)/8,2)
     CALL w3fi74(igds,icomp,gds,lengds,npts,iret)
     ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
     !  USE NCEP GRID ID TO MAKE GDS AND KGDS
  ELSEIF(iopt.GT.0.AND.iopt.LT.255) THEN
     CALL w3fi71(iopt,igds,iret)
     IF(iret.EQ.0) THEN
        icomp=mod(igds(8)/8,2)
        CALL w3fi74(igds,icomp,gds,lengds,npts,iret)
        IF(iret.EQ.0) CALL fi633(gds,kptr,kgds,iret)
     ENDIF
  ENDIF
  ! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
END SUBROUTINE makgds