#! /bin/sh
#
#    Last change: Mihai Dima 2001
#
if [ $# -lt 3 ]
then
  echo " "
  echo "            *****     APPEND     *****      "
  echo " "
  echo "      The program adds <file2> after <file1>"
  echo " "
  echo "     Syntax: append <file1> <file2> <file3>"
  echo " "
  exit 0
else
 echo " "
 echo "            *****     APPEND     *****      "
 echo " "
 echo "file1="$1
 echo "file2="$2
 echo "file3="$3
 cp $1 fort.1
 cp $2 fort.2
fi
cat > append.f <<EOF
      PROGRAM Pappend
      IMPLICIT REAL (A-H,O-Z)
      	
      PRINT *,'0**** APPEND (,OUTPUT,TAPE1,TAPE2,TAPE3)****'
      PRINT *,'   '
      NCHAN=1
      CALL PARLES (NCHAN,NDIM1,NSET)
      NCHAN=2
      CALL PARLES (NCHAN,NDIM2,NSET)
      NDIM=MAX0(NDIM1,NDIM2)
      CALL APPEND (NDIM)     
      STOP 
      END
C******************************************************************
      SUBROUTINE PARLES (NCHAN,NDIM,NSET)
      IMPLICIT REAL (A-H,O-Z)
C
      INTEGER I1,I2,I3

      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)
      NSET=NSET+1
      GOTO1
      
  101 PRINT *,' MORE HEADER THEN DATA '     
  100 PRINT 1000,NCHAN,NDIM,NSET
      RETURN    
 1000 FORMAT (1X,' TAPE=',I2,',  NDIM=',I5,',  NSET=',I5)
      END
C*****************************************************************
      SUBROUTINE APPEND(IE)
      IMPLICIT REAL (A-H,O-Z)
C
C IF YOU USE THIS PROGRAM NOT ON A CRAY-2 DELETE THE
C THE FOLLOWING TWO LINES
C      PARAMETER (IE77=9E12)
C      DIMENSION X(IE), ID(4)
C      POINTER (NPX,X)
C      ALLOCATE (NPX,IE)
C AND ACTIVATE THE FOLLOWING LINES
       PARAMETER (IE77=20000)
       DIMENSION X(IE77), ID(4)
C NOW YOU HAVE DONE ALL CHANGES FOR STANDARD FORTRAN.
      N1 = 0
      N2 = 0
      PRINT *, 'APPEND(TAPE1,TAPE2,TAPE3)'
      REWIND 1
      REWIND 2
C
    1 READ(1,END=2) ID
      READ(1) (X(I),I=1,ID(4))
      WRITE(3) ID
      WRITE(3) (X(I),I=1,ID(4))
      N1 = N1 + 1
      GOTO 1
C
    2 READ(2,END=3) ID
      READ(2) (X(I),I=1,ID(4))
      WRITE(3) ID
      WRITE(3) (X(I),I=1,ID(4))
      N2 = N2 + 1
      GOTO 2
C
    3 PRINT *, 'number of records on TAPE3: ', N1+N2
      STOP
      END
EOF
f77 append.f -o append.x
append.x $1 $2 $3
cp fort.3 $3
rm fort.1 fort.2 fort.3
rm append.x appendmore more .f
exit