C     TEST PROGRAM FOR DBGET4
C     DEMONSTRATES FORTRAN INTERFACE TO TSP DATABANKS
C
      PARAMETER (IPDATA=10000,IPDOC=1000)
      IMPLICIT INTEGER (A-Z)
      CHARACTER*128 LNAME
      CHARACTER*8 CNAME
      INTEGER*2 ECODE
      COMMON /FILINF/ PUNIT
      REAL DATA(IPDATA)
      INTEGER DOC(IPDOC)
C
      PUNIT = 1
      LNAME='DBTEST.TLB'
      CALL DBOPNR(LNAME,ECODE)
      IF (ECODE.NE.0) THEN
        WRITE(*,*) ' ECODE from DBOPNR is not zero: ',ECODE
        GOTO 900
      ENDIF
      CALL DBVERG(POS,IVERS)
      IF (IVERS.NE.2) THEN
        WRITE(*,*) ' Incompatible databank version number: ',IVERS
        GOTO 900
      ENDIF
C
      CNAME = 'X'
      LENGTH = IPDATA
      LDOC = IPDOC
      CALL DBGET4(CNAME,LENGTH,DATA,
     *                  YEAR,PERIOD,FREQ,DOC,LDOC,ECODE)
      IF (ECODE.NE.0) THEN
        WRITE(*,*) ' ECODE from DBGET4 is not zero: ',ECODE
        GOTO 900
      ENDIF
      WRITE(*,100) CNAME,LENGTH,YEAR,PERIOD,FREQ,LDOC
100   FORMAT(' series=[',A,'], length=',I6,' starting ',I4,':',I2,
     * ' freq=',I2,' length of Doc=',I4)
      IF (LDOC.NE.0) WRITE(*,101) (DOC(I),I=1,LDOC)
101   FORMAT(1X,16A4)
      WRITE(*,*) ' data=',(DATA(I),I=1,LENGTH)
C
900   STOP
      END