USGS InfoBank program -- mainatlas4kml

atlas4kml

  Program to create virtual globe presentation of 
  InfoBank Atlas areas.
    

  Clint Steele    1/11/2008


    

 open (unit = 89,
*      access = 'sequential',
*      form = 'formatted',
*      status = 'old',
*      name = '/infobank/programs/keep/gazette.keep')

 open (unit = 88,
*      access = 'sequential',
*      form = 'formatted',
*      status = 'unknown',
*      name = '/tmp/atlas4kml.gazette.temp')

 open (unit = InUnit,
*      access = 'sequential',
*      form = 'formatted',
*      status = 'old',
*      name = '/tmp/atlas4kml.temp')

 open (unit = OutUnit, 
*      access = 'sequential',
*      form = 'formatted',
*      status = 'unknown',
*      name = '/infobank/programs/kml'//
*             '/atlas4kml/atlas4kml.kml')

    

     read (89, '(q,a)', end = 98)
*        ALineLength, ALine(1:ALineLength)
    

 write (Blank, '(80x)')
     write (88, '(a)') ALine(Tab(8)+1:Tab(9)-1)//
*        char(9)//ALine(1:Tab(8)-1)//
*        ALine(Tab(9):ALineLength)//
*        char(9)//char(9)//char(9)//char(9)
 write (OutUnit, '(a)') Blank(1:Level)//
*    '<?xml version="1.0" encoding="UTF-8"?>'
 write (OutUnit, '(a)') Blank(1:Level)//
*     '<kml xmlns="http://earth.google.com/kml/2.0">'
 write (OutUnit, '(a)') Blank(1:Level)//'<Document>'
 write (OutUnit, '(a)') Blank(1:Level)//'<name>Atlas Areas</name>'
 write (OutUnit, '(a)') Blank(1:Level)//'<LookAt>'
 write (OutUnit, '(a)') 
*    Blank(1:Level)//'<longitude>-113.33</longitude>'
 write (OutUnit, '(a)') 
*    Blank(1:Level)//'<latitude>41.1</latitude>'
 write (OutUnit, '(a)') 
*    Blank(1:Level)//'<altitude>3800000.</altitude>'
 write (OutUnit, '(a)') 
*    Blank(1:Level)//
*    '<altitudeMode>relativeToGround</altitudeMode>'

 write (OutUnit, '(a)') Blank(1:Level)//'<range>10800000.</range>'
 write (OutUnit, '(a)') Blank(1:Level)//'<tilt>0</tilt>'
 write (OutUnit, '(a)') Blank(1:Level)//'<heading>0</heading>'
 write (OutUnit, '(a)') Blank(1:Level)//'</LookAt>'
 write (OutUnit, '(a)') Blank(1:Level)//'<open>1</open>'
 write (OutUnit, '(a)') Blank(1:Level)//'</Folder>'  ! limit letter
 write (OutUnit, '(a)') Blank(1:Level)//'</Folder>'  ! area limits
 write (OutUnit, '(a)') Blank(1:Level)//'</Document>'
 write (OutUnit, '(a)') Blank(1:Level)//'</kml>'
    

 call system ('/bin/sort '//
*             '/tmp/atlas4kml.gazette.temp '//
*             '> /tmp/atlas4kml.temp')

 call system ('/bin/rm '//
*             '/tmp/atlas4kml.gazette.temp')

 call includeit (OutUnit,
*    '/infobank/programs/keep/infobank_style.keep')

 call areasfolder (InUnit, OutUnit, Level)
 call limitsfolder (InUnit, OutUnit, Level)
 call kmlzip 
*  ('/infobank/programs/kml/atlas4kml/atlas4kml.kml')