if [ $# -lt 3 ]
then
  echo " "
  echo "     *****     CHANOSO     *****      "
  echo " "
  echo "  The program reverses a grid's lines order"
  echo " "
  echo " Syntax: chanoso <file1> <file2> <integer>"
  echo " "
  echo " <integer> - nr. of lines in the input file"
  echo " "
  exit 0
else
 echo " "
 echo "     *****     CHANOSO     *****      "
 echo " "
 echo "file1="$1
 echo "file2="$2
 echo "Number of lines in input file="$3
 cp $1 fort.1
fi
cat > chanoso.f <<EOF
      PROGRAM PCHANS
      IMPLICIT REAL (A-H,O-Z)
      PRINT *,'0**** CHANS (INPUT,OUTPUT,TAPE1,TAPE2)****'
      PRINT *,'   '
      NCHAN=1
      CALL PARLES (NCHAN,NDIM,NSET)
C
C* GET THE CHANS CARD.
C
      PRINT *,' *CHANS* ENTER: NROWS '
      PRINT *,'    FOR EXAMPLE:   50  '
      READ(5, *,END=12) NROWS
      IF(NROWS.EQ.0) NROWS = 1
   11 NCOL = NDIM/NROWS
      IF(NCOL*NROWS .NE. NDIM) GOTO 13
      WRITE(6,6011) ' NROWS, NCOLS '
      WRITE(6,'(8X,6(1X,I6),/)') NROWS,NCOL
C
      CALL CHANS (NROWS,NCOL)
      STOP 
   13 WRITE(6,6010) ' NRWOS * NCOL .NE. NDIM ',NROWS ,NCOL,NDIM
      STOP
   12 NROWS=32
      GOTO 11
 6010 FORMAT(' *CHANS*: ', A,/,(6I10))
 6011 FORMAT(' *CHANS*: ', A)
      END
C******************************************************************
      SUBROUTINE PARLES (NCHAN,NDIM,NSET)
      IMPLICIT REAL (A-H,O-Z)
C
      INTEGER I1,I2,I3
      REAL F

      REWIND NCHAN
      READ(NCHAN,END=10) ITIME,INAME,ILEV,ISIZE
      REWIND NCHAN

      NSET=0
      NDIM=0  
    1 READ (NCHAN,END=100) I1,I2,I3,I4
      NDIM=MAX0(NDIM,I4)
      READ (NCHAN,END=101) 
      NSET=NSET+1
      GOTO1
      
  101 PRINT *,' MORE HEADER THEN DATA '     
  100 PRINT 1000,NCHAN,NDIM,NSET
      REWIND NCHAN
      RETURN    
   10 PRINT *,'INPUT FILE ON TAPE',NCHAN,' IS EMPTY'
      STOP 'MISTAKE'
 1000 FORMAT (1X,' TAPE=',I2,',  NDIM=',I5,',  NSET=',I5)
      END

C*****************************************************************
      SUBROUTINE CHANS (IY,IX)
      IMPLICIT REAL (A-H,O-Z)
C
C      PARAMETER(IE77=9E12)
C      DIMENSION X(IX,IY)
C      POINTER (NX,X)
C      ALLOCATE (NX,IX*IY)
C AND ACTIVATE THE FOLLOWING LINES
      PARAMETER(IX77=20000,IY77=1000)
      DIMENSION X(IX77,IY77)
C NOW YOU HAVE DONE ALL CHANGES FOR STANDARD FORTRAN.
      DATA ONE/1.0/
C
C*  CHECK THE INPUT FILE.
C
      REWIND 1
      READ(1,END=10) ITIME,INAME,ILEV,ISIZE
      REWIND 1
      WRITE(6,6010) ' FIRST RECORD:',ITIME,INAME,ILEV,ISIZE
C
C*  MAIN LOOP.
C
      ICNT=0
  100 READ(1,END=9920) ITIME,INAME,ILEV,ISIZE
        READ(1,END=15) ((X(I,J),I=1,IX),J=1,IY)
        ICNT=ICNT+1
        WRITE(2) ITIME,INAME,ILEV,IX*IY
        WRITE (2) ((X(I,IY+1-J),I=1,IX),J=1,IY)
      GOTO 100
   10 WRITE(6,6010) ' INPUT FILE IS EMPTY '
         GOTO 9910
   15 WRITE(6,6010) ' UNEXPECTED *EOF* DURING READ '
 9910 CONTINUE
C
C*  ABNORMAL TERMINATION.
C
         PRINT *, '*CHANS*: FATAL END '
         CALL MINONE(ONE,XNEGO)
         ZZ=SQRT(XNEGO)
         STOP
C
C*  NORMAL TERMINATION.
C
 9920 WRITE(6,6010) '  LAST RECORD:',ITIME,INAME,ILEV,ISIZE
      WRITE(6,6060) ICNT
      PRINT *, ' *CHANS*: NORMAL END'
      RETURN
 6010 FORMAT(' *CHANS*: ', A,(6I10))
 6060 FORMAT(' *CHANS* -- READ ',I10,' RECORDS.')
      END
      SUBROUTINE MINONE(XIN,XOUT)
        XOUT=-XIN
      RETURN
      END
EOF
f77 chanoso.f -o chanoso.x
chanoso.x $1 $2 << M
$3
M
cp fort.2 $2
rm fort.1 fort.2
rm chanoso.x chanoso.f
exit