if [ $# -lt 3 ]
then
  echo " "
  echo "            *****     COARSE     *****      "
  echo " "
  echo "The program deletes every second row and column of a grid"
  echo " "
  echo "      Syntax: coarse <file1> <file2> <integer>"
  echo " "
  echo "     <integer> - number of rows in the input file"
  echo " "
  exit 0
else
 echo " "
 echo "            *****     COARSE     *****      "
 echo " "
 echo "file1="$1
 echo "file2="$2
 cp $1 fort.1
fi
cat > coarse.f <<EOF
      PROGRAM PCOARSE
      IMPLICIT REAL (A-H,O-Z)
      
      PRINT *,'0**** COARSE (INPUT,OUTPUT,TAPE1,TAPE2)****'
      PRINT *,'   '
      NCHAN=1
      CALL PARLES (NCHAN,NDIM,NSET)
      CALL COARSE (NDIM) 
      STOP 
      END
C******************************************************************
      SUBROUTINE PARLES (NCHAN,NDIM,NSET)
      IMPLICIT REAL (A-H,O-Z)
C
      INTEGER I1,I2,I3
      REAL F

      NSET=0

      REWIND NCHAN
      READ(NCHAN,END=10) ITIME,INAME,ILEV,ISIZE
      GOTO 20
   10 PRINT *,'INPUT FILE ON TAPE',NCHAN,' IS EMPTY'
      STOP 'MISTAKE'
   20 REWIND NCHAN

      NDIM=0  
    1 READ (NCHAN,END=100) I1,I2,I3,I4
      NDIM=MAX0(NDIM,I4)
      READ (NCHAN,END=101) F
      NSET=NSET+1
      GOTO1
      
  101 PRINT *,' MORE HEADER THEN DATA '     
  100 PRINT 1000,NCHAN,NDIM,NSET
      REWIND NCHAN
      RETURN    
 1000 FORMAT (1X,' TAPE=',I2,',  NDIM=',I5,',  NSET=',I5)
      END
C*****************************************************************
C
      SUBROUTINE COARSE(IE)
      IMPLICIT REAL (A-H,O-Z)
C
      PARAMETER (UNDEF=0.9E+10,UNDEFE=0.9*UNDEF)
C IF YOU USE THIS PROGRAM NOT ON A CRAY-2 ACTIVATE
C THE FOLLOWING TWO LINES
      PARAMETER (IE77=20000)
      DIMENSION X(IE77), ID(4)
C AND DELETE THE FOLLOWING LINES
C      PARAMETER (IE77=9E12)
C      DIMENSION X(IE), ID(4)
C      POINTER (NX,X)
C      ALLOCATE (NX,IE)
C NOW YOU HAVE DONE ALL CHANGES FOR STANDARD FORTRAN.
      DATA ONE /1./
      NROWS = 0
      PRINT *,' ENTER NUMBER OF ROWS: '
      READ(5,'(I10)',END=10) NROWS
   10 WRITE(6,6001) NROWS
      REWIND 1
      READ(1,END=11) ID
      WRITE(6,6010) ' FIRST RECORD:', ID
      ICNT = 1
      IF(ID(4).GT.IE) GOTO 12
      IF(NROWS.EQ.0) NROWS = 1
      NCOL = ID(4) / NROWS
      IF(NROWS*NCOL .NE. ID(4)) GOTO 14
C
C*  MAIN LOOP.
C
      REWIND 1
      ICNT=0
  100 READ(1,END=9920) ID
        ICNT=ICNT+1
        IF(ID(4).GT.IE) GOTO 12
        READ(1,END=13) (X(I),I=1,ID(4))
        IZ = 0
        DO 101 I = 1,NROWS
          IF(IFIX(I/2.)*2 .EQ. I) GOTO 101
            DO  102 J = 1, NCOL, 2
              IZ = IZ+1
  102         X(IZ) = X((I-1)*NCOL + J)
  101     CONTINUE
        ID(4) = IZ
        WRITE(2) ID
        WRITE(2) (X(I),I=1,ID(4))
      GOTO 100
C
C*  ABNORMAL TERMINATION.
C
   11 WRITE(6,6010) '<TAPE1> EMPTY '
         GOTO 9910
   12 WRITE(6,6010) 'RECORD TOO LONG: ', ID, IE, ICNT
         GOTO 9910
   13 WRITE(6,6010) ' READ ERROR: ', ID, ICNT
         GOTO 9910
   14 WRITE(6,6010) ' NROWS * NCOL .NE. ID(4) ', NROWS, NCOL, ID(4)
         GOTO 9910
 9910 CONTINUE
C
C*  CAUSE AN ABORT.
C
         PRINT *,'*COARSE*: FATAL END '
         CALL MINONE(ONE,XNEGO)
         ZZ=SQRT(XNEGO)
         STOP
C
C*  NORMAL TERMINATION.
C
 9920 WRITE(6,6010) ' LAST RECORD: ', ID, ICNT
      STOP ' REGULAR END! '
 6001 FORMAT(' NUMBER OF ROWS PER DATAFIELD',/,' NROWS ',I10)
 6010 FORMAT('0*COARSE*: ', A, 4I10,/1X,5I10)
      END
      SUBROUTINE MINONE(XIN,XOUT)
        XOUT=-XIN
      RETURN
      END
EOF
f77 coarse.f -o coarse.x
coarse.x $1 $2 $3 <<EOF
$3
EOF
cp fort.2 $2
rm fort.1 fort.2
rm coarse.x coarse.f
exit