readtest()

subroutine readtest

(

character*(*)

filename

)

Definition at line 190 of file fapidtest.f.

Referenced by fapidtest().

c
      character*(*) filename
      integer lu
      parameter(lu=10)
c time series
      integer maxsamp, nsamp
      parameter(maxsamp=100000)
      real fdata(maxsamp)
      integer idata(maxsamp)
      equivalence(fdata, idata)

c sff file header
      character code*10, timestamp*13, scs*1, date*6, time*10
      real sffversion, sc1, sc2, sc3
      character*20 source
c
c sff trace 
      logical last
      character rcs*1, wid2line*132
      real rc1, rc2, rc3, tanf, dt
      integer nstack
c 
      print *,'read test'
      print *,'---------'
      print *,'open ',filename(1:index(filename,' '))
      call sff_ropens(lu, filename, sffversion, timestamp, code, 
     &  source, scs, sc1, sc2, sc3, date, time, ierr)
      if (ierr.ne.0) then
        stop 'ERROR: opening file'
      endif
      print *,'sffversion: ', sffversion
      print *,'timestamp: ', timestamp
      print *,'code: ', code
      print *,'source: ', source
      print *,'scs, sc1, sc2, sc3: ', scs, sc1, sc2, sc3
      print *,'date time: ', date, ' ', time
      last=.false.
      do while (.not.last)
        print *,' '
        print *,'read next trace'
        nsamp=maxsamp
        call sff_rtracei(lu, tanf, dt,
     &    wid2line, nsamp, fdata,
     &    idata, code, last,
     &    rcs, rc1, rc2, rc3, nstack, ierr)
        print *,wid2line
        print *,'nsamp: ',nsamp
        print *,'rcs, rc1, rc2, rc3: ',rcs, rc1, rc2, rc3
        print *,'nstack: ',nstack
        print *,'dt: ',dt
        print *,'tanf: ',tanf
        print *,'code: ',code
        if (ierr.ne.0) then
          stop 'ERROR: reading trace'
        endif
      enddo
      return

Here is the caller graph for this function: