home *** CD-ROM | disk | FTP | other *** search
/ Graphics Plus / Graphics Plus.iso / general / hdf / unix / examples.lha / examples / ann / get1anF.f < prev    next >
Encoding:
Text File  |  1991-10-08  |  2.4 KB  |  85 lines
  1. C
  2. C $Header: /pita/work/HDF/dev/RCS/test/annotations/get1anF.f,v 1.1 90/06/26 16:49:29 mfolk beta $
  3. C
  4. C $Log:    get1anF.f,v $
  5. c Revision 1.1  90/06/26  16:49:29  mfolk
  6. c Initial revision
  7. c 
  8. C
  9.       program get1an
  10.  
  11. C Program to test routines for finding and reading label and description
  12. C
  13. C An HDF file called 'o2' must already exist with at least two SDSs
  14. C in it and both a label and description annotation corresponding 
  15. C to the second SDS.
  16. C
  17. C Mike Folk  9/12/90
  18. C
  19. C****||************************************************************
  20.  
  21.       integer DFopen, DFclose
  22.       integer dfindnr, daglab, dagdesc, dagdlen
  23.       integer desclen, ref, ret, dfile
  24.       integer DFACC_READ, DFTAG_SDG
  25.       character*20  label, filename
  26.       character*400 desc
  27.     
  28.       parameter (DFACC_READ  =   1,
  29.      $           DFTAG_SDG   = 700 )
  30.  
  31.       print *, 'Enter HDF file name.  The HDF file should contain two'
  32.       print *, 'SDSs with a label and description for the second one.'
  33.       print *, 'File name: '
  34.       read *, filename
  35.       
  36.       dfile = DFopen(filename, DFACC_READ, -1)
  37.   
  38. C****||******** find ref of second SDS in file ****************
  39.       ref = 0
  40.       ref = dfindnr(dfile, DFTAG_SDG, ref)
  41.       ref = dfindnr(dfile, DFTAG_SDG, ref)
  42.       if (ref .lt. 0) 
  43.      *  call fatalerror('Unable to find second scientific data set.')
  44.  
  45.       ret = DFclose(dfile)
  46.   
  47. C****||******** get label, then description ****************
  48.       ret = daglab(filename, DFTAG_SDG, ref, label, 11)
  49.       if (ret .lt. 0) print *, 'No label'
  50.       if (ret .ge. 0) print *,'Label: ', label
  51.  
  52.       desclen = dagdlen(filename, DFTAG_SDG, ref)
  53.       if (desclen .gt. 400) 
  54.      $    call fatalerror('Description too long. Greater than 400.')
  55.       if (desclen .lt. 0) 
  56.      $    call fatalerror('Error reading description length.')
  57.  
  58.         ret = dagdesc(filename, DFTAG_SDG, ref, desc, desclen)
  59.       if (ret .lt. 0) 
  60.      $    call fatalerror('Error reading description.')
  61.       print *,'Description: '
  62.       print *, desc
  63.       print *
  64.       print *, "+++++++++++++++++++++++++"
  65.       print *
  66.       print *
  67.       stop
  68.       end
  69.  
  70. C************************************************************
  71. * fatalerror: subroutine to report fatal error and abort
  72. *
  73. C****||***********************************************************
  74.  
  75.       subroutine fatalerror(s)
  76.       character*(*) s
  77.  
  78.       print *, s
  79.       print *, 'DFerror:', DFerrno()
  80.       print *, 'Program aborted.'
  81.       print *, ' '
  82.       stop
  83.       end
  84.  
  85.