USGS InfoBank program -- millermap0

millermap0

  Program to create GMT maps ala Greg Miller.

 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 (91, '(q,a)') WhereLength, Where(1:WhereLength)
          read (90, 2, iostat = ErrorCode2, err = 3, end = 4)
     *        north, south, east, west

              write (92, 12, iostat = ErrorCode3, err=13)
     *            west, north
              write (92, 12) east, north
              write (92, 12) east, south
              write (92, 12) west, south
              write (92, 12) west, north
                  write (93, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (93, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (93, '(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 (93, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (ALine, 32, iostat = ErrorCode5, err = 23)
     *                west, east, south, north,
     *                WIDTH, cont, smooth
                  write (93, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (ALine, 34, iostat = ErrorCode5, err = 23)
     *                west, east, south, north,
     *                WIDTH, res,
     *                glab
                  write (93, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)
                  write (ALine, 35, iostat = ErrorCode5, err = 23)
     *                WIDTH,
     *                west, east, south, north
                  write (93, '(a)', iostat = ErrorCode5, err = 23)
     *               SystemCommand(1:StringLength)

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

      call system (
     *    '/bin/pwd > /infobank/tmp/millermap.where.tem
      call openem (1, 0)
                  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)