USGS InfoBank program -- limitsfolder

limitsfolder

  Subroutine to create KML folder of atlas limits

       

  Clint Steele    1/18/2008


    

      read (InUnit, '(q,a)', end=299) 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

      write (Blank, '(80x)')
      write (OutUnit, '(a)') Blank(1:Level)//'<Folder>'  ! area 
      write (OutUnit, '(a)') Blank(1:Level)//'<name>area limits&l
          write (OutUnit, '(a)') 
     *        Blank(1:Level)//'</Folder>'  ! limit 
          write (OutUnit, '(a)') 
     *        Blank(1:Level)//'<Folder>'  ! limit 
          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)//'<visibility>0</vis
      write (OutUnit, '(a)') Blank(1:Level)//'<open>1&l
      write (OutUnit, '(a)') Blank(1:Level)//'<Styl
      write (OutUnit, '(a)') Blank(1:Level)//'<LineStyl
      write (OutUnit, '(a)') Blank(1:Level)//'<color>ffffffff<
      write (OutUnit, '(a)') Blank(1:Level)//'</LineStyl
      write (OutUnit, '(a)') Blank(1:Level)//'<PolyStyl
      write (OutUnit, '(a)') Blank(1:Level)//'<color>7f00ff00<
      write (OutUnit, '(a)') Blank(1:Level)//'</PolyStyl
      write (OutUnit, '(a)') Blank(1:Level)//'</Styl
      write (OutUnit, '(a)') Blank(1:Level)//'<LineStrin
      write (OutUnit, '(a)') Blank(1:Level)//'<extrude>1</
      write (OutUnit, '(a)') 
     *    Blank(1:Level)//'<tessellate>1</tes
      write (OutUnit, '(a)') Blank(1:Level)//'<coordinate
          write (OutUnit, '(a,f10.5,a,f9.5,a,a)') Blank(1:Level),
     *    TempWest, ',', North, ',', Z
                  write (OutUnit, '(a,f10.5,a,f9.5,a,a)') 
     *                Blank(1:Level), TempLong, ',', TempLat, ',', Z
                  write (OutUnit, '(a,f10.5,a,f9.5,a,a)') 
     *                Blank(1:Level), TempLong, ',', TempLat, ',', Z
                  write (OutUnit, '(a,f10.5,a,f9.5,a,a)') 
     *                Blank(1:Level), TempLong, ',', TempLat, ',', Z
          write (OutUnit, '(a,f10.5,a,f9.5,a,a)') Blank(1:Level),
     *        TempEast, ',', North, ',', Z
              write (OutUnit, '(a,f10.5,a,f9.5,a,a)') 
     *            Blank(1:Level), TempLong, ',', TempLat, ',', Z
                  write (OutUnit, '(a,f10.5,a,f9.5,a,a)') 
     *                Blank(1:Level), TempLong, ',', TempLat, ',', Z
                  write (OutUnit, '(a,f10.5,a,f9.5,a,a)') 
     *                Blank(1:Level), TempLong, ',', TempLat, ',', Z
          write (OutUnit, '(a,f10.5,a,f9.5,a,a)') 
     *        Blank(1:Level), TempEast, ',', South, ',', Z
                  write (OutUnit, '(a,f10.5,a,f9.5,a,a)') 
     *                Blank(1:Level), TempLong, ',', TempLat, ',', Z
                  write (OutUnit, '(a,f10.5,a,f9.5,a,a)') 
     *                Blank(1:Level), TempLong, ',', TempLat, ',', Z
                  write (OutUnit, '(a,f10.5,a,f9.5,a,a)') 
     *                Blank(1:Level), TempLong, ',', TempLat, ',', Z
          write (OutUnit, '(a,f10.5,a,f9.5,a,a)') 
     *        Blank(1:Level), TempWest, ',', South, ',', Z
              write (OutUnit, '(a,f10.5,a,f9.5,a,a)') 
     *            Blank(1:Level), TempLong, ',', TempLat, ',', Z
                  write (OutUnit, '(a,f10.5,a,f9.5,a,a)') 
     *                Blank(1:Level), TempLong, ',', TempLat, ',', Z
                  write (OutUnit, '(a,f10.5,a,f9.5,a,a)') 
     *                Blank(1:Level), TempLong, ',', TempLat, ',', Z
          write (OutUnit, '(a,f10.5,a,f9.5,a,a)') 
     *        Blank(1:Level), TempWest, ',', North, ',', Z
      write (OutUnit, '(a)') Blank(1:Level)//'</coordinate
      write (OutUnit, '(a)') Blank(1:Level)//'</LineStrin
      write (OutUnit, '(a)') Blank(1:Level)//'</Placemar

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