#! /bin/sh
#
#    Last change: Mihai Dima 2001
#
if [ $# -lt 4 ]
then
  echo " "
  echo "        *****     SPLITF     *****      "
  echo " "
  echo "  The program splits one file (file1) "
  echo "   in two files (file2) and (file3)"
  echo " Last record in <file2> corresponds to <date>"
  echo " "
  echo "  Syntax: splitf <file1> <file2> <file3> <date>"
  echo " "
  exit 0
else
 echo " "
 echo "        *****     SPLITF     *****      "
 echo " "
 echo "file1="$1
 echo "file2="$2
 echo "file3="$3
 echo "Last date in file2="$4
 cp $1 fort.1
fi
cat > splitf.f <<EOF
      PROGRAM Psplitf
      IMPLICIT REAL (A-H,O-Z)

      PRINT *,'0**** SPLIT  (INPUT,OUTPUT,TAPE1,TAPE2)****'
      PRINT *,'   '
      NCHAN=1
      CALL PARLES (NCHAN,NDIM,NSET)
      CALL SPLIT  (NDIM,NSET)     
      STOP 
      END
C******************************************************************
      SUBROUTINE PARLES (NCHAN,NDIM,NSET)
      IMPLICIT REAL (A-H,O-Z)
C
C SUBROUTINE FOR READING THE DIMENSION PARAMETER OF THE FIELDS ON
C
      INTEGER I1,I2,I3

      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 SPLIT (IE,NTS)
      IMPLICIT REAL (A-H,O-Z)
C
C      PROGRAM SPLIT(INPUT,OUTPUT,TAPE1,TAPE2,TAPE3,TAPE5=INPUT,
C                    TAPE6=OUTPUT).
C
C      *INPUT*  - NSPLIT(I.E. LAST DATE OF DATA TO BE WRITTEN ON *TAPE2*).
C      *OUTPUT* - CONTROL PRINT.
C      *TAPE1*  - T21 STANDARD INPUTFILE.
C      *TAPE2*  - T21 STANDARD OUTPUTFILE: FIRST PART OF *TAPEE1*.
C      *TAPE3*  - T21 STANDARD OUTPUTFILE: SECOND PART OF *TAPE1*.
C
C IF YOU USE THIS PROGRAM NOT ON A CRAY-2 DELETE THE
C 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.
      DATA ONE/1./

      PRINT *,' ENTER LAST DATE ON TAPE2: YYMMDD '
      READ(5,'(I10)',END=10) NSPLIT
      PRINT *, 'requested last date on tape 2: ', NSPLIT, '(YYMMDD)'
      REWIND 1
      READ(1,END=11) ID
      PRINT *, 'first record on tape1: ', ID
      ICNT = 1
      IF(ID(4).GT.IE .OR. ID(4).GT.IE77) GOTO 12
C
C*  MAIN LOOP.
C
      REWIND 1
      ICNT=0
      IT2 = 0
      IT3 = 0
      ID2A = 0
      ID3A = 0
      ID2E = 0
      ID3E = 0
  100 READ(1,END=9920) ID
        ICNT=ICNT+1
        IF(ID(4).GT.IE .OR. ID(4).GT.IE77) GOTO 12
        READ(1,END=13) (X(I),I=1,ID(4))
        IF(ID(1).LE.NSPLIT) THEN
          IF(ID2A.EQ.0) ID2A = ID(1)
          ID2E = ID(1)
          IT2 = IT2 + 1
          WRITE(2) ID
          WRITE(2) (X(I),I=1,ID(4))
        ELSE
          IF(ID3A.EQ.0) ID3A = ID(1)
          ID3E = ID(1)
          IT3 = IT3 + 1
          WRITE(3) ID
          WRITE(3) (X(I),I=1,ID(4))
        END IF
      GOTO 100
C
C*  NORMAL TERMINATION.
C
 9920 PRINT *, ' Last record on tape1: ', ID
      PRINT *, ' ', ICNT, ' records found'
      PRINT *, '  '
      PRINT *, 'First record on tape2: ', ID2A, ' last: ', ID2E,
     &         '  total of ',IT2, ' records'
      PRINT *, 'First record on tape3: ', ID3A, ' last: ', ID3E,
     &         '  total of ',IT3, ' records'
      STOP
C
C*  ABNORMAL TERMINATION.
C
   10 PRINT *, '<INPUT> EMPTY: NO DATE GIVEN '
         GOTO 9910
   11 PRINT *, '<TAPE1> EMPTY '
         GOTO 9910
   12 PRINT *, 'RECORD TOO LONG: ', ID, IE, ICNT
         GOTO 9910
   13 PRINT *, 'READ ERROR: ', ID, ICNT
         GOTO 9910
 9910 CONTINUE
C
C*  CAUSE AN ABORT.
C
      PRINT *,' *SPLIT*: FATAL END'
      CALL MINONE(ONE,XNEGO)
      ZZ=SQRT(XNEGO)
      STOP
      END
      SUBROUTINE MINONE(XIN,XOUT)
        XOUT=-XIN
      RETURN
      END
EOF
f77 splitf.f -o splitf.x
splitf.x $1 $2 $3 << M
$4
M
cp fort.2 $2
cp fort.3 $3
rm fort.1 fort.2 fort.3
rm splitf.x splitf.f
exit