c-----------------------------------------------------------------------
c
c     Test driver for    gf3_read_hdf_f77()             
c                        gf3_read_brief_hdf_f77()       and
c                        gf3_read_shapes_hdf_f77()
c
c----------------------------------------------------------------------
c     Author:       Matthew W. Choptuik
c     Institution:  The University of Texas at Austin
c     Date:         June 1994
c----------------------------------------------------------------------

      program      r0

      implicit     none

      integer      gf3_read_hdf_f77, gf3_read_brief_hdf_f77,
     &             gf3_read_shapes_hdf_f77

      integer      rc
   
      integer      maxn
      parameter  ( maxn = 10 000 )

      character*64 gf3_name
      integer      gf3_dshape(3),       gf3_cshape(3)
      real*8       gf3_c1(maxn),        gf3_c2(maxn),  
     &             gf3_c3(maxn)
      real*8       gf3_data(maxn)

      integer      nfile 
      parameter  ( nfile = 3 )
      character*8  fname(nfile)
      integer      i

      fname(1) = 'foo0.hdf'
      fname(2) = 'foo1.hdf'
      fname(3) = 'foo2.hdf'

      do i = 1 , nfile
         rc = gf3_read_hdf_f77(fname(i),gf3_name,gf3_cshape,
     &           gf3_c1,gf3_c2,gf3_c3,gf3_dshape,gf3_data)

         if( rc .eq. 0 ) then
            write(*,*) 'gf3_read_hdf_f77('//fname(i)
     &                 //',...)  succeeds'
            call ivdump(gf3_cshape,3,'cshape',6)
            call dvdump(gf3_c1,gf3_cshape(1),'c_1',6)
            call dvdump(gf3_c2,gf3_cshape(2),'c_2',6)
            call dvdump(gf3_c3,gf3_cshape(3),'c_3',6)
            call ivdump(gf3_dshape,3,'dshape',6)
            call gf3_dump(gf3_dshape,gf3_data,'data',6)
         else 
            write(*,*) 'gf3_read_hdf_f77('//fname(i)
     &                 //',...)  fails'
         end if
      end do

      rc = gf3_read_brief_hdf_f77('foo2.hdf',gf3_dshape,gf3_data)
      if( rc .eq. 0 ) then
         write(*,*) 'gf3_read_brief_hdf_f77(foo2.hdf,...) succeeds'
         call ivdump(gf3_dshape,3,'dshape',6)
         call gf3_dump(gf3_dshape,gf3_data,'data',6)
      else
         write(*,*) 'gf3_read_brief_hdf_f77(foo2.hdf,...) fails'
      end if

      rc = gf3_read_brief_hdf_f77('foo0.hdf',gf3_dshape,gf3_data)
      if( rc .eq. 0 ) then
         write(*,*) 'gf3_read_brief_hdf_f77(foo0.hdf,...) succeeds'
         call ivdump(gf3_dshape,3,'dshape',6)
         call gf3_dump(gf3_dshape,gf3_data,'data',6)
      else
         write(*,*) 'gf3_read_brief_hdf_f77(foo0.hdf,...) fails'
      end if

      rc = gf3_read_shapes_hdf_f77('foo0.hdf',gf3_cshape,gf3_dshape)
      call ivdump(gf3_cshape,3,'foo0 cshape',6)
      call ivdump(gf3_dshape,3,'foo0 dshape',6)
      rc = gf3_read_shapes_hdf_f77('foo1.hdf',gf3_cshape,gf3_dshape)
      call ivdump(gf3_cshape,3,'foo1 cshape',6)
      call ivdump(gf3_dshape,3,'foo1 dshape',6)
      rc = gf3_read_shapes_hdf_f77('foo2.hdf',gf3_cshape,gf3_dshape)
      call ivdump(gf3_cshape,3,'foo2 cshape',6)
      call ivdump(gf3_dshape,3,'foo2 dshape',6)

      stop

      end