if [ $# -lt 2 ]
then
  echo " "
  echo "            *****     TP     *****      "
  echo " "
  echo " Transposition of a data matrix given in EXTRA format"
  echo " "
  echo "          Syntax: tp <file1> <file2>"
  echo "  "
  exit 0
else
 echo " "
 echo "            *****     TP     *****      "
 echo " "
 echo "file1="$1
 echo "file2="$2
 cp $1 fort.1
fi
cat > tp.f <<EOF
      PROGRAM Ptp
      IMPLICIT REAL (A-H,O-Z)
      PRINT *,'0**** TP (INPUT,OUTPUT,TAPE1,TAPE2)****'
      PRINT *,'   '
      NCHAN=1
      CALL PARLES (NCHAN,NDIM,NSET)
      CALL TP (NDIM,NSET)     
      STOP 
      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) F
      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 TP (IE,NTS)
      IMPLICIT REAL (A-H,O-Z)
C
      PARAMETER (UNDEF=0.9E+10,UNDEFE=.9*UNDEF)
      INTEGER ID(4),IID(4)
C IF YOU USE THIS PROGRAM NOT ON A CRAY-2 DELETE THE
C FOLLOWING TWO LINES
C      PARAMETER(IE77=9E12,NTS77=9E12)
C      REAL X(IE),XS(NTS)
C      POINTER (NPX,X),(NPC,XS)
C      ALLOCATE (NPX,IE)
C      ALLOCATE (NPC,NTS)
C AND ACTIVATE THE FOLLOWING TWO LINES
      PARAMETER(IE77=3000,NTS77=3000)
      REAL X(IE77),XS(NTS77)
C NOW YOU HAVE DONE ALL CHANGES FOR STANDARD FORTRAN.
      IM = 0
C
C   check input file 1
C
      REWIND 1
      READ(1,END=401) IID
      IE4 = IID(4)
      PRINT *, 'TAPE1:vector length: ', IE4
      PRINT *, 'TAPE1:first record:  ', IID
      REWIND 1
   11 READ(1,END=200) ID
        READ(1)
        IF(ID(2).NE.IID(2)) GOTO 400
        IF(ID(3).NE.IID(3)) GOTO 400
        IF(ID(4).NE.IID(4)) GOTO 400
        IM = IM + 1
        GOTO 11
C
  200 PRINT *, 'TAPE1:last record:   ', ID
      PRINT *, 'TAPE1:number of records found: ',IM
      PRINT *,' '

      REWIND 2
      DO 100 I = 1,IE4
        REWIND 1
        DO 110 K = 1,IM
          READ(1)
          READ(1) (X(J),J=1,I)
  110     XS(K) = X(I)
        WRITE(2) I,IID(2),IID(3),IM
        IF(I.EQ.1) PRINT*, 'TAPE2:first record:  ',
     +                      I,IID(2),IID(3),IM
        IF(I.EQ.IE4) PRINT*, 'TAPE2:last record:   ',
     +                      I,IID(2),IID(3),IM
  100   WRITE(2) (XS(K),K=1,IM) 
      PRINT *, 'TAPE2:number of writen records: ',IE4
      STOP
C
C  fata ends
  400 PRINT *, ' First Record and ',IM,' Record not equal '
      PRINT *,' Record    :',IID
      PRINT *,' Record ',IM,':',ID
      STOP 'TP'
  401 PRINT *, 'Input file 1 is empty '
      STOP 'TP'
      END
EOF
f77 tp.f -o tp.x
tp.x $1 $2
cp fort.2 $2
rm fort.1 fort.2
rm tp.x tp.f
exit