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

USGS InfoBank program -- getarea

Skip navigational links
Search InfoBank
Home tab FACS tab Activities tab Atlas tab Geology School tab More tab More tab Activities tabs
   
InfoBank Programs: by Name   by Topic  
   
Expanded description
Topic Description
Name
getarea
Synopsis
/infobank/programs/nav/getarea/getarea.for
Description
  Program to generate a .area file
    
Author
  Clint Steele    10/22/97

  Clint Steele    10/29/2003  changed prompting and input handling
  Carolyn Degnan  09/08/2008  padded single point nav with 0.5 degree all around

    
Opens
      open (unit = 90,
*          access = 'sequential',
*          form = 'formatted',
*          status = 'unknown',
*          name = '/infobank/'//
*                 SmallID(1:1)//'/'//
*                 SmallID(1:IDLength)//'/nav/'//
*                 SmallID(1:1)//'-'//
*                 SmallID(2:2)//'-'//
*                 SmallID(3:4)//'-'//
*                 SmallID(5:6)//'.best_nav_file')

     open (unit = 90,
*          access = 'sequential',
*          form = 'formatted',
*          status = 'unknown',
*          name = '/infobank/'//
*                 SmallID(1:1)//'/'//
*                 SmallID(1:IDLength)//'/nav/'//
*                 SmallID(1:1)//'-'//
*                 SmallID(2:3)//'-'//
*                 SmallID(4:5)//'-'//
*                 SmallID(6:7)//'.best_nav_file')

     open (unit = 91,
*          access = 'sequential',
*          form = 'formatted',
*          status = 'unknown',
*          name = AreaPath(1:AreaPathLength))

 open (unit = 92,
*          access = 'sequential',
*          form = 'formatted',
*          status = 'unknown',
*          name = NavFile(1:NavFileLength))

     open (unit = 91,
*          access = 'sequential',
*          form = 'formatted',
*          status = 'unknown',
*          name = AreaPath(1:AreaPathLength))

    
Reads
 read (5, '(a)')  CapID
     read (91, '(t10,f9.5,t20,f9.5,t30,f10.5,t41,f10.5)')
*        OldNorth, OldSouth, OldEast, OldWest

 read (90, '(q,a)', end = 299) 
*    NavFileLength, NavFile(1:NavFileLength)

         read (92, '(q,a)', end = 99) 
*            ALineLength, ALine(1:ALineLength)

     read (ALine(1:ALineLength), '(t26,f10.5)') Long
         read (92, '(q,a)', end = 99) 
*            ALineLength, ALine(1:ALineLength)

     read (ALine(1:ALineLength), '(t16,f9.5,t26,f10.5)') Lat, Long          
    
Writes
     write (91, '(a,t10,f9.5,t20,f9.5,t30,f10.5,t41,f10.5)')
*       CapID, North, South, East, West

    
Calls
 call stringlen (CapID, IDLength)
 call changecase ('CAPITALS', CapID(1:IDLength))
 call changecase ('SMALLS', SmallID(1:IDLength))
 call stringlen (AreaPath, AreaPathLength)
     call limits (Long,Lat,West,South,East,North)
     call system ('sh /infobank/programs/shells/area2atlas '//
*        SmallID(1:IDLength))

    

Skip footer navigational links


InfoBank   Menlo Park & Santa Cruz Centers   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