USGS InfoBank program -- incgroups2kmlhelper

incgroups2kmlhelper

  Subroutine to help mainincgroups2kml to read the equipment and

  calculate the north south east west lat and long
  call getRange to calculate the range
       

  Sumy Tom        4/17/2008


    

	  read( SampsName
     * (FirstDot + 1 : LastDot - 1), '(q,a)')
     * EqTypeLength, EqTypeName
 3    read(99,'(q,a)', end = 4) 
     * ALineLength, ALine(1:ALineLength)
         read (ALine, '(t25,f9.5,t35,f10.5)', err=3) Lat, Long
         read (ALine, '(t25,f10.6,t36,f11.6)', err=3) Lat, Long

      write (82, '(6x,a)') 
     * '<name>'//EqTypeName(1:EqTypeLength) //'<
      write (82, '(a)') '  <LookA
      write (82, '(a,f11.6,a)') 
     *  '<longitude>', Long , '</lo
      write (82, '(a,f10.6,a)') '  <latitude>',Lat,'</l
      write (82, '(a,f,a)') '    <range>', Range, '<
      write (82, '(a)') '    <tilt>0&l
      write (82, '(a)') '    <heading>0</
      write (82, '(a)') '  </LookA

         open (unit = 82,
     *   access = 'sequential',
     *   form   =  'formatted',
     *   status =  'unknown',
     *   name = '/infobank/tmp/incgroups2kml.eqfolder.txt')
       open (unit = 99,
     * access = 'sequential',
     * form   =  'formatted',
     * status =  'unknown',
     * name = '/infobank/tmp/incgroups2kml.latlong2.txt')

      call system('/bin/cat '//
     * '/infobank/?/*/kml/*."'//
     *  EqTypeName(1:EqTypeLength)//
     * '".sampids.kmlinc > 
     * '/infobank/tmp/incgroups2kml.latlong1.txt')
      call system('/bin/grep Lat/Long '//
     * '/infobank/tmp/incgroups2kml.latlong1.txt '//
     * '> 
     * '/infobank/tmp/incgroups2kml.latlong2.txt' )   
      call system ('/bin/rm /infobank/tmp/incgroups2kml.latlong1.txt')
      call system ('/bin/rm /infobank/tmp/incgroups2kml.latlong2.txt')
        call  getrange( Lat, Long,
     *  North,South, East ,West, Range)