program trmmpr_vdata
  implicit none
!  
! to use read v-data
!
  integer, parameter :: READ_MODE = 1
  character :: vdata_name*128,fields*1024
  integer(4) :: hopen, hclose, vfstart, vsfgid, vsfatch
  integer(4) :: vsfinq, vsfdtch, vsfsfld, vsfread, vsffnd
  integer(4) :: file_id, vdata_id, vdata_ref
  integer(4) :: num_rec, n_scan, interlace, vdata_size
  character*(100) :: hdffile, geo_name

  integer :: istat, i

  real(8), allocatable, dimension (:) :: scan_time

  call getarg(1, hdffile)
  if(iargc() /= 1) call exit(1)

  file_id = hopen(hdffile, READ_MODE, 0)
!
  istat = vfstart(file_id)
!
  vdata_ref = vsffnd(file_id, "scan_time")
  vdata_id  = vsfatch(file_id, vdata_ref, "r")
  istat = vsfinq(vdata_id, num_rec, interlace, fields, vdata_size, vdata_name)
  n_scan = num_rec
  allocate(scan_time(n_scan))
  istat = vsfread(vdata_id, scan_time, n_scan, interlace)
! 
  istat = hclose(file_id)
!
  do i = 1, n_scan
    print "(i6,3x,2f12.4)", i, real(scan_time(i)), real(scan_time(i)/(60*60))
  enddo

end program trmmpr_vdata