USGS InfoBank program -- timelatlong

timelatlong

  Subroutine to verify that a file is in the stardard

  time/lat/long format.

       

          read (TempUnit, '(q,a)', end = 99)
     *                ALineLength, ALine(1:ALineLength)
          read (ALine(16:35), '(f9.5,1x,f10.5)',
     *                  iostat = ReadError, err = 9)
     *                  Latitude, Longitude

              write (Message, 25) InFile(1:InFileLength),
     *                    RecordCount,
     *                    'blanks are not allowed in '//
     *                    'the decimal part of the latitude field'//
     *                    ' (columns 18-24)'
              write (Message, 25) InFile(1:InFileLength),
     *                    RecordCount,
     *                    'blanks are not allowed in '//
     *                    'the decimal part of the '//
     *                    'longitude field'//
     *                    ' (columns 29-35)'
              write (Message, 25) InFile(1:InFileLength),
     *                    RecordCount,
     *                    'column 19 must be a decimal point'
              write (Message, 25) InFile(1:InFileLength),
     *                    RecordCount,
     *                    'column 25 must be blank'
              write (Message, 25) InFile(1:InFileLength),
     *                    RecordCount,
     *                    'column 30 must be a decimal point'
                  write (Message, 25) InFile(1:InFileLength),
     *                        RecordCount,
     *                        'column 36 must be blank when '//
     *                        'data exists beyond it'
          write (Message, 25) InFile(1:InFileLength),
     *                    RecordCount,
     *                    'latitude and longitude (fields 16-35) '//
     *                    'can only be decimal numbers'
              write (Message, 25) InFile(1:InFileLength),
     *                RecordCount,
     *                    'latitude field (columns 16-24) '//
     *                    'must fall within -90.0 to 90.0 range'
              write (Message, 25) InFile(1:InFileLength),
     *                RecordCount,
     *                    'longitude field (columns 26-35) '//
     *                    'must fall within -179.99999 to '//
     *                    '180.00000 range'

      open (unit   = TempUnit,
     *      access = 'sequential',
     *      form   = 'formatted',
     *      status = 'unknown',
     *      name   = InFile(1:InFileLength))

          call stringlen (ALine(1:ALineLength), 
     *                            ALineLength)