USGS InfoBank program -- area2atlas

area2atlas

  Program to generate <Activity_ID>.atlas from
  <Activity-ID>.area and gazette.keep files.
    

  Typing "area2atlas" runs "/infobank/programs/shells/area2atlas"

  it creates "/tmp/area2atlas.control.temp" (a list of activity IDs)
  and then executes this program.

  This area2atlas reads "/infobank/programs/keep/gazette.keep" 
  information.
  Loops through for each ID in the "/tmp/area2atlas.control.temp" 
  ID list.
    

  Clint Steele   3/19/2004 


    

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

 open (unit = 91,
*      access = 'sequential',
*      form = 'formatted',
*      status = 'unknown',
*      name = '/tmp/area2atlas.control.temp')

         open (unit = 93,
*            access = 'sequential',
*              form = 'formatted',
*            status = 'unknown',
*              name = '/infobank/'//
*                 ShortID(1:1)//'/'//
*                 ShortID(1:IDLength)//'/meta/'//
*                 LongID(1:IDLength+3)//'.area')

     open (unit = 94,
*          access = 'sequential',
*          form = 'formatted',
*          status = 'unknown',
*          name = '/infobank/'//
*                 ShortID(1:1)//'/'//
*                 ShortID(1:IDLength)//'/meta/'//
*                 LongID(1:IDLength+3)//'.atlas')

    

     read (90, '(q,a)', end = 99) GazLength,
*        Gaz(GazCount+1)(1:GazLength)
     read (Gaz(Loop)(             1:Tab(Loop,1)-1), '(f)') 
*         AtlasNorth(Loop)
     read (Gaz(Loop)(Tab(Loop,1)+1:Tab(Loop,2)-1), '(f)') 
*         AtlasSouth(Loop)
     read (Gaz(Loop)(Tab(Loop,2)+1:Tab(Loop,3)-1), '(f)') 
*         AtlasEast(Loop)
     read (Gaz(Loop)(Tab(Loop,3)+1:Tab(Loop,4)-1), '(f)') 
*         AtlasWest(Loop)

     read (Gaz(Loop)(Tab(Loop,4)+1:Tab(Loop,5)-1), '(a)')
*         AtlasID(Loop)
     read (Gaz(Loop)(Tab(Loop,5)+1:Tab(Loop,5)+2), '(a)')
*         RegionID(Loop)
     read (91, '(q,a)', end = 999)
*        IDLength, ShortID(1:IDLength)

         read (93, '(t10,f9.5,t20,f9.5,t30,f10.5,t41,f10.5)')
*            IDNorth, IDSouth, IDEast, IDWest 

    

         write (LongID, '(a)')
*            ShortID(1:1)//'-'//
*            ShortID(2:2)//'-'//
*            ShortID(3:4)//'-'//
*            ShortID(5:6)
         write (LongID, '(a)')
*            ShortID(1:1)//'-'//
*            ShortID(2:3)//'-'//
*            ShortID(4:5)//'-'//
*            ShortID(6:7)
     write (94, '(a,/,a,/,a,/,a,/,a,/,a)')
* '! CAUTION - THIS .KEEP FILE IS GENERATED BY AN APPLICATION',
* '!           ANY EDITING TO IT WILL BE LOST',
* '!           THE NEXT TIME THE APPLICATION IS EXECUTED',
* '! CREATED BY area2atlas',
* '! USED BY atlasgroups',
* '!' 

            write (94, '(a,t10,a)')
*           BigID(1:IDLength), BigAtlasID(Loop) 
            write (94, '(a,t10,a)')
*           BigID(1:IDLength), BigAtlasID(Loop) 
    

     call changecase ('SMALLS', AtlasID(Loop))
     call changecase ('CAPITALS', BigAtlasID(Loop))
     call changecase ('SMALLS', RegionID(Loop))
     call changecase ('SMALLS', ShortID(1:IDLength))
         call system ('/bin/rm '//
*                 '/infobank/'//
*                 ShortID(1:1)//'/'//
*                 ShortID(1:IDLength)//'/meta/'//
*                 LongID(1:IDLength+3)//'.atlas')
     call changecase ('CAPITALS', BigID)