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
      write (OutUnit, '(a)') Blank(1:Level)//'<name>area names&l
          write (OutUnit, '(a)') 
     *         Blank(1:Level)//'</Folder>' ! names l
          write (OutUnit, '(a)') 
     *         Blank(1:Level)//'<Folder>'  ! names 
          write (OutUnit, '(a)') 
     *         Blank(1:Level)//'<name>'//Current//'&l
      write (OutUnit, '(a)') Blank(1:Level)//'<Placemar
      write (OutUnit, '(a)') Blank(1:Level)//
     *    '<name>'//AtlasName(1:AtlasNameLength)//'&l
      write (OutUnit, '(a)') 
     *    Blank(1:Level)//'<styleUrl>#USGS</s
      write (OutUnit, '(a)') Blank(1:Level)//'<LookA
      write (OutUnit, '(a,f10.5,a)') Blank(1:Level)//
     *   '<longitude>', Longitude, '</lo
      write (OutUnit, '(a,f10.5,a)') Blank(1:Level)//
     *   '<latitude>', Latitude, '</l
      write (OutUnit, '(a,f,a)') Blank(1:Level)//'<range
     *    Range,'</rang
      write (OutUnit, '(a)') Blank(1:Level)//'</LookA
      write (OutUnit, '(a)') 
     *     Blank(1:Level)//'<visibility>1</vis
      write (OutUnit, '(a)') Blank(1:Level)//'<open>1&l
      write (OutUnit, '(a)') Blank(1:Level)//
     *    '<name>'//AtlasName(1:AtlasNameLength)//'&l
      write (OutUnit, '(a)') Blank(1:Level)//'<Poin
      write (OutUnit, '(a)') Blank(1:Level)//'<coordinate
      write (OutUnit, '(a,f10.5,a,f9.5,a,a)') Blank(1:Level),
     *     Longitude, ',', Latitude, ',', Z
      write (OutUnit, '(a)') Blank(1:Level)//'</coordinate
      write (OutUnit, '(a)') Blank(1:Level)//'</Poin
      write (OutUnit, '(a)') Blank(1:Level)//'<descriptio
      write (OutUnit, '(a)') Blank(1:Level)//'<![CDATA[
     *  '<a href="http://walrus.wr.usgs.gov/infobank/
     *  'gazette/htmlids/regions/'//AtlasID(1:AtlasIDLength)//
     *  '.html">Field Activities</
      write (OutUnit, '(a)') Blank(1:Level)//'</descriptio
      write (OutUnit, '(a)') 
     *    Blank(1:Level)//'<visibility>1</vis
      write (OutUnit, '(a)') Blank(1:Level)//'<open>1&l
      write (OutUnit, '(a)') Blank(1:Level)//'</Placemar
      write (OutUnit, '(a)') Blank(1:Level)//'</Folder>'  ! names 
      write (OutUnit, '(a)') Blank(1:Level)//'</Folder>'  ! area

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