USGS InfoBank program -- millermap

millermap

  Program to create GMT maps ala Greg Miller.

  
  ------------------ MAKE THE FIELD ACTIVITY MAP ---------------------
 4 6 8-1-2 4 6 8-1-2 4 6 8-1-2 4 6 8-1-2 4 6 8-1-2 4 6 8-1-2 4 6 8-1-2

       

          read (90, 10, iostat = ErrorCode2, err = 100, end = 1000)
     *        north, south, east, west
                  read (93, '(q,a)') WhereLength, Where(1:WhereLength)

              write (91, 11, iostat = ErrorCode3, err=110)
     *            west, north
              write (91, 11) east, north
              write (91, 11) east, south
              write (91, 11) west, south
              write (91, 11) west, north
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (93, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (ALine, 31, iostat = ErrorCode5, err = 23)
     *                west, east, south, north,
     *                WIDTH, res, glab
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (ALine, 32, iostat = ErrorCode5, err = 23)
     *                west, east, south, north,
     *                WIDTH, cont, smooth
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (ALine, 34, iostat = ErrorCode5, err = 23)
     *                west, east, south, north,
     *                WIDTH, res,
     *                glab
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (ALine, 35, iostat = ErrorCode5, err = 23)
     *                WIDTH,
     *                west, east, south, north
                  write (92, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)

          open (unit   = 91,
     *          access = 'sequential',
     *          form   = 'formatted',
     *          status = 'unknown',
     *          iostat = ErrorCode2,
     *          name   = '/infobank/tmp/millermap.corners.temp')
              open (unit   = 92,
     *              access = 'sequential',
     *              form   = 'formatted',
     *              status = 'unknown',
     *              iostat = ErrorCode4,
     *              name   = 'millermap.csh')
                  open (unit   = 93,
     *                access = 'sequential',
     *                form   = 'formatted',
     *                status = 'unknown',
     *                iostat = ErrorCode1,
     *                name   = '/infobank/tmp/millermap.where.temp')

      call openem (1, 0)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call system (
     *            '/bin/pwd > /infobank/tmp/millermap.where.tem
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call stringlen (SystemCommand, StringLength)
                  call removespace (ALine, ALineLength,
     *                              SystemCommand, StringLength)
                  call removespace (ALine, ALineLength,
     *                              SystemCommand, StringLength)
                  call removespace (ALine, ALineLength,
     *                              SystemCommand, StringLength)
                  call removespace (ALine, ALineLength,
     *                              SystemCommand, StringLength)
                  call system ('csh millermap.csh')
                  call forerror (ErrorCode4)
                  call forerror (ErrorCode3)
          call forerror (ErrorCode1)
              call stringlen (ALine, ALineLength)
                          call stringlen (ALine, ALineLength) 
              call stringlen (SystemCommand, StringLength)