if [ $# -lt 3 ]
then
  echo " "
  echo "            *****   MERGEF      *****      "
  echo " "
  echo "    The program merge two files, record by record"
  echo " "
  echo "               NSET1=NSET2+NSET2"
  echo "            RECORD3=RECORD1+RECORD2"
  echo " "
  echo "      Syntax: mergef <file1> <file2> <file3>"
  echo "  "
  exit 0
else
 echo " "
 echo "            *****   MERGEF      *****      "
 echo " "
 echo "file1="$1
 echo "file2="$2
 echo "file3="$3
 cp $1 fort.1
 cp $2 fort.2
fi
cat >mergef.f <<EOF
      PROGRAM Pmerge
      IMPLICIT REAL (A-H,O-Z)
      
      PRINT *,'  '
      NCHAN=1
      CALL PARLES (NCHAN,NDIM1,NSET)
      NCHAN=2
      CALL PARLES (NCHAN,NDIM2,NSET)
      NDIM=NDIM1+NDIM2
      CALL APPF (NDIM,NDIM1,NDIM2)       
      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,/)
C
      END
C*******************************************************************
      SUBROUTINE APPF (NDIM,NDA,NDB)
      IMPLICIT REAL (A-H,O-Z)
C
      DIMENSION IDATE(4), JDATE(4)
C IF YOU USE THIS PROGRAM NOT ON A CRAY-2 ACTIVATE THE
C FOLLOWING TWO LINES IN THIS PROGRAM !
      PARAMETER (NDIM77=20000) 
      DIMENSION X(NDIM77), Y(NDIM77)
C AND DELETE THE FOLLOWING LINE
C      DIMENSION X(NDIM)
C      POINTER (NPX,X)
C      ALLOCATE (NPX,NDIM)
C NOW YOU HAVE DONE ALL CHANGES FOR STANDARD FORTRAN.
      REWIND 1
      REWIND 2
      N = 0
    1 READ(1,END=10) IDATE
      READ(2,END=10) JDATE
      NFELD=IDATE(4)+JDATE(4)
      N = N + 1
      IF (N.EQ.1) PRINT 1000,N,IDATE,' FIRST RECORD ON TAPE 1 '
      IF (N.EQ.1) PRINT 1000,N,JDATE,' FIRST RECORD ON TAPE 2 '
      IF (N.EQ.1) PRINT 1000,N,(IDATE(I),I=1,3),NFELD
     +                   ,' FIRST RECORD ON TAPE 3 '
      READ(1) (X(I),I=1,IDATE(4))
      READ(2) (X(I),I=IDATE(4)+1,NFELD)
      WRITE(3) (IDATE(I),I=1,3),NFELD
      WRITE (3) (X(I), I=1,NFELD)
      GOTO 1
  10  PRINT 1000,N,IDATE,' LAST RECORD ON TAPE 1 '
      PRINT 1000,N,JDATE,' LAST RECORD ON TAPE 2 '
      PRINT 1000,N,(IDATE(I),I=1,3),NFELD
     +                  ,' LAST RECORD ON TAPE 3 '
      PRINT *, ' TOTAL NUMBER OF RECORDS FOUND: ', N
      STOP 'O.K.'
 1000 FORMAT(2I10,3I6,A30)
      END
EOF
f77 mergef.f -o mergef.x
mergef.x $1 $2 $3
cp fort.3 $3
rm fort.1 fort.2 fort.3
rm mergef.x mergef.f
exit