#! /bin/sh
#
#    Last change: Mihai Dima 2001
#
if [ $# -lt 3 ]
then
  echo " "
  echo "            *****     CDIV     *****      "
  echo " "
  echo "  The program divides the input data (file1) with"
  echo "            a constant field (file2)"
  echo " "
  echo "     Syntax: cdiv <file1> <file2> <file3>"
  echo " "
  exit 0
else
 echo " "
 echo "            *****     CDIV     *****      "
 echo " "
 echo "file1="$1
 echo "file2="$2
 echo "file3="$3
 cp $1 fort.1
 cp $2 fort.2
fi
cat > cdiv.f <<EOF
      PROGRAM Pcdiv
      IMPLICIT REAL (A-H,O-Z)
      
      PRINT *,' ****  *CDIV(,,X(T),C,X(T)/C)* '
      PRINT *,'   '
      NCHAN=1
      CALL PARLES (NCHAN,NDIM,NSET)
      CALL CDIV(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
      REWIND NCHAN
      RETURN    
 1000 FORMAT (1X,' TAPE=',I2,',  NDIM=',I5,',  NSET=',I5)
      END
C**************************************************************
      SUBROUTINE CDIV(IE)
      IMPLICIT REAL (A-H,O-Z)
C
      PARAMETER (UNDEF=0.9E+10, UNDEFE=0.9*UNDEF)
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(IE)
C      POINTER (NPX,X),(NPC,C)
C      ALLOCATE (NPX,IE)
C      ALLOCATE (NPC,IE)
C AND ACTIVATE THE FOLLOWING LINES
      PARAMETER (IE77=20000)
      DIMENSION X(IE77), ID(4), C(IE77)
C NOW YOU HAVE DONE ALL CHANGES FOR STANDARD FORTRAN.
      DATA ONE/1.0/
C
C*  INPUT OF CONSTANT FIELD.
C
      REWIND 2
      READ(2,END=19) ID
      WRITE(6,6010) ' CONSTANT FIELD ', ID
      IF(ID(4).GT.IE77.OR.ID(4).GT.IE) GOTO 20
      IEE = ID(4)
      READ(2) (C(I),I=1,IEE)
      DO 99 I=1,IEE
         IF (C(I).NE.0.0.AND.C(I).LT.UNDEFE) THEN
             C(I)=1./C(I)
         ELSE
             C(I)=UNDEF
         ENDIF
   99 CONTINUE
C
C*  MAIN LOOP.
C
      REWIND 1
      ICNT=0
  100 READ(1,END=9920) ID
        ICNT=ICNT+1
        IF(ICNT.EQ.1)THEN
          WRITE(6,6010) ' FIRST RECORD:', ID
        ENDIF
        IF(ID(4).NE.IEE) GOTO 12
        READ(1,END=13) (X(I),I=1,IEE)
        DO 101 I = 1,ID(4)
          IF(X(I).GT.UNDEFE .OR. C(I).GT.UNDEFE) THEN
              X(I) = UNDEF
            ELSE
              X(I) = X(I) * C(I)
            ENDIF
  101     CONTINUE
        WRITE(3) ID
        WRITE(3) (X(I),I=1,ID(4))
      GOTO 100
C
C*  ABNORMAL TERMINATION.
C
   12 WRITE(6,6010) 'RECORD LENGTH INSUFFICIENT: ', ID, IEE, ICNT
         GOTO 9910
   13 WRITE(6,6010) ' READ ERROR: ', ID, ICNT
         GOTO 9910
   19 WRITE(6,6010) ' NO CONSTANT FIELD FOUND '
         GOTO 9910
   20 WRITE(6,6010) ' CONSTANT FIELD TOO LONG: ', ID
         GOTO 9910
 9910 CONTINUE
C
C*  CAUSE AN ABORT.
C
         PRINT *,'*CDIV*: FATAL END '
         CALL MINONE(ONE,XNEGO)
         ZZ=SQRT(XNEGO)
         STOP
C
C*  NORMAL TERMINATION.
C
 9920 WRITE(6,6010) ' LAST RECORD: ', ID, ICNT
      STOP ' REGULAR END! '
 6010 FORMAT(' *CDIV*: ', A, 4I10,/1X,5I10)
      END
      SUBROUTINE MINONE(XIN,XOUT)
        XOUT=-XIN
      RETURN
      END
EOF
f77 cdiv.f -o cdiv.x
cdiv.x $1 $2 $3
cp fort.3 $3
rm fort.1 fort.2 fort.3
rm cdiv.x cdiv.f
exit