USGS InfoBank program -- areasfolder

areasfolder

  Subroutine to create a kml folder for atlas areas
    

  Clint Steele      1/18/2008


    

 read (InUnit, '(q,a)', end=199) ALineLength, 
*                   ALine(1:ALineLength)

 read (ALine(Tab(2)+1:Tab(3)-1), '(f16)') North
 read (ALine(Tab(3)+1:Tab(4)-1), '(f16)') South
 read (ALine(Tab(4)+1:Tab(5)-1), '(f16)') East
 read (ALine(Tab(5)+1:Tab(6)-1), '(f16)') West
         read (ALine(Tab(11)+1:Tab(12)-1), '(f)') Latitude
         read (ALine(Tab(11)+1:Tab(12)-1), '(i)') ILatitude
         read (ALine(Tab(12)+1:Tab(13)-1), '(f)') Longitude
         read (ALine(Tab(12)+1:Tab(13)-1), '(i)') ILongitude
    

 write (Blank, '(80x)')
 write (OutUnit, '(a)') Blank(1:Level)//'<Folder>'  ! area names
 write (OutUnit, '(a)') Blank(1:Level)//'<name>area names</name>'
     write (OutUnit, '(a)') 
*         Blank(1:Level)//'</Folder>' ! names letter 
     write (OutUnit, '(a)') 
*         Blank(1:Level)//'<Folder>'  ! names letter
     write (OutUnit, '(a)') 
*         Blank(1:Level)//'<name>'//Current//'</name>'

 write (OutUnit, '(a)') Blank(1:Level)//'<Placemark>'
 write (OutUnit, '(a)') Blank(1:Level)//
*    '<name>'//AtlasName(1:AtlasNameLength)//'</name>'

 write (OutUnit, '(a)') 
*    Blank(1:Level)//'<styleUrl>#USGS</styleUrl>'

 write (OutUnit, '(a)') Blank(1:Level)//'<LookAt>'
 write (OutUnit, '(a,f10.5,a)') Blank(1:Level)//
*   '<longitude>', Longitude, '</longitude>'
 write (OutUnit, '(a,f10.5,a)') Blank(1:Level)//
*   '<latitude>', Latitude, '</latitude>'

 write (OutUnit, '(a,f,a)') Blank(1:Level)//'<range>',
*    Range,'</range>'

 write (OutUnit, '(a)') Blank(1:Level)//'</LookAt>'
 write (OutUnit, '(a)') 
*     Blank(1:Level)//'<visibility>1</visibility>'
 write (OutUnit, '(a)') Blank(1:Level)//'<open>1</open>'
 write (OutUnit, '(a)') Blank(1:Level)//
*    '<name>'//AtlasName(1:AtlasNameLength)//'</name>'

 write (OutUnit, '(a)') Blank(1:Level)//'<Point>'
 write (OutUnit, '(a)') Blank(1:Level)//'<coordinates>'
 write (OutUnit, '(a,f10.5,a,f9.5,a,a)') Blank(1:Level),
*     Longitude, ',', Latitude, ',', Z

 write (OutUnit, '(a)') Blank(1:Level)//'</coordinates>'
 write (OutUnit, '(a)') Blank(1:Level)//'</Point>'
 write (OutUnit, '(a)') Blank(1:Level)//'<description>'
 write (OutUnit, '(a)') Blank(1:Level)//'<![CDATA['//
*  '<a href="http://walrus.wr.usgs.gov/infobank/'//
*  'gazette/htmlids/regions/'//AtlasID(1:AtlasIDLength)//
*  '.html">Field Activities</a>]]>'

 write (OutUnit, '(a)') Blank(1:Level)//'</description>'
 write (OutUnit, '(a)') 
*    Blank(1:Level)//'<visibility>1</visibility>'
 write (OutUnit, '(a)') Blank(1:Level)//'<open>1</open>'
 write (OutUnit, '(a)') Blank(1:Level)//'</Placemark>'
 write (OutUnit, '(a)') Blank(1:Level)//'</Folder>'  ! names letter
 write (OutUnit, '(a)') Blank(1:Level)//'</Folder>'  ! area names
    

 call changecase ('SMALLS', AtlasID(1:AtlasIDLength))
 call changecase ('CAPITALS', CapAtlasID)
 call getrange (Latitude, Longitude,
*    North, South, East, West, Range)