USGS InfoBank program -- beenthere

beenthere

  Program to determine which Activity IDs operated within a

  user-specified n/s/e/w area.
       

  Clint Steele     10/17/2006


    

      read (5, '(f)', end = 99) North
      read (5, '(f)', end = 99) South
      read (5, '(f)', end = 99) East
      read (5, '(f)', end = 99) West
      read (5, '(q,a)', end = 99) ANameLength, AName(1:ANameLength)
      read (CombinedGrids, '(i3,1x,i4,1x,q,a)', end = 98)
     *    Lat, Long, IDLength, ID(1:IDLength)
      read (UniqueIDs, '(q,a)', end=95) IDLength, ID(1:IDLength)
          read (BestFile, '(q,a)', end = 96) 
     *        InFileLength, InFile(1:InFileLength)
          read (NavFile, '(q,a)', end = 97) 
     *       ALineLength, ALine(1:ALineLength)
          read (ALine(1:ALineLength), 101) Latitude, Longitude

      write (WindowedIDs, '(a)') ID(1:IDLength)
      write (AreaKml, '(a)') '<?xml version="1.0" encoding="UTF-8"
      write (AreaKml, '(a)') 
     *     '<kml xmlns="http://earth.google.com/kml/2.0
      write (AreaKml, '(a)') '<Documen
      write (AreaKml, '(a)') '  <name>'//AName(1:ANameLength)//'&l
      write (AreaKml, '(a)') '  <ScreenOverla
      write (AreaKml, '(a)') '    <Ico
      write (AreaKml, '(a)') 
     *    '      <href>http://walrus.wr.usgs.gov/infoba
     *    'programs/html/tools/SmUSGS.png</hre
      write (AreaKml, '(a)') '    </Ico
      write (AreaKml, '(a)') '    <name>USGS Visual Identity&l
      write (AreaKml, '(a)') 
     *    '    <overlayXY x="0" y="1" xunits="fraction" 
     *    'yunits="fraction"/&g
      write (AreaKml, '(a)') 
     *    '    <screenXY x="0" y="1" xunits="fraction" 
     *    'yunits="fraction"/&g
      write (AreaKml, '(a)') '  </ScreenOverla
      write (AreaKml, '(a)') '<Placemar
      write (AreaKml, '(a)') '  <descriptio
      write (AreaKml, '(a)') '    <![CDAT
      write (AreaKml, '(a)') 
     *    '          <li><a href="http://www.googl
     *    'custom?q=usgs+br96cp+B-R-96-CP&sa=Google+Search'//
     *    '&cof=LW%3A600%3BL%3Ahttp%3A%2F%2Fwalrus.wr.usgs.gov'//
     *    '%2Finfobank%2Fprograms%2Fhtml%2Ftools%2FBanr600x72SF.jpg'//
     *    '%3BLH%3A72%3BAH%3Aleft%3BS%3Ahttp%3A%2F'//
     *    '%2Fwalrus.wr.usgs.gov%2Finfobank%3BAWFID'//
     *    '%3A6c82345480ae9708%3B&domains=walrus.wr.usgs.gov'//
     *    '&sitesearch=walrus.wr.usgs.gov">Search InfoBank&l
      write (AreaKml, '(a)') 
     *    '    <!-- Source: http://walrus.wr.usgs.gov/infobank/ --
      write (AreaKml, '(a)') '    ]]&g
      write (AreaKml, '(a)') '  </descriptio
      write (AreaKml, '(a)') '  <name>'//AName(1:ANameLength)//'&l
      write (AreaKml, '(a)') '  <visibility>1</vis
      write (AreaKml, '(a)') '  <open>0&l
      write (AreaKml, '(a)') '  <Styl
      write (AreaKml, '(a)') '    <LineStyl
      write (AreaKml, '(a)') '      <color>ff0000ff<
      write (AreaKml, '(a)') '      <width>3<
      write (AreaKml, '(a)') '    </LineStyl
      write (AreaKml, '(a)') '    <PolyStyl
      write (AreaKml, '(a)') '      <color>7f00ff00<
      write (AreaKml, '(a)') '    </PolyStyl
      write (AreaKml, '(a)') '  </Styl
      write (AreaKml, '(a)') '  <MultiGeometr
      write (AreaKml, '(a)') '    <LineStrin
      write (AreaKml, '(a)') '      <extrude>1</
      write (AreaKml, '(a)') '      <tessellate>1</tes
      write (AreaKml, '(a)') '      <coordinate
      write (AreaKml, '(a,f,a,f,a)') '        ',West,', ',North,', 0'
      write (AreaKml, '(a,f,a,f,a)') '        ',East,', ',North,', 0'
      write (AreaKml, '(a,f,a,f,a)') '        ',East,', ',South,', 0'
      write (AreaKml, '(a,f,a,f,a)') '        ',West,', ',South,', 0'
      write (AreaKml, '(a,f,a,f,a)') '        ',West,', ',North,', 0'
      write (AreaKml, '(a)') '      </coordinate
      write (AreaKml, '(a)') '    </LineStrin
      write (AreaKml, '(a)') '  </MultiGeometr
      write (AreaKml, '(a)') '</Placemar
      write (AreaKml, '(a)') '<Placemar
      write (AreaKml, '(a)') '<visibility>1</vis
      write (AreaKml, '(a)') '<name>'//AName(1:ANameLength)//'&l
      write (AreaKml, '(a)') '<Poin
      write (AreaKml, '(a)') '  <coordinate
      write (AreaKml, '(a,f,a,f,a)') '        ',West,', ',North,', 0'
      write (AreaKml, '(a)') '  </coordinate
      write (AreaKml, '(a)') '</Poin
      write (AreaKml, '(a)') '</Placemar
      write (AreaKml, '(a)') '</Documen
      write (AreaKml, '(a)') '</km

      open (unit = CombinedGrids,
     *      access = 'sequential',
     *      form = 'formatted',
     *      status = 'old',
     *      name = '/infobank/programs/dbase/grids/all.sorted')
      open (unit = WindowedIDs,
     *      access = 'sequential',
     *      form = 'formatted',
     *      status = 'unknown',
     *      name = '/infobank/tmp/beenthere.ids.temp')
      open (unit = UniqueIDs,
     *      access = 'sequential',
     *      form = 'formatted',
     *      status = 'unknown',
     *      name = '/infobank/tmp/beenthere.unique.temp')
      open (unit = BestFile,
     *      access = 'sequential',
     *      form = 'formatted',
     *      status = 'old',
     *      name = '/infobank/'//ID(1:1)//'/'//
     *             ID(1:IDLength)//'/nav/'//
     *             LongID(1:IDLength+3)//'.best_nav_file')
      open (unit = ResultsFile,
     *      access = 'sequential',
     *      form = 'formatted',
     *      status = 'unknown',
     *      name = '/infobank/tmp/beenthere.results.temp')
          open (unit = NavFile,
     *          access = 'sequential',
     *          form = 'formatted',
     *          status = 'old',
     *          name = InFile(1:InFileLength))
      open (unit = AreaKML,
     *      access = 'sequential',
     *      form = 'formatted',
     *      status = 'unknown',
     *      name = '/infobank/tmp/beenthere.'//
     *             UnderName(1:ANameLength)//
     *             '/'//UnderName(1:ANameLength)//'.kml')

      call system ('mkdir /infobank/tmp/beenthere.'//
     *             UnderName(1:ANameLength))
      call system ('sort -u /infobank/tmp/beenthere.ids.temp > 
     *             '/infobank/tmp/beenthere.unique.temp')
      call system ('rm -f /infobank/tmp/beenthere.ids.temp')
      call system ('cp -f /infobank/'//ID(1:1)//'/'//ID(1:IDLength)//
     *    '/kml/*.nav.kmz '//
     *    '/infobank/tmp/beenthere.'//UnderName(1:ANameLength)//'/')