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


    

 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 = '/tmp/beenthere.ids.temp')

 open (unit = UniqueIDs,
*      access = 'sequential',
*      form = 'formatted',
*      status = 'unknown',
*      name = '/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 = '/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 = '/tmp/beenthere.'//
*             UnderName(1:ANameLength)//
*             '/'//UnderName(1:ANameLength)//'.kml')

    

 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)') '<Document>'
 write (AreaKml, '(a)') '  <name>'//AName(1:ANameLength)//'</name>'
 write (AreaKml, '(a)') '  <ScreenOverlay>'
 write (AreaKml, '(a)') '    <Icon>'
 write (AreaKml, '(a)') 
*    '      <href>http://walrus.wr.usgs.gov/infobank/'//
*    'programs/html/tools/SmUSGS.png</href>'
 write (AreaKml, '(a)') '    </Icon>'
 write (AreaKml, '(a)') '    <name>USGS Visual Identity</name>'
 write (AreaKml, '(a)') 
*    '    <overlayXY x="0" y="1" xunits="fraction" '//
*    'yunits="fraction"/>'
 write (AreaKml, '(a)') 
*    '    <screenXY x="0" y="1" xunits="fraction" '//
*    'yunits="fraction"/>'
 write (AreaKml, '(a)') '  </ScreenOverlay>'
 write (AreaKml, '(a)') '<Placemark>'
 write (AreaKml, '(a)') '  <description>'
 write (AreaKml, '(a)') '    <![CDATA['
 write (AreaKml, '(a)') 
*    '          <li><a href="http://www.google.com/'//
*    '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</a>'
 write (AreaKml, '(a)') 
*    '    <!-- Source: http://walrus.wr.usgs.gov/infobank/ --!>'
 write (AreaKml, '(a)') '    ]]>'
 write (AreaKml, '(a)') '  </description>'
 write (AreaKml, '(a)') '  <name>'//AName(1:ANameLength)//'</name>'
 write (AreaKml, '(a)') '  <visibility>1</visibility>'
 write (AreaKml, '(a)') '  <open>0</open>'
 write (AreaKml, '(a)') '  <Style>'
 write (AreaKml, '(a)') '    <LineStyle>'
 write (AreaKml, '(a)') '      <color>ff0000ff</color>'
 write (AreaKml, '(a)') '      <width>3</width>'
 write (AreaKml, '(a)') '    </LineStyle>'
 write (AreaKml, '(a)') '    <PolyStyle>'
 write (AreaKml, '(a)') '      <color>7f00ff00</color>'
 write (AreaKml, '(a)') '    </PolyStyle>'
 write (AreaKml, '(a)') '  </Style>'
 write (AreaKml, '(a)') '  <MultiGeometry>'
 write (AreaKml, '(a)') '    <LineString>'
 write (AreaKml, '(a)') '      <extrude>1</extrude>'
 write (AreaKml, '(a)') '      <tessellate>1</tessellate>'
 write (AreaKml, '(a)') '      <coordinates>'
 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)') '      </coordinates>'
 write (AreaKml, '(a)') '    </LineString>'
 write (AreaKml, '(a)') '  </MultiGeometry>'
 write (AreaKml, '(a)') '</Placemark>'
 write (AreaKml, '(a)') '<Placemark>'
 write (AreaKml, '(a)') '<visibility>1</visibility>'
 write (AreaKml, '(a)') '<name>'//AName(1:ANameLength)//'</name>'
 write (AreaKml, '(a)') '<Point>'
 write (AreaKml, '(a)') '  <coordinates>'
 write (AreaKml, '(a,f,a,f,a)') '        ',West,', ',North,', 0'
 write (AreaKml, '(a)') '  </coordinates>'
 write (AreaKml, '(a)') '</Point>'
 write (AreaKml, '(a)') '</Placemark>'
 write (AreaKml, '(a)') '</Document>'
 write (AreaKml, '(a)') '</kml>'
    

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