USGS InfoBank program -- grids

grids

  Program to determine 1 degree grids a best_nav_file passes through.

       

  $ GRIDS

  SUBROUTINES
  getID (external shareable)
    

  Clint Steele           5/9/86         .for

  Gary Deeter         25-aug-1986
  Clint Steele          10/16/2006
  All grids in a particular run are appended to a single file.

    

      read (ControlFile, '(q,a)', end=98) IDLength, ID(1:IDLength)
          read (BestFile, '(q,a)', end = 96) 
     *        InFileLength, InFile(1:InFileLength)
          read (WCFile, '(i8)') WordCount
          read (NavFile, '(q,a)', end = 97) 
     *       ALineLength, ALine(1:ALineLength)
          read (ALine(1:ALineLength), 1) ILat, ILong

      open (unit = ControlFile,
     *      access = 'sequential',
     *      form = 'formatted',
     *      status = 'unknown',
     *      name = '/infobank/tmp/grids.control.temp')
      open (unit = BestFile,
     *      access = 'sequential',
     *      form = 'formatted',
     *      status = 'old',
     *      name = '/infobank/'//ID(1:1)//'/'//
     *             ID(1:IDLength)//'/nav/'//
     *             LongID(1:IDLength+3)//'.best_nav_file')
      open (unit = GridsFile,
     *      access = 'sequential',
     *      form = 'formatted',
     *      status = 'unknown',
     *      name = '/infobank/'//ID(1:1)//'/'//
     *             ID(1:IDLength)//'/meta/'//
     *             LongID(1:IDLength+3)//'.grids')
          open (unit = WCFile,
     *          access = 'sequential',
     *          form = 'formatted',
     *          status = 'unknown',
     *          name = '/infobank/tmp/grids.wc.temp')
          open (unit = NavFile,
     *          access = 'sequential',
     *          form = 'formatted',
     *          status = 'old',
     *          name = InFile(1:InFileLength))

      call changecase ('SMALLS', ID(1:IDLength))
      call stringlen (ID(1:IDLength), IDLength)
          call system ('wc '//InFile(1:InFileLength)//
     *                 ' > /infobank/tmp/grids.wc.tem