Link to USGS home page
USGS Home
Contact USGS
Search USGS
Coastal & Marine Geology InfoBank

USGS InfoBank program -- atlasarea

Skip navigational links
Search InfoBank
Home tab Atlas tab Activities tab FACS tab Geology School tab More tab Activities tabs
   
InfoBank Programs: by Name   by Topic  
   
Expanded description
Topic Description
Name
atlasarea
Synopsis
/infobank/programs/dbase/atlasarea/atlasarea.for
Description
  Program to determine the smallest atlas area that contains
  the navigation for a field activity.
    
Author
  Clint Steele    1/22/2008


    
Opens
 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/atlasarea.gazarea.temp')

 open (unit = 92,
*      access = 'sequential',
*      form = 'formatted',
*      status = 'unknown',
*      name = '/tmp/atlasarea.gazarea.sorted.temp')

 open (unit = 94,
*      access = 'sequential',
*      form = 'formatted',
*      status = 'old',
*      name = '/tmp/atlasarea.navareas.temp')

     open (unit = 93,
*          access = 'sequential',
*          form = 'formatted',
*          status = 'old',
*          name = IDLine(1:IDLineLength))

             open (unit = 95,
*                  access = 'sequential',
*                  form = 'formatted',
*                  status = 'unknown',
*                  name = IDLine(1:IDLineLength-8)//'.atlasarea')
    
Reads
     read (90, '(q,a)', end = 99)
*        GazLineLength, GazLine(1:GazLineLength)

         read (GazLine(1:Tab(1)-1), '(f)') GazNorth
         read (GazLine(1:Tab(1)-1), '(i)') iGazNorth
         read (GazLine(Tab(1)+1:Tab(2)-1), '(f)') GazSouth
         read (GazLine(Tab(1)+1:Tab(2)-1), '(i)') iGazSouth
         read (GazLine(Tab(2)+1:Tab(3)-1), '(f)') GazEast
         read (GazLine(Tab(2)+1:Tab(3)-1), '(i)') iGazEast
         read (GazLine(Tab(3)+1:Tab(4)-1), '(f)') GazWest
         read (GazLine(Tab(3)+1:Tab(4)-1), '(i)') iGazWest
     read (94, '(t13,q,a)', end = 999)
*        IDLineLength, IDLine(1:IDLineLength)

     read (93, '(a7, t10, f9.5, 1x, f9.5, 1x, f10.5, 1x, f10.5)')
*        ID, IDNorth, IDSouth, IDEast, IDWest
         read (92, '(q,a)', end=98) 
*            GazLineLength, GazLine(1:GazLineLength)
         read (GazLine(1:GazLineLength),
*            '(t10,f9.5,1x,f9.5,1x,f10.5,1x,f10.5)')
*            GazNorth, GazSouth, GazEast, GazWest

    
Writes
     write (91, '(f8.1,t10,f9.5,1x,f9.5,1x,f10.5,1x,f10.5,t52,a)')
*        GazArea, GazNorth, GazSouth, GazEast, GazWest,
*        GazLine(Tab(8)+1:Tab(9)-1)

             write (95, 
*                '(a7,t10,f9.5,1x,f9.5,1x,f10.5,1x,f10.5,1x,a)') 
*                ID, IDNorth, IDSouth, IDEast, IDWest, 
*                GazLine(52:GazLineLength)
    
Calls
 call system ('/bin/sort /tmp/atlasarea.gazarea.temp '//
*    ' > /tmp/atlasarea.gazarea.sorted.temp')

 call system ('/bin/rm /tmp/atlasarea.gazarea.temp')
         call inside (IDNorth, IDWest,
*            GazNorth, GazSouth, GazWest, GazEast, NW)
         call inside (IDNorth, IDEast,
*            GazNorth, GazSouth, GazWest, GazEast, NE)
         call inside (IDSouth, IDEast,
*            GazNorth, GazSouth, GazWest, GazEast, SE)
         call inside (IDSouth, IDWest,
*            GazNorth, GazSouth, GazWest, GazEast, SW)

    

Skip footer navigational links


InfoBank   Menlo Park Center   Santa Cruz Center   St. Petersburg Center   Woods Hole Center   Coastal and Marine Geology Program   Geologic Information   Ask-A-Geologist   USGS Disclaimer  

FirstGov button   Take Pride in America button