USGS InfoBank program -- whnav2daps

whnav2daps

  Program to reformat navigation files from 
  CMG Woods Hole to standard DAPS format.
    

  Clint Steele      11/15/2007


    

 open (unit = 90,
*      access = 'sequential',
*      form = 'formatted',
*      status = 'unknown',
*      name = InFile(1:InFileLength))

 open (unit = 91, 
*      access = 'sequential',
*      form = 'formatted',
*      status = 'unknown',
*      name = OutFile(1:OutFileLength))

    

 read (5, '(q,a)') InFileLength, InFile(1:InFileLength)
 read (5, '(q,a)') OutFileLength, OutFile(1:OutFileLength)
 read (5, '(i4)') Year
 read (90, '(q,a)', end = 99) ARecLength, ARec(1:ARecLength)
         read (ARec(1:LatEnd-1), '(f)') Latitude
         read (ARec(LatEnd+1:LongEnd-1), '(f)') Longitude
         read (ARec(1:LatDegEnd-1), '(i)') LatDeg
         read (ARec(LatDegEnd+1:LatMinEnd), '(f)') LatMin
         read (ARec(LatMinEnd:LongDegEnd-1), '(i)') LongDeg
         read (ARec(LongDegEnd+1:LongMinEnd-1), '(f)') LongMin
         read (ARec(1:LatEnd), '(f)') Latitude
         read (ARec(LongStart:LongEnd), '(f)') Longitude
         read (ARec(1:LatEnd), '(f)') Latitude
         read (ARec(LongStart:LongEnd-1), '(f)') Longitude
             read (ARec(ShotStart+1:ARecLength), '(i)') ShotNumber
     read (ARec(Tab2+1:ARecLength), '(i)') ShotNumber
    

 write (ARec, '(1024x)')
     write (91, '(a)') '! '//ARec(1:ARecLength)
     write (Line, '(1024x)')
 write (91, 2) Year, Day, Hour, Minute, Second, Tenths,
*     Latitude, Longitude, Line(1:LineLength), 
*     ShotNumber, Comment(1:CommentLength)
    

 call stringlen (ARec, ARecLength)
         call stringlen (Line, LineLength)